Numerical determination of stability regions for orbital motion in uniformly rotating second degree and order gravity fields

The stability of orbital motion about a uniformly rotating arbitrary second degree and order gravity field is investigated. A normalized form of the equations of motion are derived and analyzed. A numerical stability criteria is proposed and used to evaluate the stability of initially near-circular orbits in the equatorial plane of the body. Regions of stable and unstable motion are clearly delineated, and are seen to be strongly related to resonances between the mean motion and the body rotation rate.     

High-order solutions of invariant manifolds associated with libration point orbits in the elliptic restricted three-body system

High-order analytical solutions of invariant manifolds, associated with Lissajous and halo orbits in the elliptic restricted three-body problem (ERTBP), are constructed in this paper. The equations of motion of ERTBP in the pulsating synodic coordinate system have five equilibrium points, and the three collinear libration points as well as the associated center manifolds are unstable. In our calculation, the general solutions of the invariant manifolds associated with Lissajous and halo orbits around collinear libration points are expressed as power series of five parameters: the orbital eccentricity, two amplitudes corresponding to the hyperbolic manifolds, and two amplitudes corresponding to the center manifolds. The analytical solutions up to arbitrary order are constructed by means of Lindstedt–Poincaré method, and then the center and invariant manifolds, transit and non-transit trajectories in ERTBP are all parameterized. Since the circular restricted three-body problem (CRTBP) is a particular case of ERTBP when the eccentricity is zero, the general solutions constructed in this paper can be reduced to describe the dynamics around the collinear libration points in CRTBP naturally. In order to check the validity of the series expansions constructed, the practical convergence of the series expansions up to different orders is studied.     

Minimum control for spacecraft formations in a J2 perturbed environment

Large ΔV amounts are often required to maintain the relative geometry which is needed to implement a formation flying concept. A wise use of the orbital environment makes the orbit keeping phase easier, reducing the overall propellant consumption. A first important step in this direction is the selection of formation configurations and orbits which, while still satisfying the mission requirements, are less subject to perturbations resulting naturally in closed relative motion. Within this frame, a number of studies have been recently carried out in order to identify possible sets of invariant relative orbits under the effects of the Earth oblateness, a conservative force commonly referred to as J2 which is also the most important perturbation for Low Earth Orbit. These efforts clearly marked the difficulties connected with achieving genuine periodic relative motion under J2 effect, but they also showed the existence of a set of conditions on the orbital parameters which allow for quasi-periodic J2 trajectories. This paper presents these particular trajectories, by means of deeper theoretical explanations, showing the dependency of the shape of the relative configurations on the orbital inclination. Since the quasi-periodic trajectories cannot be written analytically, and moreover, they are very sensitive with respect to the initial conditions, difficulties arise when trying to exploit these paths as reference for the control of a formation. This paper proposes a novel approach to find, from the actual quasi periodic natural trajectories, minimal control periodic reference trajectories. Next, it evaluates quantitatively the amount of propellant which is needed to control a spacecraft formation along these paths. The choice of Hill’s classical circular projected configuration as a nominal trajectory is considered as a comparison, showing the clear advantages of the proposed guidance design, which assumes low-perturbed periodic reference orbits as nominal trajectories.     

Planetary dynamics in the system α Centauri: The stability diagrams

The stability of the motion of a hypothetical planet in the binary system ?? Cen A?CB has been investigated. The analysis has been performed within the framework of a planar (restricted and full) three-body problem for the case of prograde orbits. Based on a representative set of initial data, we have obtained the Lyapunov spectra of the motion of a triple system with a single planet. Chaotic domains have been identified in the pericenter distance-eccentricity plane of initial conditions for the planet through a statistical analysis of the data obtained. We have studied the correspondence of these chaotic domains to the domains of initial conditions that lead to the planet??s encounter with one of the binary??s stars or to the escape of the planet from the system. We show that the stability criterion based on the maximum Lyapunov exponent gives a more clear-cut boundary of the instability domains than does the encounterescape criterion at the same integration time. The typical Lyapunov time of chaotic motion is ??500 yr for unstable outer orbits and ??60 yr for unstable inner ones. The domain of chaos expands significantly as the initial orbital eccentricity of the planet increases. The chaos-order boundary has a fractal structure due to the presence of orbital resonances.     

The chaotic rotation of Hyperion

