On the derivation of explicit two-step peer methods

The so-called two-step peer methods for the numerical solution of Initial Value Problems (IVP) in differential systems were introduced by R. Weiner, B.A. Schmitt and coworkers as a tool to solve different types of IVPs either in sequential or parallel computers. These methods combine the advantages of Runge-Kutta (RK) and multistep methods to obtain high stage order and therefore provide in a natural way a dense output. In particular, several explicit peer methods have been proved to be competitive with standard RK methods in a wide selection of non-stiff test problems.The aim of this paper is to propose an alternative procedure to construct families of explicit two step peer methods in which the available parameters appear in a transparent way. This allows us to obtain families of fixed stepsize s stage methods with stage order 2s−1, which provide dense output without extra cost, depending on some free parameters that can be selected taking into account the stability properties and leading error terms. A study of the extension of these methods to variable stepsize is also carried out. Optimal s stage methods with s=2,3 are derived.     

Stability and convergence of boundary value methods for solving ODE †

we review some recent approaches to the numerical solution of ODE. They are based on the solution of IVP (Initial Value Problem) by means of suitable BVM (Boundary Value Methods). A discussion of their properties of stability and covergence is presented along with their practical construction. The advantages of such methods with respect to IVP methods consist in a higher accuracy in the case of stiff problems and in a more covenient implementation on parallel computers.     

A family of A-stable Runge Kutta collocation methods of higher order for initial-value problems

^{We consider the construction of a special family of Runge–Kutta(RK) collocation methods based on intra-step nodal points ofChebyshev–Gauss–Lobatto type, with A-stability andstiffly accurate characteristics. This feature with its inherentimplicitness makes them suitable for solving stiff initial-valueproblems. In fact, the two simplest cases consist in the well-knowntrapezoidal rule and the fourth-order Runge–Kutta–LobattoIIIA method. We will present here the coefficients up to eighthorder, but we provide the formulas to obtain methods of higherorder. When the number of stages is odd, we have considereda new strategy for changing the step size based on the use ofa pair of methods: the given RK method and a linear multistepone. Some numerical experiments are considered in order to checkthe behaviour of the methods when applied to a variety of initial-valueproblems.}     

Analysis of trigonometric implicit Runge–Kutta methods

Using generalized collocation techniques based on fitting functions that are trigonometric (rather than algebraic as in classical integrators), we develop a new class of multistage, one-step, variable stepsize, and variable coefficients implicit Runge–Kutta methods to solve oscillatory ODE problems. The coefficients of the methods are functions of the frequency and the stepsize. We refer to this class as trigonometric implicit Runge–Kutta (TIRK) methods. They integrate an equation exactly if its solution is a trigonometric polynomial with a known frequency. We characterize the order and A-stability of the methods and establish results similar to that of classical algebraic collocation RK methods.     

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.     

An energy-preserving exponentially-fitted continuous stage Runge–Kutta method for Hamiltonian systems

Recently, the symplectic exponentially-fitted methods for Hamiltonian systems with periodic or oscillatory solutions have been attracting a lot of interest. As an alternative to them, in this paper, we propose a class of energy-preserving exponentially-fitted methods. For this aim, we show sufficient conditions for energy-preservation in terms of the coefficients of continuous stage Runge–Kutta (RK) methods, and extend the theory of exponentially-fitted RK methods in the context of continuous stage RK methods. Then by combining these two theories, we derive second and fourth order energy-preserving exponentially-fitted schemes.     

Two-Step Runge-Kutta: Theory and Practice

Local and global error for Two-Step Runge-Kutta (TSRK) methods are analyzed using the theory of B-series. Global error bounds are derived in both constant and variable stepsize environments. An embedded TSRK pair is constructed and compared with the RK5(4)6M pair of Dormand and Prince on the DETEST set of problems. Numerical results show that the TSRK performs competitively with the RK method.     

On the stability of functionally fitted Runge–Kutta methods

Classical collocation RK methods are polynomially fitted in the sense that they integrate an ODE problem exactly if its solution is an algebraic polynomial up to some degree. Functionally fitted RK (FRK) methods are collocation techniques that generalize this principle to solve an ODE problem exactly if its solution is a linear combination of a chosen set of arbitrary basis functions. Given for example a periodic or oscillatory ODE problem with a known frequency, it might be advantageous to tune a trigonometric FRK method targeted at such a problem. However, FRK methods lead to variable coefficients that depend on the parameters of the problem, the time, the stepsize, and the basis functions in a non-trivial manner that inhibits any in-depth analysis of the behavior of the methods in general. We present the class of so-called separable basis functions and show how to characterize the stability function of the methods in this particular class. We illustrate this explicitly with an example and we provide further insight for separable methods with symmetric collocation points. AMS subject classification (2000) 65L05, 65L06, 65L20, 65L60     

