In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones. Indeed, it has already been demonstrated that such objects could stay there for hundreds of years. Our model takes into account the Earth's gravitational potential, luni-solar and planetary gravitational perturbations and direct solar radiation pressure. Based on the analysis of the energy conservation and on a comparison with a high order non-symplectic integrator, we show that our algorithm allows us to use large time steps and keep accurate results. We also propose an innovative method to model Earth's shadow crossings by means of a smooth shadow function. In the particular framework of symplectic integration, such a function needs to be included analytically in the equations of motion in order to prevent numerical drifts of the energy. For the sake of completeness, both cylindrical shadows and penumbra transitions models are considered. We show that both models are not equivalent and that big discrepancies actually appear between associated orbits, especially for high area-to-mass ratios.
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 104–105 yr.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.