The circular restricted three-body problem in curvilinear coordinates

We derive the exact equations of motion for the circular restricted three-body problem in cylindrical curvilinear coordinates together with a number of useful analytical relations linking curvilinear coordinates and classical orbital elements. The equations of motion can be seen as a generalization of Hill’s problem after including all neglected nonlinear terms. As an application of the method, we obtain a new expression for the averaged third-body disturbing function including eccentricity and inclination terms. We employ the latter to study the dynamics of the guiding center for the problem of circular coorbital motion providing an extension of some results in the literature.     

The influence of the great inequality on the secular disturbing function of the planetary system

This paper derives the contributionF ₂ ^* by the great inequality to the secular disturbing function of the principal planets. Andoyer's expansion of the planetary disturbing function and von Zeipel's method of eliminating the periodic terms is employed; thereby, the corrected secular disturbing function for the planetary system is derived. An earlier solution suggested by Hill is based on Leverrier's equations for the variation of elements of Jupiter and Saturn and on the semi-empirical adjustment of the coefficients in the secular disturbing function. Nowadays there are several modern methods of eliminating periodic terms from the Hamiltonian and deriving a purely secular disturbing function. Von Zeipel's method is especially suitable. The conclusion is drawn that the canonicity of the equations for the secular variation of the heliocentric elements can be preserved if there be retained, in the secular disturbing function, terms only of the second and fourth order relative to the eccentricity and inclinations.The Krylov-Bogolubov method is suggested for eliminating periodic terms, if it is desired to include the secular perturbations of the fifth and higher order in the heliocentric elements. The additional part of the secular disturbing functionF ₂ ^* derived in this paper can be included in existing theories of the secular effects of principal planets. A better approach would be to preserve the homogeneity of the theory and rederive all the secular perturbations of principal planets using Andoyer's symbolism, including the part produced by the great inequality.     

Some orbital characteristics of lunar artificial satellites

In this paper we present an analytical theory with numerical simulations to study the orbital motion of lunar artificial satellites. We consider the problem of an artificial satellite perturbed by the non-uniform distribution of mass of the Moon and by a third-body in elliptical orbit (Earth is considered). Legendre polynomials are expanded in powers of the eccentricity up to the degree four and are used for the disturbing potential due to the third-body. We show a new approximated equation to compute the critical semi-major axis for the orbit of the satellite. Lie-Hori perturbation method up to the second-order is applied to eliminate the terms of short-period of the disturbing potential. Coupling terms are analyzed. Emphasis is given to the case of frozen orbits and critical inclination. Numerical simulations for hypothetical lunar artificial satellites are performed, considering that the perturbations are acting together or one at a time.     

Construction of a consistent semianalytic theory of a planetary or moon orbiter perturbed by a third body

The third body perturbation of an orbiter of a planet or moon is considered. A very convenient form of the Lagrange equations is given allowing an easy derivation of the various terms of the expansion of the perturbed elements. A careful analysis of the order of magnitude of these terms indicates which ones are required for a consistent theory. It follows that in many practical cases the main term of the disturbing function has to be carried to the second order of the perturbation theory.Paper presented at the 1981 Oberwolfach Conference on Mathematical Methods in Celestial Mechanics.Dedicated to V. Szebehely on the occasion of his 60th birthday     

Nodal regression of the quadrantid meteor stream: An analytic approach

The mean orbit of the Quadrantid meteor stream has a high eccentricity and inclination with an aphelion close to the orbit of Jupiter. The nodal regression rate, a quantity which has been well determined from observations, cannot be calculated with sufficient accuracy using standard low-order expansions of the disturbing function. By using a high-order expansion of the disturbing function we show how the behavior of the longitude of ascending node of the Quadrantid stream is a result of both secular and resonant effects. Our analysis illustrates how the proximity of the stream's orbit to the 2: 1 commensurability with Jupiter dominates the short-term variations in orbital elements.     

An analytical theory for the orbit of nereid

