Several dynamical systems of interest in celestial mechanics can be written in the formqFor instance, the modified Kepler problem, the spin-orbit model and the Lane-Emden equation all belong to this class. In this work we start an investigation of these models from the point of view of contact geometry. In particular we focus on the (contact) Hamiltonisation of these models and on the construction of the corresponding geometric integrators.

Starting from a contact Hamiltonian description of Liénard systems, we introduce a new family of explicit geometric integrators for these nonlinear dynamical systems. Focusing on the paradigmatic example of the van der Pol oscillator, we demonstrate that these integrators are particularly stable and preserve the qualitative features of the dynamics, even for relatively large values of the time step and in the stiff regime.

No abstract

Contact integrators are a family of geometric numerical schemes which guarantee the conservation of the contact structure. In this work we review the construction of both the variational and Hamiltonian versions of these methods. We illustrate some of the advantages of geometric integration in the dissipative setting by focusing on models inspired by recent studies in celestial mechanics and cosmology.