Broad-search algorithms for the spacecraft trajectory design of Callisto–Ganymede–Io triple flyby sequences from 2024 to 2040, Part I: Heuristic pruning of the search space

Triple flybys of the Galilean moons of Jupiter can capture a spacecraft into orbit about Jupiter or quickly adjust the Jupiter-centered orbit of an already captured spacecraft. Because Callisto does not participate in the Laplace resonance among Ganymede, Europa, and Io, triple flyby sequences involving gravity-assists of Callisto, Ganymede, and Io occur only aperiodically for limited time windows. An exhaustive search of triple-flyby trajectories over a 16-year period from 2024 to 2040 using “blind” searching would require 8,415,358 Lambert function calls to find only 127,289 possible triple flyby trajectories. Because most of these Lambert function calls would not converge to feasible solutions, it is much more efficient to prune the solution space using a heuristic algorithm and then direct a much smaller number of Lambert function calls to find feasible triple flyby solutions. The novel “Phase Angle Pruning Heuristic” is derived and used to reduce the search space by 99%.     

面向木星卫星交会任务的探测器飞行路径规划

针对探测器在木星系统内多次借力的飞行路径和轨道优化设计问题,提出了一种基于三层优化思想的飞行路径规划方法,该方法可根据给定的任务约束和交会目标,自动搜索探测器在木星系统内的借力飞行序列,同时完成标称飞行轨道的优化设计。首先,文章在给定轨道动力学模型和木卫借力模型基础上,建立了面向木卫交会任务的两次借力飞行轨道优化设计模型和求解方法;然后,采用结合遗传算法、全局遍历和贪婪算法的三层优化设计思路,给出了一种环木飞行路径规划方法;最后,以木星四颗卫星的交会任务为例进行了仿真分析。仿真结果表明：针对木卫的交会任务,探测器速度增量需求随木卫借力次数的增多,呈现先显著减小后逐渐增大的现象;探测器采用多次木卫借力的策略,可显著降低探测器的速度增量需求;探测器速度增量达到最优之后,借力目标收敛于交会目标,且速度增量随借力次数的进一步增多而逐渐增大。     

A low-thrust transfer between the Earth–Moon and Sun–Earth systems based on invariant manifolds

A low-energy, low-thrust transfer between two halo orbits associated with two coupled three-body systems is studied in this paper. The transfer is composed of a ballistic departure, a ballistic insertion and a powered phase using low-thrust propulsion to connect these two trajectories. The ballistic departure and insertion are computed by constructing the unstable and stable invariant manifolds of the corresponding halo orbits, and a complete low-energy transfer based on the patched invariant manifolds is optimized using the particle swarm optimization (PSO) algorithm on the criterion of smallest velocity discontinuity and limited position discontinuity (less than 1 km). Then, the result is expropriated as the boundary conditions for the subsequent low-thrust trajectory design. The fuel-optimal problem is formulated using the calculus of variations and Pontryagin's Maximum Principle in a complete four-body dynamical environment. Then, a typical bang–bang control is derived and solved using the indirect method combined with a homotopic technique. The contributions of the present work mainly consist of two points. Firstly, the global search method proposed in this paper is simply handled using the PSO algorithm, a number of feasible solutions in a fairly wide range can be delivered without a priori or perfect knowledge of the transfers. Secondly, the indirect optimization method is used in the low-thrust trajectory design and the derivations of the first-order necessary conditions are simplified with a modified controlled, restricted four-body model.     

Design of transfer trajectories between resonant orbits in the Earth–Moon restricted problem

The application of dynamical systems techniques to mission design has demonstrated that employing invariant manifolds and resonant flybys enables previously unknown trajectory options and potentially reduces the ΔV$\mathrm{\Delta }V$ requirements. In this investigation, planar and three-dimensional resonant orbits are analyzed and cataloged in the Earth–Moon system and the associated invariant manifold structures are computed and visualized with the aid of higher-dimensional Poincaré maps. The relationship between the manifold trajectories associated with multiple resonant orbits is explored through the maps with the objective of constructing resonant transfer arcs. As a result, planar and three-dimensional homoclinic- and heteroclinic-type trajectories between unstable periodic resonant orbits are identified in the Earth–Moon system. To further illustrate the applicability of 2D and 3D resonant orbits in preliminary trajectory design, planar transfers to the vicinity of L₅ and an out-of-plane transfer to a 3D periodic orbit, one that tours the entire Earth–Moon system, are constructed. The design process exploits the invariant manifolds associated with orbits in resonance with the Moon as transfer mechanisms.     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.     

