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.     

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.      

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.     

On some new low storage implementations of time advancing Runge–Kutta methods

In this paper, explicit Runge–Kutta (RK) schemes with minimum storage requirements for systems with very large dimension that arise in the spatial discretization of some partial differential equations are considered. A complete study of all four stage fourth-order schemes of the minimum storage families of Williamson (1980) [2], van der Houwen (1977) [8] and Ketcheson (2010) [12] that require only two storage locations per variable is carried out. It is found that, whereas there exist no schemes of this type in the Williamson and van der Houwen families, there are two isolated schemes and a one parameter family of fourth-order schemes in four stage Ketcheson’s family. This available parameter is used to obtain the optimal scheme taking into account the _‖⋅‖2${\Vert \cdot \Vert }_2$ norm of the coefficients of the leading error term. In addition a new alternative minimum storage family to the s$s$-stage Ketcheson that depends also on 3s−3$3s-3$ free parameters is proposed. This family contains both the Williamson and van der Houwen schemes but it is not included in Ketcheson’s family. Finally, the results of some numerical experiments are presented to show the behavior of fourth-order optimal schemes for some nonlinear problems.     

Runge–Kutta methods for jump–diffusion differential equations

In this paper we consider Runge–Kutta methods for jump–diffusion differential equations. We present a study of their mean-square convergence properties for problems with multiplicative noise. We are concerned with two classes of Runge–Kutta methods. First, we analyse schemes where the drift is approximated by a Runge–Kutta ansatz and the diffusion and jump part by a Maruyama term and second we discuss improved methods where mixed stochastic integrals are incorporated in the approximation of the next time step as well as the stage values of the Runge–Kutta ansatz for the drift. The second class of methods are specifically developed to improve the accuracy behaviour of problems with small noise. We present results showing when the implicit stochastic equations defining the stage values of the Runge–Kutta methods are uniquely solvable. Finally, simulation results illustrate the theoretical findings.     

Low rank Runge–Kutta methods,symplecticity and stochastic Hamiltonian problems with additive noise

In this paper we extend the ideas of Brugnano, Iavernaro and Trigiante in their development of HBVM (s,r$s,r$) methods to construct symplectic Runge–Kutta methods for all values of s$s$ and r$r$ with s≥r$s\ge r$. However, these methods do not see the dramatic performance improvement that HBVMs can attain. Nevertheless, in the case of additive stochastic Hamiltonian problems an extension of these ideas, which requires the simulation of an independent Wiener process at each stage of a Runge–Kutta method, leads to methods that have very favourable properties. These ideas are illustrated by some simple numerical tests for the modified midpoint rule.     

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.     

A-stable Runge–Kutta methods for semilinear evolution equations

We consider semilinear evolution equations for which the linear part generates a strongly continuous semigroup and the nonlinear part is sufficiently smooth on a scale of Hilbert spaces. In this setting, we prove the existence of solutions which are temporally smooth in the norm of the lowest rung of the scale for an open set of initial data on the highest rung of the scale. Under the same assumptions, we prove that a class of implicit, A-stable Runge–Kutta semidiscretizations in time of such equations are smooth as maps from open subsets of the highest rung into the lowest rung of the scale. Under the additional assumption that the linear part of the evolution equation is normal or sectorial, we prove full order convergence of the semidiscretization in time for initial data on open sets. Our results apply, in particular, to the semilinear wave equation and to the nonlinear Schrödinger equation.     

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.     

Runge–Kutta–collocation methods for systems of functional–differential and functional equations

Systems of functional–differential and functional equations occur in many biological, control and physics problems. They also include functional–differential equations of neutral type as special cases. Based on the continuous extension of the Runge–Kutta method for delay differential equations and the collocation method for functional equations, numerical methods for solving the initial value problems of systems of functional–differential and functional equations are formulated. Comprehensive analysis of the order of approximation and the numerical stability are presented.     

Explicit Runge–Kutta pairs with lower stage-order

Explicit Runge–Kutta pairs of methods of successive orders of accuracy provide effective algorithms for approximating solutions to nonstiff initial value problems. For each explicit RK method of order of accuracy p, there is a minimum number s _p of derivative evaluations required for each step propagating the numerical solution. For p ≤ 8, Butcher has established exact values of s _p , and for p > 8, his work establishes lower bounds; otherwise, upper bounds are established by various published methods. Recently, Khashin has derived some new methods numerically, and shown that the known upper bound on s ₉ for methods of order p = 9 can be reduced from 15 to 13. His results motivate this attempt to identify parametrically exact representations for coefficients of such methods. New pairs of methods of orders 5,6 and 6,7 are characterized in terms of several arbitrary parameters. This approach, modified from an earlier one, increases the known spectrum of types of RK pairs and their derivations, may lead to the derivation of new RK pairs of higher-order, and possibly to other types of explicit algorithms within the class of general linear methods.     

Continuous two-step Runge–Kutta methods for ordinary differential equations

New classes of continuous two-step Runge-Kutta methods for the numerical solution of ordinary differential equations are derived. These methods are developed imposing some interpolation and collocation conditions, in order to obtain desirable stability properties such as A-stability and L-stability. Particular structures of the stability polynomial are also investigated.     

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.     

Numerically optimal Runge–Kutta pairs with interpolants

Explicit Runge–Kutta pairs are known to provide efficient solutions to initial value differential equations with inexpensive derivative evaluations. Two criteria for selection are proposed with a view to deriving pairs of all orders 6(5) to 9(8) which minimize computation while achieving a user-specified accuracy. Coefficients of improved pairs, their stability regions and coefficients of appended optimal interpolatory Runge–Kutta formulas are provided on the author’s website (www.math.sfu.ca/~jverner). This note reports results of tests on these pairs to illustrate their effectiveness in solving nonstiff initial value problems. These pairs and interpolants may be used for implementation, or else to provide comparison targets for other new types of methods such as explicit general linear methods.     

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...