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.     

A Hopf variables view on the libration points dynamics

This study analyzes a recently discovered class of exterior transfers to the Moon. These transfers terminate in retrograde ballistic capture orbits, i.e., orbits with negative Keplerian energy and angular momentum with respect to the Moon. Yet, their Jacobi constant is relatively low, for which no forbidden regions exist, and the trajectories do not appear to mimic the dynamics of the invariant manifolds of the Lagrange points. This paper shows that these orbits shadow instead lunar collision orbits. We investigate the dynamics of singular, lunar collision orbits in the Earth–Moon planar circular restricted three-body problem, and reveal their rich phase space structure in the medium-energy regime, where invariant manifolds of the Lagrange point orbits break up. We show that lunar retrograde ballistic capture trajectories lie inside the tube structure of collision orbits. We also develop a method to compute medium-energy transfers by patching together orbits inside the collision tube and those whose apogees are located in the appropriate quadrant in the Sun–Earth system. The method yields the novel family of transfers as well as those ending in direct capture orbits, under particular energetic and geometrical conditions.     

Earth–Moon Weak Stability Boundaries in the restricted three and four body problem

This work deals with the structure of the lunar Weak Stability Boundaries (WSB) in the framework of the restricted three and four body problem. Geometry and properties of the escape trajectories have been studied by changing the spacecraft orbital parameters around the Moon. Results obtained using the algorithm definition of the WSB have been compared with an analytical approximation based on the value of the Jacobi constant. Planar and three-dimensional cases have been studied in both three and four body models and the effects on the WSB structure, due to the presence of the gravitational force of the Sun and the Moon orbital eccentricity, have been investigated. The study of the dynamical evolution of the spacecraft after lunar capture allowed us to find regions of the WSB corresponding to stable and safe orbits, that is orbits that will not impact onto lunar surface after capture. By using a bicircular four body model, then, it has been possible to study low-energy transfer trajectories and results are given in terms of eccentricity, pericenter altitude and inclination of the capture orbit. Equatorial and polar capture orbits have been compared and differences in terms of energy between these two kinds of orbits are shown. Finally, the knowledge of the WSB geometry permitted us to modify the design of the low-energy capture trajectories in order to reach stable capture, which allows orbit circularization using low-thrust propulsion systems.     

On quasi-periodic motions around the collinear libration points in the real Earth–Moon system

Due to various perturbations, the collinear libration points of the real Earth–Moon system are not equilibrium points anymore. Under the assumption that the Moon’s motion is quasi-periodic, special quasi-periodic orbits called dynamical substitutes exist. These dynamical substitutes replace the geometrical collinear libration points as time-varying equilibrium points. In the paper, the dynamical substitutes of the three collinear libration points in the real Earth–Moon system are computed. For the points L ₁ and L ₂, linearized motions around the dynamical substitutes are described, and the variational equations of the dynamical substitutes are reduced to a form with a near constant coefficient matrix. Then higher order analytical formulae of the central manifolds are constructed. Using these analytical solutions as initial seeds, Lissajous orbits and halo orbits are computed with numerical algorithms.     

Chaotic Asteroidal Dynamics and Maximum Lyapunov Exponents

This paper describes the results of studies of dynamical chaos in the problem of the orbital dynamics of asteroids near the 3 : 1 mean-motion resonance with Jupiter. Maximum Lyapunov characteristic exponents (MLCEs) are used as an indicator and a measure of the chaoticity of motion. MLCE values are determined for trajectories calculated by the numerical integration of equations of motion in the planar elliptical restricted three-body problem. The dependence of the MLCE on the problem parameters and on the initial data is analyzed. The inference is made that the domain of chaos in the phase space of the problem considered consists of two components of different nature. The values of the MLCEs observed for one of the components (namely, for the component corresponding to low-eccentricity asteroidal orbits) are compared to the theoretical estimates obtained within the framework of model of the resonance as a perturbed nonlinear pendulum.     

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

Chaotic motion of comets in near-parabolic orbit: Mapping approaches

There exist many comets with near-parabolic orbits in the Solar System. Among various theories proposed to explain their origin, the Oort cloud hypothesis seems to be the most reasonable (Oort, 1950). The theory assumes that there is a cometary cloud at a distance 10³ – 10⁵ AU from the Sun and that perturbing forces from planets or stars make orbits of some of these comets become of near-parabolic type. Concerning the evolution of these orbits under planetary perturbations, we can raise the question: Will they stay in the Solar System forever or will they escape from it? This is an attractive dynamical problem. If we go ahead by directly solving the dynamical differential equations, we may encounter the difficulty of long-time computation. For the orbits of these comets are near-parabolic and their periods are too long to study on their long-term evolution. With mapping approaches the difficulty will be overcome. In another aspect, the study of this model has special meaning for chaotic dynamics. We know that in the neighbourhood of any separatrix i.e. the trajectory with zero frequency of the unperturbed motion of an Hamiltonian system, some chaotic motions have to be expected. Actually, the simplest example of separatrix is the parabolic trajectory of the two body problem which separates the bounded and unbounded motion. From this point of view, the dynamical study on near-parabolic motion is very important. Petrosky's elegant but more abstract deduction gives a Kepler mapping which describes the dynamics of the cometary motion (Petrosky, 1988). In this paper we derive a similar mapping directly and discuss its dynamical characters.     

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.     

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.     

