Local error estimation for multistep collocation methods

Multistep collocation methods for initial value problems in ordinary differential equations are known to be a subclass of multistep Runge-Kutta methods and a generalisation of the well-known class of one-step collocation methods as well as of the one-leg methods of Dahlquist. In this paper we derive an error estimation method of embedded type for multistep collocation methods based on perturbed multistep collocation methods. This parallels and generalizes the results for one-step collocation methods by Nørsett and Wanner. Simple numerical experiments show that this error estimator agrees well with a theoretical error estimate which is a generalisation of an error estimate first derived by Dahlquist for one-leg methods.     

A necessary condition forB-convergence of Runge-Kutta methods

Algebraic stability is a well-known property necessary forB-convergence of a Runge-Kutta method, provided its nodesc _i are such thatc _i – c _j wheneveri j. In this paper a slightly weaker property is established which becomes algebraic stability as soon asc _i – c _j , wheneveri j is excluded.Using this condition it is shown that the Lobatto IIIA methods with more than two stages cannot beB-convergent.This paper was written while the author was visiting the Rijksuniversiteit Leiden in the Netherlands, supported by the Netherlands Organization for Scientific Research (N.W.O.) and by an Erwin Schrödinger scholarship from the Fonds zur Förderung der wissenschaftlichen Forschung.     

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.     

Contractivity of Runge-Kutta methods

In this paper we present necessary and sufficient conditions for Runge-Kutta methods to be contractive. We consider not only unconditional contractivity for arbitrary dissipative initial value problems, but also conditional contractivity for initial value problems where the right hand side function satisfies a circle condition. Our results are relevant for arbitrary norms, in particular for the maximum norm.For contractive methods, we also focus on the question whether there exists a unique solution to the algebraic equations in each step. Further we show that contractive methods have a limited order of accuracy. Various optimal methods are presented, mainly of explicit type. We provide a numerical illustration to our theoretical results by applying the method of lines to a parabolic and a hyperbolic partial differential equation.Research supported by the Netherlands Organization for Scientific Research (N.W.O.) and the Royal Netherlands Academy of Arts and Sciences (K.N.A.W.)     

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.     

On the D-suitability of implicit Runge-Kutta methods

The concept of suitability means that the nonlinear equations to be solved in an implicit Runga-Kutta method have a unique solution. In this paper, we introduce the concept of D-suitability and show that previous results become special cases of ours. In addition, we also give some examples to illustrate the D-suitability of a matrixA.     

A class of pseudo Runge-Kutta methods

This paper develops a general theory for a class of Runge-Kutta methods which are based, in addition to the stages of the current step, also on the stages of the previous step. Such methods have been introduced previously for the case of one and two stages. We show that for any numbers of stages methods of orderp withs+1 p 2s can be constructed. The paper terminates with a study of step size change and stability.     

Superstable implicit Runge-Kutta methods for second order initial value problems

Using the well known properties of thes-stage implicit Runge-Kutta methods for first order differential equations, single step methods of arbitrary order can be obtained for the direct integration of the general second order initial value problemsy=f(x, y, y),y(x _o)=y _o,y(x _o)=y _o. These methods when applied to the test equationy+2y+ ² y=0, ,0, +>0, are superstable with the exception of a finite number of isolated values ofh. These methods can be successfully used for solving singular perturbation problems for which f/y and/or f/y are negative and large. Numerical results demonstrate the efficiency of these methods.     

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.     

Global error estimation and control in linearly-implicit parallel two-step peer W-methods

The class of linearly-implicit parallel two-step peer W-methods has been designed recently for efficient numerical solutions of stiff ordinary differential equations. Those schemes allow for parallelism across the method, that is an important feature for implementation on modern computational devices. Most importantly, all stage values of those methods possess the same properties in terms of stability and accuracy of numerical integration. This property results in the fact that no order reduction occurs when they are applied to very stiff problems. In this paper, we develop parallel local and global error estimation schemes that allow the numerical solution to be computed for a user-supplied accuracy requirement in automatic mode. An algorithm of such global error control and other technical particulars are also discussed here. Numerical examples confirm efficiency of the presented error estimation and stepsize control algorithm on a number of test problems with known exact solutions, including nonstiff, stiff, very stiff and large-scale differential equations. A comparison with the well-known stiff solver RODAS is also shown.     

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.     