A plot of spin rate versus orientation when Hyperion is at the pericenter of its orbit (surface of section) reveals a large chaotic zone surrounding the synchronous spin-orbit state of Hyperion, if the satellite is assumed to be rotating about a principal axis which is normal to its orbit plane. This means that Hyperion's rotation in this zone exhibits large, essentially random variations on a short time scale. The chaotic zone is so large that it surrounds the ½ and 2 states, and libration in the 3/2 state is not possible. Stability analysis shows that for libration in the synchronous and ½ states, the orientation of the spin axis normal to the orbit plane is unstable, whereas rotation in the 2 state is attitude stable. Rotation in the chaotic zone is also attitude unstable. A small deviation of the principal axis from the orbit normal leads to motion through all angles in both the chaotic zone and the attitude unstable libration regions. Measures of the exponential rate of separation of nearby trajectories in phase space (Lyapunov characteristic exponents) for these three-dimensional motions indicate the the tumbling is chaotic and not just a regular motion through large angles. As tidal dissipation drives Hyperion's spin toward a nearly synchronous value, Hyperion necessarily enters the large chaotic zone. At this point Hyperion becomes attitude unstable and begins to tumble. Capture from the chaotic state into the synchronous or ½ state is impossible since they are also attitude unstable. The 3/2 state does not exist. Capture into the stable 2 state is possible, but improbable. It is expected that Hyperion will be found tumbling chaotically.     

Lagrangian coherent structures in the planar elliptic restricted three-body problem

This study investigates Lagrangian coherent structures (LCS) in the planar elliptic restricted three-body problem (ER3BP), a generalization of the circular restricted three-body problem (CR3BP) that asks for the motion of a test particle in the presence of two elliptically orbiting point masses. Previous studies demonstrate that an understanding of transport phenomena in the CR3BP, an autonomous dynamical system (when viewed in a rotating frame), can be obtained through analysis of the stable and unstable manifolds of certain periodic solutions to the CR3BP equations of motion. These invariant manifolds form cylindrical tubes within surfaces of constant energy that act as separatrices between orbits with qualitatively different behaviors. The computation of LCS, a technique typically applied to fluid flows to identify transport barriers in the domains of time-dependent velocity fields, provides a convenient means of determining the time-dependent analogues of these invariant manifolds for the ER3BP, whose equations of motion contain an explicit dependency on the independent variable. As a direct application, this study uncovers the contribution of the planet Mercury to the Interplanetary Transport Network, a network of tubes through the solar system that can be exploited for the construction of low-fuel spacecraft mission trajectories. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Optimal transfers between unstable periodic orbits using invariant manifolds

This paper presents a method to construct optimal transfers between unstable periodic orbits of differing energies using invariant manifolds. The transfers constructed in this method asymptotically depart the initial orbit on a trajectory contained within the unstable manifold of the initial orbit and later, asymptotically arrive at the final orbit on a trajectory contained within the stable manifold of the final orbit. Primer vector theory is applied to a transfer to determine the optimal maneuvers required to create the bridging trajectory that connects the unstable and stable manifold trajectories. Transfers are constructed between unstable periodic orbits in the Sun–Earth, Earth–Moon, and Jupiter-Europa three-body systems. Multiple solutions are found between the same initial and final orbits, where certain solutions retrace interior portions of the trajectory. All transfers created satisfy the conditions for optimality. The costs of transfers constructed using manifolds are compared to the costs of transfers constructed without the use of manifolds. In all cases, the total cost of the transfer is significantly lower when invariant manifolds are used in the transfer construction. In many cases, the transfers that employ invariant manifolds are three times more efficient, in terms of fuel expenditure, than the transfer that do not. The decrease in transfer cost is accompanied by an increase in transfer time of flight.     

Invariant shape solutions of the spinning three craft Coulomb tether problem

In this paper we study shape-preserving formations of three spacecraft, where the formation keeping forces arise from the electric charges deposed on each craft. Inspired by Lagrange’s 3-body problem, the general conditions that guarantee preservation of the geometric shape of the electrically charged formation are derived. While the classical collinear configuration is a solution to the problem, the equilateral triangle configuration is found to only occur with unbounded relative motion. The three collinear spacecraft problem is analyzed and the possible solutions are categorized based on the spacecraft mass–charge ratio. Precise statements on the number of solutions associated with each category are provided. Finally, a methodology is proposed to study boundedness of the collinear solution that is inspired by past understanding and results for the 3-body problem. Given the initial position and the velocity vectors of each craft along with the charges, analytical solutions are provided describing the resulting relative motion.     

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.     

