Distant quasi-periodic orbits around Mercury

In this paper, distant quasi-periodic orbits around Mercury are studied for future Mercury missions. All of these orbits have relatively large sizes, with their altitudes near or above the Mercury sphere of influence. The research is carried out in the framework of the elliptic restricted three-body problem (ER3BP) to account for the planet’s non-negligible orbital eccentricity. Retrograde and prograde quasi-periodic trajectories in the planar ER3BP are generalized from periodic orbits in the CR3BP by the homotopy algorithm, and the shape evolution of such quasi-periodic trajectories around Mercury is investigated. Numerical simulations are performed to evaluate the stability of these distant orbits in the long term. These two classes of orbits present different characteristics: retrograde orbits can maintain shape stability with a large size, although the trajectories in some regions may oscillate with larger amplitudes; for prograde orbits, the range of existence is much smaller, and their trajectories easily move away from the vicinity of Mercury when the orbits become larger. Distant orbits can be used to explore the space environment in the vicinity of Mercury, and some orbits can be taken as transfer orbits for low-cost Mercury return missions or other programs for their high maneuverability.     

Oblateness effect of Saturn on periodic orbits in the Saturn-Titan restricted three-body problem

We explore the effect of oblateness of Saturn (more massive primary) on the periodic orbits and the regions of quasi-periodic motion around both the primaries in the Saturn-Titan system in the framework of planar circular restricted three-body problem. First order interior and exterior mean motion resonances are located. The effect of oblateness is studied on the location, nature and size of periodic and quasi-periodic orbits, using the numerical technique of Poincare surface of sections. Some of the periodic orbits change to quasi-periodic orbits due to the effect of oblateness and vice-versa. The stability of the orbits around Saturn, Titan and both varies with the inclusion of oblateness. The centers of the periodic orbits around Titan move towards Saturn, whereas those around Saturn move towards Titan. For the orbit around Titan at C=2.9992, x=0.959494, the apocenter becomes pericenter. By incorporating oblateness effect, the orbit around Titan at C=2.99345, x=0.924938 is captured by Saturn, remains in various trajectories around Saturn, and as time progresses it spirals away around both the primaries.     

Search for and study of nearly periodic orbits in the plane problem of three equal-mass bodies

We analyze nearly periodic solutions in the plane problem of three equal-mass bodies by numerically simulating the dynamics of triple systems. We identify families of orbits in which all three points are on one straight line (syzygy) at the initial time. In this case, at fixed total energy of a triple system, the set of initial conditions is a bounded region in four-dimensional parameter space. We scan this region and identify sets of trajectories in which the coordinates and velocities of all bodies are close to their initial values at certain times (which are approximately multiples of the period). We classify the nearly periodic orbits by the structure of trajectory loops over one period. We have found the families of orbits generated by von Schubart’s stable periodic orbit revealed in the rectilinear three-body problem. We have also found families of hierarchical, nearly periodic trajectories with prograde and retrograde motions. In the orbits with prograde motions, the trajectory loops of two close bodies form looplike structures. The trajectories with retrograde motions are characterized by leafed structures. Orbits with central and axial symmetries are identified among the families found.     

A novel algorithm for generating libration point orbits about the collinear points

This paper presents a numerical algorithm that can generate long-term libration points orbits (LPOs) and the transfer orbits from the parking orbits to the LPOs in the circular-restricted three-body problem (CR3BP) and the full solar system model without initial guesses. The families of the quasi-periodic LPOs in the CR3BP can also be constructed with this algorithm. By using the dynamical behavior of LPO, the transfer orbit from the parking orbit to the LPO is generated using a bisection method. At the same time, a short segment of the target LPO connected with the transfer orbit is obtained, then the short segment of LPO is extended by correcting the state towards its adjacent point on the stable manifold of the target LPO with differential evolution algorithm. By implementing the correction strategy repeatedly, the LPOs can be extended to any length as needed. Moreover, combining with the continuation procedure, this algorithm can be used to generate the families of the quasi-periodic LPOs in the CR3BP.     

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

Application of Hamiltonian structure-preserving control to formation flying on a J 2-perturbed mean circular orbit

