This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we stu...This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.展开更多
Leveraging big data signal processing offers a pathway to the development of artificial intelligencedriven equipment.The analysis of fluid flow signals and the characterization of fluid flow behavior are of critical i...Leveraging big data signal processing offers a pathway to the development of artificial intelligencedriven equipment.The analysis of fluid flow signals and the characterization of fluid flow behavior are of critical in two-phase flow studies.Significant research efforts have focused on discerning flow regimes using various signal analysis methods.In this review,recent advances in time series signals analysis algorithms for stirred tank reactors have been summarized,and the detailed methodologies are categorized into the frequency domain methods,time-frequency domain methods,and state space methods.The strengths,limitations,and notable findings of each algorithm are highlighted.Additionally,the interrelationships between these methodologies have also been discussed,as well as the present progress achieved in various applications.Future research directions and challenges are also predicted to provide an overview of current research trends in data mining of time series for analyzing flow regimes and chaotic signals.This review offers a comprehensive summary for extracting and characterizing fluid flow behavior and serves as a theoretical reference for optimizing the characterization of chaotic signals in future research endeavors.展开更多
Three functional models, polynomial, spectral analysis, and modified AR model, are studied and compared in fitting and predicting clock deviation based on the data sequence derived from two-way satellite time and freq...Three functional models, polynomial, spectral analysis, and modified AR model, are studied and compared in fitting and predicting clock deviation based on the data sequence derived from two-way satellite time and frequency transfer. A robust equivalent weight is applied, which controls the significant influence of outlying observations. Some conclusions show that the prediction precision of robust estimation is better than that of LS. The prediction precision calculated from smoothed observations is higher than that calculated from sampling observations. As a count of the obvious period variations in the clock deviation sequence, the predicted values of polynomial model are implausible. The prediction precision of spectral analysis model is very low, but the principal periods can be determined. The prediction RMS of 6-hour extrapolation interval is Ins or so, when modified AR model is used.展开更多
We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not...We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not possible.As shown in Boscarino et al.(J.Sci.Comput.68:975-1001,2016)for Runge-Kutta methods,these semi-implicit techniques give a great flexibility,and allow,in many cases,the construction of simple linearly implicit schemes with no need of iterative solvers.In this work,we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype lineal'advection-diffusion equation and in the setting of strong stability preserving(SSP)methods.Our findings are demonstrated on several examples,including nonlinear reaction-diffusion and convection-diffusion problems.展开更多
In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for c...In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.展开更多
A discrete event system is a dynamical system whose state evolves in time by the occurrence of events at possibly irregular time intervals. Timed Petri nets are a graphical and mathematical modeling tool applicable to...A discrete event system is a dynamical system whose state evolves in time by the occurrence of events at possibly irregular time intervals. Timed Petri nets are a graphical and mathematical modeling tool applicable to discrete event systems in order to represent its states evolution where the timing at which the state changes is taken into consideration. One of the most important performance issues to be considered in a discrete event system is its stability. Lyapunov theory provides the required tools needed to aboard the stability and stabilization problems for discrete event systems modeled with timed Petri nets whose mathematical model is given in terms of difference equations. By proving stability one guarantees a bound on the discrete event systems state dynamics. When the system is unstable, a sufficient condition to stabilize the system is given. It is shown that it is possible to restrict the discrete event systems state space in such a way that boundedness is achieved. However, the restriction is not numerically precisely known. This inconvenience is overcome by considering a specific recurrence equation, in the max-plus algebra, which is assigned to the timed Petri net graphical model.展开更多
The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenar...The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenarios, however, needs to be quantitatively evaluated. In this study, four compensation methods （i.e., the polynomial extrapolation, the linear acceleration extrapolation, the inverse compensation and the adaptive inverse compensation） are selected and compared experimentally using a frequency evaluation index （FEI） method. The effectiveness of the FEI method is first verified through comparison with the discrete transfer fimction approach for compensation methods assuming constant delay. Incomparable advantage is further demonstrated for the FEI method when applied to adaptive compensation methods, where the discrete transfer function approach is difficult to implement. Both numerical simulation and laboratory tests with predefined displacements are conducted using sinusoidal signals and random signals as inputs. Findings from numerical simulation and experimental results demonstrate that the FEI method is an efficient and effective approach to compare the performance of different compensation methods, especially for those requiring adaptation of compensation parameters.展开更多
We present the new predictor-corrector methods for systems of nonlinear differential equations, based on the method of exponential time differencing. We compare the present schemes with the explicit multistep exponent...We present the new predictor-corrector methods for systems of nonlinear differential equations, based on the method of exponential time differencing. We compare the present schemes with the explicit multistep exponential time differencing and Adams–Bashforth–Moulton method. The numerical results show that the schemes are more accurate and more efficient than Adams predictor-corrector method. The exponential time differencing method has been developed and perfected by the present studies.展开更多
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.展开更多
This paper discusses the modeling method of time series with neural network. In order to improve the adaptability of direct multi-step prediction models, this paper proposes a method of combining the temporal differen...This paper discusses the modeling method of time series with neural network. In order to improve the adaptability of direct multi-step prediction models, this paper proposes a method of combining the temporal differences methods with back-propagation algorithm for updating the parameters continuously on the basis of recent data. This method can make the neural network model fit the recent characteristic of the time series as close as possible, therefore improves the prediction accuracy. We built models and made predictions for the sunspot series. The prediction results of adaptive modeling method are better than that of non-adaptive modeling methods.展开更多
Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Th...Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Though, this paper shows that recent advance methods can be more favored. In this work, we have incorporated, throughout numerical comparison experiments, spectral methods, for the space discretization, in conjunction with second and fourth-order time integrating methods for approximating the solution of the reaction-diffusion differential equations. The results have revealed that these methods have advantages over the conventional methods, some of which to mention are: the ease of implementation, accuracy and CPU time.展开更多
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.展开更多
High precision position control and high speed control of the robot manipulators are fundamental and important control problems. The effectiveness of finite-time P-PI control was confirmed by end-effector position con...High precision position control and high speed control of the robot manipulators are fundamental and important control problems. The effectiveness of finite-time P-PI control was confirmed by end-effector position control of robot manipulators. However, parameter tuning method has not been proposed to finite-time P-PI control. In this paper, we propose a settling time design method and a parameter tuning method for the finite-time P-PI control. The effectiveness of the proposed parameter tuning method is confirmed by experiments of end-effcctor position control of a robot manipulator.展开更多
The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spit...The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spite of methodological progress,the issues of efficiency and ease of implementation of the solution of the associated generalized Riemann problem(GRP)remain the centre of attention in the ADER approach.In the original formulation of ADER schemes,the proposed solution procedure for the GRP was based on(i)Taylor series expansion of the solution in time right at the element interface,(ii)subsequent application of the Cauchy-Kowalewskaya procedure to convert time derivatives to functionals of space derivatives,and(iii)solution of classical Riemann problems for high-order spatial derivatives to complete the Taylor series expansion.For realistic problems the Cauchy-Kowalewskaya procedure requires the use of symbolic manipulators and being rather cumbersome its replacement or simplification is highly desirable.In this paper we propose a new class of solvers for the GRP that avoid the Cauchy-Kowalewskaya procedure and result in simpler ADER schemes.This is achieved by exploiting the history of the numerical solution that makes it possible to devise a time-reconstruction procedure at the element interface.Still relying on a time Taylor series expansion of the solution at the interface,the time derivatives are then easily calculated from the time-reconstruction polynomial.The resulting schemes are called ADER-TR.A thorough study of the linear stability properties of the linear version of the schemes is carried out using the von Neumann method,thus deducing linear stability regions.Also,via careful numerical experiments,we deduce stability regions for the corresponding non-linear schemes.Numerical examples using the present simplified schemes of fifth and seventh order of accuracy in space and time show that these compare favourably with conventional ADER methods.This paper is restricted to the one-dimensional scalar case with source term,but preliminary results for the one-dimensional Euler equations indicate that the time-reconstruction approach offers significant advantages not only in terms of ease of implementation but also in terms of efficiency for the high-order range schemes.展开更多
Based on the constructivism, using the College English（new edition） CD-ROMs for classroom teaching as main media, this paper explores the feasibility of the classroom teaching mode of appropriate combination of trad...Based on the constructivism, using the College English（new edition） CD-ROMs for classroom teaching as main media, this paper explores the feasibility of the classroom teaching mode of appropriate combination of traditional methods with the multimedia overhead projector. The new teaching mode proves to be successful in that it has greatly improved the students＇ ability to use English in an all-round way and has obvious advantages over the only traditional method.展开更多
The Burst Time Plan(BTP) generation is the key for resource allocation in Broadband Satellite Multimedia(BSM) system.The main purpose of this paper is to minimize the system response time to users' request caused ...The Burst Time Plan(BTP) generation is the key for resource allocation in Broadband Satellite Multimedia(BSM) system.The main purpose of this paper is to minimize the system response time to users' request caused by BTP generation as well as maintain the Quality of Service(QoS) and improve the channel utilization efficiency.Traditionally the BTP is generated periodically in order to simplify the implementation of the resource allocation algorithm.Based on the analysis we find that Periodical BTP Generation(P-BTPG) method cannot guarantee the delay performance,channel utilization efficiency and QoS simultaneously,especially when the capacity requests arrived randomly.The Optimized BTP Generation(O-BTPG) method is given based on the optimal scheduling period and scheduling latency without considering the signaling overhead.Finally,a novel Asynchronous BTP Generation(A-BTPG) method is proposed which is invoked according to users' requests.A BSM system application scenario is simulated.Simulation results show that A-BTPG is a trade-off between the performance and signaling overhead which can improve the system performance insensitive to the traffic pattern.This method can be used in the ATM onboard switching satellite system and further more can be expended to Digital Video Broadcasting-Return Channel Satellite(DVB-RCS) system or IP onboard routing BSM system in the future.展开更多
Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied....Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied. Furthermore, the estimation of the region-time-length (RTL) parameter has been proposed to detect seismic quiescence before the occurrence of a large earthquake. In addition, the time-to-failure method has been used to estimate the time occurrence of large earthquakes. Hence, in this study, to gain deeper insight into seismic activity in the southern Zagros region, we utilized the RTL algorithm to identify the quiescence and activation phases leading to the Fin doublet earthquakes. Temporal variations in the RTL parameter showed two significant anomalies. One corresponded to the occurrence time of the first earthquake (2017-12-12);the other anomaly was associated with the occurrence time of the second event (2021-11-14). Based on a negative value of the RTL parameter observed in the vicinity of the Fin epicenters (2021), seismic quiescence (a decrease in seismicity compared to the preceding background rate) was identified. The spatial distribution of the RTL prognostic parameters confirms the appearance of seismic quiescence surrounding the epicenter of the Fin doublet earthquakes (2021). The time-to-failure method was designed using precursory events that describe the acceleration of the seismic energy release before the mainshock. Using the time-to-failure method for the earthquake catalog, it was possible to estimate both the magnitude and time of failure of the Fin doublet. Hence, the time-tofailure technique can be a useful supplementary method to the RTL algorithm for determining the characteristics of impending earthquakes.展开更多
Time series forecasting has become an important aspect of data analysis and has many real-world applications.However,undesirable missing values are often encountered,which may adversely affect many forecasting tasks.I...Time series forecasting has become an important aspect of data analysis and has many real-world applications.However,undesirable missing values are often encountered,which may adversely affect many forecasting tasks.In this study,we evaluate and compare the effects of imputationmethods for estimating missing values in a time series.Our approach does not include a simulation to generate pseudo-missing data,but instead perform imputation on actual missing data and measure the performance of the forecasting model created therefrom.In an experiment,therefore,several time series forecasting models are trained using different training datasets prepared using each imputation method.Subsequently,the performance of the imputation methods is evaluated by comparing the accuracy of the forecasting models.The results obtained from a total of four experimental cases show that the k-nearest neighbor technique is the most effective in reconstructing missing data and contributes positively to time series forecasting compared with other imputation methods.展开更多
The Routh and Whittaker methods of reduction for Lagrange system on time scales with nabla derivatives are studied.The equations of motion for Lagrange system on time scales are established, and their cyclic integrals...The Routh and Whittaker methods of reduction for Lagrange system on time scales with nabla derivatives are studied.The equations of motion for Lagrange system on time scales are established, and their cyclic integrals and generalized energy integrals are given. The Routh functions and Whittaker functions of Lagrange system are constructed, and the order of differential equations of motion for the system are reduced by using the cyclic integrals or the generalized energy integrals with nabla derivatives. The results show that the reduced Routh equations and Whittaker equations hold the form of Lagrnage equations with nabla derivatives. Finally, two examples are given to illustrate the application of the results.展开更多
Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that ...Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that the capability to capture dynamic loading was not explored. It was found that there were different quadrature equations to predict the next step displacement increment. A modified quadrature equation of this method was derived so that the equation to determine the next step displacement was numerically equivalent to the equation used in the constant AAM. It was verified that the original form of this method, in general, had a better capability to capture dynamic loadings than the constant AAM. This excellent property, in addition to computational efficiency, will help to make this method competitive with general secondorder accurate integration methods.展开更多
文摘This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.
基金the National Natural Science Foundation of China(22078030)the National Key Research and Development Project(2019YFC1905802,2022YFB3504305)+1 种基金the Joint Funds of the National Natural Science Foundation of China(U1802255,CSTB2022NSCQ-LZX0014)the Key Project of Independent Research Project of State Key Laboratory of Coal Mine Disaster Dynamics and Control(2011DA105287-zd201902).
文摘Leveraging big data signal processing offers a pathway to the development of artificial intelligencedriven equipment.The analysis of fluid flow signals and the characterization of fluid flow behavior are of critical in two-phase flow studies.Significant research efforts have focused on discerning flow regimes using various signal analysis methods.In this review,recent advances in time series signals analysis algorithms for stirred tank reactors have been summarized,and the detailed methodologies are categorized into the frequency domain methods,time-frequency domain methods,and state space methods.The strengths,limitations,and notable findings of each algorithm are highlighted.Additionally,the interrelationships between these methodologies have also been discussed,as well as the present progress achieved in various applications.Future research directions and challenges are also predicted to provide an overview of current research trends in data mining of time series for analyzing flow regimes and chaotic signals.This review offers a comprehensive summary for extracting and characterizing fluid flow behavior and serves as a theoretical reference for optimizing the characterization of chaotic signals in future research endeavors.
基金Supported by the National Natural Science Foundations of China （No. 40474001, No. 40274002, No. 40604003）.
文摘Three functional models, polynomial, spectral analysis, and modified AR model, are studied and compared in fitting and predicting clock deviation based on the data sequence derived from two-way satellite time and frequency transfer. A robust equivalent weight is applied, which controls the significant influence of outlying observations. Some conclusions show that the prediction precision of robust estimation is better than that of LS. The prediction precision calculated from smoothed observations is higher than that calculated from sampling observations. As a count of the obvious period variations in the clock deviation sequence, the predicted values of polynomial model are implausible. The prediction precision of spectral analysis model is very low, but the principal periods can be determined. The prediction RMS of 6-hour extrapolation interval is Ins or so, when modified AR model is used.
基金Open Access funding provided by Universita degli Studi di Verona.
文摘We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not possible.As shown in Boscarino et al.(J.Sci.Comput.68:975-1001,2016)for Runge-Kutta methods,these semi-implicit techniques give a great flexibility,and allow,in many cases,the construction of simple linearly implicit schemes with no need of iterative solvers.In this work,we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype lineal'advection-diffusion equation and in the setting of strong stability preserving(SSP)methods.Our findings are demonstrated on several examples,including nonlinear reaction-diffusion and convection-diffusion problems.
基金Project supported by National Natural Science Foundation of China and China State Key project for Basic Researchcs.
文摘In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.
文摘A discrete event system is a dynamical system whose state evolves in time by the occurrence of events at possibly irregular time intervals. Timed Petri nets are a graphical and mathematical modeling tool applicable to discrete event systems in order to represent its states evolution where the timing at which the state changes is taken into consideration. One of the most important performance issues to be considered in a discrete event system is its stability. Lyapunov theory provides the required tools needed to aboard the stability and stabilization problems for discrete event systems modeled with timed Petri nets whose mathematical model is given in terms of difference equations. By proving stability one guarantees a bound on the discrete event systems state dynamics. When the system is unstable, a sufficient condition to stabilize the system is given. It is shown that it is possible to restrict the discrete event systems state space in such a way that boundedness is achieved. However, the restriction is not numerically precisely known. This inconvenience is overcome by considering a specific recurrence equation, in the max-plus algebra, which is assigned to the timed Petri net graphical model.
基金National Natural Science Foundation of China under Grant No.51378107the Fundamental Research Funds for the Central Universities and Priority Academic Program Development of Jiangsu Higher Education Institutions under Grant No.KYLX-0158the National Natural Science Foundation under Grant No.CMMI-1227962
文摘The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenarios, however, needs to be quantitatively evaluated. In this study, four compensation methods （i.e., the polynomial extrapolation, the linear acceleration extrapolation, the inverse compensation and the adaptive inverse compensation） are selected and compared experimentally using a frequency evaluation index （FEI） method. The effectiveness of the FEI method is first verified through comparison with the discrete transfer fimction approach for compensation methods assuming constant delay. Incomparable advantage is further demonstrated for the FEI method when applied to adaptive compensation methods, where the discrete transfer function approach is difficult to implement. Both numerical simulation and laboratory tests with predefined displacements are conducted using sinusoidal signals and random signals as inputs. Findings from numerical simulation and experimental results demonstrate that the FEI method is an efficient and effective approach to compare the performance of different compensation methods, especially for those requiring adaptation of compensation parameters.
基金The project supported by National Natural Science Foundation of China under Grant No.19902002
文摘We present the new predictor-corrector methods for systems of nonlinear differential equations, based on the method of exponential time differencing. We compare the present schemes with the explicit multistep exponential time differencing and Adams–Bashforth–Moulton method. The numerical results show that the schemes are more accurate and more efficient than Adams predictor-corrector method. The exponential time differencing method has been developed and perfected by the present studies.
基金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.
文摘This paper discusses the modeling method of time series with neural network. In order to improve the adaptability of direct multi-step prediction models, this paper proposes a method of combining the temporal differences methods with back-propagation algorithm for updating the parameters continuously on the basis of recent data. This method can make the neural network model fit the recent characteristic of the time series as close as possible, therefore improves the prediction accuracy. We built models and made predictions for the sunspot series. The prediction results of adaptive modeling method are better than that of non-adaptive modeling methods.
文摘Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Though, this paper shows that recent advance methods can be more favored. In this work, we have incorporated, throughout numerical comparison experiments, spectral methods, for the space discretization, in conjunction with second and fourth-order time integrating methods for approximating the solution of the reaction-diffusion differential equations. The results have revealed that these methods have advantages over the conventional methods, some of which to mention are: the ease of implementation, accuracy and CPU time.
文摘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.
文摘High precision position control and high speed control of the robot manipulators are fundamental and important control problems. The effectiveness of finite-time P-PI control was confirmed by end-effector position control of robot manipulators. However, parameter tuning method has not been proposed to finite-time P-PI control. In this paper, we propose a settling time design method and a parameter tuning method for the finite-time P-PI control. The effectiveness of the proposed parameter tuning method is confirmed by experiments of end-effcctor position control of a robot manipulator.
基金G.I.Montecinos thanks the National Chilean Fund for Scientific and Technological Development,FONDECYT(Fondo Nacional de Desarrollo Científico y Tecnológico),in the frame of the project for Initiation in Research 11180926
文摘The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spite of methodological progress,the issues of efficiency and ease of implementation of the solution of the associated generalized Riemann problem(GRP)remain the centre of attention in the ADER approach.In the original formulation of ADER schemes,the proposed solution procedure for the GRP was based on(i)Taylor series expansion of the solution in time right at the element interface,(ii)subsequent application of the Cauchy-Kowalewskaya procedure to convert time derivatives to functionals of space derivatives,and(iii)solution of classical Riemann problems for high-order spatial derivatives to complete the Taylor series expansion.For realistic problems the Cauchy-Kowalewskaya procedure requires the use of symbolic manipulators and being rather cumbersome its replacement or simplification is highly desirable.In this paper we propose a new class of solvers for the GRP that avoid the Cauchy-Kowalewskaya procedure and result in simpler ADER schemes.This is achieved by exploiting the history of the numerical solution that makes it possible to devise a time-reconstruction procedure at the element interface.Still relying on a time Taylor series expansion of the solution at the interface,the time derivatives are then easily calculated from the time-reconstruction polynomial.The resulting schemes are called ADER-TR.A thorough study of the linear stability properties of the linear version of the schemes is carried out using the von Neumann method,thus deducing linear stability regions.Also,via careful numerical experiments,we deduce stability regions for the corresponding non-linear schemes.Numerical examples using the present simplified schemes of fifth and seventh order of accuracy in space and time show that these compare favourably with conventional ADER methods.This paper is restricted to the one-dimensional scalar case with source term,but preliminary results for the one-dimensional Euler equations indicate that the time-reconstruction approach offers significant advantages not only in terms of ease of implementation but also in terms of efficiency for the high-order range schemes.
文摘Based on the constructivism, using the College English（new edition） CD-ROMs for classroom teaching as main media, this paper explores the feasibility of the classroom teaching mode of appropriate combination of traditional methods with the multimedia overhead projector. The new teaching mode proves to be successful in that it has greatly improved the students＇ ability to use English in an all-round way and has obvious advantages over the only traditional method.
基金Supported by the National Natural Science Foundation ofChina (No. 60972061,60972062,and 61032004)the Na-tional High Technology Research and Development Program of China ("863" Program) (No. 2008AA12A204,2008AA12Z307)+1 种基金Natural Science Foundation of Jiangsu Province (BK2009060)the"Triple Three" High Level Talent Development Plan of Jiangsu Province
文摘The Burst Time Plan(BTP) generation is the key for resource allocation in Broadband Satellite Multimedia(BSM) system.The main purpose of this paper is to minimize the system response time to users' request caused by BTP generation as well as maintain the Quality of Service(QoS) and improve the channel utilization efficiency.Traditionally the BTP is generated periodically in order to simplify the implementation of the resource allocation algorithm.Based on the analysis we find that Periodical BTP Generation(P-BTPG) method cannot guarantee the delay performance,channel utilization efficiency and QoS simultaneously,especially when the capacity requests arrived randomly.The Optimized BTP Generation(O-BTPG) method is given based on the optimal scheduling period and scheduling latency without considering the signaling overhead.Finally,a novel Asynchronous BTP Generation(A-BTPG) method is proposed which is invoked according to users' requests.A BSM system application scenario is simulated.Simulation results show that A-BTPG is a trade-off between the performance and signaling overhead which can improve the system performance insensitive to the traffic pattern.This method can be used in the ATM onboard switching satellite system and further more can be expended to Digital Video Broadcasting-Return Channel Satellite(DVB-RCS) system or IP onboard routing BSM system in the future.
文摘Knowledge regarding earthquake hazards and seismicity is crucial for crisis management, and the occurrence of foreshocks, seismic activity patterns, and spatiotemporal variations in seismic activity have been studied. Furthermore, the estimation of the region-time-length (RTL) parameter has been proposed to detect seismic quiescence before the occurrence of a large earthquake. In addition, the time-to-failure method has been used to estimate the time occurrence of large earthquakes. Hence, in this study, to gain deeper insight into seismic activity in the southern Zagros region, we utilized the RTL algorithm to identify the quiescence and activation phases leading to the Fin doublet earthquakes. Temporal variations in the RTL parameter showed two significant anomalies. One corresponded to the occurrence time of the first earthquake (2017-12-12);the other anomaly was associated with the occurrence time of the second event (2021-11-14). Based on a negative value of the RTL parameter observed in the vicinity of the Fin epicenters (2021), seismic quiescence (a decrease in seismicity compared to the preceding background rate) was identified. The spatial distribution of the RTL prognostic parameters confirms the appearance of seismic quiescence surrounding the epicenter of the Fin doublet earthquakes (2021). The time-to-failure method was designed using precursory events that describe the acceleration of the seismic energy release before the mainshock. Using the time-to-failure method for the earthquake catalog, it was possible to estimate both the magnitude and time of failure of the Fin doublet. Hence, the time-tofailure technique can be a useful supplementary method to the RTL algorithm for determining the characteristics of impending earthquakes.
基金This research was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF)funded by the Ministry of Education(Grant Number 2020R1A6A1A03040583).
文摘Time series forecasting has become an important aspect of data analysis and has many real-world applications.However,undesirable missing values are often encountered,which may adversely affect many forecasting tasks.In this study,we evaluate and compare the effects of imputationmethods for estimating missing values in a time series.Our approach does not include a simulation to generate pseudo-missing data,but instead perform imputation on actual missing data and measure the performance of the forecasting model created therefrom.In an experiment,therefore,several time series forecasting models are trained using different training datasets prepared using each imputation method.Subsequently,the performance of the imputation methods is evaluated by comparing the accuracy of the forecasting models.The results obtained from a total of four experimental cases show that the k-nearest neighbor technique is the most effective in reconstructing missing data and contributes positively to time series forecasting compared with other imputation methods.
基金supported by the National Natural Science Foundation of China(Grant Nos.11572212 and 11272227)the Innovation Program for Graduate Student of Jiangsu Province,China(Grant No.KYLX16-0414)
文摘The Routh and Whittaker methods of reduction for Lagrange system on time scales with nabla derivatives are studied.The equations of motion for Lagrange system on time scales are established, and their cyclic integrals and generalized energy integrals are given. The Routh functions and Whittaker functions of Lagrange system are constructed, and the order of differential equations of motion for the system are reduced by using the cyclic integrals or the generalized energy integrals with nabla derivatives. The results show that the reduced Routh equations and Whittaker equations hold the form of Lagrnage equations with nabla derivatives. Finally, two examples are given to illustrate the application of the results.
基金Science Council (NSC),Chinese Taipei Under Grant No.NSC-96-2221-E-027-030
文摘Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that the capability to capture dynamic loading was not explored. It was found that there were different quadrature equations to predict the next step displacement increment. A modified quadrature equation of this method was derived so that the equation to determine the next step displacement was numerically equivalent to the equation used in the constant AAM. It was verified that the original form of this method, in general, had a better capability to capture dynamic loadings than the constant AAM. This excellent property, in addition to computational efficiency, will help to make this method competitive with general secondorder accurate integration methods.
