Delay-dependent stability of high order Runge–Kutta methods

This paper is concerned with the study of the delay-dependent stability of Runge–Kutta methods for delay differential equations. First, a new sufficient and necessary condition is given for the asymptotic stability of analytical solution. Then, based on this condition, we establish a relationship between τ(0)-stability and the boundary locus of the stability region of numerical methods for ordinary differential equations. Consequently, a class of high order Runge–Kutta methods are proved to be τ(0)-stable. In particular, the τ(0)-stability of the Radau IIA methods is proved.     

An embedded pair of exponentially fitted explicit Runge–Kutta methods

An embedded pair of exponentially fitted explicit Runge–Kutta (RK) methods for the numerical integration of IVPs with oscillatory solutions is derived. This pair is based on the exponentially fitted explicit RK method constructed in Vanden Berghe et al., and we confirm that the methods which constitute the pair have algebraic order 4 and 3. Some numerical experiments show the efficiency of our pair when it is compared with the variable step code proposed by Vanden Berghe et al. (J. Comput. Appl. Math. 125 (2000) 107).     

High order explicit Runge–Kutta Nyström pairs

Explicit Runge–Kutta Nyström pairs provide an efficient way to find numerical solutions to second-order initial value problems when the derivative is cheap to evaluate. We present new optimal pairs of orders ten and twelve from existing families of pairs that are intended for accurate integrations in double precision arithmetic. We also present a summary of numerical comparisons between the new pairs on a set of eight problems which includes realistic models of the Solar System. Our searching for new order twelve pairs shows that there is often not quantitative agreement between the size of the principal error coefficients and the efficiency of the pairs for the tolerances we are interested in. Our numerical comparisons, as well as establishing the efficiency of the new pairs, show that the order ten pairs are more efficient than the order twelve pairs on some problems, even at limiting precision in double precision.     

Diagonally drift-implicit Runge–Kutta methods of weak order one and two for Itô SDEs and stability analysis

The class of stochastic Runge–Kutta methods for stochastic differential equations due to Rößler is considered. Coefficient families of diagonally drift-implicit stochastic Runge–Kutta (DDISRK) methods of weak order one and two are calculated. Their asymptotic stability as well as mean-square stability (MS-stability) properties are studied for a linear stochastic test equation with multiplicative noise. The stability functions for the DDISRK methods are determined and their domains of stability are compared to the corresponding domain of stability of the considered test equation. Stability regions are presented for various coefficients of the families of DDISRK methods in order to determine step size restrictions such that the numerical approximation reproduces the characteristics of the solution process.     

Approximate preservation of quadratic first integrals by explicit Runge–Kutta methods

The approximate preservation of quadratic first integrals (QFIs) of differential systems in the numerical integration with Runge–Kutta (RK) methods is studied. Conditions on the coefficients of the RK method to preserve all QFIs up to a given order are obtained, showing that the pseudo-symplectic methods studied by Aubry and Chartier (BIT 98(3):439–461, 1998) of algebraic order p preserve QFIs with order q = 2p. An expression of the error of conservation of QFIs by a RK method is given, and a new explicit six-stage formula with classical order four and seventh order of QFI-conservation is obtained by choosing their coefficients so that they minimize both local truncation and conservation errors. Several formulas with algebraic orders 3 and 4 and different orders of conservation have been tested with some problems with quadratic and general first integrals. It is shown that the new fourth-order explicit method preserves much better the qualitative properties of the flow than the standard fourth-order RK method at the price of two extra function evaluations per step and it is a practical and efficient alternative to the fully implicit methods required for a complete preservation of QFIs.     

Construction of Runge–Kutta methods of Crouch–Grossman type of high order

An approach is described to the numerical solution of order conditions for Runge–Kutta methods whose solutions evolve on a given manifold. This approach is based on least squares minimization using the Levenberg–Marquardt algorithm. Methods of order four and five are constructed and numerical experiments are presented which confirm that the derived methods have the expected order of accuracy.     

Invariants preserving schemes based on explicit Runge–Kutta methods

Numerical integration of ordinary differential equations with some invariants is considered. For such a purpose, certain projection methods have proved its high accuracy and efficiency. Unfortunately, however, sometimes they can exhibit instability. In this paper, a new, highly efficient projection method is proposed based on explicit Runge–Kutta methods. The key there is to employ the idea of the perturbed collocation method, which gives a unified way to incorporate scheme parameters for projection. Numerical experiments confirm the stability of the proposed method.     