The bounded quasi-periodic relative trajectories are investigated in this paper for on-orbit surveillance, inspection or repair, which requires rapid changes in formation configuration for full three-dimensional imaging and unpredictable evolutions of relative trajectories for non-allied spacecraft. A linearized differential equation for modeling J ₂ perturbed relative dynamics is derived without any simplified treatment of full short-period effects. The equation serves as a nominal reference model for stationkeeping controller to generate the quasi-periodic trajectories near the equilibrium, i.e., the location of the chief. The developed model exhibits good numerical accuracy and is applicable to an elliptic orbit with small eccentricity inheriting from the osculating conversion of orbital elements. A Hamiltonian structure-preserving controller is derived for the three-dimensional time-periodic system that models the J ₂-perturbed relative dynamics on a mean circular orbit. The equilibrium of the system has time-varying topological types and no fixed-dimensional unstable/stable/center manifolds, which are quite different from the two-dimensional time-independent system with a permanent pair of hyperbolic eigenvalues and fixed-dimensions of unstable/stable/ center manifolds. The unstable and stable manifolds are employed to change the hyperbolic equilibrium to elliptic one with the poles assigned on the imaginary axis. The detailed investigations are conducted on the critical controller gain for Floquet stability and the optimal gain for the fuel cost, respectively. Any initial relative position and velocity leads to a bounded trajectory around the controlled elliptic equilibrium. The numerical simulation indicates that the controller effectively stabilizes motions relative to the perturbed elliptic orbit with small eccentricity and unperturbed elliptic orbit with arbitrary eccentricity. The developed controller stabilizes the quasi-periodic relative trajectories involved in six foundational motions with different frequencies generated by the eigenvectors of the Floquet multipliers, rather than to track a reference relative configuration. Only the relative positions are employed for the feedback without the information from the direct measurement or the filter estimation of relative velocity. So the current controller has potential applications in formation flying for its less computation overload for on-board computer, less constraint on the measurements, and easily-achievable quasi-periodic relative trajectories.     

Families of periodic orbits in the planar three-body problem

A periodic orbit of the restricted circular three-body problem, selected arbitrarily, is used to generate a family of periodic motions in the general three-body problem in a rotating frame of reference, by varying the massm ₃ of the third body. This family is continued numerically up to a maximum value of the mass of the originally small body, which corresponds to a mass ratiom ₁:m ₂:m ₃?5:5:3. From that point on the family continues for decreasing massesm ₃ until this mass becomes again equal to zero. It turns out that this final orbit of the family is a periodic orbit of the elliptic restricted three body problem. These results indicate clearly that families of periodic motions of the three-body problem exist for fixed values of the three masses, since this continuation can be applied to all members of a family of periodic orbits of the restricted three-body problem. It is also indicated that the periodic orbits of the circular restricted problem can be linked with the periodic orbits of the elliptic three-body problem through periodic orbits of the general three-body problem.     

Lunar capture trajectories and homoclinic connections through isomorphic mapping

Analysis and design of low-energy transfers to the Moon has been a subject of great interest for many decades. This paper is concerned with a topological study of such transfers, with emphasis to trajectories that allow performing lunar capture and those that exhibit homoclinic connections, in the context of the circular restricted three-body problem. A fundamental theorem stated by Conley locates capture trajectories in the phase space and can be condensed in a sentence: “if a crossing asymptotic orbit exists then near any such there is a capture orbit”. In this work this fundamental theoretical assertion is used together with an original cylindrical isomorphic mapping of the phase space associated with the third body dynamics. For a given energy level, the stable and unstable invariant manifolds of the periodic Lyapunov orbit around the collinear interior Lagrange point are computed and represented in cylindrical coordinates as tubes that emanate from the transformed periodic orbit. These tubes exhibit complex geometrical features. Their intersections correspond to homoclinic orbits and determine the topological separation of long-term lunar capture orbits from short-duration capture trajectories. The isomorphic mapping is proven to allow a deep insight on the chaotic motion that characterizes the dynamics of the circular restricted three-body, and suggests an interesting interpretation, and together corroboration, of Conley’s assertion on the topological location of lunar capture orbits. Moreover, an alternative three-dimensional representation of the phase space is profitably employed to identify convenient lunar periodic orbits that can be entered with modest propellant consumption, starting from the Lyapunov orbit.     

On the possible resonant orbits of exoplanets

