On the stability properties of Brown's multistep multiderivative methods

Summary Brown introducedk-step methods usingl derivatives. We investigate for whichk andl the methods are stable or unstable. It is seen that to anyl the method becomes unstable fork large enough. All methods withk2(l+1) are stable. Fork=1,2,..., 18 there exists a _k such that the methods are stable for anyl _k and unstable for anyl < _k . The _k are given.     

On the integration of stiff systems of O.D.E.s using extended backward differentiation formulae

Summary A class of extended backward differentiation formulae suitable for the approximate numerical integration of stiff systems of first order ordinary differential equations is derived. An algorithm is described whereby the required solution is predicted using a conventional backward differentiation scheme and then corrected using an extended backward differentiation scheme of higher order. This approach allows us to developL-stable schemes of order up to 4 andL()-stable schemes of order up to 9. An algorithm based on the integration formulae derived in this paper is illustrated by some numerical examples and it is shown that it is often superior to certain existing algorithms.     

Partitioned adaptive Runge-Kutta methods and their stability

Summary This paper deals with the solution of partitioned systems of nonlinear stiff differential equations. Given a differential system, the user may specify some equations to be stiff and others to be nonstiff. For the numerical solution of such a system partitioned adaptive Runge-Kutta methods are studied. Nonstiff equations are integrated by an explicit Runge-Kutta method while an adaptive Runge-Kutta method is used for the stiff part of the system.The paper discusses numerical stability and contractivity as well as the implementation and usage of such compound methods. Test results for three partitioned stiff initial value problems for different tolerances are presented.     

On improving the absolute stability of local extrapolation

Summary A widely used technique for improving the accuracy of solutions of initial value problems in ordinary differential equations is local extrapolation. It is well known, however, that when using methods appropriate for solving stiff systems of ODES, the stability of the method can be seriously degraded if local extrapolation is employed. This is due to the fact that performing local extrapolation on a low order method is equivalent to using a higher order formula and this high order formula may not be suitable for solving stiff systems. In the present paper a general approach is proposed whereby the correction term added on in the process of local extrapolation is in a sense a rational, rather than a polynomial, function. This approach allows high order formulae with bounded growth functions to be developed. As an example we derive anA-stable rational correction algorithm based on the trapezoidal rule. This new algorithm is found to be efficient when low accuracy is requested (say a relative accuracy of about 1%) and its performance is compared with that of the more familiar Richardson extrapolation method on a large set of stiff test problems.     

On multi-grid methods for variational inequalities

Summary We consider here a general class of algorithms for the numerical solution of variational inequalities. A convergence proof is given and in particular a multi-grid method is described. Numerical results are presented for the finite-difference discretization of an obstacle problem for minimal surfaces     

A new interpolation procedure for adapting Runge-Kutta methods to delay differential equations

This paper deals with adapting Runge-Kutta methods to differential equations with a lagging argument. A new interpolation procedure is introduced which leads to numerical processes that satisfy an important asymptotic stability condition related to the class of testproblemsU(t)=U(t)+U(t–) with , C, Re()<–||, and >0. Ifc _i denotes theith abscissa of a given Runge-Kutta method, then in thenth stept _n–1t _n :=t _n–1+h of the numerical process our interpolation procedure computes an approximation toU(t _n–1+c _i h–) from approximations that have already been generated by the process at pointst _j–1+c _i h(j=1,2,3,...). For two of these new processes and a standard process we shall consider the convergence behaviour in an actual application to a given, stiff problem.     

On the order of composite multistep methods for ordinary differential equations

Summary In this paper, a general class ofk-step methods for the numerical solution of ordinary differential equations is discussed. It is shown that methods with order of consistencyq have order of convergence (q+1) if a very simple condition is satisfied. This result gives a new aspect to previous results of Spijker; it also serves as a starting point for a new theory of cyclick-step methods, completing an approach of Donelson and Hansen. It facilitates the practical determination of high-order cyclick-step methods, especially of stiffly stable,k-step methods.     

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

Asymptotic bounds on the errors of one-step methods

Summary Bulirsch and Stoer have shown how to construct asymptotic upper and lower bounds on the true (global) errors resulting from the solution by extrapolation of the initial value problem for a system of ordinary differential equations. It is shown here how to do this for any one-step method endowed with an asymptotically correct local error estimator. The one-step method can be changed at every step.This work performed at Sandia National Laboratories supported by the U.S. Department of Energy under contract number DE-AC04-76DP00789     

A 0-Stability and stiff stability of Brown's multistep multiderivative methods