Bounded relative motion under zonal harmonics perturbations

The problem of finding natural bounded relative trajectories between the different units of a distributed space system is of great interest to the astrodynamics community. This is because most popular initialization methods still fail to establish long-term bounded relative motion when gravitational perturbations are involved. Recent numerical searches based on dynamical systems theory and ergodic maps have demonstrated that bounded relative trajectories not only exist but may extend up to hundreds of kilometers, i.e., well beyond the reach of currently available techniques. To remedy this, we introduce a novel approach that relies on neither linearized equations nor mean-to-osculating orbit element mappings. The proposed algorithm applies to rotationally symmetric bodies and is based on a numerical method for computing quasi-periodic invariant tori via stroboscopic maps, including extra constraints to fix the average of the nodal period and RAAN drift between two consecutive equatorial plane crossings of the quasi-periodic solutions. In this way, bounded relative trajectories of arbitrary size can be found with great accuracy as long as these are allowed by the natural dynamics and the physical constraints of the system (e.g., the surface of the gravitational attractor). This holds under any number of zonal harmonics perturbations and for arbitrary time intervals as demonstrated by numerical simulations about an Earth-like planet and the highly oblate primary of the binary asteroid (66391) 1999 KW4.     

The Trojan problem

There has been a renewed interest in the Trojan problem in recent years. Significant progress has been made in discovering and understanding dynamical features of motion of Jupiter's Trojan asteroids. The dynamics of hypothetical Trojan-type asteroids of other major planets has also been the subject of several recent investigations. This paper offers an overview on the current status of researches on real and hypothetical Trojan asteroids of the major planets. Results of analytical and numerical works are surveyed. Questions of dynamical properties, long-term evolution of orbits, stability regions around the triangular Lagrangian points are discussed among other problems of the Trojans.     

Stochasticity of Planetary Orbits in Double Star Systems

Cosmogonical theories as well as recent observations allow us to expect the existence of planets around many stars other than the Sun. On an other hand, double and multiple star systems are established to be more numerous than single stars (such as the Sun), at least in the solar neighborhood. We are then faced to the following dynamical problem: assuming that planets can form in a binary early environment (I do not deal here with), does long-term stability for planetary orbits exist in double star systems.Although preliminary studies were rather pessimistic about the possibility of existence of stable planetary orbits in double or multiple star systems, modern computation have shown that many such stable orbits do exist (but possible chaotic behavior), either around the binary as a whole (P-type) or around one component of the binary (S-type), this latter being explored here.The dynamical model is the elliptic plane restricted three-body problem; the phase space of initial conditions is systematically explored, and limits for stability have been established. Stable S-type planetary orbits are found up to distance of their "sun" of the order of half the periastron distance of the binary; moreover, among these stable orbits, nearly-circular ones exist up to distance of their "sun" of the order of one quarter the periastron distance of the binary; finally, among the nearly-circular stable orbits, several stay inside the "habitable zone", at least for two nearby binaries which components are nearly of solar type.Nevertheless, we know that chaos may destroy this stability after a long time (sometimes several millions years). It is therefore important to compute indicators of chaos for these stable planetary orbits to investigate their actual very long-term stability. Here we give an example of such a computation for more than a billion years.     

Lunar capture in the planar restricted three-body problem

The capture dynamics is an important field in Astronomy and Astronautics. In this paper, the near-optimal lunar capture in the Earth–Moon transfer is investigated under the frame of the planar circular restricted three-body problem. We try to work out how to achieve the permanent lunar capture with the minimum maneuver consumption. This problem is decomposed into two parts: the pre-maneuver part and the post-maneuver part. In the pre-maneuver part, considering the criteria of the gravitational capture, we obtain the minimum pre-maneuver velocity via the numerical backward integration. In the post-maneuver part, using the Poincaré section and the KAM theory, we find the maximum post-maneuver velocity to achieve the permanent capture. Synthesized the results of the two parts, a new method is presented to find the near-optimal maneuver position and the minimum maneuver consumption. The method presented is simple and visible, and can provide abundant capture orbits for the design of low energy Earth–Moon transfers.     

The Two Fixed Centers: An Exceptional Integrable System

