A type of second order linear ordinary differential equations with periodic coefficients for which the characteristic exponents have exact expressions

It is shown that the characteristic exponents can be exactly expressed for a type of second order linear ordinary differential equations with periodic coefficients (Hill's equation) which appear as the variational equations of certain periodic solutions of dynamical systems. Key points are the transformation of the equation to the Gauss' hypergeometric differential equation, and evaluation of the trace of monodormy matrix in the complex plane of the independent variable. Two simple examples are given for which the stability of periodic solutions is rigorously analyzed.     

Periodic attitudes and bifurcations of a rigid spacecraft in the second degree and order gravity field of a uniformly rotating asteroid

In this work, periodic attitudes and bifurcations of periodic families are investigated for a rigid spacecraft moving on a stationary orbit around a uniformly rotating asteroid. Under the second degree and order gravity field of an asteroid, the dynamical model of attitude motion is formulated by truncating the integrals of inertia of the spacecraft at the second order. In this dynamical system, the equilibrium attitude has zero Euler angles. The linearised equations of attitude motion are utilised to study the stability of equilibrium attitude. It is found that there are three fundamental types of periodic attitude motions around a stable equilibrium attitude point. We explicitly present the linear solutions around a stable equilibrium attitude, which can be used to provide the initial guesses for computing the true periodic attitudes in the complete model. By means of a numerical approach, three fundamental families of periodic attitudes are studied, and their characteristic curves, distribution of eigenvalues, stability curves and stability distributions are determined. Interestingly, along the characteristic curves of the fundamental families, some critical points are found to exist, and these points correspond to tangent and period-doubling bifurcations. By means of a numerical approach, the bifurcated families of periodic attitudes are identified. The natural and bifurcated families constitute networks of periodic attitude families.     

Attitude stability of a spinning symmetric satellite in a planar periodic orbit

This paper contains an analysis of the attitude stability of a spinning axisymmetric satellite whose mass center moves in any known planar periodic orbit of the restricted three-body problem while the spin axis remains normal to the orbit plane. A procedure based on Floquet theory is developed for constructing attitude instability charts, and examples of these are presented for two stable periodic orbits of the Earth-Moon system—one direct and one retrograde. The physical significance of these instability predictions is then explored by means of numerical integration of the full nonlinear equations of motion. Finally, an analysis based on averaging is performed, leading to approximate instability charts and indicating a possible connection between certain orbital-attitude resonance conditions and unstable attitude motions.     

The stability of periodic orbits in the three-body problem

A method is developed to study the stability of periodic motions of the three-body problem in a rotating frame of reference, based on the notion of surface of section. The method is linear and involves the computation of a 4×4 variational matrix by integrating numerically the differential equations for time intervals of the order of a period. Several properties of this matrix are proved and also it is shown that for a symmetric periodic motion it can be computed by integrating for half the period only.This linear stability analysis is used to study the stability of a family of periodic motions of three bodies with equal masses, in a rotating frame of reference. This family represents motion such that two bodies revolve around each other and the third body revolves around this binary system in the same direction to a distance which varies along the members of the family. It was found that a large part of the family, corresponding to the case where the distance of the third body from the binary system is larger than the dimensions of the binary system, represents stable motion. The nonlinear effects to the linear stability analysis are studied by computing the intersections of several perturbed orbits with the surface of sectiony ₃=0. In some cases more than 1000 intersections are computed. These numerical results indicate that linear stability implies stability to all orders, and this is true for quite large perturbations.     

On the effect of the eccentricity of a planetary orbit on the stability of satellite orbits

The effect of the eccentricity of a planet’s orbit on the stability of the orbits of its satellites is studied. The model used is the elliptic Hill case of the planar restricted three-body problem. The linear stability of all the known families of periodic orbits of the problem is computed. No stable orbits are found, the majority of them possessing one or two pairs of real eigenvalues of the monodromy matrix, while a part of a family with complex instability is found. Two families of periodic orbits, bifurcating from the Lagrangian points L₁, L₂ of the corresponding circular case are found analytically. These orbits are very unstable and the determination of their stability coefficients is not accurate, so we compute the largest Liapunov exponent in their vicinity. In all cases these exponents are positive, indicating the existence of chaotic motions     

Parametric excitation induced by solar pressure torque on the roll-yaw attitude motion of a gravity-gradient stabilized spacecraft

The parametric excitation of a gravity gradient stabilized spacecraft induced by the periodic solar pressure torque is discussed. The solar pressure torque in the linearized equations of motion appears as linear terms with periodic coefficients. The attitude stability is analyzed numerically through the calculation of the Floquet multiplier. The perturbation method is also applied to identify the instability condition analytically. It is made clear that the periodic solar pressure torque can destabilize the coupled roll and yaw attitude motion of the spacecraft. It is also shown that the conditions of parametric resonance are included in the gravity gradient stability condition. Nonlinear simulations are also carried out to verify the effect of the parametric resonance. The numerical simulation using actual parameters shows that the spacecraft inevitably experiences a large amplitude attitude motion due to the periodic solar pressure torque even if the gravity gradient stability condition is satisfied.     

Multiple resonance in one problem of the stability of the motion of a satellite relative to the center of mass

We investigate the stability of the periodic motion of a satellite, a rigid body, relative to the center of mass in a central Newtonian gravitational field in an elliptical orbit. The orbital eccentricity is assumed to be low. In a circular orbit, this periodic motion transforms into the well-known motion called hyperboloidal precession (the symmetry axis of the satellite occupies a fixed position in the plane perpendicular to the radius vector of the center of mass relative to the attractive center and describes a hyperboloidal surface in absolute space, with the satellite rotating around the symmetry axis at a constant angular velocity). We consider the case where the parameters of the problem are close to their values at which a multiple parametric resonance takes place (the frequencies of the small oscillations of the satellite’s symmetry axis are related by several second-order resonance relations). We have found the instability and stability regions in the first (linear) approximation at low eccentricities.     

Phase-space volume of regions of trapped motion: multiple ring components and arcs

The phase-space volume of regions of regular or trapped motion, for bounded or scattering systems with two degrees of freedom respectively, displays universal properties. In particular, sudden reductions in the phase-space volume or gaps are observed at specific values of the parameter which tunes the dynamics; these locations are approximated by the stability resonances. The latter are defined by a resonant condition on the stability exponents of a central linearly stable periodic orbit. We show that, for more than two degrees of freedom, these resonances can be excited opening up gaps, which effectively separate and reduce the regions of trapped motion in phase space. Using the scattering approach to narrow rings and a billiard system as example, we demonstrate that this mechanism yields rings with two or more components. Arcs are also obtained, specifically when an additional (mean-motion) resonance condition is met. We obtain a complete representation of the phase-space volume occupied by the regions of trapped motion.     

Periodic orbits of the elliptic restricted problem for the Sun-Jupiter-Saturn system

A systematic approach to generate periodic orbits in the elliptic restricted problem of three bodies in introduced. The approach is based on (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to periodic orbits in the elliptic restricted problem. Two families of periodic orbits of the elliptic restricted problem are found by this approach. The mass ratio of the primaries of these orbits is equal to that of the Sun-Jupiter system. The sidereal mean motions between the infinitesimal body and the smaller primary are in a 2:5 resonance, so as to approximate the Sun-Jupiter-Saturn system. The linear stability of these periodic orbits are studied as functions of the eccentricities of the primaries and of the infinitesimal body. The results show that both stable and unstable periodic orbits exist in the elliptic restricted problem that are close to the actual Sun-Jupiter-Saturn system. However, the periodic orbit closest to the actual Sun-Jupiter-Saturn system is (linearly) stable.     

On the periodic solutions and resonance of spinning satellites in near-circular orbits

Attitude motion of spinning axisymmetric satellites in presence of gravity-gradient and solar radiation pressure torques is studied analytically. The approximate closed-form solution developed for the nonlinear, nonautonomous, coupled fourth-order system proves to be an excellent tool in locating periodic solutions of the system in both circular and noncircular orbits. The variational stability of the periodic motion is examined using the Floquet theory. The resonance analysis suggests the existence of critical combinations of system parameters leading to large amplitude oscillations.

     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

The Sitnikov family and the associated families of 3D periodic orbits in the photogravitational RTBP with oblateness

We consider the photogravitational restricted three-body problem with oblateness and study the Sitnikov motions. The family of straight line oscillations exists only in the case where the primaries are of equal masses as in the classical Sitnikov problem and have the same oblateness coefficients and radiation factors. A perturbation method based on Floquet theory is applied in order to study the stability of the motion and critical orbits are determined numerically at which families of three-dimensional periodic orbits of the same or double period bifurcate. Many of these families are computed.     

Properties of linearized mappings associated with periodic orbits in the three-body problem

In this paper three results on the linearized mapping associated with the plane three body problem near a periodic orbit are established. It is first shown that linear stability of such an orbit is independent of initial position on the orbit and of coordinate system. Second, the relation of Hénon connecting the rates of change of rotation angle and period on an isoenergetic family of periodic orbits is proved, together with a similar relation for families of orbits closing exactly in a rotating coordinate system. Finally, a condition for a critical orbit is given which is applicable to any family of periodic orbits.     

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.     

On the long-term motion of Pluto

A modified periodic orbit of the third kind is introduced that is closely related to periodic orbits of the third kind as defined by Poincaré. It is shown that Pluto librates about the periodic orbit with apparent stability. This further explains the librational motion of the resonant argument of Pluto and the avoidance of a Pluto-Neptune close approach as found by Cohen and Hubbard and the long-term motion of Pluto and the librational motion of the perihelion as found by Williams and Benson. With libration about a periodic orbit, the numerical solution of Williams and Benson can be extrapolated to longer times in the past and future.     

High order normal form construction near the elliptic orbit of the Sitnikov problem

We consider the Sitnikov problem; from the equations of motion we derive the approximate Hamiltonian flow. Then, we introduce suitable action–angle variables in order to construct a high order normal form of the Hamiltonian. We introduce Birkhoff Cartesian coordinates near the elliptic orbit and we analyze the behavior of the remainder of the normal form. Finally, we derive a kind of local stability estimate in the vicinity of the periodic orbit for exponentially long times using the normal form up to 40th order in Cartesian coordinates.     

Hamiltonian Stability of Spin–Orbit Resonances in Celestial Mechanics

The behaviour of ‘resonances’ in the spin-orbit coupling in celestial mechanics is investigated in a conservative setting. We consider a Hamiltonian nearly-integrable model describing an approximation of the spin-orbit interaction. The continuous system is reduced to a mapping by integrating the equations of motion through a symplectic algorithm. We study numerically the stability of periodic orbits associated to the above mapping by looking at the eigenvalues of the matrix of the linearized map over the full cycle of the periodic orbit. In particular, the value of the trace of the matrix is related to the stability character of the periodic orbit. We denote by ε_* (p/q) the value of the perturbing parameter at which a given elliptic periodic orbit with frequency p/q becomes unstable. A plot of the critical function ε_* (p/q) versus the frequency at different orbital eccentricities shows significant peaks at the synchronous resonance (for low eccentricities) and at the synchronous and 3:2 resonances (at higher eccentricities) in good agreement with astronomical observations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Periodic motion around the triangular equilibrium points of the photogravitational restricted problem of three bodies

In this article the effect of radiation pressure on the periodic motion of small particles in the vicinity of the triangular equilibrium points of the restricted three body problem is examined. Second order parametric expansions are constructed and the families of periodic orbits are determined numerically for two sets of values of the mass and radiation parameters corresponding to the non-resonant and the resonant case. The stability of each orbit is also studied.