Abstract
The paper offers the fully analytic solution to the motion of a satellite orbiting under the influence of the two major perturbations, due to the oblateness and the atmospheric drag. The solution is presented in a time-explicit form, and takes into account an exponential distribution of the atmospheric density, an assumption that is reasonably close to reality. The approach involves two essential steps. The first one concerns a new approximate mathematical model that admits a closed-form solution with respect to a set of new variables. The second step is the determination of an infinitesimal contact transformation that allows to navigate between the new and the original variables. This contact transformation is obtained in exact form, and afterwards a Taylor series approximation is proposed in order to make all the computations explicit. The aforementioned transformation accommodates both perturbations, improving the accuracy of the orbit predictions by one order of magnitude with respect to the case when the atmospheric drag is absent from the transformation. Numerical simulations are performed for a low Earth orbit starting at an altitude of 350 km, and they show that the incorporation of drag terms into the contact transformation generates an error reduction by a factor of 7 in the position vector. The proposed method aims at improving the accuracy of analytic orbit propagation and transforming it into a viable alternative to the computationally intensive numerical methods.
Similar content being viewed by others
Notes
Due to the presence of significant uncertainties in the atmosphere (influence of the solar activity, day/night density variations, the atmospheric bulge, winds), realistic long-term propagations in the presence of drag cannot be addressed with deterministic tools (see Dell’Elce and Kerschen 2014).
The critical inclination is a resonance of the dynamical system (see Lara 2015b) and therefore cannot be removed.
Deprit was not the first one to use this approach. Hori (1966) proposed a slightly different version a few years before.
They are uniquely defined up to a translation on the real numbers axis that is an integer multiple of \(2\pi \). But since only their trigonometric functions \(\sin \) and \(\cos \) are involved in computations, no restrictions need to be imposed.
A transformation is said to be canonical if its Jacobian is a symplectic matrix, i.e. it obeys the condition in Eq. (19).
Whittaker variables are singular for equatorial orbits, but this does not pose a problem for the case where only even zonal harmonics are taken into consideration, which is the case of the present paper.
Elementary transformations between different sets of coordinates are also required at some steps of the algorithm, but these coordinates are present in the analytic expressions, making this requirement obvious,.
These values were inspired by the prospective QARMAN satellite, a CubeSat mission developped by the von Karman Institute of Fluid Dynamics in Brussels and the University of Liège.
The argument of perigee drift was omitted, since it is incorporated in the argument of latitude, depicted in Fig. 2c.
Consequently, the orthogonal right-handed frame associated to the Poincaré canonical variables is found to be the same as for the equinoctial ones (see Battin 1999, Problem 10-10, p. 494).
References
Barrio, R., Palácian, J.: Lie Transforms for ordinary differential equations: taking advantage of the Hamiltonian form of terms of the perturbation. Int. J. Numer. Meth. Eng. 40(12), 2289–2300 (1997)
Barrio, R., Palacián, J.: High-order averaging of eccentric artificial satellites perturbed by the earth’potential and air-drag terms. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 459(2034), 1517–1534 (2003)
Battin, R.H.: An Introduction to the Mathematics and Methods of Astrodynamics. AIAA, Reston, VA (1999)
Brouwer, D.: Solution of the problem of artificial satellite theory without drag. Astron. J. 64, 378–396 (1959)
Brouwer, D., Hori, G.I.: Theoretical evaluation of atmospheric drag effects in the motion of an artificial satellite. Astron. J. 66, 193–225 (1961)
Celletti, A., Negrini, P.: Non-integrability of the problem of motion around an oblate planet. Celest. Mech. Dyn. Astron. 61, 253–260 (1995)
Cid, R., Lahulla, J.F.: Perturbaciones de corto periodo en el movimiento de un satélite artificial, en función de las variables de Hill. Publicaciones de la Revista de la Academia de Ciencias de Zaragoza 24, 159–165 (1969)
Dell’Elce, L., Kerschen, G.: Probabilistic assessment of the lifetime of low-earth-orbit spacecraft: uncertainty characterization. J. Guid. Control Dyn. 1, 1–13 (2014)
Deprit, A.: Canonical transformations depending on a small parameter. Celest. Mech. 1, 12–30 (1969)
Deprit, A.: The elimination of the parallax in satellite theory. Celest. Mech. 24, 111–153 (1981)
Garfinkel, B.: The orbit of a satellite of an oblate planet. Astron. J. 64, 353–366 (1959)
Goldtsein, H., Poole, ChP, Safko., J.L.: Classical Mechanics, 3rd edn. Addison-Wesley, Boston (2002)
Gurfil, P., Lara, M.: Satellite onboard orbit propagation using Deprit’s radial intermediary. Celest. Mech. Dyn. Astron. 120(2), 217–232 (2014)
Hilton, C.G., Kuhlman, J.R.: Mathematical Models for the Space Defense Center. Philco-Ford Corporation Publication U-3871, pp. 17–28 (1966)
Hori, G.: Theory of general perturbation with unspecified canonical variables. Publ. Astron. Soc. Jpn. 18, 287–296 (1966)
Irigoyen, M., Simo, C.: Non integrability of the \(J_2\) problem. Celest. Mech. Dyn. Astron. 55, 281–287 (1993)
Lane, M.H., Cranford, K.H.: American Institute of Aeronautics, Astronautics, and American Astronautical Society. In: An Improved Analytical Drag Theory for the Artificial Satellite Problem. AIAA paper. American Institute of Aeronautics and Astronautics, (1969). https://books.google.be/books?id=dxeItgAACAAJ
Lara, M.: Efficient formulation of the periodic corrections in Brouwers gravity solution. Math. Probl. Eng. 2015, 9 (2015)
Lara, M.: Leo intermediary propagation as a feasible alternative to brouwer’s gravity solution. Adv. Space Res. 56(3), 367–376 (2015b). (Advances in Asteroid and Space Debris Science and Technology—Part 1)
Lara, M., Gurfil, P.: Integrable approximation of \(J_2\)-perturbed relative orbits. Celest. Mech. Dyn. Astron. 114(3), 229–254 (2012)
Lyddane, R.H.: Small eccentricities or inclinations in the Brouwer theory of the artificial satellite. Astron. J. 68, 555–558 (1963)
Martinusi, V., Dell’Elce, L., Kerschen, G.: Analytic propagation of near-circular satellite orbits in the atmosphere of an oblate planet. Celest. Mech. Dyn. Astron. 123(1), 85–103 (2015)
Steinberg, S.: Lie series and nonlinear ordinary differential equations. J. Math. Anal. Appl. 101(1), 39–63 (1984)
Sterne, T.E.: The gravitational orbit of a satellite of an oblate planet. Astron. J. 63, 28–40 (1958)
Vallado, D.A.: Fundamentals of Astrodynamics and Applications. Microcosm, Cleveland (2001)
Vallado, D.A., Crawford, P., Hujsak, R., Kelso, T.S.: Revisiting Spacetrack Report #3: Rev 2. In: AIAA Astrodynamics Specialists Conference and Exhibit, August 2006 (2006)
Vinti, J.P.: Theory of the orbit of an artificial satellite with use of spheroidal coordinates. Astron. J. 65, 353–354 (1960)
von Zeipel, H.: Recherches sur le mouvement des petites planètes. Almqvist & Wiksells, Stockholm (1916)
Author information
Authors and Affiliations
Corresponding author
Appendix
Appendix
1.1 Lie series approximation for \(\overline{L}\)
Although the variable \(\overline{L}\) has an explicit expression, depicted in Eq. (10a), the numerical implementation of the complex error function, as well as its inverse, rise the problem of dealing with differences of very large numbers. It would be much more convenient to have an accurate approximation for \(\overline{L}\), and this is obtained by writing the solution to Eq. (8d) with the help of Lie series. For a complete description of the expansion mechanism, see Steinberg (1984).
Denote:
The IVP to solve is rewritten as:
and the change of variable:
is performed. Denote
The IVP in z is:
Consider D to be the differential operator, defined for any arbitrary function \(\phi \):
and define its formal powers as:
Consider the formal exponential operator:
Then the solution to IVP (50) is expressed as:
and explicitely:
By taking into account the expression in Eq. (51) of the operator D, the formal powers \(\left. D^{k}z\right| _{z=0}\) are expressed as:
Then \(\overline{L}=\overline{L}\left( t\right) \) is evaluated based on Eq. ( 49):
1.2 Evaluation of \(I_{1,2,3}\) integrals
In the form they are written in Eq. (13), the integrals \(I_{1,2,3}\) are not suitable for an efficient series expansion approximation, given the considerable length of the interval on which the integration is performed (i.e., for a good approximation, a large number of terms are required in the truncated series expansion). For a rapid convergence, and also to avoid differences between very large numbers (since \(L=\sqrt{\mu \alpha }\) has the same order of magnitude of the angular momentum). It would therefore be ideal to transform them to integrals from 0 to some value z that has the smallest possible absolute value, so that the accuracy of the series expansion is maximum.
First, note that \(I_{1,2,3}\) may be rewritten as:
where
To this end, denote
and perform the change of variable:
The integrals become:
By performing a Taylor series expansion for \(\exp \left( y\right) ,\) the integrals become:
and written separately:
The evaluation of \(I_{1,2,3}\) has been simplified to definite integrals that may be expressed explicitly through elementary functions.
1.3 Jacobian for coordinate transformations
The Jacobian of the transformation from Delaunay to Whittaker elements is:
The Jacobian of the transformation from Poincaré to Delaunay variables is:
The Jacobian of the transformation from Whittaker to Poincaré elements may now be explicitly derived as:
Since it involves only elementary computations, that are nevertheless leading to relatively long expressions, the full matrix form of \({\mathcal {J}} _{\mathcal {P}}^{{\mathcal {W}}}\) will not be displayed. In any case, it is obtained by the multiplication of two explicit matrices.
1.4 Transformation from Poincaré to Cartesian coordinates
The Poincaré canonical variables are mentioned and used in very few Astrodynamics/Orbital Mechanics textbooks, but they lack a comprehensive approach, and to our knowledge, a straightforward transformation from these variables to Cartesian is absent from the literature. The closest approach is found in Battin (1999), but it is made for the equinoctial variables, that are nonsingular but not canonical.
The transformation is determined here in the most direct way possible. It is known that the inertial position and velocity vectors are expressed in the perifocal frame \(\left\{ P\right\} \) as:
and their corresponding column matrices in the inertial frame originated in the attraction center are:
where the orthogonal matrices \({\mathbf {R}}_{h},\) \({\mathbf {R}}_{i},\) \({\mathbf {R}} _{g}\) are:
If a new fixed frame \(\left\{ \mathcal {P}\right\} \) is defined such that the position vector has the expression:
it will follow (by elementary computations) that the velocity in the same frame \(\left\{ \mathcal {P}\right\} \) is:
Consequently, the inertial counterparts of \({\mathbf {r}}^{\left\{ \mathcal {P} \right\} },\) \({\mathbf {v}}^{\left\{ \mathcal {P}\right\} }\) are:
The matrix \(\mathbf {A}\) is found to beFootnote 11:
Define U such that:
Then the explicit expression of matrix \(\mathbf {A,}\) in Poincaré variables, is:
Rights and permissions
About this article
Cite this article
Martinusi, V., Dell’Elce, L. & Kerschen, G. First-order analytic propagation of satellites in the exponential atmosphere of an oblate planet. Celest Mech Dyn Astr 127, 451–476 (2017). https://doi.org/10.1007/s10569-016-9734-8
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10569-016-9734-8