Runge–Kutta methods adapted to the numerical integration of oscillatory problems

New Runge–Kutta methods specially adapted to the numerical integration of IVPs with oscillatory solutions are obtained. The coefficients of these methods are frequency-dependent such that certain particular oscillatory solutions are computed exactly (without truncation errors). Based on the B-series theory and on the rooted trees we derive the necessary and sufficient order conditions for this class of RK methods. With the help of these order conditions we construct explicit methods (up to order 4) as well as pairs of embedded RK methods of orders 4 and 3. Some numerical examples show the excellent behaviour when they compete with classical RK methods.     

An efficient unified approach for the numerical solution of delay differential equations

In this paper we propose a new framework for designing a delay differential equation (DDE) solver which works with any supplied initial value problem (IVP) solver that is based on a standard step-by-step approach, such as Runge-Kutta or linear multi-step methods, and can provide dense output. This is done by treating a general DDE as a special example of a discontinuous IVP. Using this interpretation we develop an efficient technique to solve the resulting discontinuous IVP. We also give a more clear process for the numerical techniques used when solving the implicit equations that arise on a time step, such as when the underlying IVP solver is implicit or the delay vanishes. The new modular design for the resulting simulator we introduce, helps to accelerate the utilization of advances in the different components of an effective numerical method. Such components include the underlying discrete formula, the interpolant for dense output, the strategy for handling discontinuities and the iteration scheme for solving any implicit equations that arise.     

Approximate preservation of quadratic first integrals by explicit Runge–Kutta methods

The approximate preservation of quadratic first integrals (QFIs) of differential systems in the numerical integration with Runge–Kutta (RK) methods is studied. Conditions on the coefficients of the RK method to preserve all QFIs up to a given order are obtained, showing that the pseudo-symplectic methods studied by Aubry and Chartier (BIT 98(3):439–461, 1998) of algebraic order p preserve QFIs with order q = 2p. An expression of the error of conservation of QFIs by a RK method is given, and a new explicit six-stage formula with classical order four and seventh order of QFI-conservation is obtained by choosing their coefficients so that they minimize both local truncation and conservation errors. Several formulas with algebraic orders 3 and 4 and different orders of conservation have been tested with some problems with quadratic and general first integrals. It is shown that the new fourth-order explicit method preserves much better the qualitative properties of the flow than the standard fourth-order RK method at the price of two extra function evaluations per step and it is a practical and efficient alternative to the fully implicit methods required for a complete preservation of QFIs.     

Stability of Runge-Kutta methods for the generalized pantograph equation

Summary. This paper deals with stability properties of Runge-Kutta (RK) methods applied to a non-autonomous delay differential equation (DDE) with a constant delay which is obtained from the so-called generalized pantograph equation, an autonomous DDE with a variable delay by a change of the independent variable. It is shown that in the case where the RK matrix is regular stability properties of the RK method for the DDE are derived from those for a difference equation, which are examined by similar techniques to those in the case of autonomous DDEs with a constant delay. As a result, it is shown that some RK methods based on classical quadrature have a superior stability property with respect to the generalized pantograph equation. Stability of algebraically stable natural RK methods is also considered. Received May 5, 1998 / Revised version received November 17, 1998 / Published online September 24, 1999     

CHARACTERIZATIONS OF SYMMETRIC MULTISTEP RUNGE-KUTTA METHODS

Some characterizations for symmetric multistep Runge-Kutta(RK) methods are obtained. Symmetric two-step RK methods with one and two-stages are presented. Numerical examples show that symmetry of multistep RK methods alone is not sufficient for long time integration for reversible Hamiltonian systems. This is an important difference between one-step and multistep symmetric RK methods.     

Runge–Kutta interpolants for high precision computations

Runge–Kutta (RK) pairs furnish approximations of the solution of an initial value problem at discrete points in the interval of integration. Many techniques for enriching these methods with continuous approximations have been proposed. Here we construct C ¹ continuous, eighth and ninth order interpolation methods for a recently appeared RK pair of orders 9(8). These interpolants share a very small leading truncation error making them suitable for use at quadruple precision, i.e. 32–33 decimal digits of accuracy. Extended numerical results justify our effort.     

Existence and uniqueness of solutions of differential equations with respect to non-additive measures

