Linear energy-preserving integrators for Poisson systems

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 infinitely many stages.

[1]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[2]  G. Quispel,et al.  A new class of energy-preserving numerical integration methods , 2008 .

[3]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[4]  O. Gonzalez.  Time integration and discrete Hamiltonian systems , 1996 .

[5]  E. Hairer.  Energy-preserving variant of collocation methods 1 , 2010 .

[6]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

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

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

[9]  G. Quispel,et al.  Geometric integration using discrete gradients , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

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

[11]  F. Iavernaro,et al.  High-order Symmetric Schemes for the Energy Conservation of Polynomial Hamiltonian Problems 1 2 , 2009 .

[12]  J. Strelcyn,et al.  Integrals of quadratic ordinary differential equations in R3: The Lotka-Volterra system , 1990 .

[13]  E. Hairer.  Energy-Preserving Variant of Collocation Methods 12 , 2010 .

[14]  Ander Murua,et al.  An Algebraic Approach to Invariant Preserving Integators: The Case of Quadratic and Hamiltonian Invariants , 2006, Numerische Mathematik.