A note onB-stability of Runge-Kutta methods

Summary Burrage and Butcher [1, 2] and Crouzeix [4] introduced for Runge-Kutta methods the concepts ofB-stability,BN-stability and algebraic stability. In this paper we prove that for any irreducible Runge-Kutta method these three stability concepts are equivalent.Chapters 1–3 of this article have been written by the second author, whereas chapter 4 has been written by the first author     

Some stability results for explicit Runge-Kutta methods

The theory of positive real functions is used to provide bounds for the largest possible disk to be inscribed in the stability region of an explicit Runge-Kutta method. In particular, we show that the closed disk |+r| r can be contained in the stability region of an explicitm-stage Runge-Kutta method of order two if and only ifr m – 1.     

Natural Runge-Kutta and projection methods

Summary Recently the author defined the class of natural Runge-Kutta methods and observed that it includes all the collocation methods. The present paper is devoted to a complete characterization of this class and it is shown that it coincides with the class of the projection methods in some polynomial spaces.This work was supported by the Italian Ministero della Pubblica Istruzione, funds 40%     

Perturbed collocation and Runge-Kutta methods

Summary It is well known thatsome implicit Runge-Kutta methods are equivalent to collocation methods. This fact permits very short and natural proofs of order andA, B, AN, BN-stability properties for this subclass of methods (see [9] and [10]). The present paper answers the natural question, ifall RK methods can be considered as a somewhat perturbed collocation. After having introduced this notion we give a proof on the order of the method and discuss their stability properties. Much of known theory becomes simple and beautiful.     

Stability andB-convergence of linearly implicit Runge-Kutta methods

Summary In this paper we study stability and convergence properties of linearly implicit Runge-Kutta methods applied to stiff semi-linear systems of differential equations. The stability analysis includes stability with respect to internal perturbations. All results presented in this paper are independent of the stiffness of the system.     

On the relation between algebraic stability andB-convergence for Runge-Kutta methods

Summary This paper is concerned with the numerical solution of stiff initial value problems for systems of ordinary differential equations using Runge-Kutta methods. For these and other methods Frank, Schneid and Ueberhuber [7] introduced the important concept ofB-convergence, i.e. convergence with error bounds only depending on the stepsizes, the smoothness of the exact solution and the so-called one-sided Lipschitz constant . Spijker [19] proved for the case <0 thatB-convergence follows from algebraic stability, the well-known criterion for contractivity (cf. [1, 2]). We show that the order ofB-convergence in this case is generally equal to the stage-order, improving by one half the order obtained in [19]. Further it is proved that algebraic stability is not only sufficient but also necessary forB-convergence.This study was completed while this author was visiting the Oxford University Computing Laboratory with a stipend from the Netherlands Organization for Scientific Research (N.W.O.)     

A note on pseudo-symplectic Runge-Kutta methods

Aubry and Chartier introduced (1998) the concept of pseudo-symplecticness in order to construct explicit Runge-Kutta methods, which mimic symplectic ones. Of particular interest are methods of order (p, 2p), i.e., of orderp and pseudo-symplecticness order 2p, for which the growth of the global error remains linear. The aim of this note is to show that the lower bound for the minimal number of stages can be achieved forp=4 andp=5.     

Asymptotically correct error estimation for collocation methods applied to singular boundary value problems

We discuss an a posteriori error estimate for collocation methods applied to boundary value problems in ordinary differential equations with a singularity of the first kind. As an extension of previous results we show the asymptotical correctness of our error estimate for the most general class of singular problems where the coefficient matrix is allowed to have eigenvalues with positive real parts. This requires a new representation of the global error for the numerical solution obtained by piecewise polynomial collocation when applied to our problem class.     

An algebraic characterization ofB-convergent Runge-Kutta methods

Summary In the analysis of discretization methods for stiff intial value problems, stability questions have received most part of the attention in the past.B-stability and the equivalent criterion algebraic stability are well known concepts for Runge-Kutta methods applied to dissipative problems. However, for the derivation ofB-convergence results — error bounds which are not affected by stiffness — it is not sufficient in many cases to requireB-stability alone. In this paper, necessary and sufficient conditions forB-convergence are determined.This paper was written while J. Schneid was visiting the Centre for Mathematics and Computer Science with an Erwin-Schrödinger stipend from the Fonds zur Förderung der wissenschaftlichen Forschung