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 target for Venus – set oriented computation of energy efficient low thrust trajectories

Recently new techniques for the design of energy efficient trajectories for space missions have been proposed that are based on the circular restricted three body problem as the underlying mathematical model. These techniques exploit the structure and geometry of certain invariant sets and associated invariant manifolds in phase space to systematically construct energy efficient flight paths. In this paper, we extend this model in order to account for a continuously applied control force on the spacecraft as realized by certain low thrust propulsion systems. We show how the techniques for the trajectory design can be suitably augmented and compute approximations to trajectories for a mission to Venus.     

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.     

Gravity assist space pruning based on differential algebra

In this paper a differential algebra version of the gravity assist space pruning algorithm is presented. The use of differential algebraic techniques is proposed to overcome the two main drawbacks of the existing algorithm, i.e., the steep increase of the number of function evaluations with the number of planets involved in the transfer, and the use of a bounding procedure that relies on Lipschitzian tolerances. Differential algebra allows us to process boxes in place of grid points, and to substitute pointwise evaluations of the constraint functions with their Taylor expansions. Thanks to the particular instance of multi-gravity assist problems dealt with, all the planet-to-planet legs can be treated independently, and forward and backward constraining can be applied. The proposed method is applied to preprocess the search space of sample interplanetary transfers and it also serves as a stepping stone towards a fully rigorous treatment of the pruning process based on Taylor models.     

Using FLI maps for preliminary spacecraft trajectory design in multi-body environments

Fast Lyapunov Indicator (FLI) maps are presented as a tool for solving spacecraft preliminary trajectory design problems in multi-body environments with long-term stability requirements. In particular, the FLI maps are shown to provide a global overview of the dynamics in the restricted three-body problem that can guide mission designers in selecting long-term stable regions of phase space which are inherently more robust to model parameter perturbations. The FLI is also shown to numerically detect the normally hyperbolic manifolds associated with unstable periodic orbits. These, in turn, provide a global map of the principal heteroclinic connections between the various resonance regions which form the basic backbone of dynamical transfers design. Examples of maps and transfers are provided in the restricted three-body problem modeling the Jupiter–Europa system.     

A graph based methodology for mission design

A space mission design methodology is presented, where initial and final orbits are connected through segments of periodic orbits. After a discretization of the solution space, the problem of mission design is transformed into an equivalent combinatorial optimization problem. Specifically, a graph is constructed that represents periodic orbits connected by the execution of impulsive maneuvers. A low computational complexity algorithm for this transformation is introduced. An efficient combinatorial optimization algorithm that solves the shortest path problem is described. Subject to the initial discretization of the solution space, an optimal sequence of coastal arcs is determined for a low total Delta-V mission. Finally, the proposed methodology is applied to the design of a hypothetical Saturn?CTitan system mission.     

Boundary-value problem formulations for computing invariant manifolds and connecting orbits in the circular restricted three body problem

We demonstrate the remarkable effectiveness of boundary value formulations coupled to numerical continuation for the computation of stable and unstable manifolds in systems of ordinary differential equations. Specifically, we consider the circular restricted three-body problem (CR3BP), which models the motion of a satellite in an Earth–Moon-like system. The CR3BP has many well-known families of periodic orbits, such as the planar Lyapunov orbits and the non-planar vertical and halo orbits. We compute the unstable manifolds of selected vertical and halo orbits, which in several cases leads to the detection of heteroclinic connections from such a periodic orbit to invariant tori. Subsequent continuation of these connecting orbits with a suitable end point condition and allowing the energy level to vary leads to the further detection of apparent homoclinic connections from the base periodic orbit to itself, or the detection of heteroclinic connections from the base periodic orbit to other periodic orbits. Some of these connecting orbits are of potential interest in space mission design.     

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.     

Efficient design of direct low-energy transfers in multi-moon systems

In this contribution, an efficient technique to design direct (i.e., without intermediate flybys) low-energy trajectories in multi-moon systems is presented. The method relies on analytical two-body approximations of trajectories originating from the stable and unstable invariant manifolds of two coupled circular restricted three-body problems. We provide a means to perform very fast and accurate computations of the minimum-cost trajectories between two moons. Eventually, we validate the methodology by comparison with numerical integrations in the three-body problem. Motivated by the growing interest in the robotic exploration of the Jovian system, which has given rise to numerous studies and mission proposals, we apply the method to the design of minimum-cost low-energy direct trajectories between Galilean moons, and the case study is that of Ganymede and Europa.     

Cooperative evolutionary algorithm for space trajectory optimization

A hybrid evolutionary algorithm which synergistically exploits differential evolution, genetic algorithms and particle swarm optimization, has been developed and applied to spacecraft trajectory optimization. The cooperative procedure runs the three basic algorithms in parallel, while letting the best individuals migrate to the other populations at prescribed intervals. Rendezvous problems and round-trip Earth–Mars missions have been considered. The results show that the hybrid algorithm has better performance compared to the basic algorithms that are employed. In particular, for the rendezvous problem, a 100% efficiency can be obtained both by differential evolution and the genetic algorithm only when particular strategies and parameter settings are adopted. On the other hand, the hybrid algorithm always attains the global optimum, even though nonoptimal strategies and parameter settings are adopted. Also the number of function evaluations, which must be performed to attain the optimum, is reduced when the hybrid algorithm is used. In the case of Earth–Mars missions, the hybrid algorithm is successfully employed to determine mission opportunities in a large search space.     

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

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

