Intermediary LEO propagation including higher order zonal harmonics
  Two new intermediary orbits of the artificial satellite problem are proposed. The analytical solutions include higher order effects of the Geopotential, and are obtained by means of a torsion transformation applied to the quasi-Keplerian system resulting after the elimination of the parallax simplification, for the first intermediary, and after the elimination of the parallax and perigee simplifications, for the second one. The new intermediaries perform notably well for low earth orbits propagation, are free from special functions, and result advantageous, both in accuracy and efficiency, when compared to the standard Cowell integration of the J2 problem, thus providing appealing alternatives for onboard, short-term, orbit propagation under limited computational resources. 
Denis Hautesserres and Martin Lara
math.DS