Within the framework of the restricted three-body problem, the possible orbits of a small-mass exoplanet in the system with a massive exoplanet on an elliptic orbit are investigated. Possible quasi-circular orbits are sought. The dependence of the Kozdai-Lidov effect (the Kozdai resonance) on the eccentricity of the orbit of a massive planet is discussed. The effect of the commensurabilities of the mean motions on the value of the eccentricity perturbations is considered.     

The rectilinear three-body problem

The rectilinear equal-mass and unequal-mass three-body problems are considered. The first part of the paper is a review that covers the following items: regularization of the equations of motion, integrable cases, triple collisions and their vicinities, escapes, periodic orbits and their stability, chaos and regularity of motions. The second part contains the results of our numerical simulations in this problem. A classification of orbits in correspondence with the following evolution scenarios is suggested: ejections, escapes, conditional escapes (long ejections), periodic orbits, quasi-stable long-lived systems in the vicinity of stable periodic orbits, and triple collisions. Homothetic solutions ending by triple collisions and their dependence on initial parameters are found. We study how the ejection length changes in response to the variation of the triple approach parameters. Regions of initial conditions are outlined in which escapes occur after a definite number of triple approaches or a definite time. In the vicinity of a stable Schubart periodic orbit, we reveal a region of initial parameters that corresponds to trajectories with finite motions. The regular and chaotic structure of the manifold of orbits is mostly defined by this periodic orbit. We have studied the phase space structure via Poincaré sections. Using these sections and symbolic dynamics, we study the fine structure of the region of initial conditions, in particular the chaotic scattering region.     

Fast Earth–Moon transfers with ballistic capture

This contribution deals with fast Earth–Moon transfers with ballistic capture in the patched three-body model. We compute ensembles of preliminary solutions using a model that takes into account the relative inclination of the orbital planes of the primaries. The ballistic capture orbits around the Moon are obtained relying on the hyperbolic invariant structures associated to the collinear Lagrangian points of the Earth–Moon system, and the Sun–Earth system portion of the transfers are quasi-periodic orbits obtained by a genetic algorithm. The trajectories are designed to be good initial guesses to search optimal cost-efficient short-time Earth–Moon transfers with ballistic capture in more realistic models.     

Strategies for plane change of Earth orbits using lunar gravity and derived trajectories of family G

The dynamics of the circular restricted three-body Earth-Moon-particle problem predicts the existence of the retrograde periodic orbits around the Lagrangian equilibrium point L1. Such orbits belong to the so-called family G (Broucke, Periodic orbits in the restricted three-body problem with Earth-Moon masses, JPL Technical Report 32–1168, 1968) and starting from them it is possible to define a set of trajectories that form round trip links between the Earth and the Moon. These links occur even with more complex dynamical systems as the complete Sun-Earth-Moon-particle problem. One of the most remarkable properties of these trajectories, observed for the four-body problem, is a meaningful inclination gain when they penetrate into the lunar sphere of influence and accomplish a swing-by with the Moon. This way, when one of these trajectories returns to the proximities of the Earth, it will be in a different orbital plane from its initial Earth orbit. In this work, we present studies that show the possibility of using this property mainly to accomplish transfer maneuvers between two Earth orbits with different altitudes and inclinations, with low cost, taking into account the dynamics of the four-body problem and of the swing-by as well. The results show that it is possible to design a set of nominal transfer trajectories that require ΔV _Total less than conventional methods like Hohmann, bi-elliptic and bi-parabolic transfer with plane change.     

On possibility of clustering of microparticles of orbital debris in circular and elliptical orbits

Estimation is made of the possibility of clustering of debris particles in circular and elliptical orbits around the Earth due to the change in drag, caused by quasi-periodic variations of the atmospheric density in the orbit. The estimations show that the collective behavior of particles has time to be manifested in highly elliptical orbits, where the relative change in the atmospheric density along the orbital path is greater and characteristic lifetimes of the particles are longer. However, in this case the limit distributions of the particles are not realized, because the clusters form and break down several times during the lifetime of the particles in the orbit.     

A numerical investigation of the one-dimensional Newtonian three-body problem

