Explicit methods in extended phase space for inseparable Hamiltonian problems

We present a method for explicit leapfrog integration of inseparable Hamiltonian systems by means of an extended phase space. A suitably defined new Hamiltonian on the extended phase space leads to equations of motion that can be numerically integrated by standard symplectic leapfrog (splitting) methods. When the leapfrog is combined with coordinate mixing transformations, the resulting algorithm shows good long term stability and error behaviour. We extend the method to non-Hamiltonian problems as well, and investigate optimal methods of projecting the extended phase space back to original dimension. Finally, we apply the methods to a Hamiltonian problem of geodesics in a curved space, and a non-Hamiltonian problem of a forced non-linear oscillator. We compare the performance of the methods to a general purpose differential equation solver LSODE, and the implicit midpoint method, a symplectic one-step method. We find the extended phase space methods to compare favorably to both for the Hamiltonian problem, and to the implicit midpoint method in the case of the non-linear oscillator.

[1]  S. Aarseth,et al.  A Time-Transformed Leapfrog Scheme , 2002 .

[2]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[3]  Peter Deuflhard,et al.  Scientific Computing with Ordinary Differential Equations , 2002 .

[4]  William Kahan,et al.  Composition constants for raising the orders of unconventional schemes for ordinary differential equations , 1997, Math. Comput..

[5]  Alan C. Hindmarsh,et al.  Description and use of LSODE, the Livermore Solver for Ordinary Differential Equations , 1993 .

[6]  D. Merritt,et al.  IMPLEMENTING FEW-BODY ALGORITHMIC REGULARIZATION WITH POST-NEWTONIAN TERMS , 2007, 0709.3367.

[7]  S. Tremaine,et al.  A Class of Symplectic Integrators with Adaptive Time Step for Separable Hamiltonian Systems , 1999, astro-ph/9906322.

[8]  C. Lanczos.  The variational principles of mechanics , 1949 .

[9]  Parlitz,et al.  Period-doubling cascades and devil's staircases of the driven van der Pol oscillator. , 1987, Physical review. A, General physics.

[10]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[11]  S. Mikkola,et al.  Explicit Symplectic Algorithms For Time‐Transformed Hamiltonians , 1999 .

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

[13]  Balth. van der Pol,et al.  VII. Forced oscillations in a circuit with non-linear resistance. (Reception with reactive triode) , 1927 .

[14]  S. Mikkola,et al.  Algorithmic regularization of the few‐body problem , 1999 .

[15]  William B. Gragg,et al.  On Extrapolation Algorithms for Ordinary Initial Value Problems , 1965 .

[16]  G. Quispel,et al.  Acta Numerica 2002: Splitting methods , 2002 .

[17]  A. Hindmarsh.  LSODE and LSODI, two new initial value ordinary differential equation solvers , 1980, SGNM.

[18]  J. Stoer,et al.  Numerical treatment of ordinary differential equations by extrapolation methods , 1966 .

[19]  Algorithmic regularization with velocity-dependent forces , 2006, astro-ph/0605054.

[20]  S. Mikkola,et al.  Explicit algorithmic regularization in the few-body problem for velocity-dependent perturbations , 2010 .