Symplectic orbit propagation based on Deprit’s radial intermediary

The constantly challenging requirements for orbit prediction have opened the need for better onboard propagation tools. Runge-Kutta (RK) integrators have been widely used for this purpose; however RK integrators are not symplectic, which means that RK integrators may lead to incorrect global behavior and degraded accuracy. Emanating from Deprit’s radial intermediary, obtained by the elimination of the parallax transformation, we present the development of symplectic integrators of different orders for spacecraft orbit propagation. Through a set of numerical simulations, it is shown that these integrators are more accurate and substantially faster than Runge-Kutta-based methods. Moreover, it is also shown that the proposed integrators are more accurate than analytic propagation algorithms based on Deprit’s radial intermediary solution, and even other previously-developed symplectic integrators.

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

[2]  Jonathan P. How,et al.  Spacecraft Formation Flying: Dynamics, Control and Navigation , 2009 .

[3]  J. C. Simo,et al.  Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics , 1992 .

[4]  Jacques Laskar,et al.  New families of symplectic splitting methods for numerical integration in dynamical astronomy , 2012, 1208.0689.

[5]  P. Gurfil,et al.  Integrable approximation of J2-perturbed relative orbits , 2012 .

[6]  Jacques Laskar,et al.  High precision symplectic integrators for the Solar System , 2013 .

[7]  Martin Lara,et al.  LEO intermediary propagation as a feasible alternative to Brouwer’s gravity solution , 2014, 1406.2200.

[8]  Carles Simó,et al.  Non integrability of theJ2 problem , 1993 .

[9]  M. Qin,et al.  Symplectic Geometric Algorithms for Hamiltonian Systems , 2010 .

[10]  E. Imre,et al.  High-precision, symplectic numerical, relative orbit propagation , 2007 .

[11]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[12]  P. Gurfil,et al.  Variational and symplectic integrators for satellite relative orbit propagation including drag , 2018 .

[13]  Martin Lara,et al.  On Roto-Translatory Motion: Reductions and Radial Intermediaries , 2012 .

[14]  M. Suzuki,et al.  Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simulations , 1990 .

[15]  Pini Gurfil,et al.  Celestial Mechanics and Astrodynamics: Theory and Practice , 2016 .

[16]  Alessandra Celletti,et al.  Non-integrability of the problem of motion around an oblate planet , 1995 .

[17]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[18]  E. Hairer,et al.  Geometric numerical integration illustrated by the Störmer–Verlet method , 2003, Acta Numerica.

[19]  M. Irigoyen,et al.  Non integrability of the J2 problem , 1993 .

[20]  André Deprit,et al.  The elimination of the parallax in satellite theory , 1981 .

[21]  Fernando Casas,et al.  A Concise Introduction to Geometric Numerical Integration , 2016 .

[22]  Daniel J. Scheeres,et al.  Computation and Applications of an Orbital Dynamics Symplectic State Transition Matrix , 2009 .

[23]  J. Junkins,et al.  Analytical Mechanics of Space Systems , 2003 .

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

[25]  Pini Gurfil,et al.  Satellite onboard orbit propagation using Deprit’s radial intermediary , 2014 .