A complete analytical dynamic theory for the motion of Nereid has been constructed, accurate to approximately 0.01 arc second over several hundred years. The solution uses the Lie transform approach advanced by Deprit and is consistent with respect to the magnitudes of the disturbing functions, including all perturbations to an accuracy of 10^–8 relative to the two-body potential (oblateness and third-body). Multiple short-period variables in the third-body perturbations are related via the ratio of their mean motions, reducing the number of independent variables. Extensive use is made of expansions giving trigonometric functions of the true anomaly as analytical Fourier series in the mean anomaly. Initial constants and mass parameters come from the data obtained during the Voyager II encounter with Neptune in 1989.     

Minor planet short periodic perturbations: The indirect part of the disturbing function

Leverrier's development of the indirect part of the disturbing function has been extended to include terms up to degree 4 in eccentricity and inclination; the resulting series has been expressed with respect to a fixed plane, and in a computer readable form (a list of integers). Tests have been performed for the relative significance of the terms of degrees 2, 3 and 4, and estimates have been obtained for the accuracy of the short periodic perturbations of a minor planet, and of the corresponding mean orbital elements. It was found that: (i) even in extreme cases, the indirect part of the disturbing function gives rise to very small short periodic perturbations; (ii) bodies of very high eccentricity/inclination and those close to mean motion resonances are most significantly affected; (iii) indirect perturbations for minor planets can be computed up to the degree 2 terms only, without any significant loss of accuracy; and (iv) higher degree indirect perturbations appear to be important only for their contribution to the long periodic effects of higher order (with respect to the perturbing mass).     

A note on a Lagrangian formulation for motion about the collinear points

A lagrangian formulation for the three-dimensional motion of a satellite in the vicinity of the collinear points of the circular-restricted problem is reconsidered. It is shown that the influence of the primaries can be expressed in the form of two third-body disturbing functions. By use of this approach, the equations for the Lagrangian and for the motion itself are readily developed into highly compact expressions. All orders of the non-linear developments are shown to be easily obtainable using well-known recursive relationships. The resulting forms for these equations are well suited for use in the initial phase of canonical or non-canonical investigations.     

Long-term perturbations due to a disturbing body in elliptic inclined orbit

In the current study, a double-averaged analytical model including the action of the perturbing body’s inclination is developed to study third-body perturbations. The disturbing function is expanded in the form of Legendre polynomials truncated up to the second-order term, and then is averaged over the periods of the spacecraft and the perturbing body. The efficiency of the double-averaged algorithm is verified with the full elliptic restricted three-body model. Comparisons with the previous study for a lunar satellite perturbed by Earth are presented to measure the effect of the perturbing body’s inclination, and illustrate that the lunar obliquity with the value 6.68^∘ is important for the mean motion of a lunar satellite. The application to the Mars-Sun system is shown to prove the validity of the double-averaged model. It can be seen that the algorithm is effective to predict the long-term behavior of a high-altitude Martian spacecraft perturbed by Sun. The double-averaged model presented in this paper is also applicable to other celestial systems.     

Satellite motion in a non-singular gravitational potential

We study the effects of a non-singular gravitational potential on satellite orbits by deriving the corresponding time rates of change of its orbital elements. This is achieved by expanding the non-singular potential into power series up to second order. This series contains three terms, the first been the Newtonian potential and the other two, here R ₁ (first order term) and R ₂ (second order term), express deviations of the singular potential from the Newtonian. These deviations from the Newtonian potential are taken as disturbing potential terms in the Lagrange planetary equations that provide the time rates of change of the orbital elements of a satellite in a non-singular gravitational field. We split these effects into secular, low and high frequency components and we evaluate them numerically using the low Earth orbiting mission Gravity Recovery and Climate Experiment (GRACE). We show that the secular effect of the second-order disturbing term R ₂ on the perigee and the mean anomaly are 4″.307×10⁻⁹/a, and −2″.533×10⁻¹⁵/a, respectively. These effects are far too small and most likely cannot easily be observed with today’s technology. Numerical evaluation of the low and high frequency effects of the disturbing term R ₂ on low Earth orbiters like GRACE are very small and undetectable by current observational means.     

Lunar perturbations of artificial satellites of the earth

