Abstract
Low-altitude, near-polar orbits are very desirable as science orbits for missions to planetary satellites, such as the Earth's Moon. In this paper, we present an analytical theory with numerical simulations to study the orbital motion of lunar low-altitude artificial satellite. We consider the problem of an artificial satellite perturbed by the nonuniform distribution of the mass of the Moon (β, , and ). The conditions to get frozen orbits are presented. Using an approach that considers the single-averaged problem, we found families of periodic orbits for the problem of an orbiter travelling around the Moon, where frozen orbits valid for long periods of time are found. A comparison between the models for the zonal and tesseral harmonics coefficients is presented.
1. Introduction
Frozen orbits around the Moon, natural satellites, or asteroids are of current interest because several space missions have the goal of orbiting around such bodies (see, for instance, [1β5] and references therein). The dynamics of orbits around a planetary satellite, taking into account the gravitational attraction of a third-body and the nonuniform distribution of mass (, , and ) of the planetary satellite (Europa), is studied in [6β11].
This paper contains a continuation of the work developed in [1, 12, 13], where the problem of the third-body gravitational perturbation and the non-uniform distribution of mass (, ) of the planetary satellite (Moon) on a spacecraft is studied using a simplified dynamical model. Also in [14, 15], the problem of critical inclination orbits for lunar artificial satellites including the terms () zonal, () sectorial, and the lunar rotation due to the lunar potential, besides the Keplerian term is considered. In [16], a study done to control the increasing eccentricity of a polar lunar satellite for different altitudes due to the perturbations of the Earth is presented. The approach presented here is the use of low-thrust propulsion in order to keep the orbital eccentricity of the satellite at low values. A fourth-order analytical theory is presented for the accurate computation of quasiperiodic frozen orbits in [17].
In our research, we use a simplified dynamical model that considers the effects caused by the nonsphericity of the Moon (β, , and ) using single-averaged analytical model over the short satellite period. The motion of equations is obtained from the Lagrange planetary equations, and then they are integrated numerically, using the software Maple. Emphasis is given to the case of frozen orbits, defined as orbits which have constant values of eccentricity, inclination, and argument of the periapsis when the single-averaged system is considered. An approach is used for the case of frozen orbits similar to the one used in [1, 2, 5]. The basic idea is to eliminate the terms due to the short time periodic motion of the spacecraft to show smooth curves for the evolution of the mean orbital elements for a long time period. The interesting set of results allow us to find good conditions close to frozen orbits.
In other words (see [5, 11, 18, 19]), a better understanding of the physical phenomenon can be obtained and it allows the study of long-term stability of the orbits in the presence of disturbances that cause low changes in the orbital elements.
The paper has four sections. Section 2 is devoted to the terms due to the nonsphericity of the Moon and to the Hamiltonian of the system. Applications taking into account the long-period disturbing potential using the averaging method are presented in Section 3. Section 4 shows the conclusions.
2. Oblateness of the Moon
To analyze the motion of a lunar low-altitude artificial satellite, it is necessary to take into account the Moon's oblateness [1, 14, 20]. Besides that, the Moon is much less flattened than the Earth but it also causes perturbations in space vehicles. Table 1 shows the orders of magnitude for the zonal and sectorial harmonics for both celestial bodies to make a comparison. The term describes the equatorial bulge of the Moon, often referred to as oblateness. The coefficient measures the ellipticity of the equator.
The force function, the negative of the total energy as used in physics, is given by here, is the opposite of the potential energy, and is the kinetic energy. The force function can be written as or . The function , comprising all terms of except the central term, is known as the disturbing potential. The space vehicle is a point of mass in a three-dimensional orbit with orbital elements (semimajor axis), (eccentricity), (inclination), (argument of the periapsis), (longitude of the ascending node), and (mean motion) given by the third Kepler law , where is the gravitational constant and is the radius vector of the body (mass of the planetary satellite). The term due to the unperturbed potential is given by
Considering the lunar equatorial plane as the reference plane, the disturbing potential can be written in the form [21] where is the lunar gravitational constant, is the equatorial radius of the Moon (βkm), represent the Legendre polynomials, represent the associated Legendre polynomials, the angle is the latitude of the orbit with respect to the equator of the Moon, and the angle is the longitude measured from the direction of the longest axis of the Moon, where , since is the longitude measured from any fixed direction, and is the longitude of the Moon's longest meridian from the same fixed direction. However, contains explicitly the time [20]. Using spherical trigonometry, we have , where is the true anomaly of the satellite.
The following assumptions will be used in this work: (a) the motion of the Moon is uniform (librations are neglected), (b) the lunar equator lies in the ecliptic (we neglect the inclination of about 1.5Β° of the lunar equator with respect to the ecliptic and the inclination of the lunar orbit with respect to the ecliptic of about 5Β°), (c) the longitude of the lunar longest meridian is equal to the mean longitude of the Earth (librations are neglected), and (d) the mean longitude of the Earth is equal to .
Since the variables and appear only as a combination of , where , being the lunar mean motion. The degrees of freedom can be reduced by choosing as a new variable. A new term must then be added to the Hamiltonian in order to get . This term is used here in the same way it is used in [20]. The Hamiltonian is still time dependent through . Since the longest meridian is always pointing toward the Earth, it is possible to choose a rotating system whose -axis passes through this meridian.
Regarding the Earth's potential, the dominant coefficient is . The rest are of higher order [2]. Opposite to the situation of the Earth, the first harmonics of the lunar potential are almost of the same order (see Table 1). This fact complicates the choice of the harmonic where the potential can be truncated, and this fact makes its choice a little arbitrary. The influences of the Earth and of the nonsphericity of the Moon on the stability of lunar satellites were also shown by [22], but the sectorial harmonics were not considered.
In terms of the orbital elements, the Legendre-associated functions for the zonal up to and sectorial () terms can be written in the form [20, 21] where we use the shortcuts , , , and .
With this, the zonal perturbation due to the oblateness is [14]
Replacing the relation and using Cayley's tables [23] to express the true anomaly in terms of the mean anomaly, after some algebraic manipulations, we get
For the sectorial perturbation [20, 21], we get
where . With some manipulations, we get where the disturbing potential is written in the form .
2.1. The Hamiltonian System
The Hamiltonian of the dynamical system associated with the problem of the orbital motion of the satellite around the Moon taking into account the non-uniform distribution of the mass of the planetary satellite (β, , and ) can be written in the following form: where
Now, using , and , we get where , , , , , .
We arrange the Hamiltonian as follows where
The motion of the spacecraft is studied under the single-averaged analytical model. The average is taken over the mean motion of the satellite. The standard definition for average used in this research is [24] where is the mean anomaly of the satellite. The single-average method is applied in (2.13) to eliminate the terms of short period, to analyze the effect of the disturbing potential on the orbital elements. Periodic terms will be calculated, replacing the result in the Lagrange planetary equations [25] and numerically integrated using the Maple software. The long-period disturbing potential is given in the appendix.
3. Applications: The Long-Period Disturbing Potential Using the Averaging Method
With a simplified model for the disturbing potential, it is possible to analyse the orbital motion of the satellite.
To study the lifetimes of low-altitude artificial satellites of the Moon [26] we took into account the zonal (β) and sectorial (, ) terms in the disturbing potential. However, [5] shows that the effect of the term is an order of magnitude larger than the β terms (see Table 2, model given by [27]), and therefore it should be considered in the potential. In [5], an analytic model is presented to find frozen orbits for lunar satellites, where the potential only due to and zonal harmonics is taken into account. In this paper, we present an approach taking into account the β, , and terms.
3.1. Frozen Orbits Solutions
The definition of frozen orbit is written in the form
From this, we can obtain frozen orbits, say, for or , and or (see [1, 13]).
Replacing by the potential in (A.1) in the Lagrange's planetary equations [25] and solving the equation , we get an equation that depends on three variables, , , and , that can be represented as a 3-D surface as shown in Figures 1, 2, 3, and 4. Points on this surface correspond to equilibria of the system in equations and , that is, frozen orbits. Figures 1 and 3 take into account the coefficients β and , while Figures 2 and 4 take into account the β, , and terms. Comparing the models given by [27] (see Table 2) and [28] (see Table 3), we can analyze that the results are very close and the model given by [28] presents better results to obtain a family of frozen orbits for inclinations above 69Β°. In [30], an analytical theory based on Lie-Deprit transformation is used to locate family of frozen orbits for asteroids and natural satellites.
For inclinations larger than 69Β°, we found a family of frozen orbit, as we can see in Figures 1, 2, 3, and 4, mainly considering the model LP165P for the coefficients of the lunar gravity field [28]. We also observed that the term contributed efficiently to obtain of such orbits.
The curves obtained for constant-a sections are very similar; thus, in order to go into more detail, we define a constant value for the semimajor axis βkm, which indeed corresponds to low orbits. For the chosen value of , we plot the curve versus (see Figures 5 and 6). Figure 5 does not take into account the term in the potential, while Figure 6 considers the sectorial term. We observed an increase of the eccentricity for inclinations between 28.65Β° and 45.84Β° (comparing Figures 5 and 6). For orbits with inclination larger than 69Β° degree, we find a family of frozen orbits (see Figure 6). The region where the eccentricity is frozen is for small values of the initial eccentricity, when compared with Figure 5.
Taking into account that Figure 6 represents a family of frozen orbits for inclinations of 90Β°, an example to illustrate the behavior of the frozen orbits for a case where the semimajor axis is 1861βkm is shown in Figures 7, 9, and 10. Figures 7β10 are obtained by numerical integration (using the software Maple) of the Lagrange planetary equations for different values of the eccentricity. Figure 8 shows frozen orbits for a case where the semimajor axis is 1838βkm.
We found frozen orbits with initial eccentricity equal to 0.02 that can be applied for missions of long period without the needs for doing orbit corrections for a period of very long time (larger than 4000 days). The initial orbit with appropriate characteristics can extend the time life and reduce maintenance costs. The set of initial conditions for frozen orbits with long lifetime is given by βkm, , Β°, Β°, and Β°.
In the literature, in general, frozen polar orbits are found with periapsis Β° [3, 31β33], however, for orbits with inclinations near 90Β°, the argument of periapsis moves to Β° [2, 31, 34]. These authors take into account only the zonal terms. In [34], there is an approach that compares the disturbing potential when a model of the gravitational field of high order taking into account the terms , , and is considered. For small inclinations, the three models show similar results for low-altitude orbits (see [34]), while for high inclination orbits the truncation only approximates roughly the case, giving frozen orbits that, in general, show lower eccentricities and shifting the zero eccentricity solutions to lower inclinations by several degrees. In this case, the potential truncated in the term provides a qualitative description only correct in the case of low inclinations. In [35], a numerical study considering models for the lunar gravity field of high order is applied to near-circular, low-altitude orbits. The behavior of the periapsis altitude is analyzed for some harmonics, for example, and .
Here, we truncate the potential of the lunar gravity field up to the term, to analyze the effect caused by this harmonic in particular [5]. When we consider the potential taking into account the terms β and (zonal and sectorial), the frozen orbit (Β°, see Figure 10) librates around the equilibrium point for Β° which is in agreement with the literature (see [3, 31, 32]). But, when we take into account the potential up to the term, we found frozen orbits that librate around the equilibrium point for Β°, which is in agreement with [5]. In [34], a frozen orbit considering the following initial conditions: βkm, , Β°, and Β°, is found, and it corresponds to a near-polar orbit.
Then, we develop the disturbing potential up to the term (β, ) to analyze what happens with the value of for a low-altitude, near-circular, near-polar orbit. Figure 11 shows the behavior for the diagram versus for a polar orbit which librates around the equilibrium point for the value of Β° (see [33] for more details with respect to the potential truncated up to the term.). In this case, the initial eccentricity varies from 0.001 to 0.004. The numerical values used for the harmonics up to is given in Table 3. Thus, we conclude that, because of the characteristics of the lunar gravity field, the potential truncated up to the produces the effect of the frozen orbit obtained that librates around Β° and not Β°.
The idea here is to consider a simplified potential to analyze the effect of the harmonics on the orbital elements and . Therefore, this study did not realize the full potential integrations, which can be viewed in [33, 34].
Another important factor that we must analyze is the value of the argument of the periapsis, that strongly affects the determination of the frozen orbit. As we previously showed, the values of the periapsis to assure the condition of frozen orbit are Β° or Β°, equally for the longitude of the ascending node that is Β° or Β°. Since that, in general, we find frozen orbits for the value of Β°. Figures 12β15 show the behavior of the eccentricity as a function of time for different values of the inclinations. Note that Figure 12 presents an orbit with lesser variation of the eccentricity for inclination of Β°, considering Β°. When we compare Figure 12 with Figure 13, we get that, for the value of Β°, the orbit presents great variations of the eccentricity. We also analyze the case for the initial eccentricity of the order of 10β3 and verify that the behavior of the orbit does not suffer the same effect as occurred for the case of the eccentricity of order 10β2. Here, the value of Β° or Β° has smaller effect on the amplitude of the eccentricity (see Figures 14 and 15).
4. Conclusions
In this paper, the dynamics of orbits around a planetary satellite (Moon), taking into account the non-uniform distribution of mass (β and ), was studied. Several analyses were performed, using the long-period disturbing potential obtained through the single-averaged model.
Considering the long-period disturbing potential, we found a family of frozen orbits for inclinations larger than 69Β° considering the model LP165P for the coefficients of the lunar gravity field. It was observed that the term contributed efficiently to obtain such orbits.
Low-altitude, near-polar orbits are very desirable for scientific missions to study natural satellites, such as the Moon. Due to the characteristics of the lunar gravity field, the potential truncated up to the produces the effect of the frozen orbit that librates around Β° and not Β°. Considering the long-period disturbing potential using the single-averaged method, we found near-circular, near-polar frozen orbits where this choice reduces the maintenance cost considerably.
In general, we conclude that, for the case of the Moon, it is necessary to consider a gravity model of high order, including the zonal and sectorial terms. Because the question of the coefficients are almost of the same order, this fact complicates the choice of the harmonic where the potential can be truncated and fact makes its choice a little arbitrary.
To study the lifetime of a lunar artificial satellite considering an analytical model, the zonal (β) and sectorial () terms of the harmonic coefficients should be taken into account. We observe that the frozen orbits could be strongly affected by terms of the potential containing the coefficient due to the equatorial ellipticity of the Moon () and by terms containing the argument of the periapsis and the longitude of the ascending node.
Appendix
The long-period disturbing potential up to the term is
where the values for the zonal (β) and tesseral () harmonics coefficients taking into account the two models are given in Table 2 [27] and Table 3 [28].
Acknowledgment
This work was accomplished with the support of the FAPESP under the contract no. 2007/04413-7 and 2006/04997-6, SP-Brazil, and CNPQ (300952/2008-2).