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.     

A Class of Explicit Exponential General Linear Methods

In this paper, we consider a class of explicit exponential integrators that includes as special cases the explicit exponential Runge–Kutta and exponential Adams–Bashforth methods. The additional freedom in the choice of the numerical schemes allows, in an easy manner, the construction of methods of arbitrarily high order with good stability properties. We provide a convergence analysis for abstract evolution equations in Banach spaces including semilinear parabolic initial-boundary value problems and spatial discretizations thereof. From this analysis, we deduce order conditions which in turn form the basis for the construction of new schemes. Our convergence results are illustrated by numerical examples. AMS subject classification (2000) 65L05, 65L06, 65M12, 65J10     

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 …     

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.     

Volterra方程数值解中带自适应步长控制的Runge-Kutta方法(英文)

在本文中我们构造了解第二类Volterra方程的一般Runge-Kutta方法,并且研究了第二类Volterra方程数值解法的自适应步长控制。     

Geometric Integration Methods that Preserve Lyapunov Functions

We consider ordinary differential equations (ODEs) with a known Lyapunov function V. To ensure that a numerical integrator reflects the correct dynamical behaviour of the system, the numerical integrator should have V as a discrete Lyapunov function. Only second-order geometric integrators of this type are known for arbitrary Lyapunov functions. In this paper we describe projection-based methods of arbitrary order that preserve any given Lyapunov function. AMS subject classification (2000) 65L05, 65L06, 65L20, 65P40     

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.     

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

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

STABILITY ANALYSIS OF RUNGE-KUTTA METHODS FOR NONLINEAR SYSTEMS OF PANTOGRAPH EQUATIONS

This paper is concerned with numerical stability of nonlinear systems of pantograph equations. Numerical methods based on （k, l）-algebraically stable Runge-Kutta methods are suggested. Global and asymptotic stability conditions for the presented methods are derived.     

Dissipativity of Runge-Kutta methods for dynamical systems with delays

This paper is concerned with the numerical solution of dissipativeinitial value problems with delays by Runge-Kutta methods. Asufficient condition for the dissipativity of the systems isgiven. The concepts of D(l)-dissipativity and GD(l)-dissipativityare introduced. We investigate the dissipativity propertiesof (k,l)-algebraically stable Runge-Kutta methods with piecewiseconstant or linear interpolation procedures for finite-dimensionaland infinite-dimensional dynamical systems with delays.     

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.     

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

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

Maximal order for second derivative general linear methods with Runge-Kutta stability

An extension of general linear methods (GLMs), so-called SGLMs (GLMs with second derivative), was introduced to the case in which second derivatives, as well as first derivatives, can be calculated. SGLMs are divided into four types, depending on the nature of the differential system to be solved and the computer architecture that is used to implement these methods. In this paper, we obtain maximal order for two types of SGLMs with Runge-Kutta stability (RKS) property. Also, we construct methods of these types which possess RKS property and A-stability. Efficiency of the constructed methods is shown by numerical experiments.     

Runge-Kutta Methods for the Numerical Solution of Stiff Semilinear Systems

This paper studies the stability and convergence properties of general Runge-Kutta methods when they are applied to stiff semilinear systems y(t) = J(t)y(t) + g(t, y(t)) with the stiffness contained in the variable coefficient linear part.We consider two assumptions on the relative variation of the matrix J(t) and show that for each of them there is a family of implicit Runge-Kutta methods that is suitable for the numerical integration of the corresponding stiff semilinear systems, i.e. the methods of the family are stable, convergent and the stage equations possess a unique solution. The conditions on the coefficients of a method to belong to these families turn out to be essentially weaker than the usual algebraic stability condition which appears in connection with the B-stability and convergence for stiff nonlinear systems. Thus there are important RK methods which are not algebraically stable but, according to our theory, they are suitable for the numerical integration of semilinear problems.This paper also extends previous results of Burrage, Hundsdorfer and Verwer on the optimal convergence of implicit Runge-Kutta methods for stiff semilinear systems with a constant coefficients linear part.     

On the Explicit Two-Stage Fourth-Order Accurate Time Discretizations

This paper continues to study the explicit two-stage fourth-order accurate time discretizations [5-7]. By introducing variable weights, we propose a class of more general explicit one-step two-stage time discretizations, which are different from the existing methods, e.g. the Euler methods, Runge-Kutta methods, and multistage multiderivative methods etc. We study the absolute stability, the stability interval, and the intersection between the imaginary axis and the absolute stability region. Our results show that our two-stage time discretizations can be fourth-order accurate conditionally, the absolute stability region of the proposed methods with some special choices of the variable weights can be larger than that of the classical explicit fourth- or fifth-order Runge-Kutta method, and the interval of absolute stability can be almost twice as much as the latter. Several numerical experiments are carried out to demonstrate the performance and accuracy as well as the stability of our proposed methods.     

A modified Runge-Kutta method for the numerical solution of ODE''s with oscillation solutions

A modified Runge-Kutta method with minimal phase-lag is developed for the numerical solution of Ordinary Differential Equations with oscillating solutions. The method is based on the accurate Runge-Kutta method of Sharp and Smart RK4SS(5) (see [1]) of order five. Numerical and theoretical results show that this new approach is more efficient, compared with the fifth order Runge-Kutta Sharp and Smart method.     

ORDER PROPERTIES AND CONSTRUCTION OF SYMPLECTIC RUNGE-KUTTA METHODS

1. IntroductionFOr a given s stage Runge-Kutta methodwith A = [ail], p = [pl, PZt... 5 P.]T and ac = [afl, ry23... ) %]T / 0, we introduce thefollowing simplifying conditions as in Butcher [1]and make the notational convensionwhere 1 5 m? pi(x), i ~ 1, 2, 3,' ? are arbitrarily given i--th polynomials with the property that pi(0) = 0,Note that B(P), C(P) and D(P) are equivalent to BI,. = 0, CI,P = 0 and DI,. = 0respectively. We shall always denote BI,., CI,., DI,. and VI,. by B, …     

Two-Step Runge-Kutta: Theory and Practice

Local and global error for Two-Step Runge-Kutta (TSRK) methods are analyzed using the theory of B-series. Global error bounds are derived in both constant and variable stepsize environments. An embedded TSRK pair is constructed and compared with the RK5(4)6M pair of Dormand and Prince on the DETEST set of problems. Numerical results show that the TSRK performs competitively with the RK method.     

求解多延迟微分方程的Runge-Kutta方法的收缩性

该文涉及多延迟微分方程MDDEs系统的理论解与数值解的收缩性．为此,一些新的稳定性概念诸如：BN_f^(m)-稳定性及GRN_m-稳定性稳定性被引入．该探讨得出：Ｒｕｎｇｅ Ｋｕｔｔａ（ＲＫ）方法及相应的连续插值的BN^(m)-稳定性导致求解MDDEs的方法的收缩性（GRN_m-稳定性）．     

Construction of two-step Runge-Kutta methods of high order for ordinary differential equations

The construction of two-step Runge-Kutta methods of order p and stage order q=p with stability polynomial given in advance is described. This polynomial is chosen to have a large interval of absolute stability for explicit methods and to be A-stable and L-stable for implicit methods. After satisfying the order and stage order conditions the remaining free parameters are computed by minimizing the sum of squares of the difference between the stability function of the method and a given polynomial at a sufficiently large number of points in the complex plane. This revised version was published online in August 2006 with corrections to the Cover Date.