Non-smooth data error estimates for linearly implicit Runge-Kutta methods

Linearly implicit time discretizations of semilinear parabolicequations with non-smooth initial data are studied. The analysisuses the framework of analytic semigroups which includes reaction-diffusionequations and the incompressible Navier-Stokes equations. Itis shown that the order of convergence on finite time intervalsis essentially one. Applications to the long-term behaviourof linearly implicit Runge-Kutta methods are given.     

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.     

On the implementation of implicit Runge-Kutta methods

The modified Newton iterations in the implementation of ans stage implicit Runge-Kutta method for ann dimensional differential equation system require 2s ³ n ³/3+O(n ²) operations for theLU factorisations and 2s ² n ²+O(n) operations for the back substitutions. This paper describes a method for transforming the linear system so as to reduce these operation counts.     

On global error control in nested implicit Runge-Kutta methods of the gauss type

Automatic global error control based on a combined control of step size and order presented by Kulikov and Khrustaleva in 2008 is investigated. Special attention is given to the efficiency of computation, because the implicit extrapolation based on multistage implicit Runge-Kutta schemes may be expensive. Specifically, we discuss a technique of global error estimation and control in order to compute a numerical solution satisfying the user-supplied accuracy conditions (in exact arithmetic) automatically. The theoretical results of this paper are confirmed by numerical experiments on test problems.     

On implicit Runge-Kutta methods with higher derivatives

Two families of implicit Runge-Kutta methods with higher derivatives are (re-)considered generalizing classical Runge-Kutta methods of Butcher type and f Ehle type. For generalized Butcher methods the characteristic functionG() is represented by means of the node polynomial directly, thereby showing that in methods of maximum order,G() is connected withs-orthogonal polynomials in exactly the same way as Padé approximations in the classical case.     

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.     

OnA-stable implicit Runge-Kutta methods

An elementary proof is given of theA-stability of implicit Runge-Kutta methods for which the corresponding rational function is on the diagonal or one of the first two subdiagonals of the Padé table for the exponential function. The result is extended to give necessary and sufficient conditions for theA-stability ofn-stage methods of order greater than or equal to 2n–2.     

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.     

Diagonally implicit Runge-Kutta methods for stiff problems

Diagonally implicit Runge-Kutta methods are examined. It is shown that, for stiff problems, the methods based on the minimization of certain error functions have advantages over other methods; these functions are determined in terms of the errors for simplest model equations. Methods of orders three, four, five, and six are considered.     

Implications of order reduction for implicit Runge-Kutta methods

Stability analysis of Runge-Kutta (RK) formulas was originally limited to linear ordinary differential equations (ODEs). More recently such analysis has been extended to include the behaviour of solutions to nonlinear problems. This extension led to additional stability requirements for RK methods. Although the class of problems has been widened, the analysis is still restricted to a fixed stepsize. In the case of differential algebraic equations (DAEs), additional order conditions must be satisfied [6] to achieve full classical ODE order and avoid possible order reduction. In this case too, a fixed stepsize analysis is employed. Such analysis may be of only limited use in quantifying the effectiveness of adaptive methods on stiff problems.In this paper we examine the phenomenon of order reduction and its implications on variable-step algorithms. We introduce a global measure of order referred to here as the observed order which is based on the average stepsize over the region of integration. This measure may be better suited to the study of stiff systems, where the stepsize selection algorithm will vary the stepsize considerably over the interval of integration. Observed order gives a better indication of the relationship between accuracy and cost. Using this measure, the observed order reduction will be seen to be less severe than that predicated by fixed stepsize order analysis.Supported by the Information Technology Research Centre of Ontario, and the Natural Science and Engineering Research Council of Canada.     

Optimal order diagonally implicit Runge-Kutta methods