Design of Interplanetary Transfers with Passive Gravity Assists and Deep Space Maneuvers

In the paper, the problem of designing interplanetary trajectories with several swing-bys and deep-space maneuvers is solved using the method of virtual trajectories developed by the authors. The algorithms for the calculation of both heliocentric and planetocentric trajectory arcs are presented, including the case of resonant trajectories. The results of applying the method of virtual trajectories to the problem of designing an interplanetary transfer to Jupiter are given and compared with the baseline trajectories for the Juno, Europa Clipper, and Laplace missions.     

Halo orbits in the vicinity of the L 2 libration point in the Sun-Earth system

Methods are proposed for constructing the orbits of spacecraft remaining for long periods of time in the vicinity of the L ₂ libration point in the Sun-Earth system (so-called halo orbits), and the trajectories of uncontrolled flights from low near-Earth orbits to halo orbits. Halo orbits and flight trajectories are constructed in two stages: A suitable solution to a circular restricted three-body problem is first constructed and then transformed into the solution for a restricted four-body problem in view of the real motions of the Sun, Earth, and Moon. For a halo orbit, its prototype in the first stage is a combination of a periodic Lyapunov solution in the vicinity of the L ₂ point and lying in the plane of large-body motion, with the solution for the linear second-order system describing small deviations of the spacecraft from this plane along the periodic solution. The desired orbit is found as the solution to the three-body problem best approximating the prototype in the mean square. The constructed orbit serves as a similar prototype in the second stage. In both stages, the approximating solution is constructed by continuation along a parameter that is the length of the approximation interval. Flight trajectories are constructed in a similar manner. The prototype orbit in the first stage is a combination of a solution lying in the plane of large-body motion and a solution for a linear second-order system describing small deviations of the spacecraft from this plane. The planar solution begins near the Earth and over time tends toward the Lyapunov solution existing in the vicinity of the L ₂ point. The initial conditions of both prototypes and the approximating solutions correspond to the spacecraft’s departure from a low near-Earth orbit at a given distance, perigee, and inclination.     

Optimization of low-thrust orbit transfers with constraints

Approximate numerical methods of optimization are presented for multi-orbit noncoplanar orbit transfers of low-thrust spacecraft. The linear representation of derivatives of boundary parameters is used in the vicinity of a reference trajectory with its discretization into small segments. For each segment a set of pseudo-impulses is introduced, representing possible directions of the thrust vector. In order to solve essentially nonlinear problems, the iterative process of upgrading the reference trajectory is used. At each iteration the linear programming problem of high dimensionality is solved, its boundary conditions being represented in the form of a linear matrix equation. Interval constraints are considered in the form of blocking the propulsion system operation on shadow segments of the orbit, as well as cycling constraints, and constraints on total duration of maneuvers at certain trajectory segments. The results of comparison with solutions obtained by other methods are presented together with examples illustrating the convergence of iterative processes. Optimizations of the trajectories for launching geosynchronous satellites to their orbits and of the trajectories of a noncoplanar transfer from low to high-elliptic Molniya orbit are considered under these constraints.     

Broad-search algorithms for the spacecraft trajectory design of Callisto–Ganymede–Io triple flyby sequences from 2024 to 2040, Part II: Lambert pathfinding and trajectory solutions

Triple-satellite-aided capture employs gravity-assist flybys of three of the Galilean moons of Jupiter in order to decrease the amount of ΔV$\mathrm{\Delta }V$ required to capture a spacecraft into Jupiter orbit. Similarly, triple flybys can be used within a Jupiter satellite tour to rapidly modify the orbital parameters of a Jovicentric orbit, or to increase the number of science flybys. In order to provide a nearly comprehensive search of the solution space of Callisto–Ganymede–Io triple flybys from 2024 to 2040, a third-order, Chebyshev's method variant of the p-iteration solution to Lambert's problem is paired with a second-order, Newton–Raphson method, time of flight iteration solution to the _V∞$V_{\infty }$-matching problem. The iterative solutions of these problems provide the orbital parameters of the Callisto–Ganymede transfer, the Ganymede flyby, and the Ganymede–Io transfer, but the characteristics of the Callisto and Io flybys are unconstrained, so they are permitted to vary in order to produce an even larger number of trajectory solutions. The vast amount of solution data is searched to find the best triple-satellite-aided capture window between 2024 and 2040.     

Optimization of low-thrust transfers in the three body problem

