For Hamiltonian systems with non-canonical structure matrix a new class of numerical integrators is proposed. The methods exactly preserve energy, are invariant with respect to linear transformations, and have arbitrarily high order. Those of optimal order also preserve quadratic Casimir functions. The discussion of the order is based on an interpretation as partitioned Runge–Kutta method with ...

Article history: Received 27 April 2010 Received in revised form 2 August 2010 Accepted 19 August 2010 Available online 26 August 2010

In this work, we derive the order conditions for fourth order time splitting schemes in the case of the 1D Vlasov-Poisson system. Computations to obtain such conditions are motivated by the specific Poisson structure of the Vlasov-Poisson system : this structure is similar to Runge-Kutta-Nyström systems. The obtained conditions are proved to be the same as RKN conditions derived for ODE up to t...

Spectral deferred correction is a flexible technique for constructing high-order, stiffly-stable time integrators using a low order method as a base scheme. Here we examine their use in conjunction with splitting methods to solve initial-boundary value problems for partial differential equations. We exploit their close connection with implicit Runge–Kutta methods to prove that up to the full ac...

In this paper we analyze the consistency and stability properties of Runge-Kutta discrete adjoints. Discrete adjoints are very popular in optimization and control since they can be constructed automatically by reverse mode automatic differentiation. The consistency analysis uses the concept of elementary differentials and reveals that the discrete Runge-Kutta adjoint method has the same order o...

A connection between the algebra of rooted trees used in renormalization theory and Runge-Kutta methods is pointed out. Butcher’s group and B-series are shown to provide a suitable framework for renormalizing a toy model of field theory, following Kreimer’s approach. Finally B-series are used to solve a class of non-linear partial differential equations.

Recently, the Runge-Kutta methods, obtained via Taylor’s expansion is exist in the literature. In this study, we have derived explicit methods for problems of the form y′ = f(y) including second and third derivatives , by considering available Two-Derivative Runge-Kutta methods (TDRK). The methods use one evaluation of first derivative, one evaluation of second derivative and many evaluations o...

The aim of this paper is to design a new family of numerical methods of arbitrarily high order for systems of rst-order diierential equations which are to be termed pseudo two-step Runge-Kutta methods. By using collocation techniques, we can obtain an arbitrarily high-order stable pseudo two-step Runge-Kutta method with any desired number of implicit stages in retaining the two-step nature. In ...

Abstract. Space discretization of some time-dependent partial differential equations gives rise to stiff systems of ordinary differential equations. In this case, implicit methods should be used and therefore, in general, nonlinear systems must be solved. The solutions to these systems are approximated by iterative schemes and, in order to obtain an efficient code, good initializers should be u...