Extended two-body problem for rotating rigid bodies

Given the interest in future space missions devoted to the exploration of key moons in the solar system and that may involve libration point orbits, an efficient design strategy for transfers between moons is introduced that leverages the dynamics in these multi-body systems. The moon-to-moon analytical transfer (MMAT) method is introduced, comprised of a general methodology for transfer design between the vicinities of the moons in any given system within the context of the circular restricted three-body problem, useful regardless of the orbital planes in which the moons reside. A simplified model enables analytical constraints to efficiently determine the feasibility of a transfer between two different moons moving in the vicinity of a common planet. In particular, connections between the periodic orbits of such two different moons are achieved. The strategy is applicable for any type of direct transfers that satisfy the analytical constraints. Case studies are presented for the Jovian and Uranian systems. The transition of the transfers into higher-fidelity ephemeris models confirms the validity of the MMAT method as a fast tool to provide possible transfer options between two consecutive moons.

     

Escape from planetary neighbourhoods

In this paper we use recently developed phase-space transport theory coupled with a so-called classical spectral theorem to develop a dynamically exact and computationally efficient procedure for studying escape from a planetary neighbourhood. The 'planetary neighbourhood' is a bounded region of phase space where entrance and escape are only possible by entering or exiting narrow 'bottlenecks' created by the influence of a saddle point. The method therefore immediately applies to, for example, the circular restricted three-body problem and Hill's lunar problem (which we use to illustrate the results), but it also applies to more complex, and higher-dimensional, systems possessing the relevant phase-space structure. It is shown how one can efficiently compute the mean passage time through the planetary neighbourhood, the phase-space flux in, and out, of the planetary neighbourhood, the phase-space volume of initial conditions corresponding to trajectories that escape from the planetary neighbourhood, and the fraction of initial conditions in the planetary neighbourhood corresponding to bound trajectories. These quantities are computed for Hill's problem. We study the dependence of the proportions of these quantities on energy and dimensionality (two-dimensional planar and three-dimensional spatial Hill's problem). The methods and quantities presented are of central interest for many celestial and stellar dynamical applications such as, for example, the capture and escape of moons near giant planets, the formation of binaries in the Kuiper belt and the escape of stars from star clusters orbiting about a galaxy.     

Numerical study of the three-dimensional transit orbits in the circular restricted three-body problem

Transit orbits are defined as the trajectories that can pass through the neck region of the zero velocity surface in the circular restricted three-body problem (CR3BP). The low-energy transfers in the CR3BP or between two CR3BPs are always through the instrumentality of the transit orbits. In this paper, the distribution of the transit orbits in the six-dimensional phase space is explored by using numerical methods. The necessary and sufficient condition of transition is introduced, which defines the distribution of the transit orbits by using the manifolds of the vertical and horizontal Lyapunov orbits and the transit cones. The relationship between the manifolds of the libration point orbits and the boundary of the transit orbits is discovered. By using this relationship, a fast algorithm for detecting the boundary of the transit orbits is developed. Moreover, this boundary is parametrized by using Fourier series, which makes easy to use the conclusions of this paper in future trajectory optimization and mission design. All the analyses in this paper are based on the Sun?CEarth CR3BP, but the methods introduced here can be extended to any CR3BPs.     

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.     

Dynamical model of binary asteroid systems through patched three-body problems

The paper presents a strategy for trajectory design in the proximity of a binary asteroid pair. A novel patched approach has been used to design trajectories in the binary system, which is modeled by means of two different three-body systems. The model introduces some degrees of freedom with respect to a classical two-body approach and it is intended to model to higher accuracy the peculiar dynamical properties of such irregular and low gravity field bodies, while keeping the advantages of having a full analytical formulation and low computational cost required. The neighborhood of the asteroid couple is split into two regions of influence where two different three-body problems describe the dynamics of the spacecraft. These regions have been identified by introducing the concept of surface of equivalence (SOE), a three-dimensional surface that serves as boundary between the regions of influence of each dynamical model. A case of study is presented, in terms of potential scenario that may benefit of such an approach in solving its mission analysis. Cost-effective solutions to land a vehicle on the surface of a low gravity body are selected by generating Poincaré maps on the SOE, seeking intersections between stable and unstable manifolds of the two patched three-body systems.     

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.     

A fast tour design method using non-tangent v-infinity leveraging transfer

The announced missions to the Saturn and Jupiter systems renewed the space community interest in simple design methods for gravity assist tours at planetary moons. A key element in such trajectories are the V-Infinity Leveraging Transfers (VILT) which link simple impulsive maneuvers with two consecutive gravity assists at the same moon. VILTs typically include a tangent impulsive maneuver close to an apse location, yielding to a desired change in the excess velocity relative to the moon. In this paper we study the VILT solution space and derive a linear approximation which greatly simplifies the computation of the transfers, and is amenable to broad global searches. Using this approximation, Tisserand graphs, and heuristic optimization procedure we introduce a fast design method for multiple-VILT tours. We use this method to design a trajectory from a highly eccentric orbit around Saturn to a 200-km science orbit at Enceladus. The trajectory is then recomputed removing the linear approximation, showing a Δv change of <4%. The trajectory is 2.7 years long and comprises 52 gravity assists at Titan, Rhea, Dione, Tethys, and Enceladus, and several deterministic maneuvers. Total Δv is only 445 m/s, including the Enceladus orbit insertion, almost 10 times better then the 3.9 km/s of the Enceladus orbit insertion from the Titan–Enceladus Hohmann transfer. The new method and demonstrated results enable a new class of missions that tour and ultimately orbit small mass moons. Such missions were previously considered infeasible due to flight time and Δv constraints.     

Ballistic capture into distant retrograde orbits around Phobos: an approach to entering orbit around Phobos without a critical maneuver near the moon

This paper discusses an approach for designing missions to Phobos that do not require a critical maneuver in proximity of the moon. A low-energy transfer is designed that utilizes the aspherical mass distribution of Phobos to capture a spacecraft into a distant retrograde orbit (DRO) for the mission duration. The process for generating stable DROs in the Mars–Phobos system is discussed along with the method for generating and surveying a set of ballistic capture trajectories (BCTs) for DROs with altitudes between 0.5 and 14 km above Phobos. Statistical analysis of the BCT set reveals options for designing a mission to the desired DRO. This approach can be used in any three-body system when a significant perturbation is present, such as Phobos’ aspherical co-rotating gravity field.     

Detailed survey of the phase space around Nix and Hydra

We present a detailed survey of the dynamical structure of the phase space around the new moons of the Pluto–Charon system. The spatial elliptic restricted three-body problem was used as model and stability maps were created by chaos indicators. The orbital elements of the moons are in the stable domain on the semimajor axis, eccentricity and inclination spaces. The structures related to the 4:1 and 6:1 mean motion resonances are clearly visible on the maps. They do not contain the positions of the moons, confirming previous studies. We showed the possibility that Nix might be in the 4:1 resonance if its argument of pericentre or longitude of node falls in a certain range. The results strongly suggest that Hydra is not in the 6:1 resonance for arbitrary values of the argument of pericentre or longitude of node.     

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.     

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

Numerical orbit integrations have been conducted to characterize the types of trajectories in the one-dimensional Newtonian three-body problem with equal masses and negative energy. Essentially three different types of motions were found to exist. They may be classified according to the duration of the bound three-body state. There are zero-lifetime predictable trajectories, finite lifetime apparently chaotic orbits, and infinite lifetime quasi-periodic motions. The quasi-periodic orbits are confined to the neighbourhood of Schubart's stable periodic orbit. For all other trajectories the final state is of the type binary + single particle in both directions of time. The boundaries of the different orbit-type regions seem to be sharp. We present statistical results for the binding energies and for the duration of the bound three-body state. Properties of individual orbits are also summarized in the form of various graphical maps in a two-dimensional grid of parameters defining the orbit. Supported by the Academy of Finland.     

Trajectory design strategies applied to temporary comet capture including Poincaré maps and invariant manifolds

Temporary satellite capture (TSC) of Jupiter-family comets has been a focus of investigation within the astronomy community for decades. More recently, TSC has been approached from the perspective of dynamical systems theory, within the context of the circular restricted three-body problem (CR3BP). Thus, this problem serves as a testbed for exploring techniques that support trajectory design in similar dynamical regimes. In particular, an association between the invariant manifolds of libration point orbits and the paths of comets that experience TSC has been explored. In this investigation, TSC is further examined from the perspective of transit, that is, transition through the gateways associated with the collinear libration points, in the three-body problem. Periapsis Poincaré maps, previously employed for trajectory design in several investigations, are used to deliver insight into the nature of transit trajectories for energy levels near those associated with several Jupiter-family comets. The evolution of transit trajectories with increasing energy is explored, and the existence of solutions with similar characteristics to the paths of comets P/1996 R2, 82P/Gehrels 3, and 147P/Kushida–Muramatsu is demonstrated within the context of the planar CR3BP using planar periapsis maps. During TSC, the path of comet 111P/Helin–Roman–Crockett is highly inclined with respect to Jupiter; the motion of this comet is examined relative to invariant manifolds in the spatial CR3BP. A method to display the information contained in higher-dimensional Poincaré maps is also demonstrated, and is employed to locate a trajectory possessing the same qualitative characteristics as the path of 111P/Helin–Roman–Crockett.     

Free Return Trajectories in Lunar Missions

Four types of symmetric free return trajectories in the planar circular restricted three-body problem are computed and compared with each other. One of these four types is most applicable in practice. Concentrating ourselves on this special type of free return trajectory, the corresponding planar asymmetric cases are studied. Then the studies are generalized to the three-dimensional case. The restrictions on the inclination angle of the probe at the perilune are discussed. It is found that the maximum inclination at the perilune between the probe's orbit plane and the Moon's orbit plane is restricted by the heights of the perigee and the perilune of the free return trajectory. However, the inclination at the perigee is nearly not affected by them. At last, a strategy to design free return trajectories in the real Earth-Moon system is proposed. Some numerical simulations are done to show the feasibility of this orbit design strategy. The discussions from the planar case to the spatial case and then to the real force model can also be applied to the other three types of planar free return trajectories.     

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

Generalization of Levi-Civita regularization in the restricted three-body problem

A family of polynomial coupled function of n degree is proposed, in order to generalize the Levi-Civita regularization method, in the restricted three-body problem. Analytical relationship between polar radii in the physical plane and in the regularized plane are established; similar for polar angles. As a numerical application, trajectories of the test particle using polynomial functions of 2,3,…,8 degree are obtained. For the polynomial of second degree, the Levi-Civita regularization method is found.     

Evolution of the Janus-Epimetheus coorbital resonance due to torques from Saturn's ring

We analyze the interactions between Saturn's coorbital satellites, Janus and Epimetheus, and the outer edge of the A ring, which is presumably maintained by these moons at their 7:6 resonance. Using two distinct but conceptually related methods, we show that ring torques are driving these satellites into a tighter lock. Unless there is a counterbalancing force which we have neglected, their orbital configuration will evolve from the current horseshoe-type lock to one of tadpole orbits around a single Lagrange point in ～20 myr. This finding adds an additional member to the list of short time scale problems associated with the interactions between Saturn's rings and its inner moons     

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.     

Multi-SunSynchronous Orbits in the Solar System

Performances of a planetary observation system are strongly related to the choice of the orbit used. Trajectories with characteristics of periodicity are very useful for the assessment of time-varying phenomena and thus Periodic SunSynchronous and Periodic Multi-SunSynchronous Orbits are particularly suitable to this end. In this paper, the research into these kinds of orbits, previously proposed for the Earth and Mars, has been extended to planets of the Solar System and to their principal moons. In general, these trajectories are typically obtained under the hypothesis that the J₂ harmonic is predominant with respect to the other orbital perturbations, since this allows an analytical solution. However, the hypothesis of J₂ predominant is not always verified in the Solar System and so analytical techniques must be replaced by numerical simulations. Interesting results have been obtained for the planets Mars and Jupiter and for the moons Europa, Callisto and Titan, where periodic trajectories with reduced revisit times and low altitudes have been found. These solutions allow the observation of time-varying phenomena with high spatial and temporal resolution.