Extension of the King-Hele orbit contraction method for accurate, semi-analytical propagation of non-circular orbits

Abstract

Numerical integration of orbit trajectories for a large number of initial conditions and for long time spans is computationally expensive. Semi-analytical methods were developed to reduce the computational burden. An elegant and widely used method of semi-analytically integrating trajectories of objects subject to atmospheric drag was proposed by King-Hele (KH). However, the analytical KH contraction method relies on the assumption that the atmosphere density decays strictly exponentially with altitude. If the actual density profile does not satisfy the assumption of a fixed scale height, as is the case for Earth's atmosphere, the KH method introduces potentially large errors for non-circular orbit configurations.In this work, the KH method is extended to account for such errors by using a newly introduced atmosphere model derivative. By superimposing exponentially decaying partial atmospheres, the superimposed KH method can be applied accurately while considering more complex density profiles. The KH method is further refined by deriving higher order terms during the series expansion. A variable boundary condition to choose the appropriate eccentricity regime, based on the series truncation errors, is introduced. The accuracy of the extended analytical contraction method is shown to be comparable to numerical Gauss-Legendre quadrature. Propagation using the proposed method compares well against non-averaged integration of the dynamics, while the computational load remains very low.

Analytical Propagation of Space Debris Density for Collisions near Sun-Synchronous Orbits

“…This effect is neglected here. For extensions to modeling the evolution of the distribution of the eccentricity, one can refer to the work in [36]. In the follow-ing parts, we validate the analytical propagation method for evolution of debris clouds in the SSO region.…”

A Machine Learning and Feature Engineering Approach for the Prediction of the Uncontrolled Re-Entry of Space Objects

“…The output hidden state is computed in two steps. First, the candidate memory cell c ∈ R k is calculated according to Equation (11). It takes as input the previous hidden state, which is multiplied by the reset gate, and the current input.…”

“…However, in this work, we focused on the latest stage of the re-entry prediction of low eccentricity orbits; therefore, the dominant force is atmospheric drag and the main uncertainties are related to the atmospheric density. It is modelled through empirical models which describe the density variations on spatial and temporal scales [7][8][9][10][11]. However, they are affected by two types of uncertainties [10]: the simplified physical modelling on which they are based, and the uncertainties and complexities in predicting the space weather, in particular the solar index.…”

Short and long-term reconstruction of in-orbit fragmentation events

“…Frey et al [4] developed a method to search for a fragmentation event in the long-term, i.e. in the order of years, exploiting a continuum approach for fragmentation modelling in LEO, based on the density of objects.…”