Lunisolar perturbations of an artificial satellite for general terms of the disturbing function were derived by Kaula (1962). However, his formulas use equatorial elements for the Moon and do not give a definite algorithm for computational procedures. As Kozai (1966, 1973) noted, both inclination and node of the Moon's orbit with respect to the equator of the Earth are not simple functions of time, while the same elements with respect to the ecliptic are well approximated by a constant and a linear function of time, respectively. In the present work, we obtain the disturbing function for the Lunar perturbations using ecliptic elements for the Moon and equatorial elements for the satellite. Secular, long-period, and short-period perturbations are then computed, with the expressions kept in closed form in both inclination and eccentricity of the satellite. Alternative expressions for short-period perturbations of high satellites are also given, assuming small values of the eccentricity. The Moon's position is specified by the inclination, node, argument of perigee, true (or mean) longitude, and its radius vector from the center of the Earth. We can then apply the results to numerical integration by using coordinates of the Moon from ephemeris tapes or to analytical representation by using results from lunar theory, with the Moon's motion represented by a precessing and rotating elliptical orbit.     

Lunar frozen orbits revisited

Lunar frozen orbits, characterized by constant orbital elements on average, have been previously found using various dynamical models, incorporating the gravitational field of the Moon and the third-body perturbation exerted by the Earth. The resulting mean orbital elements must be converted to osculating elements to initialize the orbiter position and velocity in the lunar frame. Thus far, however, there has not been an explicit transformation from mean to osculating elements, which includes the zonal harmonic \(J_2\), the sectorial harmonic \(C_{22}\), and the Earth third-body effect. In the current paper, we derive the dynamics of a lunar orbiter under the mentioned perturbations, which are shown to be dominant for the evolution of circumlunar orbits, and use von Zeipel’s method to obtain a transformation between mean and osculating elements. Whereas the dynamics of the mean elements do not include \(C_{22}\), and hence does not affect the equilibria leading to frozen orbits, \(C_{22}\) is present in the mean-to-osculating transformation, hence affecting the initialization of the physical circumlunar orbit. Simulations show that by using the newly-derived transformation, frozen orbits exhibit better behavior in terms of long-term stability about the mean values of eccentricity and argument of periapsis, especially for high orbits.     

Explicit Semianalytical Theory of Asteroid Motion

A new semianalytical theory of asteroid motion is presented. The theory is developed on the basis of Kaula's expansion of the disturbing function including terms up to the second order with respect to the masses of disturbing bodies. The theory is constructed in explicit form that gives the possibility to study separately the influence of different perturbations in the dynamics of minor planets. The mean-motion resonances with major planets as well as mixed three-body resonances can also be taken into account. For the non-resonant case the formulas obtained can be used for deriving the second transformation to calculate the proper elements of an asteroid orbit in closed form with respect to inclinations and eccentricities. This revised version was published online in July 2006 with corrections to the Cover Date.     

Averaging the Equations of a Planetary Problem in an Astrocentric Reference Frame

A system of averaged equations of planetary motion around a central star is constructed. An astrocentric coordinate system is used. The two-planet problem is considered, but all constructions are easily generalized to an arbitrary number N of planets. The motion is investigated in modified (complex) Poincarécanonical elements. The averaging is performed by the Hori–Deprit method over the fast mean longitudes to the second order relative to the planetary masses. An expansion of the disturbing function is constructed using the Laplace coefficients. Some terms of the expansion of the disturbing function and the first terms of the expansion of the averaged Hamiltonian are given. The results of this paper can be used to investigate the evolution of orbits with moderate eccentricities and inclinations in various planetary systems.     

Kerr场中缓慢粒子的运动

本文用后－后牛顿近似讨论Ｋｅｒｒ场中缓慢粒子的运动，我们用Ｂｏｙｅｒ－Ｌｉｎｄｑｕｉｓｔ坐标，导出试验粒子的运动方程，把它与有心力场中粒子作二体运动之球坐标形式下的运动方程对比，得出由于Ｋｅｒｒ场的作用而引起的试验粒子的等效摄动加速度，利用球面三角公式把它换算到行星运动摄动方程的形状，对摄动方程进行积分，我们得出了试验粒子绕中心天体运动一周后粒子轨道根数的变化以及单位时间中轨道根数的平均变化，运用     

