Fourth-order theories for orbit predictions for low- and high-eccentricity orbits in an oblate atmosphere with scale height dependent on altitude

Two new fourth-order non-singular analytical theories for the motion of near-Earth satellite orbits with air drag are developed for low- and high-eccentricity orbits in an oblate atmosphere with variation of density scale height with altitude. Uniformly regular Kustaanheimo–Stiefel (KS) canonical elements are utilized for low-eccentricity orbits and KS element equations are employed for high-eccentricity orbits. Only two of the nine equations are solved analytically to compute the state vector and change in energy at the end of each revolution, due to symmetry in the equations of motion. The analytical solutions are compared with the numerically integrated values up to 100 revolutions, and found to be quite accurate over a wide range of eccentricity, perigee height and inclination.