This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering a...This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a...Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.展开更多
The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this ...The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.展开更多
Landfalling typhoons can cause disasters over large regions.The government and emergency responders need to take measures to mitigate disasters according to the forecast of landfall position,while slight timing error ...Landfalling typhoons can cause disasters over large regions.The government and emergency responders need to take measures to mitigate disasters according to the forecast of landfall position,while slight timing error can be ignored.The reliability of operational model forecasts of typhoon landfall position needs to be evaluated beforehand,according to the forecasts and observation of historical cases.In the evaluation of landfalling typhoon track,the traditional method based on point-to-point matching methods could be influenced by the predicted typhoon translation speed.Consequently,the traditional track evaluation method may result in a large track error even if the predicted landfall position is close to observation.The purpose of this paper is to address the above issue using a simple evaluation method of landfalling typhoon track forecast based on the time neighborhood approach.In this new method,the timing error was lessened to highlight the importance of the position error during the landfall of typhoon.The properties of the time neighborhood method are compared with the traditional method based on numerical forecast results of 12 landfalling typhoon cases.Results demonstrated that the new method is not sensitive to the sampling frequency,and that the difference between the time neighborhood and traditional method will be more obvious when the moving speed of typhoon is moderate(between 15−30 km h^(−1)).The time neighborhood concept can be easily extended to a broader context when one attempts to examine the position error more than the timing error.展开更多
A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the c...A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.展开更多
In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its ...In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.展开更多
In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finit...In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.展开更多
We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatia...We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.展开更多
In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The ob...In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.展开更多
The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the co...The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.展开更多
The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experi...The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.展开更多
The particle residence time distribution(RTD)and axial dispersion coefficient are key parameters for the design and operation of a pressurized circulating fluidized bed(PCFB).In this study,the effects of pressure(0.1-...The particle residence time distribution(RTD)and axial dispersion coefficient are key parameters for the design and operation of a pressurized circulating fluidized bed(PCFB).In this study,the effects of pressure(0.1-0.6 MPa),fluidizing gas velocity(2-7 m·s^(-1)),and solid circulation rate(10-90 kg·m^(-2)·s^(-1))on particle RTD and axial dispersion coefficient in a PCFB are numerically investigated based on the multiphase particle-in-cell(MP-PIC)method.The details of the gas-solid flow behaviors of PCFB are revealed.Based on the gas-solid flow pattern,the particles tend to move more orderly under elevated pressures.With an increase in either fluidizing gas velocity or solid circulation rate,the mean residence time of particles decreases while the axial dispersion coefficient increases.With an increase in pressure,the core-annulus flow is strengthened,which leads to a wider shape of the particle RTD curve and a larger mean particle residence time.The back-mixing of particles increases with increasing pressure,resulting in an increase in the axial dispersion coefficient.展开更多
An effective Dynamic Moire method was presented to determine the critical time of crack instable propagation in rock dynamic fracture. Two pieces of grating were installed near the notch of Short Rod specimen to form ...An effective Dynamic Moire method was presented to determine the critical time of crack instable propagation in rock dynamic fracture. Two pieces of grating were installed near the notch of Short Rod specimen to form the Moire fringes, then the COD versus time could be monitored from the movement of the Moire fringes, and finally the critical time could be determined from the velocity of COD. This method was also compared with another one.It could be concluded that the critical time determined by Dynamic Moire method corresponds with that of the Transmitted Wave method at the loading rates from 103 to 104 MPa·m(1/2)·S(-1).展开更多
In order to evaluate the efficiency of the automated storage/retrieval system(AS/RS)accurately,and compare different layouts of the AS/RS using mean travel time,under randomized storage conditions,an exact,geometry-ba...In order to evaluate the efficiency of the automated storage/retrieval system(AS/RS)accurately,and compare different layouts of the AS/RS using mean travel time,under randomized storage conditions,an exact,geometry-based analytical model is presented.The model can be used to compute the expected single-command and dual-command travel time for a storage/retrieval(S/R)machine which can travel simultaneously horizontally and vertically as it moves along a storage aisle.The rack may be either square in time or non square in time.Additionally,the alternative layouts of the AS/RS and travel-time models are examined.Comparing with setting the I/O point at the left-lower corner of the rack,setting the I/O point at any point at the vertical edge can help enhance the efficiency of the AS/RS.展开更多
To study the effects of different storage methods and time on content of nutrients in biogas slurry of straw, two storage methods were carried out on biogas slurry between open storage and airtight storage conditions ...To study the effects of different storage methods and time on content of nutrients in biogas slurry of straw, two storage methods were carried out on biogas slurry between open storage and airtight storage conditions at normal atmospheric temperature. The contents of N, P, K, and organic matter in biogas slurry of straw were determined in different storage times. The results showed that: during the pro-cess of biogas slurry storage, little change occurred in the content of the organic matter while the total content of N, P, K significantly declined; up to 50 days, the total content of N, P, K reduced to nearly 80%-90%. Because the contents of N, P, K in biogas slurry reduced less in airtight storage conditions so that a better re-sult was found on airtight storage methods than open storage methods in fertilizer field of biogas slurry of straw.展开更多
The effects of the calorimetric buffer solutions were investigated while the two colorimetric reactions of AI-ferron complex and Fe-ferron complex occurred individually, and the effects of the testing wavelength and t...The effects of the calorimetric buffer solutions were investigated while the two colorimetric reactions of AI-ferron complex and Fe-ferron complex occurred individually, and the effects of the testing wavelength and the pH of the solutions were also investigated. A timed complexatian colorimetric analysis method of Al-Fe-ferron in view of the total concentration of {AI + Fe} was then established to determine the species distribution of polymeric Al-Fe. The testing wavelength was recommended at 362 net and the testing pH value was 5. With a comparison of the ratios of n(Al)/n(Fe), the standard adsorption curves of the polymeric Al-Fe solutions were derived from the experimental results. Furthermore, the solutions' composition were carious in both the molar n(Al)/n(Fe) ratios, i.e. 0/0, 5/5, 9/1 and 0/10, and the concentrations associated with the total ( Al + Fe which ranged from 10(-5) to 10(-4) mol/L..展开更多
The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
In endurance time(ET) method structures are subjected to a set of predesigned intensifying excitations. These excitations are produced in a way that their response spectrum, while complying with a specifi ed spectrum,...In endurance time(ET) method structures are subjected to a set of predesigned intensifying excitations. These excitations are produced in a way that their response spectrum, while complying with a specifi ed spectrum, intensifi es with time so they can be used approximately to simulate the average effects of several ground motions scaled to different intensities. In this paper applicability of the ET method for evaluating collapse potential of buildings is investigated. A set of four steel moment frames is used for collapse assessment. The process of using ET method in collapse evaluation is explained and the results are compared with incremental dynamic analysis(IDA) results. It is shown that although the computational effort using the ET method is much less than the IDA analysis, the results of both methods are consistent. Finally collapse fragility curves using ET and IDA methods are produced and it is shown that the probabilities of collapse in different hazard levels are also consistent.展开更多
基金funded by the National Key Research and Development Program of China(No.2021YFB2600704)the National Natural Science Foundation of China(No.52171272)the Significant Science and Technology Project of the Ministry of Water Resources of China(No.SKS-2022112).
文摘This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
文摘Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.
文摘The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.
基金supported by the National Natural Science Foundation of China(Grant Nos.U1811464,U2142213)Guangdong Basic and Applied Basic Research Foundation(Grants Nos.2020A1515110275,2020A1515110040,2022A1515011870)the Special program for innovation and development of China Meteorological Administration(CXFZ2021Z006,CXFZ2022P026).
文摘Landfalling typhoons can cause disasters over large regions.The government and emergency responders need to take measures to mitigate disasters according to the forecast of landfall position,while slight timing error can be ignored.The reliability of operational model forecasts of typhoon landfall position needs to be evaluated beforehand,according to the forecasts and observation of historical cases.In the evaluation of landfalling typhoon track,the traditional method based on point-to-point matching methods could be influenced by the predicted typhoon translation speed.Consequently,the traditional track evaluation method may result in a large track error even if the predicted landfall position is close to observation.The purpose of this paper is to address the above issue using a simple evaluation method of landfalling typhoon track forecast based on the time neighborhood approach.In this new method,the timing error was lessened to highlight the importance of the position error during the landfall of typhoon.The properties of the time neighborhood method are compared with the traditional method based on numerical forecast results of 12 landfalling typhoon cases.Results demonstrated that the new method is not sensitive to the sampling frequency,and that the difference between the time neighborhood and traditional method will be more obvious when the moving speed of typhoon is moderate(between 15−30 km h^(−1)).The time neighborhood concept can be easily extended to a broader context when one attempts to examine the position error more than the timing error.
基金the National Natural Science Foundation of China(No.11701103)the Young Top-notch Talent Program of Guangdong Province of China(No.2017GC010379)+4 种基金the Natural Science Foundation of Guangdong Province of China(No.2022A1515012147)the Project of Science and Technology of Guangzhou of China(No.202102020704)the Opening Project of Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University of China(2021023)the Science and Technology Development Fund,Macao SAR(File No.0005/2019/A)the University of Macao of China(File Nos.MYRG2020-00035-FST,MYRG2018-00047-FST).
文摘A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.
文摘In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.
基金supported by the State Key Program of National Natural Science Foundation of China(11931003)the National Natural Science Foundation of China(41974133)。
文摘In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.
文摘We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.
文摘In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.
文摘The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.
基金supported by National Natural Science Foundation of China(Nos.92066108 and 51277061)。
文摘The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.
基金Financial support of this work by National Natural Science Foundation of China(51976037)。
文摘The particle residence time distribution(RTD)and axial dispersion coefficient are key parameters for the design and operation of a pressurized circulating fluidized bed(PCFB).In this study,the effects of pressure(0.1-0.6 MPa),fluidizing gas velocity(2-7 m·s^(-1)),and solid circulation rate(10-90 kg·m^(-2)·s^(-1))on particle RTD and axial dispersion coefficient in a PCFB are numerically investigated based on the multiphase particle-in-cell(MP-PIC)method.The details of the gas-solid flow behaviors of PCFB are revealed.Based on the gas-solid flow pattern,the particles tend to move more orderly under elevated pressures.With an increase in either fluidizing gas velocity or solid circulation rate,the mean residence time of particles decreases while the axial dispersion coefficient increases.With an increase in pressure,the core-annulus flow is strengthened,which leads to a wider shape of the particle RTD curve and a larger mean particle residence time.The back-mixing of particles increases with increasing pressure,resulting in an increase in the axial dispersion coefficient.
文摘An effective Dynamic Moire method was presented to determine the critical time of crack instable propagation in rock dynamic fracture. Two pieces of grating were installed near the notch of Short Rod specimen to form the Moire fringes, then the COD versus time could be monitored from the movement of the Moire fringes, and finally the critical time could be determined from the velocity of COD. This method was also compared with another one.It could be concluded that the critical time determined by Dynamic Moire method corresponds with that of the Transmitted Wave method at the loading rates from 103 to 104 MPa·m(1/2)·S(-1).
基金The National Key Technology R&D Program of China during the 11th Five-Year Plan Period(No.2006BAH02A06)
文摘In order to evaluate the efficiency of the automated storage/retrieval system(AS/RS)accurately,and compare different layouts of the AS/RS using mean travel time,under randomized storage conditions,an exact,geometry-based analytical model is presented.The model can be used to compute the expected single-command and dual-command travel time for a storage/retrieval(S/R)machine which can travel simultaneously horizontally and vertically as it moves along a storage aisle.The rack may be either square in time or non square in time.Additionally,the alternative layouts of the AS/RS and travel-time models are examined.Comparing with setting the I/O point at the left-lower corner of the rack,setting the I/O point at any point at the vertical edge can help enhance the efficiency of the AS/RS.
基金Supported by Key Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period(2012BAJ21B04-04)Shandong Environment Bottlenecks Analysis Project(SDHBPJ-ZB-05)~~
文摘To study the effects of different storage methods and time on content of nutrients in biogas slurry of straw, two storage methods were carried out on biogas slurry between open storage and airtight storage conditions at normal atmospheric temperature. The contents of N, P, K, and organic matter in biogas slurry of straw were determined in different storage times. The results showed that: during the pro-cess of biogas slurry storage, little change occurred in the content of the organic matter while the total content of N, P, K significantly declined; up to 50 days, the total content of N, P, K reduced to nearly 80%-90%. Because the contents of N, P, K in biogas slurry reduced less in airtight storage conditions so that a better re-sult was found on airtight storage methods than open storage methods in fertilizer field of biogas slurry of straw.
基金TheNationalNaturalScienceFoundationofChina (No .2 96 770 0 4)
文摘The effects of the calorimetric buffer solutions were investigated while the two colorimetric reactions of AI-ferron complex and Fe-ferron complex occurred individually, and the effects of the testing wavelength and the pH of the solutions were also investigated. A timed complexatian colorimetric analysis method of Al-Fe-ferron in view of the total concentration of {AI + Fe} was then established to determine the species distribution of polymeric Al-Fe. The testing wavelength was recommended at 362 net and the testing pH value was 5. With a comparison of the ratios of n(Al)/n(Fe), the standard adsorption curves of the polymeric Al-Fe solutions were derived from the experimental results. Furthermore, the solutions' composition were carious in both the molar n(Al)/n(Fe) ratios, i.e. 0/0, 5/5, 9/1 and 0/10, and the concentrations associated with the total ( Al + Fe which ranged from 10(-5) to 10(-4) mol/L..
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
文摘In endurance time(ET) method structures are subjected to a set of predesigned intensifying excitations. These excitations are produced in a way that their response spectrum, while complying with a specifi ed spectrum, intensifi es with time so they can be used approximately to simulate the average effects of several ground motions scaled to different intensities. In this paper applicability of the ET method for evaluating collapse potential of buildings is investigated. A set of four steel moment frames is used for collapse assessment. The process of using ET method in collapse evaluation is explained and the results are compared with incremental dynamic analysis(IDA) results. It is shown that although the computational effort using the ET method is much less than the IDA analysis, the results of both methods are consistent. Finally collapse fragility curves using ET and IDA methods are produced and it is shown that the probabilities of collapse in different hazard levels are also consistent.