A method for averaging the perturbing function in Hill’s problem

A modified method for averaging the perturbing function in Hill’s problem is suggested. The averaging is performed in the revolution period of the satellite over the mean anomaly of its motion with a full allowance for a variation in the position of the perturbing body. At its fixed position, the semimajor axis of the satellite orbit during the revolution of the satellite is constant in view of the evolution equations, while the remaining orbital elements undergo secular and long-period perturbations. Therefore, when the motion of the perturbing body is taken into account, the semimajor axis of the satellite orbit undergoes the strongest perturbations. The suggested approach generalizes the averaging method in which only the linear (in time) term is included in the perturbing function. This method requires no expansion in powers of time. The described method is illustrated by calculating the perturbations of the semimajor axes for two distant satellites of Saturn, S/2000 S 1 and S/2000 S5. An approximate analytic solution is compared with the results of numerical integration of the averaged system of equations of motion for these satellites.     

Using wavelet approximation to study the influence of Jupiter on the orbital evolution of a distant satellite of Saturn

We suggest a nonstandard methodology for studying the influence of Jupiter on the secular orbital evolution of a distant satellite of Saturn. This influence is tangible only in short time spans near the times of the smallest separation between Jupiter and Saturn, i.e., when the heliocentric longitudes of the two planets coincide. These times are spaced about 20 years apart. To describe the jumplike behavior of perturbations, we suggest approximating the principal part of the perturbing function averaged over the satellite’s motion by a two-parameter exponential wavelet-type (burst) function. The subsequent averaging (smoothing) of the perturbing function allows us to eliminate the 20-year-period terms and obtain an approximate analytical solution in a special case of the problem. The results are illustrated by plots of the variations in the averaged perturbing function and the orbital eccentricity of Saturn’s outer satellite S/2000 S1, which is most strongly perturbed by Jupiter.     

The Translational–Rotational Motion of an Earth Spheroid Satellite

In this paper we consider the translational–rotational motion of a spheroid satellite in the gravitational field, taking into account the asphericity of the earth. The harmonic coefficients of the earth’s gravitational field are taken up to J ₄. The equations of motion are obtained in terms of the canonical elements of Delaunay-Andoyer. A first order solution is obtained using the perturbing technique of Lie series.     

On the motion of a mass point inside a rotating inhomogeneous ellipsoidal body

The plane motion of a mass point inside an inhomogeneous rotating ellipsoidal body with a homothetic density distribution is considered. The force function of the problem is expanded in terms of the ellipsoid's second eccentricities up to the fourth order, which are taken as small parameters. We derive an expression for the perturbing function and solve the equations of perturbed motion in orbital elements.     

Resonant and non-resonant gravity-gradient perturbations of a tumbling tri-axial satellite

Gravity-gradient perturbations of the attitude motion of a tumbling tri-axial satellite are investigated. The satellite center of mass is considered to be in an elliptical orbit about a spherical planet and to be tumbling at a frequency much greater than orbital rate. In determining the unperturbed (free) motion of the satellite, a canonical form for the solution of the torque-free motion of a rigid body is obtained. By casting the gravity-gradient perturbing torque in terms of a perturbing Hamiltonian, the long-term changes in the rotational motion are derived. In particular, far from resonance, there are no long-period changes in the magnitude of the rotational angular momentum and rotational energy, and the rotational angular momentum vector precesses abound the orbital angular momentum vector.At resonance, a low-order commensurability exists between the polhode frequency and tumbling frequency. Near resonance, there may be small long-period fluctuations in the rotational energy and angular momentum magnitude. Moreover, the precession of the rotational angular momentum vector about the orbital angular momentum vector now contains substantial long-period contributions superimposed on the non-resonant precession rate. By averaging certain long-period elliptic functions, the mean value near resonance for the precession of the rotational angular momentum vector is obtained in terms of initial conditions.     

Satellite Relative Motion Using Differential Equinoctial Elements

