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.     

Interplanetary CubeSats system for space weather evaluations and technology demonstration

The paper deals with the mission analysis and conceptual design of an interplanetary 6U CubeSats system to be implemented in the L₁ Earth–Sun Lagrangian Point mission for solar observation and in-situ space weather measurements.     

Single-impulse transfers from orbits of artificial satellites to orbits around <Emphasis Type="Italic">L</Emphasis><Subscript>1</Subscript> or <Emphasis Type="Italic">L</Emphasis><Subscript>2</Subscript> libration points

In the context of the restricted circular three-body problem a method for constructing families of periodic orbits is described. Each orbit contains a segment of transfer from artificial satellite orbit of a smaller body to an orbit around L ₁ or L ₂ points of the Sun-Earth and Earth-Moon systems, a segment of multiple flyby of this libration point, and a segment of return to the artificial satellite orbit. Dependences of velocities at the pericenter on the pericenter radius are given.     

Periodic solutions in three-dimensional restricted three-body problem: 2. Loss of symmetry at the 1:1 resonance

Within the framework of the circular restricted three-body problem a family of inverse periodic orbits around the two attracting bodies (the Egorov’s family) and families generated by it at the 1:1, 2:1, and 3:1 resonances for three-dimensional orbits in the Sun-Earth and Earth-Moon systems are considered. Their relationship with families generated by orbits around the libration points L ₁, L ₂ and L ₃ is investigated. One of the families contains periodic solutions that seem promising as possible orbits for the space radio telescope of the Millimetron project.     

Stable spatial orbits around collinear libration points

A technique of generation of spatial periodic solutions to the restricted circular three-body problem from periodic orbits of the planar problem has been used for the families of orbits around collinear libration points L ₁ and L ₂. Developing the families obtained at the 1: 1 resonance, we have obtained stable solutions both in the Earth-Moon system and in the Sun-Earth system. Of course, the term “around the libration point” is rather conventional; the obtained orbits become more similar to the orbits around the smaller attracting body. The further development of the family of orbits “around” the libration point L ₂ in the Sun-Earth system made it possible to find the orbits satisfying the new, much more rigorous constraints on cooling the spacecraft of the Millimetron project.     

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.     

Lagrangian coherent structures in various map representations for application to multi-body gravitational regimes

The Finite-Time Lyapunov Exponent (FTLE) has been demonstrated as an effective metric for revealing distinct, bounded regions within a flow. The dynamical differential equations derived in multi-body gravitational environments model a flow that governs the motion of a spacecraft. Specific features emerge in an FTLE map, denoted Lagrangian Coherent Structures (LCS), that define the extent of regions that bound qualitatively different types of behavior. Consequently, LCS supply effective barriers to transport in a generic system, similar to the notion of invariant manifolds in autonomous systems. Unlike traditional invariant manifolds associated with solutions in an autonomous system, LCS evolve with the flow in time-dependent systems while continuing to bound distinct regions of behavior. Moreover, in general, FTLE values supply information describing the relative sensitivity in the neighborhood of a trajectory. Here, different models and variable representations are used to generate maps of FTLE, and the resulting structures are applied to design and analysis within an astrodynamical context. Application of FTLE and LCS to transfers from LEO to the L₁ region in the Earth–Moon system are presented and discussed. In an additional example, an FTLE analysis is offered of a few stationkeeping maneuvers from the Earth–Moon mission ARTEMIS (Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun).     

Planetary Defense From Space: Part 2 (Simple) Asteroid Deflection Law

A system of two space bases housing missiles for an efficient Planetary Defense of the Earth from asteroids and comets was firstly proposed by this author in 2002. It was then shown that the five Lagrangian points of the Earth–Moon system lead naturally to only two unmistakable locations of these two space bases within the sphere of influence of the Earth. These locations are the two Lagrangian points L1 (in between the Earth and the Moon) and L3 (in the direction opposite to the Moon from the Earth). In fact, placing missiles based at L1 and L3 would enable the missiles to deflect the trajectory of incoming asteroids by hitting them orthogonally to their impact trajectory toward the Earth, thus maximizing the deflection at best. It was also shown that confocal conics are the only class of missile trajectories fulfilling this “best orthogonal deflection” requirement.The mathematical theory developed by the author in the years 2002–2004 was just the beginning of a more expanded research program about the Planetary Defense. In fact, while those papers developed the formal Keplerian theory of the Optimal Planetary Defense achievable from the Earth–Moon Lagrangian points L1 and L3, this paper is devoted to the proof of a simple “(small) asteroid deflection law” relating directly the following variables to each other:

