Stability of Runge-Kutta methods for delay differential systems with multiple delays

The stability of Runge-Kutta methods for systems of delay differentialequations (DDEs) with multiple delays is considered. The stabilityregions of explicit and implicit Runge-Kutta methods are discussedwhen they are applied to asymptotically stable linear DDEs withmultiple delays. A simple estimate on the stability regionsof explicit Runge-Kutta methods is presented. It is shown thatthe stable step-size for numerical integration of DDEs withmultiple delays can be easily selected by means of the estimate.     

D-CONVERGENCE OF RUNGE-KUTTA METHODS FOR STIFF DELAY DIFFERENTIAL EQUATIONS

1. IntroductionWhen considering the applicability of numerical methods for the solution of the delay differential equation (DDE) y'(t) = f(t, y(t), y(t - T)), it is necessary to analyze the error behaviourof the methods. In fact, many papers have investigated the local and global error behaviour ofDDE solvers (cL[1,2,14]). These error analyses are based on the assumption that the fUnctionf(t,y,z) satisfies Lipschitz conditions in both the last two variables. They are suitable fornonstiff …     

A stability property ofA-stable natural Runge-Kutta methods for systems of delay differential equations

A natural Runge-Kutta method is a special type of Runge-Kutta method for delay differential equations (DDEs); it is known that any one-step collocation method is equivalent to one of such methods. In this paper, we consider a linear constant-coefficient system of DDEs with a constant delay, and discuss the application of natural Runge-Kutta methods to the system. We show that anA-stable method preserves the asymptotic stability property of the analytical solutions of the system.     

Convergence Aspects of Step-Parallel Iteration of Runge-Kutta Methods for Delay Differential Equations

Implicit Runge-Kutta methods are known as highly accurate and stable methods for solving differential equations. However, the iteration technique used to solve implicit Runge-Kutta methods requires a lot of computational efforts. To lessen the computational effort, one can iterate simultaneously at a number of points along the t-axis. In this paper, we extend the PDIRK (Parallel Diagonal Iterated Runge-Kutta) methods to delay differential equations (DDEs). We give the region of convergence and analyze the speed of convergence in three parts for the P-stability region of the Runge-Kutta corrector. It is proved that PDIRK methods to DDEs are efficient, and the diagonal matrix D of the PDIRK methods for DDES can be selected in the same way as for ordinary differential equations (ODEs).     

On the asymptotic stability properties of Runge-Kutta methods for delay differential equations

Summary. In this paper asymptotic stability properties of Runge-Kutta (R-K) methods for delay differential equations (DDEs) are considered with respect to the following test equation: where and is a continuous real-valued function. In the last few years, stability properties of R-K methods applied to DDEs have been studied by numerous authors who have considered regions of asymptotic stability for “any positive delay” (and thus independent of the specific value of ). In this work we direct attention at the dependence of stability regions on a fixed delay . In particular, natural Runge-Kutta methods for DDEs are extensively examined. Received April 15, 1996 / Revised version received August 8, 1996     

Stability Analysis of Runge-Kutta Methods for Non-Linear Delay Differential Equations

This paper is concerned with the numerical solution of delay differential equations(DDEs). We focus on the stability behaviour of Runge-Kutta methods for nonlinear DDEs. The new concepts of GR(l)-stability, GAR(l)-stability and weak GAR(l)-stability are further introduced. We investigate these stability properties for (k, l)-algebraically stable Runge-Kutta methods with a piecewise constant or linear interpolation procedure.     

CONVERGENCE OF PARALLEL DIAGONAL ITERATION OF RUNGE-KUTTA METHODS FOR DELAY DIFFERENTIAL EQUATIONS

Implicit Runge-Kutta method is highly accurate and stable for stiff initial value prob-lem.But the iteration technique used to solve implicit Runge-Kutta method requires lotsof computational efforts.In this paper,we extend the Parallel Diagonal Iterated Runge-Kutta(PDIRK)methods to delay differential equations(DDEs).We give the convergenceregion of PDIRK methods,and analyze the speed of convergence in three parts for theP-stability region of the Runge-Kutta corrector method.Finally,we analysis the speed-upfactor through a numerical experiment.The results show that the PDIRK methods toDDEs are efficient.     

B-Theory of Runge-Kutta methods for stiff Volterra functional differential equations

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 equatioons(IDEs)and VFDEs of     

多步Runge-Kutta方法关于一类多延迟系统的GPk-稳定性

本文多步Ｒｕｎｇｅ－Ｋｕｔｔａ方法关于延迟微分方程系统的渐近稳定性,在本文中我们证明了在适当条件下常微多步Ｒｕｎｇｅ－Ｋｕｔｔａ方法的Ａ－稳定性等价于相应求解多延迟微分方程系统的ＧＰｋ－稳定性。     

B-Theory of Runge-Kutta methods for stiff volterra functional differential equations

B-stability andB-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra functional differential equations (VFDEs) are established which provide unified theoretical foundation for the study of Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differential equations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs of other type which appear in practice.     

