Approximating Runge-Kutta matrices by triangular matrices

The implementation of implicit Runge-Kutta methods requires the solution of large systems of non-linear equations. Normally these equations are solved by a modified Newton process, which can be very expensive for problems of high dimension. The recently proposed triangularly implicit iteration methods for ODE-IVP solvers [5] substitute the Runge-Kutta matrixA in the Newton process for a triangular matrixT that approximatesA, hereby making the method suitable for parallel implementation. The matrixT is constructed according to a simple procedure, such that the stiff error components in the numerical solution are strongly damped. In this paper we prove for a large class of Runge-Kutta methods that this procedure can be carried out and that the diagnoal entries ofT are positive. This means that the linear systems that are to be solved have a non-singular matrix. The research reported in this paper was supported by STW (Dutch Foundation for Technical Sciences).     

High orderP-stable formulae for the numerical integration of periodic initial value problems

Summary Recently there has been considerable interest in the approximate numerical integration of the special initial value problemy=f(x, y) for cases where it is known in advance that the required solution is periodic. The well known class of Störmer-Cowell methods with stepnumber greater than 2 exhibit orbital instability and so are often unsuitable for the integration of such problems. An appropriate stability requirement for the numerical integration of periodic problems is that ofP-stability. However Lambert and Watson have shown that aP-stable linear multistep method cannot have an order of accuracy greater than 2. In the present paper a class of 2-step methods of Runge-Kutta type is discussed for the numerical solution of periodic initial value problems.P-stable formulae with orders up to 6 are derived and these are shown to compare favourably with existing methods.     

Total variation diminishing Runge-Kutta schemes

In this paper we further explore a class of high order TVD (total variation diminishing) Runge-Kutta time discretization initialized in a paper by Shu and Osher, suitable for solving hyperbolic conservation laws with stable spatial discretizations. We illustrate with numerical examples that non-TVD but linearly stable Runge-Kutta time discretization can generate oscillations even for TVD (total variation diminishing) spatial discretization, verifying the claim that TVD Runge-Kutta methods are important for such applications. We then explore the issue of optimal TVD Runge-Kutta methods for second, third and fourth order, and for low storage Runge-Kutta methods.

     

A New Integrator for Special Third Order Differential Equations With Application to Thin Film Flow Problem

In recent time, Runge-Kutta methods that integrate special third order ordinary differential equations (ODEs) directly are proposed to address efficiency issues associated with classical Runge-Kutta methods. Albeit, the methods require evaluation of three set of equations to proceed with the numerical integration. In this paper, we propose a class of multistep-like Runge-Kutta methods (hybrid methods), which integrates special third order ODEs directly. The method is completely derivative-free. Algebraic order conditions of the method are derived. Using the order conditions, a four-stage method is presented. Numerical experiment is conducted on some test problems. The method is also applied to a practical problem in Physics and engineering to ascertain its validity. Results from the experiment show that the new method is more accurate and efficient than the classical Runge-Kutta methods and a class of direct Runge-Kutta methods recently designed for special third order ODEs.     

The efficiency of Singly-implicit Runge-Kutta methods for stiff differential equations

Singly-implicit Runge-Kutta methods are considered to be good candidates for stiff problems because of their good stability and high accuracy. The existing methods, SIRK (Singly-implicit Runge-Kutta), DESI (Diagonally Extendable Singly-implicit Runge-Kutta), ESIRK (Effective order Singly-implicit Rung-Kutta) and DESIRE (Diagonally Extended Singly-implicit Runge-Kutta Effective order) methods have been shown to be efficient for stiff differential equations, especially for high dimensional stiff problems. In this paper, we measure the efficiency for the family of singly-implicit Runge-Kutta methods using the local truncation error produced within one single step and the count of number of operations. Verification of the error and the computational costs for these methods using variable stepsize scheme are presented. We show how the numerical results are effected by the designed factors: additional diagonal-implicit stages and effective order.     

Multiple invariants conserving Runge-Kutta type methods for Hamiltonian problems