Stability Domain and Invariant Manifolds of 2d Area-Preserving Diffeomorphisms

We study the stability domain of generic 2D area-preserving polynomialdiffeomorphisms. The starting point of our analysis is the study of thedistribution of stable and unstable fixed points. We show that the locationof fixed points and their stability type are linked to the degree of thepolynomial map. These results are based on a classification Theorem forplane automorphisms by Friedland and Milnor. Then we discuss the problem ofdetermining the domain in phase space where stable motion occurs. We showthat the boundary of the stability domain is given by the invariantmanifolds emanating from the outermost unstable fixed point of low period(one or two). This fact extends previous results obtained for reversiblearea-preserving polynomial maps of the plane. This analysis is based onanalytical arguments and is supported by the results of numericalsimulations.     

Out-of-plane equilibrium points and their stability in the Robe’s problem with oblateness and triaxiality

This paper examines the existence and stability of the out-of-plane equilibrium points of a third body of infinitesimal mass when the equations of motion are written in the three dimensional form under the set up of the Robe’s circular restricted three-body problem, in which the hydrostatic equilibrium figure of the first primary is an oblate spheroid and the second one is a triaxial rigid body under the full buoyancy force of the fluid. The existence of the out of orbital plane equilibrium points lying on the xz-plane is noticed. These points are however unstable in the linear sense.     

Role of the Rossiter-McLaughlin effect in the study of close binaries

This paper is devoted to binary stars belonging to the class of eclipsing-variable systems.Photometric and spectroscopic analysis of eclipses allows us to determine geometric parameters of the orbit and physical characteristics of stellar components as well as inclinations of stellar equators to the orbital plane. Estimations of inclinations can be obtained from measurement of the Rossiter-McLaughlin effect, which is discussed using examples of some eccentric binaries with an anomalous apsidal effect. Our task is to find the complete spectrum of solutions of the equation of apsidal motion, depending on the inclinations of the polar axes of the components to the orbital one for these systems, based on their individual spectroscopic and photometric observational data. The matrix of solutions allows us to select those pairs of polar inclinations that provide agreement with the observational apsidal period.     

Cylindrical isomorphic mapping applied to invariant manifold dynamics for Earth–Moon Missions

Several families of periodic orbits exist in the context of the circular restricted three-body problem. This work studies orbital motion of a spacecraft among these periodic orbits in the Earth–Moon system, using the planar circular restricted three-body problem model. A new cylindrical representation of the spacecraft phase space (i.e., position and velocity) is described, and allows representing periodic orbits and the related invariant manifolds. In the proximity of the libration points, the manifolds form a four-fold surface, if the cylindrical coordinates are employed. Orbits departing from the Earth and transiting toward the Moon correspond to the trajectories located inside this four-fold surface. The isomorphic mapping under consideration is also useful for describing the topology of the invariant manifolds, which exhibit a complex geometrical stretch-and-folding behavior as the associated trajectories reach increasing distances from the libration orbit. Moreover, the cylindrical representation reveals extremely useful for detecting periodic orbits around the primaries and the libration points, as well as the possible existence of heteroclinic connections. These are asymptotic trajectories that are ideally traveled at zero-propellant cost. This circumstance implies the possibility of performing concretely a variety of complex Earth–Moon missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining a suitable, convenient end-of-life strategy for spacecraft placed in any of the unstable orbits. The final disposal orbit is an externally confined trajectory, never approaching the Earth or the Moon, and can be entered by means of a single velocity impulse (of modest magnitude) along the right unstable manifold that emanates from the Lyapunov orbit at \(L_2\) .     

Interchange stability of a rapidly rotating magnetosphere

A rotation-dominated magnetosphere is unstable to magnetic flux-tube interchange motions if and only if the plasma content of a unit magnetic flux tube is a decreasing function of distance from the spin axis. For a spin-aligned dipole field the marginally stable distribution is approximately ρr^9/2 = constant, where ρ is the plasma mass density at the radial distance r in the equatorial plane. Plasma filling the Jovian magnetosphere from internal sources would initially violate this stability criterion so that interchange motions would act to establish the marginally stable distribution.     

Classes of orbits in the main problem of satellite theory