For precise control, to minimize the fuel consumption, and to maximize the lifetime of satellite formations a precise analytic solution is needed for the relative motion of satellites. Based on the relationship between the relative states and the differential orbital elements, the state transition matrix for the linearized relative motion that includes the effects due to the reference orbit eccentricity and the gravitational perturbations is derived. This method is called the Geometric Method. To avoid any singularities at zero eccentricity and zero inclination, equinoctial variables are used to derive the relative motion state transition matrices for both mean and osculating elements. This approach can be extended easily to include other perturbing forces.     

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.     

A Perturbative Treatment of The Co-Orbital Motion

We develop a formalism of the non-singular evaluation of the disturbing function and its derivatives with respect to the canonical variables. We apply this formalism to the case of the perturbed motion of a massless body orbiting the central body (Sun) with a period equal to that of the perturbing (planetary) body. This situation is known as the co-orbital motion, or equivalently, as the 1/1 mean motion commensurability. Jupiter's Trojan asteroids, Earth's co-orbital asteroids (e.g., (3753) Cruithne, (3362) Khufu), Mars' co-orbital asteroids (e.g., (5261) Eureka), and some Jupiter-family comets are examples of the co-orbital bodies in our solar system. Other examples are known in the satellite systems of the giant planets. Unlike the classical expansions of the disturbing function, our formalism is valid for any values of eccentricities and inclinations of the perturbed and perturbing body. The perturbation theory is used to compute the main features of the co-orbital dynamics in three approximations of the general three-body model: the planar-circular, planar-elliptic, and spatial-circular models. We develop a new perturbation scheme, which allows us to treat cases where the classical perturbation treatment fails. We show how the families of the tadpole, horseshoe, retrograde satellite and compound orbits vary with the eccentricity and inclination of the small body, and compute them also for the eccentricity of the perturbing body corresponding to a largely eccentric exoplanet's orbit.This revised version was published online in October 2005 with corrections to the Cover Date.     

On the analytic lunar and solar perturbations of a near Earth satellite

The disturbing function of the Moon (Sun) is expanded as a sum of products of two harmonic functions, one depending on the position of the satellite and the other on the position of the Moon (Sun). The harmonic functions depending on the position of the perturbing body are developed into trigonometric series with the ecliptic elementsl, l′, F, D and Γ of the lunar theory which are nearly linear with respect to time. Perturbation of elements are in the form of trigonometric series with the ecliptic lunar elements and the equatorial elements ω and Ω of the satellite so that analytic integration is simple and the results accurate over a long period of time.     

Particular solutions of the singly averaged Hill problem

We consider the case of averaging the perturbing function of the Hill problem over the fastest variable, the mean anomaly of the satellite. In integrable special cases, we found solutions to the evolutionary system of equations in elements.     

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.     

On the true circular orbit of a satellite

The osculating orbit of a planetary satellite moving in the equatorial plane of the central body under the influence of a rotational symmetric perturbation force is elliptical in first order approximation even if the true orbit is always circular. The satellite motion is influenced by a resonance effect due to this perturbing force. An inclined true satellite orbit cannot be circular.     

任意倾角的共轨运动研究

共轨运动天体与摄动天体的半长径相同,处于1:1平运动共振中.太阳系内多个行星的特洛伊天体即为处于蝌蚪形轨道的共轨运动天体,其中一些高轨道倾角特洛伊天体的轨道运动与来源仍未被完全理解.利用一个新发展的适用于处理1:1平运动共振的摄动函数展开方式,对三维空间中的共轨运动进行考察,计算不同初始轨道根数情况下共轨轨道的共振中心、共振宽度,分析轨道类型与初始轨道根数的关系.并将分析方法所得结果与数值方法的结果相互比较验证,得到了广阔初始轨道根数空间内共轨运动的全局图景.     

Secular evolution of the orbits of hypothetical satellites of Uranus

The problem of the secular perturbations of the orbit of a test satellite with a negligible mass caused by the joint influence of the oblateness of the central planet and the attraction by its most massive (or main) satellites and the Sun is considered. In contrast to the previous studies of this problem, an analytical expression for the full averaged perturbing function has been derived for an arbitrary orbital inclination of the test satellite. A numerical method has been used to solve the evolution system at arbitrary values of the constant parameters and initial elements. The behavior of some set of orbits in the region of an approximately equal influence of the perturbing factors under consideration has been studied for the satellite system of Uranus on time scales of the order of tens of thousands of years. The key role of the Lidov–Kozai effect for a qualitative explanation of the absence of small bodies in nearly circular equatorial orbits with semimajor axes exceeding ~1.8 million km has been revealed.     

