General presentation of secular part of the perturbing function of mutual attraction in the satellite system of a planet

A special presentation of secular part of the perturbing function of mutual attraction in the satellite system is presented. In contrast to known ones, it is given in general analytical form for any ratio between semimajor axes for perturbed and perturbing satellites.     

Constructive analytical solution of the evolution hill problem

A new analytical solution of the system of differential equations describing secular perturbations and long-period solar perturbations of mean orbits of outer satellites of giant planets was obtained. As distinct from other solutions, the solution constructed using von Zeipel’s method approximately takes into account, in the secular part of the perturbing function, the totality of fourth order with respect to the small parameter m of the ratio of the mean motions of the primary planet and the satellite. This enables us to describe more accurately the evolution of satellite orbits with large apocentric distances, which in the course of evolution may exceed the halved radius of the Hill sphere of the planet with respect to the Sun. Among these are the orbits of the two outermost Neptunian satellites N10 (Psamathe) and N13 (Neso). For these satellites, the parameter m amounts to 0.152 and 0.165, respectively. Different from a purely analytical solution, the proposed solution requires preliminary calculations for each satellite. More precisely, in doing so, we need to construct some simple functions to approximate more complex ones. This is why we use the phrase “constructive analytical.” To illustrate the solution, we compare it with the results of the numerical integration of the strict motion equations of the satellites N10 and N13 over time intervals 5–15 thousand years.     

On the classical expansion of the perturbing function in individual orbital elements

Starting from the classical expansion of the perturbing function in the three‐body problem, the transformation to individual orbital elements is performed in principle up to any degree in small parameters. Some corrections to the results presented in the well‐known article by Yuasa on secular perturbations of asteroids are given. Consequences for the expansion of the indirect part of the perturbing function are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

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.     

第三体摄动分析解的一种表达式

在太阳系中,大行星、小行星和卫星（包括自然卫星和人造卫星）等对应的运动问题,都可以处理成受摄二体问题,而摄动源又多为第三体,作为第三体的摄动天体,有的比运动天体离中心天体近,有的则相反,前者称为内摄内体,全者则称为外摄天体,对一个具体的运动天体,可以同时出现这两个摄动天体,但是,只要运动天体与摄动天体的轨道都建立在以中心天体（质心）为坐标原点的同一坐标系内,那么在一定条件下（即除运动天体与摄动天体     

A secular theory of coplanar, non-resonant planetary system

We present the secular theory of coplanar N -planet system, in the absence of mean motion resonances between the planets. This theory relies on the averaging of a perturbation to the two-body problem over the mean longitudes. We expand the perturbing Hamiltonian in Taylor series with respect to the ratios of semimajor axes which are considered as small parameters, without direct restrictions on the eccentricities. Next, we average out the resulting series term by term. This is possible thanks to a particular but in fact quite elementary choice of the integration variables. It makes it possible to avoid Fourier expansions of the perturbing Hamiltonian. We derive high-order expansions of the averaged secular Hamiltonian (here, up to the order of 24) with respect to the semimajor axes ratio. The resulting secular theory is a generalization of the octupole theory. The analytical results are compared with the results of numerical (i.e. practically exact) averaging. We estimate the convergence radius of the derived expansions, and we propose a further improvement of the algorithm. As a particular application of the method, we consider the secular dynamics of three-planet coplanar system. We focus on stationary solutions in the HD 37124 planetary system.     

On the stability of circular ‘asteroid’ orbits in anN-planetary system

The restricted problem of the motion of a point of negligible mass (asteroid) in anN-planetary system is considered. It is assumed that all the planets move about the central body (Sun) along circular orbits in the same plane and the mean motions of the asteroid and the planets are incommensurable. The asteroid orbit evolution is described as a first approximation by secular equations with the perturbing function averaged by the mean longitudes of the asteroid and the planets. For small values of the asteroid orbit eccentricity an expression for the secular part of the perturbing function has been obtained. This expression holds for the arbitrary values of the asteroid orbit semiaxis which are different from those of the planet orbit radii. The stability of the asteroid circular orbits in a linear approximation with respect to the eccentricity is studied. The critical inclinations for a Solar system model are calculated.     

The construction of a fifth-order Jupiter-Saturn motion: Analytical theory

We construct a fifth-order with respect to masses Jupiter-Saturn secular theory by Hori-Lie canonical technique. The J-S Hamiltonian includes both parts of the perturbing function. The influence of the 2:5 critical terms is taken into consideration. The Jacobi-Radau system of origins is adopted and the theory is expressed in terms of the Poincaré canonical variables.     

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.     

Free and forced obliquities of the Galilean satellites of Jupiter

The obliquity, or angular separation between orbit normal and spin pole, is an important parameter for the geodynamics of most Solar System bodies. Tidal dissipation has driven the obliquities of the Galilean satellites of Jupiter to small, but non-zero values. We present estimates of the free and forced obliquities of these satellites using a simple secular variation model for the orbits, and spin pole precession rate estimates based on gravity field parameters derived from Galileo spacecraft encounters. The free obliquity values are not well constrained by observations, but are presumed to be very small. The forced obliquity variations depend only on the orbital variations and the spin pole precession rate parameters, which are quite well known. These variations are large enough to influence spatial and temporal patterns of tidal dissipation and tidal stress.     

