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.     

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

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.     

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 two-level perturbation method for connecting unstable periodic orbits with low fuel cost and short time of flight: application to a lunar observation mission

The proposed method connects two unstable periodic orbits by employing trajectories of their associated invariant manifolds that are perturbed in two levels. A first level of velocity perturbations is applied on the trajectories of the discretized manifolds at the points where they approach the nominal unstable periodic orbit in order to accelerate them. A second level of structured velocity perturbations is applied to trajectories that have already been subjected to first level perturbations in order to approximately meet the necessary conditions for a low \(\varDelta \text {V}\) transfer. Due to this two-level perturbation approach, the number of the trajectories obtained is significantly larger compared with approaches that employ traditional invariant manifolds. For this reason, the problem of connecting two unstable periodic orbits through perturbed trajectories of their manifolds is transformed into an equivalent discrete optimization problem that is solved with a very low computational complexity algorithm that is proposed in this paper. Finally, the method is applied to a lunar observation mission of practical interest and is found to perform considerably better in terms of \(\varDelta \text {V}\) cost and time of flight when compared with previous techniques applied to the same project.     

On the construction of low-energy cislunar and translunar transfers based on the libration points

There exist cislunar and translunar libration points near the Moon, which are referred to as the LL ₁ and LL ₂ points, respectively. They can generate the different types of low-energy trajectories transferring from Earth to Moon. The time-dependent analytic model including the gravitational forces from the Sun, Earth, and Moon is employed to investigate the energy-minimal and practical transfer trajectories. However, different from the circular restricted three-body problem, the equivalent gravitational equilibria are defined according to the geometry of the instantaneous Hill boundary due to the gravitational perturbation from the Sun. The relationship between the altitudes of periapsis and eccentricities is achieved from the Poincaré mapping for all the captured lunar trajectories, which presents the statistical feature of the fuel cost and captured orbital elements rather than generating a specified Moon-captured segment. The minimum energy required by the captured trajectory on a lunar circular orbit is deduced in the spatial bi-circular model. The idea is presented that the asymptotical behaviors of invariant manifolds approaching to/traveling from the libration points or halo orbits are destroyed by the solar perturbation. In fact, the energy-minimal cislunar transfer trajectory is acquired by transiting the LL ₁ point, while the energy-minimal translunar transfer trajectory is obtained by transiting the LL ₂ point. Finally, the transfer opportunities for the practical trajectories that have escaped from the Earth and have been captured by the Moon are yielded by the transiting halo orbits near the LL ₁ and LL ₂ points, which can be used to generate the whole of the trajectories.     

On some applications of invariant manifolds

Taking transfer orbits of a collinear libration point probe, a lunar probe and an interplanetary probe as examples, some applications of stable and unstable invariant manifolds of the restricted three-body problem are discussed. Research shows that transfer energy is not necessarily conserved when invariant manifolds are used. For the cases in which the transfer energy is conserved, the cost is a much longer transfer time.     

A motivating exploration on lunar craters and low-energy dynamics in the Earth–Moon system

It is known that most of the craters on the surface of the Moon were created by the collision of minor bodies of the Solar System. Main Belt Asteroids, which can approach the terrestrial planets as a consequence of different types of resonance, are actually the main responsible for this phenomenon. Our aim is to investigate the impact distributions on the lunar surface that low-energy dynamics can provide. As a first approximation, we exploit the hyberbolic invariant manifolds associated with the central invariant manifold around the equilibrium point L ₂ of the Earth–Moon system within the framework of the Circular Restricted Three-Body Problem. Taking transit trajectories at several energy levels, we look for orbits intersecting the surface of the Moon and we attempt to define a relationship between longitude and latitude of arrival and lunar craters density. Then, we add the gravitational effect of the Sun by considering the Bicircular Restricted Four-Body Problem. In the former case, as main outcome, we observe a more relevant bombardment at the apex of the lunar surface, and a percentage of impact which is almost constant and whose value depends on the assumed Earth–Moon distance d_EM. In the latter, it seems that the Earth–Moon and Earth–Moon–Sun relative distances and the initial phase of the Sun θ ₀ play a crucial role on the impact distribution. The leading side focusing becomes more and more evident as d_EM decreases and there seems to exist values of θ ₀ more favorable to produce impacts with the Moon. Moreover, the presence of the Sun makes some trajectories to collide with the Earth. The corresponding quantity floats between 1 and 5 percent. As further exploration, we assume an uniform density of impact on the lunar surface, looking for the regions in the Earth–Moon neighbourhood these colliding trajectories have to come from. It turns out that low-energy ejecta originated from high-energy impacts are also responsible of the phenomenon we are considering.     

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.     

满足一定约束条件的登月飞行轨道的设计

讨论满足约束条件的登月飞行轨道的设计问题，将约束条件分类为只与太阳，月球，地球，飞行器和观测站之间的相对位置有关的运动学约束条件以及小及到飞行器轨道云动的动力学约束条件，在考虑登月飞行轨道的特征后，给出一种设计满足约束条件的标准飞行轨道的方法，并将方法应用于不同约束条件下的我国登月飞行以及月球卫星的轨道预测计。     

A greedy global search algorithm for connecting unstable periodic orbits with low energy cost.

A method for space mission trajectory design is presented in the form of a greedy global search algorithm. It uses invariant manifolds of unstable periodic orbits and its main advantage is that it performs a global search for the suitable legs of the invariant manifolds to be connected for a preliminary transfer design, as well as the appropriate points of the legs for maneuver application. The designed indirect algorithm bases the greedy choice on the optimality conditions that are assumed for the theoretical minimum transfer cost of a spacecraft when using invariant manifolds. The method is applied to a test case space mission design project in the Earth–Moon system and is found to compare favorably with previous techniques applied to the same project.     

