Symplectic integration of space debris motion considering several Earth's shadowing models

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

Influence of Earth’s shadowing effects on space debris stability

International audienc

On the Dynamics of Space Debris: 1:1 and 2:1 Resonances

We study the dynamics of the space debris in the 1:1 and 2:1 resonances, where geosynchronous and GPS satellites are located. By using Hamiltonian formalism, we consider a model including the geopotential contribution for which we compute the secular and resonant expansions of the Hamiltonian. Within such model we are able to detect the equilibria and to study the main features of the resonances in a very effective way. In particular, we analyze the regular and chaotic behavior of the 1:1 and 2:1 resonant regions by analytical methods and by computing the Fast Lyapunov Indicators, which provide a cartography of the resonances. This approach allows us to detect easily the location of the equilibria, the amplitudes of the libration islands and the main dynamical stability features of the resonances, thus providing an overview of the 1:1 and 2:1 resonant domains under the effect of Earth's oblateness. The results are validated by a comparison with a model developed in Cartesian coordinates, including the geopotential, the gravitational attraction of Sun and Moon and the solar radiation pressure