On the perturbations of the five outer planets by the four inner ones

The expressions given by Clemence were checked by the comparison of numerical values of their second differentials with the numerical values of the perturbing forces. Agreement was good in most cases, save that the use of the second differentials unduly magnified some terms rightly neglected by Clemence. In the appendix the work of Clemence is compared with that of Carpenter and the differences are found to be all smaller than the contributions from the secular variation of the eccentricity, which are so small that they were neglected by Clemence. The differences can be removed by considering the secular variation of the eccentricity and by making a small adjustment in the integration constants.     

Dynamics of close binary systems with mass exchange

The effects of non-isotropic ejection of mass from either component of a binary system on the orbital elements are studied, for the case of a small initial eccentricity of the relative orbit, when all the ejected mass falls on the other component. The problem is transformed to an equivalent two-body problem with isotropic variation of mass, plus a perturbing force which is a function of the intial conditions of ejection of the particles and their final, positions and velocities when they fall on the surface of the other star. The variation of the orbital elements are derived. It is shown that, to first-order terms in the eccentricity, the secular change of the semimajor axis is equal to the one corresponding to the case of zero initial eccentricity. On the contrary, the secular change of the eccentricity is smaller and it depends on the variations of mass ejection due to the finite eccentricity.     

Nonlinear theory of secular perturbations of satellites of an oblate planet

Anonlinear analytical theory of secular perturbations in the problem of the motion of a systemof small bodies around a major attractive center has been developed. Themutual perturbations of the satellites and the influence of the oblateness of the central body are taken into account in the model. In contrast to the classical Laplace-Lagrange theory based on linear equations for Lagrange elements, the third-degree terms in orbital eccentricities and inclinations are taken into account in the equations. The corresponding improvement of the solution turns out to be essential in studying the evolution of orbits over long time intervals. A program inC has been written to calculate the corrections to the fundamental frequencies of the solution and the third-degree secular perturbations in orbital eccentricities and inclinations. The proposed method has been applied to investigate the motion of the major Uranian satellites. Over time intervals longer than 100 years, allowance for the nonlinear terms in the equations is shown to give corrections to the coordinates of Miranda on the order of the orbital eccentricity, which is several thousand kilometers in linear measure. For other satellites, the effect of allowance for the nonlinear terms turns out to be smaller. Obviously, when a general analytical theory of motion for the major Uranian satellites is constructed, the nonlinear terms in the equations for the secular perturbations should be taken into account.     

Possible Effect of the Earth’s Inertial Induction on the Orbital Decay of LAGEOS

The theory of velocity dependent inertial induction, based upon extended Mach’s principle, has been able to generate many interesting results related to celestial mechanics and cosmological problems. Because of the extremely minute magnitude of the effect its presence can be detected through the motion of accurately observed bodies like Earth satellites. LAGEOS I and II are medium altitude satellites with nearly circular orbits. The motions of these satellites are accurately recorded and the past data of a few decades help to test many theories including the general theory of relativity. Therefore, it is hoped that the effect of the Earth’s inertial induction can have any detectable effect on the motion of these satellites. It is established that the semi-major axis of LAGEOS I is decreasing at the rate of 1.3 mm/d. As the atmospheric drag is negligible at that altitude, a proper explanation of the secular change has been wanting, and, therefore, this paper examines the effect of the Earth’s inertial induction effect on LAGEOS I. Past researches have established that Yarkovsky thermal drag, charged and neutral particle drag might be the possible mechanisms for this orbital decay. Inertial induction is found to generate a perturbing force that results in 0.33 mm/d decay of the semi major axis. Some other changes are also predicted and the phenomenon also helps to explain the observed changes in the orbits of a few other satellites. The results indicate the feasibility of the theory of inertial induction i.e. the dynamic gravitation phenomenon of the Earth on its satellites as a possible partial cause for orbital decay.     

A semi-analytical model of the Galilean satellites’ dynamics

The Galilean satellites’ dynamics has been studied extensively during the last century. In the past it was common to use analytical expansions in order to get simple models to integrate, but with the new generation of computers it became prevalent the numerical integration of very sophisticated and almost complete equations of motion. In this article we aim to describe the resonant and secular motion of the Galilean satellites through a Hamiltonian, depending on the slow angles only, obtained with an analytical expansion of the perturbing functions and an averaging operation. In order to have a model as near as possible to the actual dynamics, we added perturbations and we considered terms that in similar studies of the past were neglected, such as the terms involving the inclinations and the Sun’s perturbation. Moreover, we added the tidal dissipation into the equations, in order to investigate how well the model captures the evolution of the system.     

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).     

Secular love and tidal numbers of synchronously orbiting satellites

The secular Love and the secular tidal numbers have been computed for eight synchronously orbiting satellites in the solar system for which the triaxiality parameters and satellite-centric gravitational constant are available. Excepting Deimos the total mass of which should be first refined, the secular Love and tidal numbers are rather close to unity, as a rule. That is why, the centrifugal and tidal distortions can be assumed responsible for the actual figures of the synchronously orbiting satellites resulting from the primordial spheres, as well as, their static equilibrium nearly satisfied. The hypothesis of the origin of synchronously orbiting satellites by accretion in orbits is supported by the results obtained.     

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.     

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.