It is usually believed that we know everything to be known for any separable Hamiltonian system, i.e. an integrable system in which we can separate the variables in some coordinate system (e.g. see Lichtenberg and Lieberman 1992, Regular and Chaotic Dynamics, Springer). However this is not always true, since through the separation the solutions may be found only up to quadratures, a form that might not be particularly useful. A good example is the two-fixed-centers problem. Although its integrability was discovered by Euler in the 18th century, the problem was far from being considered as completely understood. This apparent contradiction stems from the fact that the solutions of the equations of motion in the confocal ellipsoidal coordinates, in which the variables separate, are written in terms of elliptic integrals, so that their properties are not obvious at first sight. In this paper we classify the trajectories according to an exhaustive scheme, comprising both periodic and quasi-periodic ones. We identify the collision orbits (both direct and asymptotic) and find that collision orbits are of complete measure in a 3-D submanifold of the phase space while asymptotically collision orbits are of complete measure in the 4-D phase space. We use a transformation, which regularizes the close approaches and, therefore, enables the numerical integration of collision trajectories (both direct and asymptotic). Finally we give the ratio of oscillation period along the two axes (the ‘rotation number’) as a function of the two integrals of motion. This revised version was published online in July 2006 with corrections to the Cover Date.     

New contributions to the problem of capture

We present a treatment of libration-point capture in the restricted three-body problem. Examples of capture are given, and a long-term numerical integration is presented, to illustrate major features of orbits arising from capture. A theory of lifetimes is given, providing order-of-magnitude (though rather conservative) estimates of the time a body remains captured. A general capture criterion, giving bounds on admissible values of the postcapture semimajor axis, for given values of eccentricity and inclination. This criterion is used to demonstrate that, in general, direct postcapture orbits lie outside retrograde ones. We also emphasize the importance of mass-change, of one or both primaries, in producing capture. This phenomenon is shown to give rise to a new type of capture, “pull-down capture,” which produces retrograde orbits. The effects of nebular drag also are noted.These results suggest the improbability of a capture origin for Jupiter's outer satellites within the last 4+ billion years, or since the solar system reached its present dynamical configuration. Computations indicate, however, that either mass-change or nebular drag could have been effective in producing capture. The outer satellite groups are shown to resemble Hirayama families physically, thus supporting a hypothesis of capture followed by collisional fragmentation.     

Modeling the 3-D secular planetary three-body problem: Discussion on the outer υ Andromedae planetary system

The three-dimensional secular behavior of a system composed of a central star and two massive planets is modeled semi-analytically in the frame of the general three-body problem. The main dynamical features of the system are presented in geometrical pictures allowing us to investigate a large domain of the phase space of this problem without time-expensive numerical integrations of the equations of motion and without any restriction on the magnitude of the planetary eccentricities, inclinations and mutual distance. Several regimes of motion of the system are observed. With respect to the secular angle Δ?, possible motions are circulations, oscillations (around 0° and 180°), and high-eccentricity/inclination librations in secular resonances. With respect to the arguments of pericenter, ω₁ and ω₂, possible motions are direct circulation and high-inclination libration around ±90° in the Lidov-Kozai resonance. The regions of transition between domains of different regimes of motion are characterized by chaotic behavior. We apply the analysis to the case of the two outer planets of the υ Andromedae system, observed edge-on. The topology of the 3-D phase space of this system is investigated in detail by means of surfaces of section, periodic orbits and dynamical spectra, mapping techniques and numerical simulations. We obtain the general structure of the phase space, and the boundaries of the spatial secular stability. We find that this system is secularly stable in a large domain of eccentricities and inclinations.     

The dynamics of the elliptic Hill problem: periodic orbits and stability regions

The motion of a satellite around a planet can be studied by the Hill model, which is a modification of the restricted three body problem pertaining to motion of a satellite around a planet. Although the dynamics of the circular Hill model has been extensively studied in the literature, only few results about the dynamics of the elliptic model were known up to now, namely the equations of motion and few unstable families of periodic orbits. In the present study we extend these results by computing a large set of families of periodic orbits and their linear stability and classify them according to their resonance condition. Although most of them are unstable, we were able to find a considerable number of stable ones. By computing appropriate maps of dynamical stability, we study the effect of the planetary eccentricity on the stability of satellite orbits. We see that, even for large values of the planetary eccentricity, regular orbits can be found in the vicinity of stable periodic orbits. The majority of irregular orbits are escape orbits.     

A keplerian method to estimate perturbations in the restricted three-body problem

We develop a new and fast method to estimate perturbations by a planet on cometary orbits. This method allows us to identify accurately the cases of large perturbations in a set of fictitious orbits. Hence, it can be used in constructing perturbation samples for Monte Carlo simulations in order to maximize the amount of information. Furthermore, the estimated perturbations are found to yield a good approximation to the real perturbation sample. This is shown by a comparison of the perturbations obtained by the new estimator with the results of numerical integration of regularized equations of motion for the same orbits in the same dynamical model: the three-dimensional elliptic restricted three-body problem (Sun-Jupiter-comet).