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

Abstract 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.

[1]  P. Farinella,et al.  Solar radiation pressure perturbations for Earth satellites. 1: A complete theory including penumbra transitions , 1993 .

[2]  Leland E. Cunningham,et al.  On the computation of the spherical harmonic terms needed during the numerical integration of the orbital motion of an artificial satellite , 1970 .

[3]  Giovanni F. Gronchi,et al.  An Algebraic Method to Compute the Critical Points of the Distance Function Between Two Keplerian Orbits , 2005 .

[4]  S. F. Mello.  Analytical Study of the Earth's Shadowing Effects on Satellite Orbits , 1972 .

[5]  Seppo Mikkola.  Efficient Symplectic Integration of Satellite Orbits , 1999 .

[6]  P. Farinella,et al.  Solar radiation pressure perturbations for Earth satellites II. an approximate method to model penumbra transitions and their long-term orbital effects on LAGEOS , 1994 .

[7]  Nicolas Bourbaki,et al.  Eléments de mathématique : groupes et algèbres de Lie , 1972 .

[8]  J. Laskar,et al.  High order symplectic integrators for perturbed Hamiltonian systems , 2000 .

[9]  A. Milani,et al.  Theory of Orbit Determination , 2009 .

[10]  C. Chao,et al.  Applied Orbit Perturbation and Maintenance , 2005 .

[11]  J. Liou,et al.  Orbital Dynamics of High Area-To-Mass Ratio Debris and Their Distribution in the Geosynchronous Region , 2005 .

[12]  C. Pardini,et al.  Orbital Evolution of Geosynchronous Objects with High Area-To-Mass Ratios , 2005 .

[13]  Robert Bryant,et al.  THE EFFECT OF SOLAR RADIATION PRESSURE ON THE MOTION OF AN ARTIFICIAL SATELLITE , 1961 .

[14]  A. Compère,et al.  NIMASTEP: a software to modelize, study, and analyze the dynamics of various small objects orbiting specific bodies , 2011, 1112.1304.

[15]  N. Delsate,et al.  Global dynamics of high area-to-mass ratios GEO space debris by means of the MEGNO indicator , 2008, 0810.1859.

[16]  Francois Mignard,et al.  SOLAR RADIATION PRESSURE PERTURBATIONS FOR EARTH SATELLITES. IV. EFFECTS OF THE EARTH'S POLAR FLATTENING ON THE SHADOW STRUCTURE AND THE PENUMBRA TRAN SITIONS , 1996 .

[17]  Carles Simó,et al.  Phase space structure of multi-dimensional systems by means of the mean exponential growth factor of nearby orbits , 2003 .

[18]  Nicolas Bourbaki,et al.  Elements de Mathematiques , 1954, The Mathematical Gazette.

[19]  Phil Palmer,et al.  An Implementation of the Logarithmic Hamiltonian Method for Artificial Satellite Orbit Determination , 2002 .

[20]  P. R. Escobal,et al.  Methods of orbit determination , 1976 .

[21]  C. Murray,et al.  Solar System Dynamics: Expansion of the Disturbing Function , 1999 .

[22]  Anne Lemaitre,et al.  Semi-analytical investigations of high area-to-mass ratio geosynchronous space debris including Earth’s shadowing effects , 2008 .

[23]  Johannes Herzog,et al.  AIUB Efforts to Survey, Track, and Characterize Small-size Objects at High Altitudes , 2012 .

[24]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[25]  Luciano Anselmo,et al.  Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios , 2008 .

[26]  S. Breiter.  Explicit Symplectic Integrator for Highly Eccentric Orbits , 1998 .

[27]  B. Melendo,et al.  Long-term predictability of orbits around the geosynchronous altitude , 2005 .

[28]  P. Farinella,et al.  Solar radiation pressure perturbations for Earth satellites. III. Global atmospheric phenomena and the albedo effect , 1994 .

[29]  Oliver Montenbruck,et al.  Satellite Orbits: Models, Methods and Applications , 2000 .

[30]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[31]  I. M. Levitt.  Advances in the astronautical sciences: Vol. 6, edited by Horace Jacobs and Eric Burgess. 898 pages, diagrams, illustrations, 612 × 978 in. New York, The Macmillan Co., 1961. Price, $25.00 , 1961 .

[32]  A web of secondary resonances for large A/m geostationary debris , 2009 .

[33]  P. Lala,et al.  The Earth's Shadowing Effects in the Short-Periodic Perturbations of Satellite Orbits , 1970 .

[34]  K. Aksnes,et al.  Short-period and long-period perturbations of a spherical satellite due to direct solar radiation , 1976 .

[35]  N. K. Pavlis,et al.  The Development of the Joint NASA GSFC and the National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96 , 1998 .

[36]  F. Ryland.  The effects of solar radiation pressure on the motion of an artificial satellite , 1979 .