(1) the speed of the arriving asteroid with respect to the Earth (known from the astrometric observations);

(2) the asteroid's size and density (also supposed to be known from astronomical observations of various types);

(3) the “security radius” of the Earth, that is, the minimal sphere around the Earth outside which we must force the asteroid to fly if we want to be safe on Earth. Typically, we assume the security radius to equal about 10,000 km from the Earth center, but this number might be changed by more refined analyses, especially in the case of “rubble pile” asteroids;

(4) the distance from the Earth of the two Lagrangian points L1 and L3 where the defense missiles are to be housed;

(5) the deflecting missile's data, namely its mass and especially its “extra-boost”, that is, the extra-energy by which the missile must hit the asteroid to achieve the requested minimal deflection outside the security radius around the Earth.

This discovery of the simple “asteroid deflection law” presented in this paper was possible because:

(1) In the vicinity of the Earth, the hyperbola of the arriving asteroid is nearly the same as its own asymptote, namely, the asteroid's hyperbola is very much like a straight line. We call this approximation the line/circle approximation. Although “rough” compared to the ordinary Keplerian theory, this approximation simplifies the mathematical problem to such an extent that two simple, final equations can be derived.

(2) The confocal missile trajectory, orthogonal to this straight line, ceases then to be an ellipse to become just a circle centered at the Earth. This fact also simplifies things greatly. Our results are thus to be regarded as a good engineering approximation, valid for a preliminary astronautical design of the missiles and bases at L1 and L3.

Still, many more sophisticated refinements would be needed for a complete Planetary Defense System:

(1) taking into account many perturbation forces of all kinds acting on both the asteroids and missiles shot from L1 and L3;

(2) adding more (non-optimal) trajectories of missiles shot from either the Lagrangian points L4 and L5 of the Earth–Moon system or from the surface of the Moon itself;

(3) encompassing the full range of missiles currently available to the USA (and possibly other countries) so as to really see “which missiles could divert which asteroids”, even just within the very simplified scheme proposed in this paper.

In summary: outlined for the first time in February 2002, our Confocal Planetary Defense concept is a simplified Keplerian Theory that already proved simple enough to catch the attention of scholars, popular writers, and representatives of the US Military. These developments would hopefully mark the beginning of a general mathematical vision for building an efficient Planetary Defense System in space and in the vicinity of the Earth, although not on the surface of the Earth itself!We must make a real progress beyond academic papers, Hollywood movies and secret military plans, before asteroids like 99942 Apophis get close enough to destroy us in 2029 or a little later.     

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.     

Gravitation Maneuver Using the Families of Super-Unstable Orbits around the Libration Points

The families of periodic solutions to an autonomous Hamiltonian system in that part where the solutions are unstable have their specific field of influence. Under strong instability, the orbits that have fallen in such a field of influence are drawn into the family as in a whirlpool and then are thrown away from it. In the particular case of the restricted three-body problem, the orbits around the libration points L ₁ and L ₂ are super-unstable and the Keplerian elements in motion over these orbits change drastically. When the orbits fall into such a domain, for some time the motion is close to the motion around the libration point, and after being thrown out of this domain, the Keplerian orbital elements also change essentially.     

Regularization of the Perturbed Spatial Restricted Three-Body Problem by <Emphasis Type="Italic">L</Emphasis>-Transformations

Equations of motion for the perturbed circular restricted three-body problem have been regularized in canonical variables in a moving coordinate system. Two different L-matrices of the fourth order are used in the regularization. Conditions for generalized symplecticity of the constructed transform have been checked. In the unperturbed case, the regular equations have a polynomial structure. The regular equations have been numerically integrated using the Runge–Kutta–Fehlberg method. The results of numerical experiments are given for the Earth–Moon system parameters taking into account the perturbation of the Sun for different L-matrices.     