Numerical solution of stochastic differential equations by second order Runge–Kutta methods

In this paper we propose the numerical solutions of stochastic initial value problems via random Runge–Kutta methods of the second order and mean square convergence of these methods is proved. A random mean value theorem is required and established. The concept of mean square modulus of continuity is also introduced. Expectation and variance of the approximating process are computed. Numerical examples show that the approximate solutions have a good degree of accuracy.     

Efficient exponential Runge–Kutta methods of high order: construction and implementation

BIT Numerical Mathematics - Exponential Runge–Kutta methods have shown to be competitive for the time integration of stiff semilinear parabolic PDEs. The current construction of stiffly...     

Second derivative general linear methods with inherent Runge–Kutta stability

In this paper, we find some relationships among the coefficients matrices of second derivative general linear methods (SGLMs) which are sufficient conditions, but not necessary, to ensure the methods have Runge–Kutta stability (RKS) property. Considering these conditions, we construct some A– and L–stable SGLMs with inherent RKS of orders up to five. Also, some numerical experiments for the constructed methods in variable stepsize environment are given.     

Runge–Kutta methods and Banach algebras

The equations defining both the exact and the computed solution to an initial value problem are related to a single functional equation, which can be regarded as prototypical. The functional equation can be solved in terms of a formal Taylor series, which can also be generated using an iteration process. This leads to the formal Taylor expansions of the solution and approximate solutions to initial value problems. The usual formulation, using rooted trees, can be modified to allow for linear combinations of trees, and this gives an insight into the nature of order conditions for explicit Runge–Kutta methods. A short derivation of the family of fourth order methods with four stages is given.     

A class of symplectic partitioned Runge–Kutta methods

This paper deals with some relevant properties of Runge–Kutta (RK) methods and symplectic partitioned Runge–Kutta (PRK) methods. First, it is shown that the arithmetic mean of a RK method and its adjoint counterpart is symmetric. Second, the symplectic adjoint method is introduced and a simple way to construct symplectic PRK methods via the symplectic adjoint method is provided. Some relevant properties of the adjoint method and the symplectic adjoint method are discussed. Third, a class of symplectic PRK methods are proposed based on Radau IA, Radau IIA and their adjoint methods. The structure of the PRK methods is similar to that of Lobatto IIIA–IIIB pairs and is of block forms. Finally, some examples of symplectic partitioned Runge–Kutta methods are presented.     

Practical symplectic partitioned Runge–Kutta and Runge–Kutta–Nyström methods

We present new symmetric fourth and sixth-order symplectic partitioned Runge–Kutta and Runge–Kutta–Nyström methods. We studied compositions using several extra stages, optimising the efficiency. An effective error, E_f, is defined and an extensive search is carried out using the extra parameters. The new methods have smaller values of E_f than other methods found in the literature. When applied to several examples they perform up to two orders of magnitude better than previously known method, which is in very good agreement with the values of E_f.     

Runge–Kutta methods and viscous wave equations

We study the numerical time integration of a class of viscous wave equations by means of Runge–Kutta methods. The viscous wave equation is an extension of the standard second-order wave equation including advection–diffusion terms differentiated in time. The viscous wave equation can be very stiff so that for time integration traditional explicit methods are no longer efficient. A-Stable Runge–Kutta methods are then very good candidates for time integration, in particular diagonally implicit ones. Special attention is paid to the question how the A-Stability property can be translated to this non-standard class of viscous wave equations.      

Optimal order a posteriori error estimates for a class of Runge–Kutta and Galerkin methods

We derive a posteriori error estimates, which exhibit optimal global order, for a class of time stepping methods of any order that include Runge–Kutta Collocation (RK-C) methods and the continuous Galerkin (cG) method for linear and nonlinear stiff ODEs and parabolic PDEs. The key ingredients in deriving these bounds are appropriate one-degree higher continuous reconstructions of the approximate solutions and pointwise error representations. The reconstructions are based on rather general orthogonality properties and lead to upper and lower bounds for the error regardless of the time-step; they do not hinge on asymptotics.