On the perturbations of a close-Earth satellite due to lunar inequalities

The lunar disturbing function for a close-Earth satellite is expressed as a sum of products of harmonics of the satellite's position and harmonics of the Moon's position, and the latter are expanded about a rotating and precessing elliptic orbit inclined to the ecliptic. The deviations of the Moon from this approximate orbit are computed from Brown's lunar theory andthe perturbations in satellite orbital elements due to these inequalities are derived. Numerical calculations indicate that several perturbations in the position of the satellite's node and perigee have magnitudes on the order of one meter.The author is supported in part by a National Science Foundation Graduate Fellowship.     

Analytical methods for the orbits of artificial satellites of the Moon

The motion of a close artificial satellite of the Moon is considered. The principal perturbations taken into account are caused by the nonsphericity of the Moon and the attraction of the Earth and the Sun. To begin with, the expansions of the disturbing functions due to the nonsphericity of the primary body and the action of the disturbing mass-point body have been derived. The second expansion is produced in terms of the Keplerian elements of a satellite and the spherical coordinates of the disturbing body. Both expansions are valid for an arbitrary reference plane. The motion of a satellite of the Moon is studied in the selenocentric coordinate system referred to the Lunar equator and rotating with respect to the fixed ecliptic system. However, the coordinate exes in the equatorial plane are chosen so that the angular speed of rotation of the system is small. The motion of the satellite is described by means of the contact elements which enable one to utilize the conventional Lagrange's planetary equations and may be regarded as the generalization of the notion of the osculating elements to the case of the disturbing function depending not only o the coordinates and the time but on the velocities as well. Two methods are proposed to represent the motion of Lunar satellites over long intervals of time: the von Zeipel method and the Euler method of analytical integration with application of the variation-of-elements technique at every step of integration. The second method is exposed in great detail.Presented at the Meeting of Commission 7 of the IAU on Analytical Methods for the Orbits of Artificial Celestial Objects 14-th General Assembly of the IAU, Brighton, 1970.     

A numerical investigation of secular terms of the planetary disturbing function

The secular terms of the planetary disturbing function are given, after elimination of short period terms by von Zeipel's transformation. The adequacy of this expansion up to terms of eighth order in the inclination and eccentricity is investigated by numerical processes, as a function of the Keplerian elementsa, e andi. The eccentricityé of the outer planet, is taken equal to zero. It is concluded that for values ofi which are not small the inclusion of additional terms in the expression for the disturbing function, results to drastic changes in its values, while larger values ofe do not have an equaly large effect on the disturbing function.     

The impact of Earth’s shadow on the long-term evolution of space debris

Direct solar radiation pressure and Earth’s shadow crossings are known to be responsible for short-term variations of space debris orbital elements, the higher the area-to-mass ratio the larger the perturbation. Nevertheless, existing studies have always been performed on periods of time shorter than 150 years. Considering longer time scales of the order of a 1000 years, this paper focuses on the long-term periodic evolution of space debris trajectories caused by successive Earth’s shadow crossings. Other perturbations as the geopotential and third-body gravitational attractions obviously play a role and compete with the one which is described in this paper. Symplectic numerical propagations and new (semi-)analytical models are developed to identify a frequency associated to shadow entry and exit eccentric anomalies. It is shown that Earth’s shadow is responsible for large deviations from the initial orbital elements, even on shorter period of times, and that this effect increases along with the area-to-mass ratio.     

A set modified equinoctial orbit elements

Modified equinoctial elements are introduced which are suitable for perturbation analysis of all kinds of orbit. Equations of motion in Lagrangian and Gaussian forms are derived. Identities connecting the partial derivatives of the disturbing function with respect to equinoctial elements are established. Numerical comparisons of the evolution of a perturbed, highly eccentric, elliptic orbit analysed in equinoctial elements and by Cowell's method show satisfactory agreement.