A Retrospective Geometric Analysis of the Long-Term Evolution of Orbits and the Ballistic Lifetime for Satellites of the Prognoz Series

We consider the applied aspects of the geometrical analysis of solutions in a restricted circular double-averaged three-body problem that are concerned with the design of high-apogee satellite orbits. Based on the analysis of the long-term evolution and the ballistic lifetime for orbit families of the Prognoz satellites launched into their orbits in the period 1972–1995, we suggest some practical suggestions for choosing long-lived high-apogee orbits with taking into account various requirements for the domain of evolution of the orbital elements.     

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.     

Artificial equilibrium points for a solar balloon in the α Centauri system

This paper discusses the generation, stability, and control of artificial equilibrium points for a solar balloon spacecraft in the α Centauri A and B binary star system. The continuous propulsive acceleration provided by a solar balloon is shown to be able to modify the position of the (classical) Lagrangian equilibrium points of the three-body system on a locus whose geometrical form is known analytically. A linear stability analysis reveals that the new generated equilibrium points are usually unstable, but part of them can be stabilized with a simple feedback control logic.     

Horseshoe orbits in the Earth–Moon system

Horseshoe orbits in the restricted three-body problem have been mostly considered in the Sun–Jupiter system and, in recent years, in the Sun–Earth system. Here, these orbits have been used to find asteroids that have orbits of this kind. We have built a planar family of horseshoe orbits in the Earth–Moon system and determined the points of planar and 1/1 vertical resonances on this family. We have presented examples of orbits generated by these spatial families.     

Periodic orbits in the restricted four-body problem

One solution to the restricted three-body problem is that of three-dimensional, periodic halo orbits. These orbits emanate from the three collinear libration points and exist at all possible mass ratios of the primaries, with infinite possible trajectories. Their existence motivated this study of the effect of an additional gravitational influence on the motion. The approach is to seek a solution in the restricted four-body problem using halo solutions as an initial approximation. The method first solves for periodic, coplanar motion of the three primaries under their mutual gravitational attractions and represents them as trigonometric series. Then, under the modified gravity force model, a three-dimensional solution is obtained in the problem of four bodies, periodic with respect to the synodic system of the three primaries. The additional primary remains relatively far removed and acts as a perturbing influence on the original motion. Some shape and stability characteristics are presented for three such solutions.     

An Empirical Model of the Radiation Belt of Helium Nuclei

Based on satellite data, we present the results of modeling the spatial and energy distributions of integral fluxes of He nuclei (α particles) with E > 1, 2, 4, and 7 MeV at L = 1.1–6.6 in a broad range of B/B ₀ (E is the kinetic energy of particles, L is the drift shell parameter, and B/B ₀ is the magnetic field ratio). Some ways of practically applying the model are considered. The results of calculation of α-particle fluxes for a circular orbit with a height of 300 km and an inclination of 50° are presented.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 243–247.Original Russian Text Copyright © 2005 by Getselev, Sosnovets, Kovtyukh, Dmitriev, Podzolko, Vlasova, Reizman.     

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.     

Two-craft tether formation relative equilibria about circular orbits and libration points

The relative equilibria of a two spacecraft tether formation connected by line-of-sight elastic forces moving in the context of a restricted two-body system and a circularly restricted three-body system are investigated. For a two spacecraft formation moving in a central gravitational field, a common assumption is that the center of the circular orbit is located at the primary mass and the center of mass of the formation orbits around the primary in a great-circle orbit. The relative equilibrium is called great-circle if the center of mass of the formation moves on the plane with the center of the gravitational field residing on it; otherwise, it is called a nongreat-circle orbit. Previous research shows that nongreat-circle equilibria in low Earth orbits exhibit a deflection of about a degree from the great-circle equilibria when spacecraft with unequal masses are separated by 350 km. This paper studies these equilibria (radial, along-track and orbit-normal in circular Earth orbit and Earth–Moon Libration points) for a range of inter-craft distances and semi-major axes of the formation center of mass. In the context of a two-spacecraft Coulomb formation with separation distances on the order of dozens of meters, this paper shows that the equilibria deflections are negligible (less than 10^?6°) even for very heterogeneous mass distributions. Furthermore, the nongreat-circle equilibria conditions for a two spacecraft tether structure at the Lagrangian libration points are developed.