Optimal low-thrust transfers between libration point orbits

Over the past three decades, ballistic and impulsive trajectories between libration point orbits (LPOs) in the Sun–Earth–Moon system have been investigated to a large extent. It is known that coupling invariant manifolds of LPOs of two different circular restricted three-body problems (i.e., the Sun–Earth and the Earth–Moon systems) can lead to significant mass savings in specific transfers, such as from a low Earth orbit to the Moon’s vicinity. Previous investigations on this issue mainly considered the use of impulsive maneuvers along the trajectory. Here we investigate the dynamical effects of replacing impulsive ΔV’s with low-thrust trajectory arcs to connect LPOs using invariant manifold dynamics. Our investigation shows that the use of low-thrust propulsion in a particular phase of the transfer and the adoption of a more realistic Sun–Earth–Moon four-body model can provide better and more propellant-efficient solution. For this purpose, methods have been developed to compute the invariant tori and their manifolds in this dynamical model.     

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.     

Low-energy Lunar Trajectories with Lunar Flybystwo

The low-energy lunar trajectories with lunar flybys are investigated based on the Sun-Earth-Moon bicircular problem (BCP). The characteristics of the distribution of trajectories in the phase space are summarized. Using the invariant manifolds in the BCP system, the low-energy lunar trajectories with lunar flybys are sought. Then, take time as an augmented dimension in the phase space of a nonautonomous system, we present the state space map and reveal the distribution of these lunar trajectories in the phase space. Consequently, we find that the low-energy lunar trajectories exist as families, and that the every moment in the Sun-Earth-Moon synodic period can be the departure date. Finally, we analyse the velocity increment, transfer duration, and system energy for the different trajectory families, and obtain the velocity-impulse optimal family and the transfer-duration optimal family, respectively.     

Study of the transfer between libration point orbits and lunar orbits in Earth–Moon system

This paper is devoted to the study of the transfer problem from a libration point orbit of the Earth–Moon system to an orbit around the Moon. The transfer procedure analysed has two legs: the first one is an orbit of the unstable manifold of the libration orbit and the second one is a transfer orbit between a certain point on the manifold and the final lunar orbit. There are only two manoeuvres involved in the method and they are applied at the beginning and at the end of the second leg. Although the numerical results given in this paper correspond to transfers between halo orbits around the \(L_1\) point (of several amplitudes) and lunar polar orbits with altitudes varying between 100 and 500 km, the procedure we develop can be applied to any kind of lunar orbits, libration orbits around the \(L_1\) or \(L_2\) points of the Earth–Moon system, or to other similar cases with different values of the mass ratio.     

Three-body problem, its Lagrangian points and how to exploit them using an alternative transfer to L4 and L5

An alternative transfer strategy to send spacecraft to stable orbits around the Lagrangian equilibrium points L4 and L5 based in trajectories derived from the periodic orbits around L1 is presented in this work. The trajectories derived, called Trajectories G, are described and studied in terms of the initial generation requirements and their energy variations relative to the Earth through the passage by the lunar sphere of influence. Missions for insertion of spacecraft in elliptic orbits around L4 and L5 are analysed considering the restricted three-body problem Earth–Moon-particle and the results are discussed starting from the thrust, time of flight and energy variation relative to the Earth.     

Design of lunar landing trajectory under constraints

The design of a lunar landing trajectory which satisfies certain constraints is considered and discussed. The constraints are of two kinds, kinetic constraints, which deal with the relative positions among the Sun, the Moon, the Earth, the spacecraft and tracking stations, and dynamic constraints, which deal with the orbital motion of the spacecraft. After a discussion of the characteristics of lunar flight trajectory, a method of designing standard flight trajectory is suggested that satisfies the constraints. This method is applied to the Chinese lunar landing flight and to the pre-design of the orbit of a lunar satellite.     

Transport orbits in an equilateral restricted four-body problem

In this paper we consider a restricted equilateral four-body problem where a particle of negligible mass is moving under the Newtonian gravitational attraction of three masses (called primaries) which move on circular orbits around their center of masses such that their configuration is always an equilateral triangle (Lagrangian configuration). We consider the case of two bodies of equal masses, which in adimensional units is the parameter of the problem. We study numerically the existence of families of unstable periodic orbits, whose invariant stable and unstable manifolds are responsible for the existence of homoclinic and heteroclinic connections, as well as of transit orbits traveling from and to different regions. We explore, for three different values of the mass parameter, what kind of transits and energy levels exist for which there are orbits with prescribed itineraries visiting the neighborhood of different primaries.     

The use of invariant manifolds for transfers between unstable periodic orbits of different energies

Techniques from dynamical systems theory have been applied to the construction of transfers between unstable periodic orbits that have different energies. Invariant manifolds, trajectories that asymptotically depart or approach unstable periodic orbits, are used to connect the initial and final orbits. The transfer asymptotically departs the initial orbit on a trajectory contained within the initial orbit’s unstable manifold and later asymptotically approaches the final orbit on a trajectory contained within the stable manifold of the final orbit. The manifold trajectories are connected by the execution of impulsive maneuvers. Two-body parameters dictate the selection of the individual manifold trajectories used to construct efficient transfers. A bounding sphere centered on the secondary, with a radius less than the sphere of influence of the secondary, is used to study the manifold trajectories. A two-body parameter, κ, is computed within the bounding sphere, where the gravitational effects of the secondary dominate. The parameter κ is defined as the sum of two quantities: the difference in the normalized angular momentum vectors and eccentricity vectors between a point on the unstable manifold and a point on the stable manifold. It is numerically demonstrated that as the κ parameter decreases, the total cost to complete the transfer decreases. Preliminary results indicate that this method of constructing transfers produces a significant cost savings over methods that do not employ the use of invariant manifolds.