Abstract:For systems of oscillatory second-order differential equations y″+My=f with M∈Rm×m, a symmetric positive semi-definite matrix, X. Wu et al. have proposed the multidimensional ARKN methods [X. Wu, X. You, J. Xia, Order conditions for ARKN methods solving oscillatory systems, Comput. Phys. Comm. 180 (2009) 2250–2257], which are an essential generalization of J.M. Franco's ARKN methods for one-dimensional problems or for systems with a diagonal matrix M=w2I [J.M. Franco, Runge–Kutta–Nyström methods adapted to the numerical integration of perturbed oscillators, Comput. Phys. Comm. 147 (2002) 770–787]. One of the merits of these methods is that they integrate exactly the unperturbed oscillators y″+My=0. Regretfully, even for the unperturbed oscillators the internal stages Yi of an ARKN method fail to equal the values of the exact solution y(t) at tn+cih, respectively. Recently H. Yang et al. proposed the ERKN methods to overcome this drawback [H.L. Yang, X.Y. Wu, Xiong You, Yonglei Fang, Extended RKN-type methods for numerical integration of perturbed oscillators, Comput. Phys. Comm. 180 (2009) 1777–1794]. However, the ERKN methods in that paper are only considered for the special case where M is a diagonal matrix with nonnegative entries. The purpose of this paper is to extend the ERKN methods to the general case with M∈Rm×m, and the perturbing function f depends only on y. Numerical experiments accompanied demonstrates that the ERKN methods are more efficient than the existing methods for the computation of oscillatory systems. In particular, if M∈Rm×m is a symmetric positive semi-definite matrix, it is highly important for the new ERKN integrators to show the energy conservation in the numerical experiments for problems with Hamiltonian H(p,q)=12pTp+12qTMq+V(q) in comparison with the well-known methods in the scientific literature. Those so called separable Hamiltonians arise in many areas of physical sciences, e.g., macromolecular dynamics, astronomy, and classical mechanics.

Novel Improved Multidimensional Störmer–Verlet Formulas with Applications to Four Aspects in Scientific Computation

High-oder Symplectic FDTD Scheme for Solving Time-Dependent Schrdinger Equation

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

ERKN methods for long-term integration of multidimensional orbital problems

Three kinds of novel multi-symplectic methods for stochastic Hamiltonian partial differential equations

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.

A Symplectic Approximation with Nonlinear Stability and Convergence Analysis for Efficiently Solving Semi-Linear Klein–Gordon Equations

Symplectic Exponential Runge-Kutta Methods for Solving Nonlinear Hamiltonian Systems

A highly accurate explicit symplectic ERKN method for multi-frequency and multidimensional oscillatory Hamiltonian systems

High-Order Symplectic and Symmetric Composition Methods for Multi-Frequency and Multi-Dimensional Oscillatory Hamiltonian Systems

An Efficient High-Order Explicit Scheme for Solving Hamiltonian Nonlinear Wave Equations

On Extended RKN Integrators for Multidimensional Perturbed Oscillators with Applications

On Energy Stable Runge-Kutta Methods for the Water Wave Equation and Its Simplified Non-Local Hyperbolic Model

Embedded Exponential Runge–Kutta–Nyström Methods for Highly Oscillatory Hamiltonian Systems

Two Second-Order and Linear Numerical Schemes for the Multi-Dimensional Nonlinear Time-Fractional Schrödinger Equation

ERKN integrators for systems of oscillatory second-order differential equations

Explicit Multi-Symplectic Extended Leap-Frog Methods for Hamiltonian Wave Equations

A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations: I. Scheme and Theoretical Analysis

An improved Störmer-Verlet method based on exact discretization for nonlinear oscillators

Explicit Symplectic Multidimensional Exponential Fitting Modified Runge-Kutta-Nyström Methods

Multi-symplectic exponentially fitted RKN methods for wave equations