Optimisation of Low-Thrust and Hybrid Earth-Moon Transfers

This paper presents an optimization procedure to generate fast and low-∆v Earth-Moon transfer trajectories, by exploiting the multi-body dynamics of the Sun-Earth-Moon system. Ideal (first-guess) trajectories are generated at first, using two coupled planar circular restricted three-body problems, one representing the Earth-Moon system, and one representing the Sun-Earth. The trajectories consist of a first ballistic arc in the Sun-Earth system, and a second ballistic arc in the Earth-Moon system. The two are connected at a patching point at one end (with an instantaneous ∆v), and they are bounded at Earth and Moon respectively at the other end. Families of these trajectories are found by means of an evolutionary optimization method. Subsequently, they are used as first-guess for solving an optimal control problem, in which the full three-dimensional 4-body problem is introduced and the patching point is set free. The objective of the optimisation is to reduce the total ∆v, and the time of flight, together with introducing the constraints on the transfer boundary conditions and of the considered propulsion technology. Sets of different optimal trajectories are presented, which represents trade-off options between ∆v and time of flight. These optimal transfers include conventional solar-electric low-thrust and hybrid chemical/solar-electric high/low-thrust, envisaging future spacecraft that can carry both systems. A final comparison is made between the optimal transfers found and only chemical high-thrust optimal solutions retrieved from literature

A Heuristic Strategy to Compute Ensemble of Trajectories for 3D Low Cost Earth-Moon Transfers

The problem of finding optimal trajectories is essential for modern space mission design. When considering multibody gravitational dynamics and exploiting both low-thrust and high-thrust and alternative forms of propulsion such as solar sailing, sets of good initial guesses are fundamental for the convergence to local or global optimal solutions, using both direct or indirect methods available to solve the optimal control problem. This paper deals with obtaining preliminary trajectories that are designed to be good initial guesses as input to search optimal low-energy short-time Earth-Moon transfers with ballistic capture. A more realistic modelling is introduced, in which the restricted four-body system Sun-Earth-Moon-Spacecraft is decoupled in two patched planar Circular Restricted Three-Body Problems, taking into account the inclination of the orbital plane of the Moon with respect to the ecliptic. We present a heuristic strategy based on the hyperbolic invariant manifolds of the Lyapunov orbits around the Lagrangian points of the Earth- Moon system to obtain ballistic capture orbits around the Moon that fulfill specific mission requirements. Moreover, quasi-periodic orbits of the Sun-Earth system are exploited using a genetic algorithm to find optimal solutions with respect to total Dv, time of flight and altitude at departure. Finally, the procedure is illustrated and the full transfer trajectories assessed in view of relevant properties. The proposed methodology provides sets of low-cost and shorttime initial guesses to serve as inputs to compute fully optimized three-dimensional solutions considering different propulsion technologies, such as low, high, and hybrid thrust, and/or using more realistic models

Trajectory Design Combining Invariant Manifolds with Discrete Mechanics and Optimal Control

A mission design technique that combines invariant manifold techniques, discrete mechanics, and optimal control produces locally optimal low-energy trajectories. Previously, invariant manifolds of the planar circular restricted three-body problem have been used to design trajectories with relatively small midcourse change in velocity ΔV. A different method of using invariant manifolds is explored to design trajectories directly in the four-body problem. Then, using the local optimal control method DMOC (Discrete Mechanics and Optimal Control), it is possible to reduce the midcourse ΔV to zero. The influence of different boundary conditions on the optimal trajectory is also demonstrated. These methods are tested on a trajectory that begins in Earth orbit and ends in ballistic capture at the moon. Impulsive DMOC trajectories require up to 19% less ΔV than trajectories using a Hohmann transfer. When applied to low-thrust trajectories, DMOC produces an improvement of up to 59% in the mass fraction and 22% in travel time when compared with results from shooting methods

RAPID SPACE TRAJECTORY GENERATION USING A FOURIER SERIES SHAPE-BASED APPROACH

With the insatiable curiosity of human beings to explore the universe and our solar system, it is essential to benefit from larger propulsion capabilities to execute efficient transfers and carry more scientific equipment. In the field of space trajectory optimization the fundamental advances in using low-thrust propulsion and exploiting the multi-body dynamics has played pivotal role in designing efficient space mission trajectories. The former provides larger cumulative momentum change in comparison with the conventional chemical propulsion whereas the latter results in almost ballistic trajectories with negligible amount of propellant. However, the problem of space trajectory design translates into an optimal control problem which is, in general, time-consuming and very difficult to solve. Therefore, the goal of the thesis is to address the above problem by developing a methodology to simplify and facilitate the process of finding initial low-thrust trajectories in both two-body and multi-body environments. This initial solution will not only provide mission designers with a better understanding of the problem and solution but also serves as a good initial guess for high-fidelity optimal control solvers and increases their convergence rate. Almost all of the high-fidelity solvers enjoy the existence of an initial guess that already satisfies the equations of motion and some of the most important constraints. Despite the nonlinear nature of the problem, it is sought to find a robust technique for a wide range of typical low-thrust transfers with reduced computational intensity. Another important aspect of our developed methodology is the representation of low-thrust trajectories by Fourier series with which the number of design variables reduces significantly. Emphasis is given on simplifying the equations of motion to the possible extent and avoid approximating the controls. These facts contribute to speeding up the solution finding procedure. Several example applications of two and three-dimensional two-body low-thrust transfers are considered. In addition, in the multi-body dynamic, and in particular the restricted-three-body dynamic, several Earth-to-Moon low-thrust transfers are investigated