In a recent series of papers, the class of energy-conserving Runge-Kutta methods named Hamiltonian BVMs (HBVMs) has been defined and studied. Such methods have been further generalized for the efficient solution of general conservative problems, thus providing the class of Line Integral Methods (LIMs). In this paper we derive a further extension, which we name Enhanced Line Integral Methods (ELIMs), more tailored for Hamiltonian problems, allowing for the conservation of multiple invariants of the continuous dynamical system. The analysis of the methods is fully carried out and some numerical tests are reported, in order to confirm the theoretical achievements.     

A note on error expressions for reflected and averaged implicit Runge-Kutta methods

In addition to their usefulness in the numerical solution of initial value ODE's, the implicit Runge-Kutta (IRK) methods are also important for the solution of two-point boundary value problems. Recently, several classes of modified IRK methods which improve significantly on the efficiency of the standard IRK methods in this application have been presented. One such class is the Averaged IRK methods; a member of the class is obtained by applying an averaging operation to a non-symmetric IRK method and its reflection. In this paper we investigate the forms of the error expressions for reflected and averaged IRK methods. Our first result relates the expression for the local error of the reflected method to that of the original method. The main result of this paper relates the error expression of an averaged method to that of the method upon which it is based. We apply these results to show that for each member of the class of the averaged methods, there exists an embedded lower order method which can be used for error estimation, in a formula-pair fashion.This work was supported by the Natural Science and Engineering Research Council of Canada.     

A bound on the maximum strong order of stochastic Runge-Kutta methods for stochastic ordinary differential equations

In Burrage and Burrage [1] it was shown that by introducing a very general formulation for stochastic Runge-Kutta methods, the previous strong order barrier of order one could be broken without having to use higher derivative terms. In particular, methods of strong order 1.5 were developed in which a Stratonovich integral of order one and one of order two were present in the formulation. In this present paper, general order results are proven about the maximum attainable strong order of these stochastic Runge-Kutta methods (SRKs) in terms of the order of the Stratonovich integrals appearing in the Runge-Kutta formulation. In particular, it will be shown that if ans-stage SRK contains Stratonovich integrals up to orderp then the strong order of the SRK cannot exceed min{(p+1)/2, (s−1)/2},p≥2,s≥3 or 1 ifp=1.     

A special family of Runge-Kutta methods for solving stiff differential equations

An efficient way of implementing Implicit Runge-Kutta Methods was proposed by Butcher [3]. He showed that the most efficient methods when using this implementation are those whose characteristic polynomial of the Runge-Kutta matrix has a single reals-fold zero. In this paper we will construct such a family of methods and give some results concerning their maximum attainable order and stability properties. Some consideration is also given to showing how these methods can be efficiently implemented and, in particular, how local error estimates can be obtained by the use of embedding techniques.     

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.     

Trigonometrically fitted multi-step Runge-Kutta methods for solving oscillatory initial value problems

In this paper, trigonometrically fitted multi-step Runge-Kutta (TFMSRK) methods for the numerical integration of oscillatory initial value problems are proposed and studied. TFMSRK methods inherit the frame of multi-step Runge-Kutta (MSRK) methods and integrate exactly the problem whose solutions can be expressed as the linear combinations of functions from the set of \(\{\exp (\mathrm {i}wt),\exp (-\mathrm {i}wt)\},\) or equivalently the set \(\{\cos (wt),\sin (wt)\}\), where w represents an approximation of the main frequency of the problem. The general order conditions are given and four new explicit TFMSRK methods with order three and four, respectively, are constructed. Stability of the new methods is examined and the corresponding regions of stability are depicted. Numerical results show that our new methods are more efficient in comparison with other well-known high quality methods proposed in the scientific literature.     

High order algebraically stable Runge-Kutta methods

In Burrage and Butcher [3] the concept of Algebraic Stability was introduced in the study of Runge-Kutta methods. In this paper an analysis is made of the family ofs-stage Runge-Kutta methods of order 2s—2 or more which possesses this property.     

On implicit Runge-Kutta methods with high stage order