Equations of motion of the restricted problem of three bodies with variable mass

In this paper we have deduced the differential equations of motion of the restricted problem of three bodies with decreasing mass, under the assumption that the mass of the satellite varies with respect to time. We have applied Jeans law and the space time transformation contrast to the transformation of Meshcherskii. The space time transformation is applicable only in the special casen=1,k=0,q=1/2. The equations of motion of our problem differ from the equations of motion of the restricted three body problem with constant mass only by small perturbing forces.     

A semi-analytical method to study perturbed rotational motion

A semi-analytical method is presented to study the system of differential equations governing the rotational motion of an artificial satellite. Gravity gradient and non gravitational torques are considered. Operations with trigonometric series were performed using an algebraic manipulator. Andoyer's variables are used to describe the rotational motion. The osculating elements are transformed analytically into a mean set of elements. As the differential equations in the mean elements are free of fast frequency terms, their numerical integration can be performed using a large step size.     

General relativity and satellite orbits: The motion of a test particle in the Schwarzschild metric

The motion of a satellite with negligible mass in the Schwarzschild metric is treated as a problem in Newtonian physics. The relativistic equations of motion are formally identical with those of the Newtonian case of a particle moving in the ordinary inverse-square law field acted upon by a disturbing function which varies asr ^?3. Accordingly, the relativistic motion is treated with the methods of celestial mechanics. The disturbing function is expressed in terms of the Keplerian elements of the orbit and substituted into Lagrange's planetary equations. Integration of the equations shows that a typical Earth satellite with small orbital eccentricity is displaced by about 17 cm from its unperturbed position after a single orbit, while the periodic displacement over the orbit reaches a maximum of about 3 cm. Application of the equations to the planet Mercury gives the advance of the perihelion and a total displacement of about 85 km after one orbit, with a maximum periodic displacement of about 13 km.     

Periodic solutions of the singly averaged hill problem

Based on the ideas of Lyapunov’s method, we construct a family of symmetric periodic solutions of the Hill problem averaged over the motion of a zero-mass point (a satellite). The low eccentricity of the satellite orbit and the sine of its inclination to the plane of motion of the perturbing body are parameters of the family. We compare the analytical solution with numerical solutions of the averaged evolutionary system and the rigorous (nonaveraged) equations of the restricted circular three-body problem.     

Perturbations in close binary systems produced by the tides lagging in latitude

The aim of this investigation is to present the periodic and secular perturbations of the orbital elements of close binary systems due to tidal lag in latitude. The variational equations of the problem of plane motion will be set up in terms of the rectengular componentsR, S, andW of the disturbing accelerations. These equations are highly nonlinear with respect to the orbital elements and we present analytic approximations to the effects produced by the perturbing acceleration due to dynamical tides lagging in latitude. The perturbed elements of the orbit have been expressed by means of Hansen coefficients in the compact form of summations.     

Expansion of the perturbing function of a satellite restricted four-body problem

A satellite four-body problem is the problem of motion of an artificial satellite of a planet in a region of the space where perturbations due to the gravitational field of the planet are of the same order as perturbations due to influences of two perturbing bodies. In this paper an expansion of the perturbing function into a Fourier series in terms of angular Keplerian elements ( _j , _j ,M _j :j=0,1,2) (designations are standard) is obtained taking into account a sharp commensurability of the typen/ ₀=(p+q)/p (n is the mean motion of the artificial satellite and ₀ is the angular velocity of rotation of the planet,p andq are integers).The coefficients of the Fourier series are the functions of the positional Keplerian elements (a _j ,e _j ,i _j ;j=0, 1, 2) (designations are standard) and, in particular, are series in terms ofe _j that, generally speaking, can be written out to an accuracy ofe _j ¹⁹ .The expansion obtained can be used for the construction of a semianalytical theory of motion of resonant satellites on the basis of conditionally periodic solutions of the restricted four-body problem.