We consider the main problem in satellite theory restricted to the polar plane. For suitable values of the energy the system has two unstable periodic orbits. We classify the trajectories in terms of their ultimate behavior with respect these periodic orbits in: oscillating, asymptotic and capture orbits. We study the energy level set and the existence and properties of the mentioned types of motion.     

The eleventh motion constant of the two-body problem

The two-body problem is a twelfth-order time-invariant dynamic system, and therefore has eleven mutually-independent time-independent integrals, here referred to as motion constants. Some of these motion constants are related to the ten mutually-independent algebraic integrals of the n-body problem, whereas some are particular to the two-body problem. The problem can be decomposed into mass-center and relative-motion subsystems, each being sixth-order and each having five mutually-independent motion constants. This paper presents solutions for the eleventh motion constant, which relates the behavior of the two subsystems. The complete set of mutually-independent motion constants describes the shape of the state-space trajectories. The use of the eleventh motion constant is demonstrated in computing a solution to a two-point boundary-value problem.     

An efficient method for studying the stability and dynamics of the rotational motions of celestial bodies

We use the alternative MEGNO (Mean Exponential Growth of Nearby Orbits) technique developed by Cincotta and Simó to study the stability of orbital—rotational motions for plane oscillations and three-dimensional rotations. We present a detailed numerical—analytical study of a rigid body in the case where the proper rotation of the body is synchronized with its orbital motion as 3: 2(Mercurian—type synchronism). For plane rotations, the loss of stability of the periodic solution that corresponds to a 3: 2 resonance is shown to be soft, which should be taken into account to estimate the upper limit for the ellipticity of Mercury. In studying stable and chaotic translational—rotational motions, we point out that the MEGNO criterion can be effectively used. This criterion gives a clear picture of the resonant structures and allows the calculations to be conveniently presented in the form of the corresponding MEGNO stability maps for multidimensional systems. We developed an appropriate software package.     

Analysis of the neighborhood of the 2: 1 resonance in the equal-mass three-body problem

We consider the trajectories in the neighborhood of a 2: 1 resonance (in periods of osculating motions of the outer and inner binaries) in the plane equal-mass three-body problem. We identified the zones of motions that are stable on limited time intervals. All of them correspond to the retrograde motions of the outer and inner subsystems. The prograde motions are unstable: the triple system breaks up into a final binary and an escaping component. In the barycentric nonrotating coordinate system, the trajectories occasionally form symmetric structures composed of several leaves. These structures persist for a long time, and, subsequently, the trajectories of the bodies fill compact regions in coordinate space.     

On the relation between resonance and instability in planetary systems

The factors which affect the linear stability of a periodic planetary orbit in the plane are studied. It is proved that planetary systems with two planets describing nearly circular orbits in the same direction are linearly stable and no perturbation exists which destroys the stability, unless a resonance of the form 1/3, 3/5, 5/7, ... among the orbits of the planets occurs. This latter resonant case is always unstable. Retrograde motion is always linearly stable. Planetary systems with three or more planets in nearly circular orbits in the same direction are proved to be unstable, in the sense that a Hamiltonian perturbation always exists which destroys the stability. The generation of instability in the case of three or more planets is not only due to the existence of resonances, as in the case of two planets, but also to the nonexistence of integrals of motion, apart from the energy and angular momentum integrals. It is also proved that planetary systems with nearly elliptic orbits of the planets are unstable.     

Long-term evaluation of three-dimensional heliocentric solar sail trajectories with arbitrary fixed sail setting

The solar radiation effects upon the orbital behaviour of an arbitrarily shaped spacecraft (or a solar sail in particular) in a general fixed orientation with respect to the local coordinate frame are investigated. Through introduction of a quasi-angle in the osculating plane, the motion of the orbital plane becomes uncoupled from the in-plane perturbations. Exact solutions in the form of conic sections and logarithmic spirals can readily be formulated for certain specific initial conditions. An effective out-of-plane spiral transfer trajectory is obtained by reversing the force component normal to the orbital plane at specified positions in the orbit. By choosing the appropriate control angles for the sail orientation, any point in space can be reached eventually. In the case of general initial conditions, the long-term orbital behaviour is assessed asymptotically by means of the two-variable expansion procedure. An implicit expression for the eccentricity is derived and explicit results are established by an iteration scheme. The other orbital elements can be expressed in terms of the eccentricity and their asymptotic series for near-circular initial orbits are also obtained. While equations for the higher-order contributions as well as the periodic parts of their solutions can be formulated readily, their secular terms are determined only for a circular initial orbit.