Low Energy Transfer to the Moon

In 1991, the Japanese Hiten mission used a low energy transfer with a ballistic capture at the Moon which required less Vthan a standard Hohmann transfer. In this paper, we apply the dynamical systems techniques developed in our earlier work to reproduce systematically a Hiten-like mission. We approximate the Sun–Earth–Moon-spacecraft 4-body system as two 3-body systems. Using the invariant manifold structures of the Lagrange points of the 3-body systems, we are able to construct low energy transfer trajectories from the Earth which execute ballistic capture at the Moon. The techniques used in the design and construction of this trajectory may be applied in many situations.     

Digital complex correlator for a C-band polarimetry survey

The international Galactic Emission Mapping project aims to map and characterize the polarization field of the Milky Way. In Portugal it will map the sky polarized emission of the Northern Hemisphere in C-band and provide templates for map calibration and foreground control of microwave space data to be provided by ESA Planck Surveyor mission and later missions. The receiver system is equipped with a novel receiver with a full digital back-end using a low-cost Field Programmable Gate Array without compromising its performance relation. This new digital backend comprises a base-band complex cross-correlator outputting the four Stokes parameters of the incoming polarized radiation. In this document we describe the design and implementation of the complex correlator using the FPGA and the dedicated digitizers at each receiver arm, detailing the method applied at the several algorithm stages. This correlator is suitable for large sky area polarization continuum surveys.     

Quasi-terminator orbits near primitive bodies

Quasi-terminator orbits are introduced as a class of quasi-periodic trajectories in the solar radiation pressure (SRP) perturbed Hill dynamics. These orbits offer significant displacements along the Sun-direction without the need for station-keeping maneuvers. Thus, quasi-terminator orbits have application to primitive-body mapping missions, where a variety of observation geometries relative to the Sun (or other directions) can be achieved. This paper describes the characteristics of these orbits as a function of normalized SRP strength and invariant torus frequencies and presents a discussion of mission design considerations for a global surface mapping orbit design.     

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.     

Comparison of low-energy lunar transfer trajectories to invariant manifolds

In this study, transfer trajectories from the Earth to the Moon that encounter the Moon at various flight path angles are examined, and lunar approach trajectories are compared to the invariant manifolds of selected unstable orbits in the circular restricted three-body problem. Previous work focused on lunar impact and landing trajectories encountering the Moon normal to the surface, and this research extends the problem with different flight path angles in three dimensions. The lunar landing geometry for a range of Jacobi constants is computed, and approaches to the Moon via invariant manifolds from unstable orbits are analyzed for different energy levels.     

Analytical design methods for transfer trajectories between the Earth and the Lunar Orbital Station

Lunar Orbital Station (LOS) is proposed as support of manned lunar exploration missions. A fast-converging iteration method for determining the initial conditions of two-impulse transfer trajectories between the Earth and the LOS is proposed based on the patched conic approach. In the Earth phase, near Earth state is connected with the state at the lunar sphere of influence (LSOI) based on the relationship between the initial and terminal orbital state. Then, an analytical algorithm is proposed to find the state vector at LSOI, such to satisfy the LOS orbital constraint. An iterative process is finally adopted to generate favorable initial solutions that satisfy the constraint near the Earth and at the perilune. The algorithm convergence is investigated, and two types of transfer trajectories are found for both Earth-LOS and LOS-Earth transfer. Based on the algorithm, orbital transfer windows, velocity impulse and time of flight are analyzed in the typical years 2025 and 2034. At last, the initial solution is corrected with a high fidelity model based on the active-set method, which shows the precision of this algorithm. The novel procedure for the transfer trajectories design and the analytic result can be used as a basis for rapid mission evaluation and design for future manned lunar missions based on the LOS.     

New Challenges to Trajectory Design by the Use of Electric Propulsion and Other New Means of Wandering in the Solar System

The design of spacecraft trajectories is a crucial part of a space mission design. Often the mission goal is tightly related to the spacecraft trajectory. A geostationary orbit is indeed mandatory for a stationary equatorial position. Visiting a solar system planet implies that a proper trajectory is used to bring the spacecraft from Earth to the vicinity of the planet. The first planetary missions were based on conventional trajectories obtained with chemical engine rockets. The manoeuvres could be considered 'impulsive' and clear limitations to the possible missions were set by the energy required to reach certain orbits. The gravity-assist trajectories opened a new way of wandering through the solar system, by exploiting the gravitational field of some planets. The advent of other propulsion techniques, as electric or ion propulsion and solar sail, opened a new dimension to the planetary trajectory, while at the same time posing new challenges. These 'low thrust' propulsion techniques cannot be considered 'impulsive' anymore and require for their study mathematical techniques which are substantially different from before. The optimisation of such trajectories is also a new field of flight dynamics, which involves complex treatments especially in multi-revolution cases as in a lunar transfer trajectory. One advantage of these trajectories is that they allow to explore regions of space where different bodies gravitationally compete with each other. We can exploit therefore these gravitational perturbations to save fuel or reduce time of flight. The SMART-1 spacecraft, first European mission to the Moon, will test for the first time all these techniques. The paper is a summary report on various activities conducted by the project team in these areas.     

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.