Open AccessThe Global Symplectic Integrator: an efficient tool for stability studies of dynamical systems. Application to the Kozai resonance in the restricted three‐body problem
Author(s) -
Libert A.S.,
Hubaux Ch.,
Carletti T.
Publication year - 2011
Publication title -
monthly notices of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.058
H-Index - 383
eISSN - 1365-2966
pISSN - 0035-8711
DOI - 10.1111/j.1365-2966.2011.18431.x
Subject(s) - symplectic integrator , symplectic geometry , integrator , physics , hamiltonian system , variational integrator , semi implicit euler method , orbit (dynamics) , hamiltonian (control theory) , stability (learning theory) , symplectomorphism , mathematics , classical mechanics , symplectic manifold , mathematical analysis , computer science , mathematical optimization , quantum mechanics , aerospace engineering , voltage , machine learning , backward euler method , engineering , euler equations
Following the discovery of extrasolar systems, the study of long‐term evolution and stability of planetary systems is enjoying a renewed interest. While non‐symplectic integrators are very time‐consuming because of the very long time‐scales and the small integration steps required to have a good energy preservation, symplectic integrators are well suited for the study of such orbits on long time‐spans. However, stability studies of dynamical systems generally rely on non‐symplectic integrations of deviation vectors. In this work we propose a numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, hereby called Global Symplectic Integrator . It consists of the simultaneous integration of the orbit and the deviation vectors using a symplectic scheme of any order. In particular, due to its symplectic properties, the proposed method allows us to recover the correct orbit characteristics using very large integration time‐steps, fluctuations of energy around a constant value and short CPU times. It proves to be more efficient than non‐symplectic schemes to correctly identify the behaviour of a given orbit, especially on dynamics acting on long time‐scales. To illustrate the numerical performances of the global symplectic integrator, we will apply it to the well‐known toy problem of Hénon–Heiles and the challenging problem of the Kozai resonance in the restricted three‐body problem, whose secular effects have periods of the order of 10 4 –10 5 yr.