Summary Brown [1] introducedk-step methods usingl derivatives. Necessary and sufficient conditions forA ₀-stability and stiff stability of these methods are given. These conditions are used to investigate for whichk andl the methods areA ₀-stable. It is seen that for allk andl withk1.5 (l+1) the methods areA ₀-stable and stiffly stable. This result is conservative and can be improved forl sufficiently large. For smallk andl A ₀-stability has been determined numerically by implementing the necessary and sufficient condition.     

On the numerical treatment of a small divisor problem

Summary Quasiperiodic solutions of perturbed integrable Hamiltonian equations such as weakly coupled harmonic oscillators can be found by constructing an appropriate coordinate transformation which leads to a small divisor problem. However the numerical difficulties are not merely caused by the small divisors but rather by the appearence of ghost solutions, which appear in any reasonable discretization of the problem. Our numerical treatment, based on a Newton-type iteration, guarantees an approximation of the relevant solution of the nonlinear problem. Numerical solutions are found up to a critical value of the coupling constant, which is much larger than the coupling constants allowed by the existence theory available so far.     

A-priori error estimates of Galerkin backward differentiation methods in time-inhomogeneous parabolic problems

Summary Backward differentiation methods up to orderk=5 are applied to solve linear ordinary and partial (parabolic) differential equations where in the second case the space variables are discretized by Galerkin procedures. Using a mean square norm over all considered time levels a-priori error estimates are derived. The emphasis of the results lies on the fact that the obtained error bounds do not depend on a Lipschitz constant and the dimension of the basic system of ordinary differential equations even though this system is allowed to have time-varying coefficients. It is therefore possible to use the bounds to estimate the error of systems with arbitrary varying dimension as they arise in the finite element regression of parabolic problems.     

Unconditionally stable explicit methods for parabolic equations

Summary This paper discussesrational Runge-Kutta methods for stiff differential equations of high dimensions. These methods are explicit and in addition do not require the computation or storage of the Jacobian. A stability analysis (based onn-dimensional linear equations) is given. A second orderA ₀-stable method with embedded error control is constructed and numerical results of stiff problems originating from linear and nonlinear parabolic equations are presented.     

On the stability of bilinear-constant velocity-pressure finite elements

Summary An analysis of the Babuka stability of bilinear/constant finite element pairs for viscous flow calculations is given. An unstable mode not of the checkerboard type is given for which the stability constant turns out to beO(h). Thus, the indicated spaces are not stable in general for numerical calculation.Work supported by U.S. Air Force Office of Scientific Research under grant AF-AFOSR-82-0213     

High order methods for the numerical integration of ordinary differential equations

Summary High order implicit integration formulae with a large region of absolute stability are developed for the approximate numerical integration of both stiff and non-stiff systems of ordinary differential equations. The algorithms derived behave essentially like one step methods and are demonstrated by direct application to certain particular examples.     

Stability properties of variable stepsize variable formula methods

Summary When variable stepsize variable formula methods (VSVFM's) are used in the solution of systems of first order differential equations instability arises sometimes. Therefore it is important to find VSVFM's whose zerostability properties are not affected by the choice of both the stepsize and the formula. The Adams VSVFM's are such methods. In this work a more general class of methods which contains the Adams VSVFM's is discussed and it is proved that the zero-stability of the class is not affected by the choice of the stepsize and of the formula.     

A study of Rosenbrock-type methods of high order

Summary This paper deals with the solution of nonlinear stiff ordinary differential equations. The methods derived here are of Rosenbrock-type. This has the advantage that they areA-stable (or stiffly stable) and nevertheless do not require the solution of nonlinear systems of equations. We derive methods of orders 5 and 6 which require one evaluation of the Jacobian and oneLU decomposition per step. We have written programs for these methods which use Richardson extrapolation for the step size control and give numerical results.     

The stability function for multistep collocation methods

Summary C-polynomials for rational approximation to the exponential function was introduced by Nørsett [7] to study stability properties of one-step methods. For one-step collocation methods theC-polynomial has a very simple form. In this paper we studyC-polynomials for multistep collocation methods and obtain results that generalize those in the one-step case, and provide a way to analyze linear stability of such methods.     

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.     

Higher order methods for a singularly perturbed problem

Summary A finite element method using piecewise polynomials of degree k is used to approximate the problem u+u=f, >0 a small parameter. A very irregular mesh is used. On this mesh error estimates of order0(h ^k+1) are obtained uniformly in ,h the maximum stepsize, fork2. The condition number of the system of linear equations one has to solve in order to get the approximation is estimated. Extension of the results to more complicated problems is briefly indicated. Finally, a numerical example is given.Work performed while visiting the IBM Thomas J. Watson Research Center, Yorktown Heights, N.Y.