This paper deals with the stability of linear multistep methods for multidimensional differential systems with distributed delays. The delay-dependent stability of linear multistep methods with compound quadrature rul...This paper deals with the stability of linear multistep methods for multidimensional differential systems with distributed delays. The delay-dependent stability of linear multistep methods with compound quadrature rules is studied. Several new sufficient criteria of delay-dependent stability are obtained by means of the argument principle. An algorithm is provided to check delay-dependent stability. An example that illustrates the effectiveness of the derived theoretical results is given.展开更多
The stability analysis of linear multistep methods for the numerical solutions of the systems of generalized neutral delay differential equations is discussed. The stability behaviour of linear multistep methods was a...The stability analysis of linear multistep methods for the numerical solutions of the systems of generalized neutral delay differential equations is discussed. The stability behaviour of linear multistep methods was analysed for the solution of the generalized system of linear neutral test equations, After the establishment of a sufficient condition for asymptotic stability of the solutions of the generalized system, it is shown that a linear multistep method is NGP(G)-stable if and only if it is A-stable.展开更多
This paper deals with the stability analysis of the linear multistep (LM) methods in the numerical solution of delay differential equations. Here we provide a qualitative stability estimates, pertiment to the classica...This paper deals with the stability analysis of the linear multistep (LM) methods in the numerical solution of delay differential equations. Here we provide a qualitative stability estimates, pertiment to the classical scalar test problem of the form y′(t)=λy(t)+μy(t-τ) with τ>0 and λ,μ are complex, by using (vartiant to) the resolvent condition of Kreiss. We prove that for A stable LM methods the upper bound for the norm of the n th power of square matrix grows linearly with the order of the matrix.展开更多
Studies the numerical stability region of linear multistep(LM) methods applied to linear test equation of the form y′(t)=ay(t)+by(t-1), t>0, y(t)=g(t)-1≤t≤0, a,b∈R, proves through delay dependent stability ...Studies the numerical stability region of linear multistep(LM) methods applied to linear test equation of the form y′(t)=ay(t)+by(t-1), t>0, y(t)=g(t)-1≤t≤0, a,b∈R, proves through delay dependent stability analysis that the intersection of stability regions of the equation and the method is not empty, in addition to approaches to the boundary of the delay differential equation(DDEs) in the limiting case of step size boundary of the stability region of linear multistep methods.展开更多
The numerical stability of linear neutral Volterra delay integral differential equations is dealt with.Under a sufficient condition such that this system with a lagging argument is asymptotically stable,it is proved t...The numerical stability of linear neutral Volterra delay integral differential equations is dealt with.Under a sufficient condition such that this system with a lagging argument is asymptotically stable,it is proved that every A-stable linear multistep method preserves the delay-independent stabil-ity of its exact solutions.Finally,some numerical experiments are given to demonstrate the main conclu-sion.展开更多
The stability analysis of linear multistep (LM) methods is carried out under Kreiss resolvent condition when they are applied to neutral delay differential equations of the form y′(t)=ay(t)+by(t-τ)+ cy′(t- τ) y(t)...The stability analysis of linear multistep (LM) methods is carried out under Kreiss resolvent condition when they are applied to neutral delay differential equations of the form y′(t)=ay(t)+by(t-τ)+ cy′(t- τ) y(t)=g(t) -τ≤t≤0 with τ>0 and a, b and c∈, and it is proved that the ‖B n‖ is suitably bounded, where B is the companion matrix.展开更多
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.展开更多
Forced and damped oscillators appear in the mathematical modelling of many problems in pure and applied sciences such as physics, engineering and celestial mechanics among others. Although the accuracy of the T-functi...Forced and damped oscillators appear in the mathematical modelling of many problems in pure and applied sciences such as physics, engineering and celestial mechanics among others. Although the accuracy of the T-functions series method is high, the calculus of their coefficients needs specific recurrences in each case. To avoid this inconvenience, the T-functions series method is transformed into a multistep method whose coefficients are calculated using recurrence procedures. These methods are convergent and have the same properties to the T-functions series method. Numerical examples already used by other authors are presented, such as a stiff problem, a Duffing oscillator and an equatorial satellite problem when the perturbation comes from zonal harmonics J2.展开更多
Stability and global error bounds are studied for a class of stepsize-dependent linear multistep methods for nonlinear evolution equations governed by ω-dissipative vector fields in Banach space.To break through the ...Stability and global error bounds are studied for a class of stepsize-dependent linear multistep methods for nonlinear evolution equations governed by ω-dissipative vector fields in Banach space.To break through the order barrier p≤1 of unconditionally contractive linear multistep methods for dissipative systems,strongly dissipative systems are introduced.By employing the error growth function of the methods,new contractivity and convergence results of stepsize-dependent linear multistep methods on infinite integration intervals are provided for strictly dissipative systems(ω<0)and strongly dissipative systems.Some applications of the main results to several linear multistep methods,including the trapezoidal rule,are supplied.The theoretical results are also illustrated by a set of numerical experiments.展开更多
This paper deals with a delay-dependent treatment of linear multistep methods for neutral delay differential equations y'(t) = ay(t) + by(t - τ) + cy'(t - τ), t > 0, y(t) = g(t), -τ≤ t ≤ 0, a,b andc ∈...This paper deals with a delay-dependent treatment of linear multistep methods for neutral delay differential equations y'(t) = ay(t) + by(t - τ) + cy'(t - τ), t > 0, y(t) = g(t), -τ≤ t ≤ 0, a,b andc ∈ R. The necessary condition for linear multistep methods to be Nτ(0)-stable is given. It is shown that the trapezoidal rule is Nτ(0)-compatible. Figures of stability region for some linear multistep methods are depicted.展开更多
Presents information on a study which focused on the numerical solution of initial value problems for systems of neutral differential equations. Adaptations of linear multistep methods; Linear stability of linear mult...Presents information on a study which focused on the numerical solution of initial value problems for systems of neutral differential equations. Adaptations of linear multistep methods; Linear stability of linear multistep method; Presentation of numerical equations.展开更多
Presents information on a study which dealt with the error behavior and the stability analysis of a class of linear multistep methods with the Lagrangian interpolation as applied to the nonlinear delay differential eq...Presents information on a study which dealt with the error behavior and the stability analysis of a class of linear multistep methods with the Lagrangian interpolation as applied to the nonlinear delay differential equations. Methods and the basic lemmas; Analysis of convergence and stability.展开更多
Implicit-explicit （IMEX） linear multistep methods are popular techniques for solving partial differential equations （PDEs） with terms of different types. While fixed timestep versions of such schemes have been dev...Implicit-explicit （IMEX） linear multistep methods are popular techniques for solving partial differential equations （PDEs） with terms of different types. While fixed timestep versions of such schemes have been developed and studied, implicit-explicit schemes also naturally arise in general situations where the temporal smoothness of the solution changes. In this paper we consider easily implementable variable step-size implicit-explicit （VSIMEX） linear multistep methods for time-dependent PDEs. Families of order-p, pstep VSIMEX schemes are constructed and analyzed, where p ranges from 1 to 4. The corresponding schemes are simple to implement and have the property that they reduce to the classical IMEX schemes whenever constant time step-sizes are imposed. The methods are validated on the Burgers＇ equation. These results demonstrate that by varying the time step-size, VSIMEX methods can outperform their fixed time step counterparts while still maintaining good numerical behavior.展开更多
For the numerical treatment of Hamiltonian differential equations, symplectic integrators are the most suitable choice, and methods that are conjugate to a symplectic integrator share the same good long-time behavior....For the numerical treatment of Hamiltonian differential equations, symplectic integrators are the most suitable choice, and methods that are conjugate to a symplectic integrator share the same good long-time behavior. This note characterizes linear multistep methods whose underlying one-step method is conjugate to a symplectic integrator. The bounded- hess of parasitic solution components is not addressed.展开更多
The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equ...The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equation is discretized by a mixed finite element method. The electron and hole density equations are treated by implicit-explicit multistep finite element methods. The schemes are very efficient. The optimal order error estimates both in time and space are derived.展开更多
This paper considers the asymptotic stability of linear multistep(LM)methods for neutral systems with distributed delays.In particular,several sufficient conditions for delay-dependent stability of numerical solutions...This paper considers the asymptotic stability of linear multistep(LM)methods for neutral systems with distributed delays.In particular,several sufficient conditions for delay-dependent stability of numerical solutions are obtained based on the argument principle.Compound quadrature formulae are used to compute the integrals.An algorithm is proposed to examine the delay-dependent stability of numerical solutions.Several numerical examples are performed to verify the theoretical results.展开更多
A series of sufficient and necessary conditions for the algebraic stability of multistepRunge-Kutta methods is obtained, most of which can be regarded as extension of the relevant results available for Runge-Kutta met...A series of sufficient and necessary conditions for the algebraic stability of multistepRunge-Kutta methods is obtained, most of which can be regarded as extension of the relevant results available for Runge-Kutta methods, especially, for Radau Ⅰ A, Radau Ⅱ A and Gaussian Runge-Kutta methods.展开更多
This paper is concerned with the numerical dissipativity of multistep Runge-Kutta methods for nonlinear neutral delay-integro-differential equations.We investigate the dissipativity properties of-algebraically stable ...This paper is concerned with the numerical dissipativity of multistep Runge-Kutta methods for nonlinear neutral delay-integro-differential equations.We investigate the dissipativity properties of-algebraically stable multistep Runge-Kutta methods with constrained grid.The finite-dimensional and infinite-dimensional dissipativity results of-algebraically stable multistep Runge-Kutta methods are obtained.展开更多
The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the soluti...The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the solution,characterized by||u_(n)||L^(2)(Ω)=O(t^(−αs)) as t→∞,is determined by the minimum order α_(s) of the time-fractional derivatives.Building on this foundational result,this article pursues two primary objectives.First,we introduce a strongly A-stable fractional linear multistep method and derive the numerical stability region for the governing equation.Second,we rigorously prove the long-term decay rate of the numerical solution through a detailed singularity analysis of its generating function.Notably,the numerical decay rate||u_(n)||L^(2)(Ω)=O(t_(n)^(−α_(s)) as t_(n)→∞aligns precisely with the continuous case.Theoretical findings are further validated through comprehensive numerical simulations,underscoring the robustness of our proposed method.展开更多
基金Project supported by the National Natural Science Foundation of China(No.11471217)
文摘This paper deals with the stability of linear multistep methods for multidimensional differential systems with distributed delays. The delay-dependent stability of linear multistep methods with compound quadrature rules is studied. Several new sufficient criteria of delay-dependent stability are obtained by means of the argument principle. An algorithm is provided to check delay-dependent stability. An example that illustrates the effectiveness of the derived theoretical results is given.
文摘The stability analysis of linear multistep methods for the numerical solutions of the systems of generalized neutral delay differential equations is discussed. The stability behaviour of linear multistep methods was analysed for the solution of the generalized system of linear neutral test equations, After the establishment of a sufficient condition for asymptotic stability of the solutions of the generalized system, it is shown that a linear multistep method is NGP(G)-stable if and only if it is A-stable.
文摘This paper deals with the stability analysis of the linear multistep (LM) methods in the numerical solution of delay differential equations. Here we provide a qualitative stability estimates, pertiment to the classical scalar test problem of the form y′(t)=λy(t)+μy(t-τ) with τ>0 and λ,μ are complex, by using (vartiant to) the resolvent condition of Kreiss. We prove that for A stable LM methods the upper bound for the norm of the n th power of square matrix grows linearly with the order of the matrix.
文摘Studies the numerical stability region of linear multistep(LM) methods applied to linear test equation of the form y′(t)=ay(t)+by(t-1), t>0, y(t)=g(t)-1≤t≤0, a,b∈R, proves through delay dependent stability analysis that the intersection of stability regions of the equation and the method is not empty, in addition to approaches to the boundary of the delay differential equation(DDEs) in the limiting case of step size boundary of the stability region of linear multistep methods.
基金Supported by the Natural Science Foundation of Heilongjiang Province(A200602)the Project of Science Research Foundation(HITC200710)the Project of Development Program for Outstanding Young Teachers in Harbin Institute of Technology(HITQNJS.2006.053)
文摘The numerical stability of linear neutral Volterra delay integral differential equations is dealt with.Under a sufficient condition such that this system with a lagging argument is asymptotically stable,it is proved that every A-stable linear multistep method preserves the delay-independent stabil-ity of its exact solutions.Finally,some numerical experiments are given to demonstrate the main conclu-sion.
文摘The stability analysis of linear multistep (LM) methods is carried out under Kreiss resolvent condition when they are applied to neutral delay differential equations of the form y′(t)=ay(t)+by(t-τ)+ cy′(t- τ) y(t)=g(t) -τ≤t≤0 with τ>0 and a, b and c∈, and it is proved that the ‖B n‖ is suitably bounded, where B is the companion matrix.
基金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.
文摘Forced and damped oscillators appear in the mathematical modelling of many problems in pure and applied sciences such as physics, engineering and celestial mechanics among others. Although the accuracy of the T-functions series method is high, the calculus of their coefficients needs specific recurrences in each case. To avoid this inconvenience, the T-functions series method is transformed into a multistep method whose coefficients are calculated using recurrence procedures. These methods are convergent and have the same properties to the T-functions series method. Numerical examples already used by other authors are presented, such as a stiff problem, a Duffing oscillator and an equatorial satellite problem when the perturbation comes from zonal harmonics J2.
基金supported by the Natural Science Foundation of China(Grant Nos.12271367,11771060)by the Science and Technology Innovation Plan of Shanghai,China(Grant No.20JC1414200)sponsored by the Natural Science Foundation of Shanghai,China(Grant No.20ZR1441200).
文摘Stability and global error bounds are studied for a class of stepsize-dependent linear multistep methods for nonlinear evolution equations governed by ω-dissipative vector fields in Banach space.To break through the order barrier p≤1 of unconditionally contractive linear multistep methods for dissipative systems,strongly dissipative systems are introduced.By employing the error growth function of the methods,new contractivity and convergence results of stepsize-dependent linear multistep methods on infinite integration intervals are provided for strictly dissipative systems(ω<0)and strongly dissipative systems.Some applications of the main results to several linear multistep methods,including the trapezoidal rule,are supplied.The theoretical results are also illustrated by a set of numerical experiments.
基金This work was supported by the NSF of P.R.of China(10271036)
文摘This paper deals with a delay-dependent treatment of linear multistep methods for neutral delay differential equations y'(t) = ay(t) + by(t - τ) + cy'(t - τ), t > 0, y(t) = g(t), -τ≤ t ≤ 0, a,b andc ∈ R. The necessary condition for linear multistep methods to be Nτ(0)-stable is given. It is shown that the trapezoidal rule is Nτ(0)-compatible. Figures of stability region for some linear multistep methods are depicted.
文摘Presents information on a study which focused on the numerical solution of initial value problems for systems of neutral differential equations. Adaptations of linear multistep methods; Linear stability of linear multistep method; Presentation of numerical equations.
基金National Natural Science Foundation of China!No.69974018 Postdoctoral Science Foundation of China.
文摘Presents information on a study which dealt with the error behavior and the stability analysis of a class of linear multistep methods with the Lagrangian interpolation as applied to the nonlinear delay differential equations. Methods and the basic lemmas; Analysis of convergence and stability.
基金supported by an NSERC Canada Postgraduate Scholarshipsupported by a grant from NSERC Canada
文摘Implicit-explicit （IMEX） linear multistep methods are popular techniques for solving partial differential equations （PDEs） with terms of different types. While fixed timestep versions of such schemes have been developed and studied, implicit-explicit schemes also naturally arise in general situations where the temporal smoothness of the solution changes. In this paper we consider easily implementable variable step-size implicit-explicit （VSIMEX） linear multistep methods for time-dependent PDEs. Families of order-p, pstep VSIMEX schemes are constructed and analyzed, where p ranges from 1 to 4. The corresponding schemes are simple to implement and have the property that they reduce to the classical IMEX schemes whenever constant time step-sizes are imposed. The methods are validated on the Burgers＇ equation. These results demonstrate that by varying the time step-size, VSIMEX methods can outperform their fixed time step counterparts while still maintaining good numerical behavior.
基金the Swiss National Science Foundation, project No.200020-121561
文摘For the numerical treatment of Hamiltonian differential equations, symplectic integrators are the most suitable choice, and methods that are conjugate to a symplectic integrator share the same good long-time behavior. This note characterizes linear multistep methods whose underlying one-step method is conjugate to a symplectic integrator. The bounded- hess of parasitic solution components is not addressed.
文摘The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equation is discretized by a mixed finite element method. The electron and hole density equations are treated by implicit-explicit multistep finite element methods. The schemes are very efficient. The optimal order error estimates both in time and space are derived.
基金supported by the National Natural Science Foundation of China(No.11971303).
文摘This paper considers the asymptotic stability of linear multistep(LM)methods for neutral systems with distributed delays.In particular,several sufficient conditions for delay-dependent stability of numerical solutions are obtained based on the argument principle.Compound quadrature formulae are used to compute the integrals.An algorithm is proposed to examine the delay-dependent stability of numerical solutions.Several numerical examples are performed to verify the theoretical results.
文摘A series of sufficient and necessary conditions for the algebraic stability of multistepRunge-Kutta methods is obtained, most of which can be regarded as extension of the relevant results available for Runge-Kutta methods, especially, for Radau Ⅰ A, Radau Ⅱ A and Gaussian Runge-Kutta methods.
基金Inner Mongolia University 2020 undergraduate teaching reform research and construction project-NDJG2094。
文摘This paper is concerned with the numerical dissipativity of multistep Runge-Kutta methods for nonlinear neutral delay-integro-differential equations.We investigate the dissipativity properties of-algebraically stable multistep Runge-Kutta methods with constrained grid.The finite-dimensional and infinite-dimensional dissipativity results of-algebraically stable multistep Runge-Kutta methods are obtained.
文摘The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the solution,characterized by||u_(n)||L^(2)(Ω)=O(t^(−αs)) as t→∞,is determined by the minimum order α_(s) of the time-fractional derivatives.Building on this foundational result,this article pursues two primary objectives.First,we introduce a strongly A-stable fractional linear multistep method and derive the numerical stability region for the governing equation.Second,we rigorously prove the long-term decay rate of the numerical solution through a detailed singularity analysis of its generating function.Notably,the numerical decay rate||u_(n)||L^(2)(Ω)=O(t_(n)^(−α_(s)) as t_(n)→∞aligns precisely with the continuous case.Theoretical findings are further validated through comprehensive numerical simulations,underscoring the robustness of our proposed method.
