Variable steps for reversible integration methods

Conventional variable-step implementation of symplectic or reversible integration methods destroy the symplectic or reversible structure of the system. We show that to preserve the symplectic structure of a method the step size has to be kept almost constant. For reversible methods variable steps are possible but the step size has to be equal for “reflected” steps. We demonstrate possible ways to construct reversible variable step size methods. Numerical experiments show that for the Kepler problem the new methods perform better than conventional variable step size methods or symplectic constant step size methods. In particular they exhibit linear growth of the global error (as symplectic methods with constant step size).ZusammenfassungWerden symplektische oder reversible Integrationsverfahren mit herkömmlichen Schrittweitensteuerungen verwendet, so geht die symplektische, bzw. die reversible Struktur des Problems verloren. In dieser Arbeit wird gezeigt, dass die symplektische Struktur des Verfahrens nur dann erhalten wird, wenn die Schrittweite fast konstant bleibt. Für reversible Verfahren sind echt variable Schrittweiten möglich, die Schrittweite muss jedoch für “gespiegelte” Schritte gleich sein. Es werden verschiedene Wege aufgezeigt, um reversible, variable Schrittweiten zu konstruieren. Numerische Experimente zeigen, dass für das Keplerproblem die neuen Methoden den herkömmlichen Schrittweitensteuerungen oder den symplektischen Verfahren mit konstanter Schrittweite überlegen sind. Insbesondere wächst der globale Fehler linear.

[1]  J. Moser Stable and Random Motions in Dynamical Systems: With Special Emphasis on Celestial Mechanics. , 1973 .

[2]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

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

[4]  C. Scovel,et al.  Symplectic integration of Hamiltonian systems , 1990 .

[5]  Haruo Yoshida Recent progress in the theory and application of symplectic integrators , 1993 .

[6]  Yifa Tang The necessary condition for a Runge-Kutta scheme to be symplectic for Hamiltonian systems , 1993 .

[7]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[8]  Robert D. Skeel,et al.  An explicit Runge-Kutta-Nystro¨m method is canonical if and only if its adjoint is explicit , 1992 .

[9]  C. Scovel,et al.  Symplectic integration of Hamiltonian systems , 1990 .

[10]  Urs Kirchgraber,et al.  Multi-step methods are essentially one-step methods , 1986 .

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

[12]  K. Nipp,et al.  Invariant curves for variable step size integrators , 1991 .

[13]  Daniel Stoffer,et al.  General linear methods: connection to one step methods and invariant curves , 1993 .

[14]  J. M. Sanz-Serna,et al.  Symplectic integrators for Hamiltonian problems: an overview , 1992, Acta Numerica.

[15]  C. Scovel,et al.  A survey of open problems in symplectic integration , 1993 .

[16]  F. Lasagni Canonical Runge-Kutta methods , 1988 .