The minimal stage, energy preserving Runge–Kutta method for polynomial Hamiltonian systems is the averaged vector field method

Celledoni, Elena; Owren, Brynjulf; Sun, Yajuan

doi:10.1090/S0025-5718-2014-02805-6

The minimal stage, energy preserving Runge–Kutta method for polynomial Hamiltonian systems is the averaged vector field method
HTML articles powered by AMS MathViewer

by Elena Celledoni, Brynjulf Owren and Yajuan Sun PDF

Math. Comp. 83 (2014), 1689-1700 Request permission

Abstract:

No Runge–Kutta method can be energy preserving for all Hamiltonian systems. But for problems in which the Hamiltonian is a polynomial, the averaged vector field (AVF) method can be interpreted as a Runge–Kutta method whose weights $b_i$ and abscissae $c_i$ represent a quadrature rule of degree at least that of the Hamiltonian. We prove that when the number of stages is minimal, the Runge–Kutta scheme must in fact be identical to the AVF scheme.

References

M. P. Calvo and J. M. Sanz-Serna, Canonical B-series, Numer. Math. 67 (1994), no. 2, 161–175. MR 1262779, DOI 10.1007/s002110050022
Elena Celledoni, Robert I. McLachlan, David I. McLaren, Brynjulf Owren, G. Reinout W. Quispel, and William M. Wright, Energy-preserving Runge-Kutta methods, M2AN Math. Model. Numer. Anal. 43 (2009), no. 4, 645–649. MR 2542869, DOI 10.1051/m2an/2009020
Elena Celledoni, Robert I. McLachlan, Brynjulf Owren, and G. R. W. Quispel, Energy-preserving integrators and the structure of B-series, Found. Comput. Math. 10 (2010), no. 6, 673–693. MR 2728426, DOI 10.1007/s10208-010-9073-1
E. Celledoni, B. Owren, and Y. Sun, The minimal stage, energy preserving Runge–Kutta method for polynomial Hamiltonian systems is the Averaged Vector Field method, arXiv:1203.3252.
Philippe Chartier, Erwan Faou, and Ander Murua, An algebraic approach to invariant preserving integrators: the case of quadratic and Hamiltonian invariants, Numer. Math. 103 (2006), no. 4, 575–590. MR 2221062, DOI 10.1007/s00211-006-0003-8
Erwan Faou, Ernst Hairer, and Truong-Linh Pham, Energy conservation with non-symplectic methods: examples and counter-examples, BIT 44 (2004), no. 4, 699–709. MR 2211040, DOI 10.1007/s10543-004-5240-6
E. Hairer, Energy-preserving variant of collocation methods, JNAIAM. J. Numer. Anal. Ind. Appl. Math. 5 (2010), no. 1-2, 73–84. MR 2833602
Ernst Hairer, Christian Lubich, and Gerhard Wanner, Geometric numerical integration, 2nd ed., Springer Series in Computational Mathematics, vol. 31, Springer-Verlag, Berlin, 2006. Structure-preserving algorithms for ordinary differential equations. MR 2221614
Ernst Hairer and Christophe J. Zbinden, On conjugate symplecticity of B-series integrators, IMA J. Numer. Anal. 33 (2013), no. 1, 57–79. MR 3020950, DOI 10.1093/imanum/drs010
F. Iavernaro and D. Trigiante, High-order symmetric schemes for the energy conservation of polynomial Hamiltonian problems, JNAIAM J. Numer. Anal. Ind. Appl. Math. 4 (2009), no. 1-2, 87–101. MR 2598786
Benedict Leimkuhler and Sebastian Reich, Simulating Hamiltonian dynamics, Cambridge Monographs on Applied and Computational Mathematics, vol. 14, Cambridge University Press, Cambridge, 2004. MR 2132573
A. Murua, Formal series and numerical integrators. I. Systems of ODEs and symplectic integrators, Appl. Numer. Math. 29 (1999), no. 2, 221–251. MR 1666537, DOI 10.1016/S0168-9274(98)00064-6
Matteo Petrera, Andreas Pfadler, and Yuri B. Suris, On integrability of Hirota-Kimura type discretizations, Regul. Chaotic Dyn. 16 (2011), no. 3-4, 245–289. MR 2810980, DOI 10.1134/S1560354711030051
G. R. W. Quispel and D. I. McLaren, A new class of energy-preserving numerical integration methods, J. Phys. A 41 (2008), no. 4, 045206, 7. MR 2451073, DOI 10.1088/1751-8113/41/4/045206
Ge Zhong and Jerrold E. Marsden, Lie-Poisson Hamilton-Jacobi theory and Lie-Poisson integrators, Phys. Lett. A 133 (1988), no. 3, 134–139. MR 967725, DOI 10.1016/0375-9601(88)90773-6