Order results for Rosenbrock type methods on classes of stiff equations

Summary In this paper we give conditions for theB-convergence of Rosenbrock type methods when applied to stiff semi-linear systems. The convergence results are extended to stiff nonlinear systems in singular perturbation form. As a special case partitioned methods are considered. A third order method is constructed.Dedicated to the memory of Professor Lothar Collatz     

An inverse eigenvalue problem: ComputingB-stable Runge-Kutta methods having real poles

The implementation of implicit Runge-Kutta methods requires the solution of large sets of nonlinear equations. It is known that on serial machines these costs can be reduced if the stability function of ans-stage method has only ans-fold real pole. Here these so-called singly-implicit Runge-Kutta methods (SIRKs) are constructed utilizing a recent result on eigenvalue assignment by state feedback and a new tridiagonalization, which preserves the entries required by theW-transformation. These two algorithms in conjunction with an unconstrained minimization allow the numerical treatment of a difficult inverse eigenvalue problem. In particular we compute an 8-stage SIRK which is of order 8 andB-stable. This solves a problem posed by Hairer and Wanner a decade ago. Furthermore, we finds-stageB-stable SIRKs (s=6,8) of orders, which are evenL-stable.     

B-convergence properties of defect correction methods. I

Summary B-convergence properties of defect correction methods based on the implicit Euler and midpoint schemes are discussed. The property ofB-convergence means that there exist global error bounds for nonlinear stiff problems independent of their stiffness. It turns out that the orders ofB-convergence of these methods coincide with the conventional orders of convergence of these methods derived under the assumption that.hL is small (whereL is a Lipschitz constant of the right-hand side). In Part I these assertions are reduced to the validity of the so-called Hypothesis A which is discussed in greater detail in Part II. Numerical experiments confirming the theoretical analysis are also given in Part II.     

B-convergence properties of defect correction methods. II

Summary B-convergence properties of defect correction methods based on the implicit Euler and midpoint schemes are discussed. The property ofB-convergence means that there exist global error bounds for nonlinear stiff problems independent of their stiffness. It turns out that the orders ofB-convergence of these methods coincide with the conventional orders of small (whereL is a Lipschitz constant of the right-hand side). In Part I these assertions are reduced to the validity of the so-called HypothesisA which is discussed in greater detail in Part II. Numerical experiments confirming the theoretical analysis are also given in Part II.     

The order ofB-convergence of algebraically stable Runge-Kutta methods

In a previous paper it was shown that for a class of semi-linear problems many high order Runge-Kutta methods have order of optimalB-convergence one higher than the stage order. In this paper we show that for the more general class of nonlinear dissipative problems such as result holds only for a small class of Runge-Kutta methods and that such methods have at most classical order 3.     

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     

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

The asymptotic stability of one-parameter methods for neutral differential equations

This paper deals with the asymptotic stability of theoretical solutions and numerical methods for systems of neutral differential equationsx=Ax(t–)+Bx(t)+Cx(t–), whereA, B, andC are constant complexN ×N matrices, and >0. A necessary and sufficient condition such that the differential equations are asymptotically stable is derived. We also focus on the numerical stability properties of adaptations of one-parameter methods. Further, we investigate carefully the characterization of the stability region.     

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 the relations between B2V Ms and Runge-Kutta collocation methods

The principal aim of this paper is a rigorous analysis of the relations between Block Boundary Value Methods (B₂V Ms) with minimal blocksize defined over a suitable nonuniform finer mesh and well-known Runge-Kutta collocation methods. Moreover, a further aspect that will be briefly investigated is the construction of an extended finer mesh for building B₂V Ms with nonminimal blocksize. Some advantages that may arise from the use of the so-obtained methods will be also discussed.     

A criterion forP-stability properties of Runge-Kutta methods

P-stability is an analogous stability property toA-stability with respect to delay differential equations. It is defined by using a scalar test equation similar to the usual test equation ofA-stability. EveryP-stable method isA-stable, but anA-stable method is not necessarilyP-stable. We considerP-stability of Runge-Kutta (RK) methods and its variation which was originally introduced for multistep methods by Bickart, and derive a sufficient condition for an RK method to have the stability properties on the basis of an algebraic characterization ofA-stable RK methods recently obtained by Schere and Müller. By making use of the condition we clarify stability properties of some SIRK and SDIRK methods, which are easier to implement than fully implicit methods, applied to delay differential equations.     

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.     

NP-stability of Runge-Kutta methods based on classical quadrature

Recently Bellen, Jackiewicz and Zennaro have studied stability properties of Runge-Kutta (RK) methods for neutral delay differential equations using a scalar test equation. In particular, they have shown that everyA-stable collocation method isNP-stable, i.e., the method has an analogous stability property toA-stability with respect to the test equation. Consequently, the Gauss, Radau IIA and Lobatto IIIA methods areNP-stable. In this paper, we examine the stability of RK methods based on classical quadrature by a slightly different approach from theirs. As a result, we prove that the Radau IA and Lobatto IIIC methods equipped with suitable continuous extensions are alsoNP-stable by virtue of fundamental notions related to those methods such as simplifying conditions, algebraic stability, and theW-transformation.     

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.     

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.     

B-Convergence of Lobatto IIIC formulas

Summary A completion ofB-convergence results of Lobatto IIIC schemes is presented. In particular, it is shown that Lobatto IIIC schemes with more than two stages areB-convergent when applied to IVPs with a negative one-sided Lipschitz constantm; they are notB-convergent, however, for IVPs with a non-negativem.     

Error growth analysis via stability regions for discretizations of initial value problems

This paper deals with numerical methods for the solution of linear initial value problems. Two main theorems are presented on the stability of these methods. Both theorems give conditions guaranteeing a mild error growth, for one-step methods characterized by a rational function ϕ(z). The conditions are related to the stability regionS={z:z∈ℂ with |ϕ(z)|≤1}, and can be viewed as variants to the resolvent condition occurring in the reputed Kreiss matrix theorem. Stability estimates are presented in terms of the number of time stepsn and the dimensions of the space. The first theorem gives a stability estimate which implies that errors in the numerical process cannot grow faster than linearly withs orn. It improves previous results in the literature where various restrictions were imposed onS and ϕ(z), including ϕ′(z)≠0 forz∈σS andS be bounded. The new theorem is not subject to any of these restrictions. The second theorem gives a sharper stability result under additional assumptions regarding the differential equation. This result implies that errors cannot grow faster thann ^β, with fixed β<1. The theory is illustrated in the numerical solution of an initial-boundary value problem for a partial differential equation, where the error growth is measured in the maximum norm.     

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     

Sur laB-stabilité des méthodes de Runge-Kutta

Résumé Dans cet article, nous modifions légèrement la définition de laB-stabilité donnée par J.C. Butcher [1] afin qu'elle s'applique à une plus large classe d'équations différentielles et nous donnons des caractérisations simples de cette propriété.

---

OnB-stability of the methods of Runge Kutta

Summary In this paper, we slightly modify the definition ofB-stability of Butcher [1], so as to cover a wider class of differential equations, and we give simple characterizations of this property.