Explicit symplectic multidimensional exponential fitting modified Runge-Kutta-Nyström methods

This paper is concerned with multidimensional exponential fitting modified Runge-Kutta-Nyström (MEFMRKN) methods for the system of oscillatory second-order differential equations q″(t)+Mq(t)=f(q(t)), where M is a d×d symmetric and positive semi-definite matrix and f(q) is the negative gradient of a potential scalar U(q). We formulate MEFMRKN methods and show clearly the relationship between MEFMRKN methods and multidimensional extended Runge-Kutta-Nyström (ERKN) methods proposed by Wu et al. (Comput. Phys. Comm. 181:1955–1962, 2010). Taking into account the fact that the oscillatory system is a separable Hamiltonian system with Hamiltonian $H(p,q)=\frac{1}{2}p^{T}p+ \frac{1}{2}q^{T}Mq+U(q)$, we derive the symplecticity conditions for the MEFMRKN methods. Two explicit symplectic MEFMRKN methods are proposed. Numerical experiments accompanied demonstrate that our explicit symplectic MEFMRKN methods are more efficient than some well-known numerical methods appeared in the scientific literature.

[1]  Bin Wang,et al.  Multidimensional adapted Runge-Kutta-Nyström methods for oscillatory systems , 2010, Comput. Phys. Commun..

[2]  M. J,et al.  RUNGE-KUTTA SCHEMES FOR HAMILTONIAN SYSTEMS , 2005 .

[3]  R. Ruth.  A Can0nical Integrati0n Technique , 1983, IEEE Transactions on Nuclear Science.

[4]  Hans Van de Vyver,et al.  A symplectic exponentially fitted modified Runge–Kutta–Nyström method for the numerical integration of orbital problems , 2005 .

[5]  Ernst Hairer,et al.  Long-Time Energy Conservation of Numerical Methods for Oscillatory Differential Equations , 2000, SIAM J. Numer. Anal..

[6]  Jesús María Sanz-Serna,et al.  Mollified Impulse Methods for Highly Oscillatory Differential Equations , 2008, SIAM J. Numer. Anal..

[7]  Manuel Calvo,et al.  Symmetric and symplectic exponentially fitted Runge-Kutta methods of high order , 2010, Comput. Phys. Commun..

[8]  E. Hairer,et al.  Numerical Energy Conservation for Multi-Frequency Oscillatory Differential Equations , 2005 .

[9]  J. M. Franco.  Exponentially fitted explicit Runge-Kutta-Nyström methods , 2004 .

[10]  Ernst Hairer,et al.  Explicit, Time Reversible, Adaptive Step Size Control , 2005, SIAM J. Sci. Comput..

[11]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[12]  Marlis Hochbruck,et al.  A Gautschi-type method for oscillatory second-order differential equations , 1999, Numerische Mathematik.

[13]  Jianlin Xia,et al.  Order conditions for ARKN methods solving oscillatory systems , 2009, Comput. Phys. Commun..

[14]  R. D. Vogelaere,et al.  Methods of Integration which Preserve the Contact Transformation Property of the Hamilton Equations , 1956 .

[15]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[16]  T. E. Simos,et al.  Exponentially fitted symplectic integrator. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[18]  J. M. Franco.  New methods for oscillatory systems based on ARKN methods , 2006 .

[19]  Xinyuan Wu,et al.  Note on derivation of order conditions for ARKN methods for perturbed oscillators , 2009, Comput. Phys. Commun..

[20]  On Difference Schemes and Symplectic Geometry ? X1 Introductory Remarks , 2022 .

[21]  Robert D. Skeel,et al.  Long-Time-Step Methods for Oscillatory Differential Equations , 1998, SIAM J. Sci. Comput..

[22]  G. Rowlands.  A numerical algorithm for Hamiltonian systems , 1991 .

[23]  Bin Wang,et al.  ERKN integrators for systems of oscillatory second-order differential equations , 2010, Comput. Phys. Commun..

[24]  Marlis Hochbruck,et al.  Error analysis of exponential integrators for oscillatory second-order differential equations , 2006 .

[25]  Wojciech Rozmus,et al.  A symplectic integration algorithm for separable Hamiltonian functions , 1990 .

[26]  Xinyuan Wu,et al.  Trigonometrically-fitted ARKN methods for perturbed oscillators , 2008 .

[27]  Jesús Vigo-Aguiar,et al.  Symplectic conditions for exponential fitting Runge-Kutta-Nyström methods , 2005, Math. Comput. Model..

[28]  Guido Vanden Berghe,et al.  Symplectic exponentially-fitted four-stage Runge–Kutta methods of the Gauss type , 2011, Numerical Algorithms.

[29]  Bin Wang,et al.  Two-step extended RKN methods for oscillatory systems , 2011, Comput. Phys. Commun..