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.

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

Frey, S.;Colombo, C.;
2019-01-01

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.
2019
Orbit decay, Atmospheric drag, Semi-analytical propagation, King-Hele
File in questo prodotto:
File Dimensione Formato  
FREYS02-19.pdf

Accesso riservato

Descrizione: Paper
: Publisher’s version
Dimensione 1.57 MB
Formato Adobe PDF
1.57 MB Adobe PDF   Visualizza/Apri
FREYS_OA_02-19.pdf

accesso aperto

Descrizione: Paper open access
: Publisher’s version
Dimensione 1.59 MB
Formato Adobe PDF
1.59 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11311/1079624
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 12
  • ???jsp.display-item.citation.isi??? 3
social impact