An idea of relaxing the effect of delay when computing the Runge-Kutta stages in the current step and a class of two-step continuity Runge-Kutta methods (TSCRK) is presented. Their construction, their order conditio...An idea of relaxing the effect of delay when computing the Runge-Kutta stages in the current step and a class of two-step continuity Runge-Kutta methods (TSCRK) is presented. Their construction, their order conditions and their convergence are studied. The two-step continuity Runge-Kutta methods possess good numerical stability properties and higher stage-order, and keep the explicit process of computing the Runge-Kutta stages. The numerical experiments show that the TSCRK methods are efficient.展开更多
In this paper we discuss diagonally implicit and semi-implicit methods based on the three-stage stiffly accurate Runge-Kutta methods for solving Stratonovich stochastic differential equations(SDEs).Two methods,a thr...In this paper we discuss diagonally implicit and semi-implicit methods based on the three-stage stiffly accurate Runge-Kutta methods for solving Stratonovich stochastic differential equations(SDEs).Two methods,a three-stage stiffly accurate semi-implicit(SASI3) method and a three-stage stiffly accurate diagonally implicit (SADI3) method,are constructed in this paper.In particular,the truncated random variable is used in the implicit method.The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of stiff SDEs.展开更多
In this paper, the Ito-Taylor expansion of stochastic differential equation is briefly introduced. The colored rooted tree theory is applied to derive strong order 1.0 implicit stochastic Runge-Kutta method(SRK). Two ...In this paper, the Ito-Taylor expansion of stochastic differential equation is briefly introduced. The colored rooted tree theory is applied to derive strong order 1.0 implicit stochastic Runge-Kutta method(SRK). Two fully implicit schemes are presented and their stability qualities are discussed. And the numerical report illustrates the better numerical behavior.展开更多
This paper deals with H-stability of the Runge-Kutta methods with a general variable stepsize for the system of pantograph equations with two delay terms. It is shown that the Runge-Kutta methods with a regular matrix...This paper deals with H-stability of the Runge-Kutta methods with a general variable stepsize for the system of pantograph equations with two delay terms. It is shown that the Runge-Kutta methods with a regular matrix A are H-stable if and only if the modulus of the stability function at infinity is less than 1.展开更多
Differential equation has widely applied in science and engineering calculation. Runge Kutta method is a main method for solving differential equations. In this paper, the numerical properties of Runge-Kutta methods f...Differential equation has widely applied in science and engineering calculation. Runge Kutta method is a main method for solving differential equations. In this paper, the numerical properties of Runge-Kutta methods for the equation u′(t) = au(t)+bu([K/N* t]) is dealed with, where K and N is relatively prime and K < N,K,N∈ Z+. The conditions are obtained under which the numerical solutions preserve the analytical stability properties of the analytic ones and some numerical experiments are given.展开更多
In this paper, a class of two-step continuity Runge-Kutta(TSCRK) methods for solving singular delay differential equations(DDEs) is presented. Analysis of numerical stability of this methods is given. We consider ...In this paper, a class of two-step continuity Runge-Kutta(TSCRK) methods for solving singular delay differential equations(DDEs) is presented. Analysis of numerical stability of this methods is given. We consider the two distinct cases: (i)τ≥ h, (ii)τ 〈 h, where the delay τ and step size h of the two-step continuity Runge-Kutta methods are both constant. The absolute stability regions of some methods are plotted and numerical examples show the efficiency of the method.展开更多
Presents the stability analysis of theoretical solutions for a class of nonlinear neutral delay-differential equations (NDDE). Discussion on the numerical analogous results of the natural Runge-Kutta (NRK) methods for...Presents the stability analysis of theoretical solutions for a class of nonlinear neutral delay-differential equations (NDDE). Discussion on the numerical analogous results of the natural Runge-Kutta (NRK) methods for the same class of nonlinear NDDE; Review of the related concepts and results on RK methods; Information on the asymptotic stability and global stability of the induced NRK method.展开更多
B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations (VFDEs) are established which provide unified theoretical foundation for the studyof Runge-...B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations (VFDEs) are established which provide unified theoretical foundation for the studyof Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differentialequations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs ofother type which appear in practice.展开更多
Focuses on a study which examined the numerical solution of delay differential equations (DDE). Information on the Runge-Kutta methods for DDE; Results of the D-convergence analysis of Runge-Kutta methods; Details on ...Focuses on a study which examined the numerical solution of delay differential equations (DDE). Information on the Runge-Kutta methods for DDE; Results of the D-convergence analysis of Runge-Kutta methods; Details on the equations with several delays.展开更多
This paper is concerned with the numerical stability of implicit Runge-Kutta methods for nonlinear neutral Volterra delay-integro-differential equations with constant delay.Using a Halanay inequality generalized by Li...This paper is concerned with the numerical stability of implicit Runge-Kutta methods for nonlinear neutral Volterra delay-integro-differential equations with constant delay.Using a Halanay inequality generalized by Liz and Trofimchuk,we give two sufficient conditions for the stability of the true solution to this class of equations.Runge-Kutta methods with compound quadrature rule are considered.Nonlinear stability conditions for the proposed methods are derived.As an illustration of the application of these investigations,the asymptotic stability of the presented methods for Volterra delay-integro-differential equations are proved under some weaker conditions than those in the literature.An extension of the stability results to such equations with weakly singular kernel is also discussed.展开更多
Presents a study which dealt with stability of the implicit Runge-Kutta methods for the numerical solutions of the systems of delay differential equations. Stability behavior of the methods; Exponential solutions of t...Presents a study which dealt with stability of the implicit Runge-Kutta methods for the numerical solutions of the systems of delay differential equations. Stability behavior of the methods; Exponential solutions of the equations.展开更多
The adapting Runge-Kutta methods with a new interpolation procedure to delay differential equations was introduced by K.J. in't Hout in 1992[1], he proved that the numerical process, satisfies an important asympto...The adapting Runge-Kutta methods with a new interpolation procedure to delay differential equations was introduced by K.J. in't Hout in 1992[1], he proved that the numerical process, satisfies an important asymptotic stability condition. In this paper the convergence of this method under the asymptotic stability and other conditions in theorem 3 is proved.展开更多
The aim of this paper is to obtain numerical solutions of the one-dimensional,two-dimensional and coupled Burgers' equations through the generalized differential quadrature method(GDQM).The polynomial-based differ...The aim of this paper is to obtain numerical solutions of the one-dimensional,two-dimensional and coupled Burgers' equations through the generalized differential quadrature method(GDQM).The polynomial-based differential quadrature(PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta(TVD-RK) method.The numerical solutions are satisfactorily coincident with the exact solutions.The method can compete against the methods applied in the literature.展开更多
This paper is describing in detail the way we define the equations which give the formulas in the methods Runge-Kutta 6<sup>th</sup> order 7 stages with the incorporated control step size in the numerical ...This paper is describing in detail the way we define the equations which give the formulas in the methods Runge-Kutta 6<sup>th</sup> order 7 stages with the incorporated control step size in the numerical solution of Ordinary Differential Equations (ODE). The purpose of the present work is to construct a system of nonlinear equations and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (6,7) method (6<sup>th</sup> order and 7 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is complicated, all coefficients are found with respect to 7 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically. Some examples for five different choices of the arbitrary values of the systems are presented in this paper.展开更多
The purpose of the present work is to construct a nonlinear equation system (85 × 53) using Butcher’s Table and then by solving this system to find the values of all set parameters and finally the reduction form...The purpose of the present work is to construct a nonlinear equation system (85 × 53) using Butcher’s Table and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (7,9) method (7<sup>th</sup> order and 9 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is too complicated, we introduce a subsystem from the original system where all coefficients are found with respect to 9 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically.展开更多
In this paper the Modified Equations of Emden type (MEE), χ+αχχ+βχ 3 is solved numerically by the differential transform method. This technique doesn’t require any discretization, linearization or small perturb...In this paper the Modified Equations of Emden type (MEE), χ+αχχ+βχ 3 is solved numerically by the differential transform method. This technique doesn’t require any discretization, linearization or small perturbations and therefore it reduces significantly the numerical computation. The current results of this paper are in excellent agreement with those provided by Chandrasekar et al. [1] and thereby illustrate the reliability and the performance of the differential transform method. We have also compared the results with the classical Runge-Kutta 4 (RK4) Method.展开更多
It is well known that mono-implicit Runge-Kutta methods have been applied in the efficient numerical solution of initial or boundary value problems of ordinary differential equations. Burrage (1994) has shown that the...It is well known that mono-implicit Runge-Kutta methods have been applied in the efficient numerical solution of initial or boundary value problems of ordinary differential equations. Burrage (1994) has shown that the order of an s-stage monoimplicit Runge-Kutta method is at most s+1 and the stage order is at most 3. In this paper, it is shown that the order of an s-stage mono-implicit Runge-Kutta method being algebraically stable is at most min((s) over tilde, 4), and the stage order together with the optimal B-convergence order is at most min(s, 2), where [GRAPHICS]展开更多
We study the construction of symplectic Runge-Kutta methods for stochastic Hamiltonian systems(SHS).Three types of systems,SHS with multiplicative noise,special separable Hamiltonians and multiple additive noise,respe...We study the construction of symplectic Runge-Kutta methods for stochastic Hamiltonian systems(SHS).Three types of systems,SHS with multiplicative noise,special separable Hamiltonians and multiple additive noise,respectively,are considered in this paper.Stochastic Runge-Kutta(SRK)methods for these systems are investigated,and the corresponding conditions for SRK methods to preserve the symplectic property are given.Based on the weak/strong order and symplectic conditions,some effective schemes are derived.In particular,using the algebraic computation,we obtained two classes of high weak order symplectic Runge-Kutta methods for SHS with a single multiplicative noise,and two classes of high strong order symplectic Runge-Kutta methods for SHS with multiple multiplicative and additive noise,respectively.The numerical case studies confirm that the symplectic methods are efficient computational tools for long-term simulations.展开更多
In [4] we proved that all Gauss methods areNtau(0)-compatible for neutral delay differential equations (NDDEs) of the form y'(t) = ay(t) + by(t-tau) + cy'(t-tau), t >0, (0.1) y(t) = g(t), -tau less than or ...In [4] we proved that all Gauss methods areNtau(0)-compatible for neutral delay differential equations (NDDEs) of the form y'(t) = ay(t) + by(t-tau) + cy'(t-tau), t >0, (0.1) y(t) = g(t), -tau less than or equal to t less than or equal to 0, where a, b, c are real, tau > 0, g(t) is a continuous real valued function. In this paper we are going to use the theory of order stars to characterize the asymptotic stability properties of Gauss methods for NDDEs. And then proved that all Gauss methods are Ntau(0)-stable.展开更多
文摘An idea of relaxing the effect of delay when computing the Runge-Kutta stages in the current step and a class of two-step continuity Runge-Kutta methods (TSCRK) is presented. Their construction, their order conditions and their convergence are studied. The two-step continuity Runge-Kutta methods possess good numerical stability properties and higher stage-order, and keep the explicit process of computing the Runge-Kutta stages. The numerical experiments show that the TSCRK methods are efficient.
基金supported by the NSF(10926158) of ChinaDoctoral Fund(20090061120038) of Ministry of Education of ChinaBasic Scientific Research Foundation(200903287) of Jilin University
文摘In this paper we discuss diagonally implicit and semi-implicit methods based on the three-stage stiffly accurate Runge-Kutta methods for solving Stratonovich stochastic differential equations(SDEs).Two methods,a three-stage stiffly accurate semi-implicit(SASI3) method and a three-stage stiffly accurate diagonally implicit (SADI3) method,are constructed in this paper.In particular,the truncated random variable is used in the implicit method.The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of stiff SDEs.
文摘In this paper, the Ito-Taylor expansion of stochastic differential equation is briefly introduced. The colored rooted tree theory is applied to derive strong order 1.0 implicit stochastic Runge-Kutta method(SRK). Two fully implicit schemes are presented and their stability qualities are discussed. And the numerical report illustrates the better numerical behavior.
文摘This paper deals with H-stability of the Runge-Kutta methods with a general variable stepsize for the system of pantograph equations with two delay terms. It is shown that the Runge-Kutta methods with a regular matrix A are H-stable if and only if the modulus of the stability function at infinity is less than 1.
基金This work is supported by the Research Fund of the Natural Science Foundation of Heilongjiang Province (No. A201214) and the National Natural Science Foundation of China(61501148).
文摘Differential equation has widely applied in science and engineering calculation. Runge Kutta method is a main method for solving differential equations. In this paper, the numerical properties of Runge-Kutta methods for the equation u′(t) = au(t)+bu([K/N* t]) is dealed with, where K and N is relatively prime and K < N,K,N∈ Z+. The conditions are obtained under which the numerical solutions preserve the analytical stability properties of the analytic ones and some numerical experiments are given.
文摘In this paper, a class of two-step continuity Runge-Kutta(TSCRK) methods for solving singular delay differential equations(DDEs) is presented. Analysis of numerical stability of this methods is given. We consider the two distinct cases: (i)τ≥ h, (ii)τ 〈 h, where the delay τ and step size h of the two-step continuity Runge-Kutta methods are both constant. The absolute stability regions of some methods are plotted and numerical examples show the efficiency of the method.
文摘Presents the stability analysis of theoretical solutions for a class of nonlinear neutral delay-differential equations (NDDE). Discussion on the numerical analogous results of the natural Runge-Kutta (NRK) methods for the same class of nonlinear NDDE; Review of the related concepts and results on RK methods; Information on the asymptotic stability and global stability of the induced NRK method.
文摘B-stability and B-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra func-tional differential equations (VFDEs) are established which provide unified theoretical foundation for the studyof Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differentialequations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs ofother type which appear in practice.
基金NSF of China (No.19871070) and China Postdoctoral Science Foundation.
文摘Focuses on a study which examined the numerical solution of delay differential equations (DDE). Information on the Runge-Kutta methods for DDE; Results of the D-convergence analysis of Runge-Kutta methods; Details on the equations with several delays.
基金supported by NSF of China(Grant No.11001033)Natural Science Foundation of Hunan Province(Grant No.10JJ4003)Chinese Society for Electrical Engineering,and Graduates’innovation fund of HUST(No.HF-08-02-2011-011).
文摘This paper is concerned with the numerical stability of implicit Runge-Kutta methods for nonlinear neutral Volterra delay-integro-differential equations with constant delay.Using a Halanay inequality generalized by Liz and Trofimchuk,we give two sufficient conditions for the stability of the true solution to this class of equations.Runge-Kutta methods with compound quadrature rule are considered.Nonlinear stability conditions for the proposed methods are derived.As an illustration of the application of these investigations,the asymptotic stability of the presented methods for Volterra delay-integro-differential equations are proved under some weaker conditions than those in the literature.An extension of the stability results to such equations with weakly singular kernel is also discussed.
文摘Presents a study which dealt with stability of the implicit Runge-Kutta methods for the numerical solutions of the systems of delay differential equations. Stability behavior of the methods; Exponential solutions of the equations.
文摘The adapting Runge-Kutta methods with a new interpolation procedure to delay differential equations was introduced by K.J. in't Hout in 1992[1], he proved that the numerical process, satisfies an important asymptotic stability condition. In this paper the convergence of this method under the asymptotic stability and other conditions in theorem 3 is proved.
文摘The aim of this paper is to obtain numerical solutions of the one-dimensional,two-dimensional and coupled Burgers' equations through the generalized differential quadrature method(GDQM).The polynomial-based differential quadrature(PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta(TVD-RK) method.The numerical solutions are satisfactorily coincident with the exact solutions.The method can compete against the methods applied in the literature.
文摘This paper is describing in detail the way we define the equations which give the formulas in the methods Runge-Kutta 6<sup>th</sup> order 7 stages with the incorporated control step size in the numerical solution of Ordinary Differential Equations (ODE). The purpose of the present work is to construct a system of nonlinear equations and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (6,7) method (6<sup>th</sup> order and 7 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is complicated, all coefficients are found with respect to 7 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically. Some examples for five different choices of the arbitrary values of the systems are presented in this paper.
文摘The purpose of the present work is to construct a nonlinear equation system (85 × 53) using Butcher’s Table and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (7,9) method (7<sup>th</sup> order and 9 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is too complicated, we introduce a subsystem from the original system where all coefficients are found with respect to 9 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically.
文摘In this paper the Modified Equations of Emden type (MEE), χ+αχχ+βχ 3 is solved numerically by the differential transform method. This technique doesn’t require any discretization, linearization or small perturbations and therefore it reduces significantly the numerical computation. The current results of this paper are in excellent agreement with those provided by Chandrasekar et al. [1] and thereby illustrate the reliability and the performance of the differential transform method. We have also compared the results with the classical Runge-Kutta 4 (RK4) Method.
文摘It is well known that mono-implicit Runge-Kutta methods have been applied in the efficient numerical solution of initial or boundary value problems of ordinary differential equations. Burrage (1994) has shown that the order of an s-stage monoimplicit Runge-Kutta method is at most s+1 and the stage order is at most 3. In this paper, it is shown that the order of an s-stage mono-implicit Runge-Kutta method being algebraically stable is at most min((s) over tilde, 4), and the stage order together with the optimal B-convergence order is at most min(s, 2), where [GRAPHICS]
基金This work was supported by NSFC(91130003)The first authors is also supported by NSFC(11101184,11271151)+1 种基金the Science Foundation for Young Scientists of Jilin Province(20130522101JH)The second and third authors are also supported by NSFC(11021101,11290142).The authors would like to thank anonymous reviewers for careful reading and invaluable suggestions,which greatly improved the presentation of the paper.
文摘We study the construction of symplectic Runge-Kutta methods for stochastic Hamiltonian systems(SHS).Three types of systems,SHS with multiplicative noise,special separable Hamiltonians and multiple additive noise,respectively,are considered in this paper.Stochastic Runge-Kutta(SRK)methods for these systems are investigated,and the corresponding conditions for SRK methods to preserve the symplectic property are given.Based on the weak/strong order and symplectic conditions,some effective schemes are derived.In particular,using the algebraic computation,we obtained two classes of high weak order symplectic Runge-Kutta methods for SHS with a single multiplicative noise,and two classes of high strong order symplectic Runge-Kutta methods for SHS with multiple multiplicative and additive noise,respectively.The numerical case studies confirm that the symplectic methods are efficient computational tools for long-term simulations.
基金the National Natural Science Foundation of China.
文摘In [4] we proved that all Gauss methods areNtau(0)-compatible for neutral delay differential equations (NDDEs) of the form y'(t) = ay(t) + by(t-tau) + cy'(t-tau), t >0, (0.1) y(t) = g(t), -tau less than or equal to t less than or equal to 0, where a, b, c are real, tau > 0, g(t) is a continuous real valued function. In this paper we are going to use the theory of order stars to characterize the asymptotic stability properties of Gauss methods for NDDEs. And then proved that all Gauss methods are Ntau(0)-stable.