We consider transfers with low thrust in an arbitrary field of forces. The modified method of transporting trajectory [1–4] is used for optimization of the transfers. The complexity of finding the transporting trajectory of a preset type can be the main obstacle to application of this method. This challenge is solved for the three-body problem in the Hill motion model. Numerical analysis of the method is performed using an example of the transfers to halo-orbits around the solar-terrestrial libration points.     

Earth–Moon libration point orbit stationkeeping: Theory,modeling, and operations

Collinear Earth–Moon libration points have emerged as locations with immediate applications. These libration point orbits are inherently unstable and must be maintained regularly which constrains operations and maneuver locations. Stationkeeping is challenging due to relatively short time scales for divergence, effects of large orbital eccentricity of the secondary body, and third-body perturbations. Using the Acceleration Reconnection and Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) mission orbit as a platform, the fundamental behavior of the trajectories is explored using Poincaré maps in the circular restricted three-body problem. Operational stationkeeping results obtained using the Optimal Continuation Strategy are presented and compared to orbit stability information generated from mode analysis based in dynamical systems theory.     

Venus transfer design by combining invariant manifolds and low-thrust arcs

The design of interplanetary trajectories based on patched circular restricted three body models is gradually becoming a valuable alternative to the classical patched conic approach. The main advantage offered by such a model is the possibility to exploit the manifold dynamics to move naturally far from or toward a body. Generally, propulsive maneuvers are required to match these structures. Low-thrust arcs offer the possibility to have a significant propellant mass reduction when moving from manifold to manifold. The aim of this paper is to present a methodology to design low-thrust trajectories between two planetary orbits connecting the manifolds of two circular three body systems. The approach is based on a grid search on the main parameters governing the solution to identify those trajectories moving within the manifold images on given Poincarè sections. The value of the Jacoby constant of the target libration point periodic orbit is chosen as stop condition for the thrusting phases. Ballistic arcs follow up to the proper Poincarè section intersection. A grid search for an Earth to Venus transfer is presented as test case.     

Low-energy Earth–Moon transfers involving manifolds through isomorphic mapping

Analysis and design of low-energy transfers to the Moon has been a subject of great interest for decades. Exterior and interior transfers, based on the transit through the regions where the collinear libration points are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is concerned with a geometrical approach for low-energy Earth-to-Moon mission analysis, based on isomorphic mapping. The isomorphic mapping of trajectories allows a visual, intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Two types of Earth-to-Moon missions are considered. The first mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a periodic orbit around the Moon. The second mission is composed of the following arcs: (i) transfer trajectory from a circular low Earth orbit to the stable invariant manifold associated with the Lyapunov orbit at L₁ (corresponding to a specified energy level) and (ii) transfer trajectory along the unstable manifold associated with the Lyapunov orbit at L₁, with final injection in a capture (non-periodic) orbit around the Moon. In both cases three velocity impulses are needed to perform the transfer: the first at an unknown initial point along the low Earth orbit, the second at injection on the stable manifold, the third at injection in the final (periodic or capture) orbit. The final goal is in finding the optimization parameters, which are represented by the locations, directions, and magnitudes of the velocity impulses such that the overall delta-v of the transfer is minimized. This work proves how isomorphic mapping (in two distinct forms) can be profitably employed to optimize such transfers, by determining in a geometrical fashion the desired optimization parameters that minimize the delta-v budget required to perform the transfer.     

On solving the problems of optimization of trajectories of many-revolution orbit transfers of spacecraft

The problem of optimal control over many-revolution spacecraft orbit transfers between circular coplanar orbits of satellites is considered. The spacecraft flight is controlled by a thrust vector of a jet engine with restricted thrust (JERT). The mass expenditure is minimized at a limited time of flight. The optimal control problem is solved based on the maximum principle. The boundary value problem of the maximum principle is solved numerically using the shooting method. A modified computation scheme of the shooting method is suggested (multi-point shooting), as well as a method (correlated with the scheme) of choosing the initial approximation with the use of a solution to the optimization problem in the impulse formulation. The scheme and method allow one to construct many-revolution spacecraft orbit transfers.     

Control of satellite imaging formations in multi-body regimes