It is well known that high stage order is a desirable property for implicit Runge-Kutta methods. In this paper it is shown that it is always possible to construct ans-stage IRK method with a given stability function and stage orders−1 if the stability function is an approximation to the exponential function of at least orders. It is further indicated how to construct such methods as well as in which cases the constructed methods will be stiffly accurate.     

On projected Runge-Kutta methods for differential-algebraic equations

Ascher and Petzold recently introducedprojected Runge-Kutta methods for the numerical solution of semi-explicit differential-algebraic systems of index 2. Here it is shown that such a method can be regarded as the limiting case of a standard application of a Runge-Kutta method with a very small implicit Euler step added to it. This interpretation allows a direct derivation of the order conditions and of superconvergence results for the projected methods from known results for standard Runge-Kutta methods for index-2 differential-algebraic systems, and an extension to linearly implicit differential-algebraic systems.     

Essentially Symplectic Boundary Value Methods for Linear Hamiltonian Systems

1.IlltroductiollInmanyareasofphysics,mechanics,etc.,HamiltoniansystemsofODEsplayaveryimportantrole.Suchsystemshavethefollowinggeneralform:where,bydenotingwithOfandimthenullmatrixandtheidentitymatrixofordermarespectively,SimplepropertiesofthematrixJZmarethefollowingones:Inequation(1)AH(~,t)isthegradientofascalarfunctionH(y,t),usuallycalledHamiltonian.InthecasewhereH(y,t)=H(y),thenthevalueofthisfunctionremainsconstantalongt.hesollltion7/(t),t,hatis'*ReceivedFebruaryI3,1995.l)Worksupporte…     

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.     

General linear methods with external stages of different orders

The approach introduced recently by Albrecht to derive order conditions for Runge-Kutta formulas based on the theory of A-methods is also very powerful for the general linear methods. In this paper, using Albrecht's approach, we formulate the general theory of order conditions for a class of general linear methods where the components of the propagating vector of approximations to the solution have different orders. Using this theory we derive a class of diagonally implicit multistage integration methods (DIMSIMs) for which the global order is equal to the local order. We also derive a class of general linear methods with two nodal approximations of different orders which facilitate local error estimation. Our theory also applies to the class of two-step Runge-Kutta introduced recently by Jackiewicz and Tracogna.The work of the first author was supported by the National Science Foundation under grant NSF DMS-9208048. The work of the second author was supported by the Italian Consiglio Nazionale delle Richerche.     

On the relation between stability and contractivity

This paper concerns the rate of growth of numerical approximations obtained by one-step methods for solving linear stiff initial value problems. For some of these methods weak stability with respect to arbitrary norms is shown to be equivalent to contractivity. This kind of stability is also proved to entail a barrierp1 for the order of accuracyp within a broad class of methods, including general Runge-Kutta methods withm1 stages.Dedicated to Professor Germund Dahlquist on the occasion of his sixtieth birthday.     

Equilibria of Runge-Kutta methods

Summary It is known that certain Runge-Kutta methods share the property that, in a constant-step implementation, if a solution trajectory converges to a bounded limit then it must be a fixed point of the underlying differential system. Such methods are calledregular. In the present paper we provide a recursive test to check whether given method is regular. Moreover, by examining solution trajectories of linear equations, we prove that the order of ans-stage regular method may not exceed 2[(s+2)/2] and that the maximal order of regular Runge-Kutta method with an irreducible stability function is 4.     

Pseudo-symplectic Runge-Kutta methods

Apart from specific methods amenable to specific problems, symplectic Runge-Kutta methods are necessarily implicit. The aim of this paper is to construct explicit Runge-Kutta methods which mimic symplectic ones as far as the linear growth of the global error is concerned. Such method of orderp have to bepseudo-symplectic of pseudosymplecticness order2p, i.e. to preserve the symplectic form to within ⊗(h ^2p )-terms. Pseudo-symplecticness conditions are then derived and the effective construction of methods discussed. Finally, the performances of the new methods are illustrated on several test problems.