A CLASS OF TWO-STEP CONTINUITY RUNGE-KUTTA METHODS FOR SOLVING SINGULAR DELAY DIFFERENTIAL EQUATIONS AND ITS STABILITY ANALYSIS

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.     

Dissipativity of multistep runge-kutta methods for nonlinear volterra delay-integro-differential equations

This paper is concerned with the numerical dissipativity of multistep Runge-Kutta methods for nonlinear Volterra delay-integro-differential equations.We investigate the dissipativity properties of (k,l)algebraically stable multistep Runge-Kutta methods with constrained grid and an uniform grid.The finitedimensional and infinite-dimensional dissipativity results of (k,l)-algebraically stable Runge-Kutta methods are obtained.     

含分布时滞的时滞微分系统多步龙格-库塔方法的时滞相关稳定性

本文研究了一类含分布时滞的时滞微分系统的多步龙格-库塔方法的稳定性.基于辐角原理，本文给出了多步龙格-库塔方法弱时滞相关稳定性的充分条件，并通过数值算例验证了理论结果的有效性.     

Numerical solution by LMMs of stiff delay differential systems modelling an immune response

Summary. We consider the application of linear multistep methods (LMMs) for the numerical solution of initial value problem for stiff delay differential equations (DDEs) with several constant delays, which are used in mathematical modelling of immune response. For the approximation of delayed variables the Nordsieck's interpolation technique, providing an interpolation procedure consistent with the underlying linear multistep formula, is used. An analysis of the convergence for a variable-stepsize and structure of the asymptotic expansion of global error for a fixed-stepsize is presented. Some absolute stability characteristics of the method are examined. Implementation details of the code DIFSUB-DDE, being a modification of the Gear's DIFSUB, are given. Finally, an efficiency of the code developed for solution of stiff DDEs over a wide range of tolerances is illustrated on biomedical application model. Received March 23, 1994 / Revised version received March 13, 1995     

Non-linear stability of a general class of differential equation methods

For a class of methods sufficiently general as to include linear multistep and Runge-Kutta methods as special cases, a concept known as algebraic stability is defined. This property is based on a non-linear test problem and extends existing results on Runge-Kutta methods and on linear multistep and one-leg methods. The algebraic stability properties of a number of particular methods in these families are studied and a generalization is made which enables estimates of error growth to be provided for certain classes of methods.     

求解延迟微分代数方程的多步Runge-Kutta方法的渐近稳定性

延迟微分代数方程(DDAEs)广泛出现于科学与工程应用领域．本文将多步Runge-Kutta方法应用于求解线性常系数延迟微分代数方程，讨论了该方法的渐近稳定性．数值试验表明该方法对求解DDAEs是有效的．     

CHARACTERIZATIONS OF SYMMETRIC MULTISTEP RUNGE-KUTTA METHODS

Some characterizations for symmetric multistep Runge-Kutta(RK) methods are obtained. Symmetric two-step RK methods with one and two-stages are presented. Numerical examples show that symmetry of multistep RK methods alone is not sufficient for long time integration for reversible Hamiltonian systems. This is an important difference between one-step and multistep symmetric RK methods.     

Asymptotic stability analysis of Runge-Kutta methods for nonlinear systems of delay differential equations

Summary. We consider systems of delay differential equations (DDEs) of the form with the initial condition . Recently, Torelli [10] introduced a concept of stability for numerical methods applied to dissipative nonlinear systems of DDEs (in some inner product norm), namely RN-stability, which is the straighforward generalization of the wellknown concept of BN-stability of numerical methods with respect to dissipative systems of ODEs. Dissipativity means that the solutions and corresponding to different initial functions and , respectively, satisfy the inequality , and is guaranteed by suitable conditions on the Lipschitz constants of the right-hand side function . A numerical method is said to be RN-stable if it preserves this contractivity property. After showing that, under slightly more stringent hypotheses on the Lipschitz constants and on the delay function , the solutions and are such that , in this paper we prove that RN-stable continuous Runge-Kutta methods preserve also this asymptotic stability property. Received March 29, 1996 / Revised version received August 12, 1996     

MULTISTEP DISCRETIZATION OF INDEX 3 DAES

1. IntroductionIn this paper, we will consider the multistep discrezations of the differential--algebraicequations (DAEs) in Hessenberg formwhere F e AN M M R", K E AN M L - AM, G E AN - RL, the initial value(yo, ic, no) at xo are assumed to be consistent, i.e., they satisfyWe supposes, F, G and K are sufficiently differentiable, and thatin a neighbourhood of the solution. Such problems often appear in the simulation ofmechanical systems with constraints and the singularly perturbed…     

Dissipativity of Runge-Kutta methods for Volterra functional differential equations

This paper is concerned with the numerical dissipativity of nonlinear Volterra functional differential equations (VFDEs). We give some dissipativity results of Runge-Kutta methods when they are applied to VFDEs. These results provide unified theoretical foundation for the numerical dissipativity analysis of systems in ordinary differential equations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs), Volterra delay integro-differential equations (VDIDEs) and VFDEs of other type which appear in practice. Numerical examples are given to confirm our theoretical results.