In this paper, the optimal order of non-confluent Diagonally Implicit Runge-Kutta (DIRK) methods with non-zero weights is examined. It is shown that the order of aq-stage non-confluent DIRK method with non-zero weights cannot exceedq+1. In particular the optimal order of aq stage non-confluent DIRK method with non-zero weights isq+1 for 1q5. DIRK methods of orders five and six in four and five stages respectively are constructed. It is further shown that the optimal order of a non-confluentq stage DIRK method with non-zero weights isq, forq6.     

An efficient scheme for the implementation of implicit Runge-Kutta methods

A scheme is proposed for solving nonlinear algebraic equations arising in the implementation of the implicit Runge-Kutta methods. In contrast to the available schemes, not only the starting values of the variables but also those of the derivatives are predicted. This makes it possible to reduce the number of evaluations of the function (the right-hand side) at each implicit stage without significantly reducing the accuracy of integration.     

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.     

A stability property of implicit Runge-Kutta methods

A class of implicit Runge-Kutta methods is shown to possess a stability property which is a natural extension of the notion ofA-stability for non-linear systems.     

On the iterative solution of the algebraic equations in fully implicit Runge-Kutta methods

This paper deals with the iterative solution of stage equations which arise when some fully implicit Runge-Kutta methods, in particular those based on Gauss, Radau and Lobatto points, are applied to stiff ordinary differential equations. The error behaviour in the iterates generated by Newton-type and, particularly, by single-Newton schemes which are proposed for the solution of stage equations is studied. We consider stiff systems y'(t) = f(t,y(t)) which are dissipative with respect to a scalar product and satisfy a condition on the relative variation of the Jacobian of f(t,y) with respect to y, similar to the condition considered by van Dorsselaer and Spijker in [7] and [17]. We prove new convergence results for the single-Newton iteration and derive estimates of the iteration error that are independent of the stiffness. Finally, some numerical experiments which confirm the theoretical results are presented. This revised version was published online in June 2006 with corrections to the Cover Date.     

On the solvability of the systems of equations arising in implicit Runge-Kutta methods

In this paper we present a new condition under which the systems of equations arising in the application of an implicit Runge-Kutta method to a stiff initial value problem, has unique solutions. We show that our condition is weaker than related conditions presented previously. It is proved that the Lobatto IIIC methods fulfil the new condition.     

The order of some implicit Runge-Kutta methods

The orders of some single step methods for the solution of a general system of differential equations are established. Leading error terms are given.     

On theA-stability of implicit Runge-Kutta processes

A new technique to calculate the characteristic functions and to examine theA-stability of implicit Runge-Kutta processes is presented. This technique is based on a direct algebraic approach and an application of theC-polynomial theory of Nørsett. New processes are suggested. These processes can be exponentially fitted in anA-stable manner.     

Monotonicity and boundedness in implicit Runge-Kutta methods

Summary This paper concerns the analysis of implicit Runge-Kutta methods for approximating the solutions to stiff initial value problems. The analysis includes the case of (nonlinear) systems of differential equations that are essentially more general than the classical test equationU=U (with a complex constant). The properties of monotonicity and boundedness of a method refer to specific moderate rates of growth of the approximations during the numerical calculations. This paper provides necessary conditions for these properties by using the important concept of algebraic stability (introduced by Burrage, Butcher and by Crouzeix). These properties will also be related to the concept of contractivity (B-stability) and to a weakened version of contractivity.     

On linearly implicit IMEX Runge-Kutta methods for degenerate convection-diffusion problems modeling polydisperse sedimentation

Implicit-explicit (IMEX) Runge-Kutta (RK) methods are suitable for the solution of nonlinear, possibly strongly degenerate, convection-diffusion problems, since the stability restrictions, coming from the explicitly treated convective part, are much less severe than those that would be deduced from an explicit treatment of the diffusive term. A particularly efficient variant of these schemes, so-called linearly implicit IMEX-RK schemes, arise from discretizing the diffusion terms in a way that more carefully distinguishes between stiff and nonstiff dependence, such that in each time step only a linear system needs to be solved. These schemes provide an efficient tool for the numerical exploration of sediment formation and composition under a strongly degenerate polydisperse sedimentation model.