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.

Nonlinear Stability and Convergence of ERKN Integrators for Solving Nonlinear Multi-Frequency Highly Oscillatory Second-Order ODEs with Applications to Semi-Linear Wave Equations

AnEssential Extension of the Finite-Energy Condition for ExtendedRunge-Kutta-Nyström Integratorswhen Applied to Nonlinear Wave Equations

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

Effective Integrators for Nonlinear Second-Order Oscillatory Systems with a Time-Dependent Frequency Matrix

Error bounds for explicit ERKN integrators for systems of multi-frequency oscillatory second-order differential equations

ERKN Integrators Solving Multi-Frequency Highly Oscillatory Systems with Applications

Explicit Multi-Frequency Symmetric Extended RKN Integrators for Solving Multi-Frequency and Multidimensional Oscillatory Reversible Systems

ERKN integrators for systems of oscillatory second-order differential equations

Efficient Implementation of the ARKN and ERKN Integrators for Multi-Frequency Oscillatory Systems with Multiple Time Scales

ERKN methods for long-term integration of multidimensional orbital problems

Numerical Conservation of Energy, Momentum and Actions for Extended RKN Methods when Applied to Nonlinear Wave Equations Via Spatial Spectral Semi-Discretizations

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

Cost-reduction implicit exponential Runge-Kutta methods for highly oscillatory systems

Long-term Analysis of Symplectic or Symmetric Extended RKN Methods for Nonlinear Wave Equations

Error Analysis of Explicit TSERKN Methods for Highly Oscillatory Systems

Extended Phase Properties and Stability Analysis of RKN-type Integrators for Solving General Oscillatory Second-Order Initial Value Problems

Global Error Bounds of One-Stage Extended RKN Integrators for Semilinear Wave Equations

Symplectic Exponential Runge-Kutta Methods for Solving Nonlinear Hamiltonian Systems

Long-time Analysis of an Extended RKN Integrator for Hamiltonian Systems with a Solution-Dependent High Frequency

ENERGY-PRESERVING CONTINUOUS-STAGE EXPONENTIAL RUNGE--KUTTA INTEGRATORS FOR EFFICIENTLY SOLVING HAMILTONIAN SYSTEMS

Improved Uniform Error Bounds of Exponential Wave Integrator Method for Long-Time Dynamics of the Space Fractional Klein-Gordon Equation with Weak Nonlinearity