The one-dimensional Newtonian three-body problem is known to have stable (quasi-)periodic orbits when the masses are equal. The existence and size of the stable region is discussed here in the case where the three masses are arbitrary. We consider only the stability of the periodic (generalized) Schubart's (1956) orbit. If this orbit is linearly stable it is almost always surrounded by a region of stable quasi-periodic orbits and the size and shape of this stable region depends on the masses. The three-dimensional linear stability of the periodic orbits is also determined. Final results show that the region of stability has a complicated shape and some of the stable regions in the mass-plane are quite narrow. The non-linear three-dimensional stability is studied independently by extensive numerical integrations and the results are found to be in agreement with the linear stability analysis. The boundaries of stable region in the mass-plane are given in terms of polynomial approximations. The results are compared with a similar work by Héenon (1977).We thank the referee for pointing out this reference to us.     

Numerical models for the study of motion of lunar satellites

The present study deals with numerical modeling of the elliptic restricted three-body problem as well as of the perturbed elliptic restricted three-body (Earth-Moon-Satellite) problem by a fourth body (Sun). Two numerical algorithms are established and investigated. The first is based on the method of the series solution of the differential equations and the second is based on a 5th-order Runge-Kutta method. The applications concern the solution of the equations and integrals of motion of the circular and elliptical restricted three-body problem as well as the search for periodic orbits of the natural satellites of the Moon in the Earth-Moon system in both cases in which the Moon describes circular or elliptical orbit around the Earth before the perturbations induced by the Sun. After the introduction of the perturbations in the Earth-Moon-Satellite system the motions of the Moon and the Satellite are studied with the same initial conditions which give periodic orbits for the unperturbed elliptic problem.     

Existence of quasi-periodic solutions to the three-body problem

A rigorous proof is given for the existence of quasi-periodic solutions with only two degrees of freedom to a planar three-body problem. The solution corresponds physically to the small bodies moving on different, nearly elliptical orbits about a large mass located at a focus. The perihelia of the two orbits are locked in such a way that the difference of the two perihelia has mean value zero.     

Lambert problem solution in the hill model of motion

The goal of this paper is obtaining a solution of the Lambert problem in the restricted three-body problem described by the Hill equations. This solution is based on the use of pre determinate reference orbits of different types giving the first guess and defining the sought-for transfer type. A mathematical procedure giving the Lambert problem solution is described. This procedure provides step-by-step transformation of the reference orbit to the sought-for transfer orbit. Numerical examples of the procedure application to the transfers in the Sun–Earth system are considered. These examples include transfer between two specified positions in a given time, a periodic orbit design, a halo orbit design, halo-to-halo transfers, LEO-to-halo transfer, analysis of a family of the halo-to-halo transfer orbits. The proposed method of the Lambert problem solution can be used for the two-point boundary value problem solution in any model of motion if a set of typical reference orbits can be found.     

The Restricted Hill Four-Body Problem with Applications to the Earth–Moon–Sun System

Starting from the four-body problem a generalization of both the restricted three-body problem and the Hill three-body problem is derived. The model is time periodic and contains two parameters: the mass ratio ν of the restricted three-body problem and the period parameter m of the Hill Variation orbit. In the proper coordinate frames the restricted three-body problem is recovered as m → 0 and the classical Hill three-body problem is recovered as ν → 0. This model also predicts motions described by earlier researchers using specific models of the Earth–Moon–Sun system. An application of the current model to the motion of a spacecraft in the Sun perturbed Earth–Moon system is made using Hill's Variation orbit for the motion of the Earth–Moon system. The model is general enough to apply to the motion of an infinitesimal mass under the influence of any two primaries which orbit a larger mass. Using the model, numerical investigations of the structure of motions around the geometric position of the triangular Lagrange points are performed. Values of the parameter ν range in the neighborhood of the Earth–Moon value as the parameter m increases from 0 to 0.195 at which point the Hill Variation orbit becomes unstable. Two families of planar periodic orbits are studied in detail as the parameters m and ν vary. These families contain stable and unstable members in the plane and all have the out-of-plane stability. The stable and unstable manifolds of the unstable periodic orbits are computed and found to be trapped in a geometric area of phase space over long periods of time for ranges of the parameter values including the Earth–Moon–Sun system. This model is derived from the general four-body problem by rigorous application of the Hill and restricted approximations. The validity of the Hill approximation is discussed in light of the actual geometry of the Earth–Moon–Sun system. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.