Libration point orbits may be ideal locations for satellite imaging formations. Therefore, control of these arrays in multi-body regimes is critical. A continuous feedback control algorithm is developed that maintains a formation of satellites in motion that is bounded relative to a halo orbit. This algorithm is derived based on the dynamic characteristics of the phase space near periodic orbits in the circular restricted three-body problem (CR3BP). By adjusting parameters of the control algorithm appropriately, satellites in the formation follow trajectories that are particularly advantageous to imaging arrays. Image reconstruction and coverage of the (u, v) plane are simulated for interferometric satellite configurations, demonstrating potential applications of the algorithm and the resulting motion.     

Solving Optimization Problems for the Flight Trajectories of a Spacecraft with a High-Thrust Jet Engine in Pulse Formulation for an Arbitrary Gravitational Field in a Vacuum

A mathematically well-posed technique is suggested to obtain first-order necessary conditions of local optimality for the problems of optimization to be solved in a pulse formulation for flight trajectories of a spacecraft with a high-thrust jet engine (HTJE) in an arbitrary gravitational field in vacuum. The technique is based on the Lagrange principle of derestriction for conditional extremum problems in a function space. It allows one to formalize an algorithm of change from the problems of optimization to a boundary-value problem for a system of ordinary differential equations in the case of any optimization problem for which the pulse formulation makes sense. In this work, such a change is made for the case of optimizing the flight trajectories of a spacecraft with a HTJE when terminal and intermediate conditions (like equalities, inequalities, and the terminal functional of minimization) are taken in a general form. As an example of the application of the suggested technique, we consider in this work, within the framework of a bounded circular three-point problem in pulse formulation, the problem of constructing the flight trajectories of a spacecraft with a HTJE through one or several libration points (including the case of going through all libration points) of the Earth–Moon system. The spacecraft is launched from a circular orbit of an Earth's artificial satellite and, upon passing through a point (or points) of libration, returns to the initial orbit. The expenditure of mass (characteristic velocity) is minimized at a restricted time of transfer.     

Electric transfer optimization for mars sample return mission

The aim of this paper is to analyse an alternative scenario for Mars Sample Return Orbiter mission, where electric propulsion is used for Earth-Mars and Mars-Earth heliocentric cruises and for Mars orbit insertion / escape transfers, whereas chemical propulsion is used for final Mars rendezvous. The problem consists in minimizing the initial vehicle mass to obtain a specific final dry mass in reasonable time. The planetocentric phases correspond to continuous low-thrust trajectories, spiraling around Mars between a low orbit and the influence sphere altitude. The heliocentric phases consist of a succession of low-thrust and coasting arcs with specific departure and arrival conditions at the Earth. For these two types of transfer, efficient optimal control tools exist based on Pontryagin's maximum principle. Thanks to the coordination between planetocentric and heliocentric phases, the solution obtained with these two separate tools gives a good upper bound of the optimal solution in terms of propellant consumption and duration. This optimization procedure is described and finally applied to the proposed mission. The numerical results are presented and compared with the baseline chemical mission solution. The electric option could allow to decrease the spacecraft departure mass but may lead to rather long mission duration.     

Aladdin mission concept

Aladdin, one of five Concept Study winners for NASA's Discovery AO98-OSS-04, was a mission to obtain samples from the two Martian moons using several unique mission design and sample collection techniques. The mission design enabled sample return from two bodies at the relatively low cost of a Discovery-class mission. It featured a phasing orbit, multiple flybys of the Martian moons, and a short overall mission duration. The phasing orbit greatly reduced the post-launch Δv requirement, thus permitting the use of a Delta II launch vehicle. Multiple moon flybys provided ample opportunities for sample collection and science observations. The short overall mission duration reduced program costs. Aladdin's sample collection, unlike traditional sample collection methods, used a “launch-and-catch” technique to obtain samples. Projectiles would be launched to the moon's surface during a close flyby and the ejected particles gathered for Earth return and analysis. This innovative technique, the Aladdin mission, and the possible extension of the technique to other bodies are described.     

从地球到日-火系平动点轨道的转移

平动点轨道特殊的空间位置及动力学特征,使其在深空探测中具有重要的应用。以日-火系平动点轨道(Lissajous与Halo轨道)任务为目标,结合平动点轨道的不变流形理论,研究了小推力转移问题。首先给出了圆型限制性三体动力学模型下平动点附近不变流形(稳定和不稳定流形)高阶分析解以及相应的计算实例。接着以流形分析解为基础,建立了初始小推力轨道优化模型,并利用改进的协作进化算法求解初始小推力轨道。最后将初始轨道离散,采用多点打靶法将最优控制问题转化为参数优化问题,并用序列二次规划方法(SQP)求解。仿真结果证明轨道设计方法的有效性。