By taking Sugeno-derivative into account, first, we investigate the existence of solutions to the initial value problems (IVP) of first-order differential equations with respect to non-additive measure (more precisely, distorted Lebesgue measure). It particularly occurs in the mathematical modeling of biology. We begin by expressing the differential equation in terms of ordinary derivative and the derivative with respect to the distorted Lebesgue measure. Then, by using the fixed point theorem on cones, we construct an operator and prove the existence of positive non-decreasing solutions on cones in semi-order Banach spaces. In addition, we also use Picard–Lindelöf theorem to prove the existence and uniqueness of the solution of the equation. Second, we investigate the existence of a solution to the boundary value problem (BVP) on cones with integral boundary conditions of a mix-order differential equation with respect to non-additive measures. Moreover, the Krasnoselskii fixed point theorem is also applied to both BVP and IVP and obtains at least one positive non-decreasing solution. Examples with graphs are provided to validate the results.     

Symplectic RK Methods and Symplectic PRK Methods with Real Eigenvalues

Properties of symplectic Runge-Kutta (RK) methods and symplectic partitioned Runge-Kutta (PRK) methods with real eigenvalues are discussed in this paper. It is shown that an s stage such method can‘t reach order more than s 1. Particularly, we prove that no symplectic RK method with real eigenvalues exists in stage s of order s 1 when s is even. But an example constructed by using the W-transformation shows that PRK method of this type does not necessarily meet this order barrier. Another useful way other than W-transformation to construct symplectic PRK method with real eigenvalues is then presented. Finally, a class of efficient symplectic methods is recommended.     

Towards optimal Runge-Kutta smoothers for multigrid methods for unsteady flow problems

We consider Runge-Kutta smoothers in a dual time stepping multigrid method for unsteady flow problems. These smoothers are easily parallelizable and Jacobian-free, making them very attractive for 3D calculations. Existing methods have been designed for steady flows, leading to slow convergence for unsteady problems. Here we determine the free parameters of the smoother to provide optimal damping for high frequency components for the unsteady linear advection equation. This is compared with an RK smoother designed for steady state problems, as commonly used in CFD codes. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new class of efficient one-step contractivity preserving high-order time discretization methods of order 5 to 14

A family of one-step, explicit, contractivity preserving, multi-stage, multi-derivative, Hermite–Birkhoff–Taylor methods of order p =?5,6,…,14, that we denote by CPHBT_RK4(d,s,p), with nonnegative coefficients are constructed by casting s-stage Runge–Kutta methods of order 4 with Taylor methods of order d. The constructed CPHBT_RK4 methods are implemented using efficient variable step control and are compared to other well-known methods on a variety of initial value problems. A selected method: CP 6-stages 9-derivative HBT method of order 12, denoted by CPHBTRK412, has larger region of absolute stability than Dormand–Prince DP(8,7)13M and Taylor method T(12) of order 12. It is superior to DP(8,7)13M and T(12) methods on the basis the number of steps, CPU time, and maximum global error on several problems often used to test higher-order ODE solvers. Also, we show that the contractivity preserving property of CPHBTRK412is very efficient in suppressing the effect of the propagation of discretization errors and the new method compares positively with explicit 17 stages Runge-Kutta-Nyström pair of order 12 by Sharp et al. on a long-term integration of a standard N-body problem. The selected CPHBTRK412is listed in the Appendix.     

An embedded pair of exponentially fitted explicit Runge–Kutta methods

An embedded pair of exponentially fitted explicit Runge–Kutta (RK) methods for the numerical integration of IVPs with oscillatory solutions is derived. This pair is based on the exponentially fitted explicit RK method constructed in Vanden Berghe et al., and we confirm that the methods which constitute the pair have algebraic order 4 and 3. Some numerical experiments show the efficiency of our pair when it is compared with the variable step code proposed by Vanden Berghe et al. (J. Comput. Appl. Math. 125 (2000) 107).     

Convergence of a class of runge-kutta methods for differential-algebraic systems of index 2

This paper deals with convergence results for a special class of Runge-Kutta (RK) methods as applied to differential-algebraic equations (DAE's) of index 2 in Hessenberg form. The considered methods are stiffly accurate, with a singular RK matrix whose first row vanishes, but which possesses a nonsingular submatrix. Under certain hypotheses, global superconvergence for the differential components is shown, so that a conjecture related to the Lobatto IIIA schemes is proved. Extensions of the presented results to projected RK methods are discussed. Some numerical examples in line with the theoretical results are included.