Low-thrust trajectories design for the European Student Moon Orbiter mission

The following paper presents the mission analysis studies performed for the phase A of the solar electric propulsion option of the European Student Moon Orbiter (ESMO) mission. ESMO is scheduled to be launched in 2011, as an auxiliary payload on board of Ariane 5. Hence the launch date will be imposed by the primary payload. A method to efficiently assess wide launch windows for the Earth-Moon transfer is presented here. Sets of spirals starting from the GTO were propagated forward with a continuous tangential thrust until reaching an apogee of 280,000 km. Concurrently, sets of potential Moon spirals were propagated backwards from the lunar orbit injection. The method consists of ranking all the admissible lunar spiral-down orbits that arrive to the target orbit with a simple tangential thrust profile after a capture through the L1 Lagrange point. The 'best' lunar spiral is selected for each Earth spiral. Finally,comparing the value of the ranking function for each launch date, the favourable and unfavourable launch windows are identified

Low-Thrust Lyapunov to Lyapunov and Halo to Halo with L2L^2-Minimization

In this work, we develop a new method to design energy minimum low-thrust missions (L2-minimization). In the Circular Restricted Three Body Problem, the knowledge of invariant manifolds helps us initialize an indirect method solving a transfer mission between periodic Lyapunov orbits. Indeed, using the PMP, the optimal control problem is solved using Newton-like algorithms finding the zero of a shooting function. To compute a Lyapunov to Lyapunov mission, we first compute an admissible trajectory using a heteroclinic orbit between the two periodic orbits. It is then used to initialize a multiple shooting method in order to release the constraint. We finally optimize the terminal points on the periodic orbits. Moreover, we use continuation methods on position and on thrust, in order to gain robustness. A more general Halo to Halo mission, with different energies, is computed in the last section without heteroclinic orbits but using invariant manifolds to initialize shooting methods with a similar approach

Design of a Multi-Moon Orbiter

The Multi-Moon Orbiter concept is introduced, wherein a single spacecraft orbits several moons of Jupiter, allowing long duration observations. The ΔV requirements for this mission can be low if ballistic captures and resonant gravity assists by Jupiter’s moons are used. For example, using only 22 m/s, a spacecraft initially injected in a jovian orbit can be directed into a capture orbit around Europa, orbiting both Callisto and Ganymede enroute. The time of flight for this preliminary trajectory is four years, but may be reduced by striking a compromise between fuel and time optimization during the inter-moon transfer phases

Construction and Verification of a Solution of the 8th Global Trajectory Optimization Competition Problem. TEAM 13: GlasgowJena+

This paper describes the methodology to find and verify the solution to the 8th Global Trajectory Optimization Competition (GTOC) problem, developed by Team 13, GlasgowJena+. We chose a stochastic approach to quickly assess a large number (about 1010) of 3-spacecraft formations. A threshold was used to select promising solutions for further optimization. Our search algorithm (implemented in Java) is based on three C++ algorithms called via Java native interface (JNI). A great deal was given to the verification process, which became a core part of our solution. Our final solution has a performance index of 75.9710kmJ=×, 40 distinct observations, and the sum of the final masses of the three spacecraft is 5846.57 kg

A Homotopy-Based Method for Optimization of Hybrid High-Low Thrust Trajectories

Space missions require increasingly more efficient trajectories to provide payload transport and mission goals by means of lowest fuel consumption, a strategic mission design key-point. Recent works demonstrated that the combined (or hybrid) use of chemical and electrical propulsion can give important advantages in terms of fuel consumption, without losing the ability to reach other mission objectives: as an example the Hohmann Spiral Transfer, applied in the case of a transfer to GEO orbit, demonstrated a fuel mass saving between 5-10% of the spacecraft wet mass, whilst satisfying a pre-set boundary constraint for the time of flight. Nevertheless, methods specifically developed for optimizing space trajectories considering the use of hybrid high-low thrust propulsion systems have not been extensively developed, basically because of the intrinsic complexity in the solution of optimal problem equations with existent numerical methods. The study undertaken and presented in this paper develops a numerical strategy for the optimization of hybrid high-low thrust space trajectories. An indirect optimization method has been developed, which makes use of a homotopic approach for numerical convergence improvement. The adoption of a homotopic approach provides a relaxation to the optimal problem, transforming it into a simplest problem to solve in which the optimal problem presents smoother equations and the shooting function acquires an increased convergence radius: the original optimal problem is then reached through a homotopy parameter continuation. Moreover, the use of homotopy can make possible to include a high thrust impulse (treated as velocity discontinuity) to the low thrust optimal control obtained from the indirect method. The impulse magnitude, location and direction are obtained following from a numerical continuation in order to minimize the problem cost function. The initial study carried out in this paper is finally correlated with particular test cases, in order to validate the work developed and to start investigating in which cases the effectiveness of hybrid-thrust propulsion subsists