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.     

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.     

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.     

Optimal aeroassisted return from high earth orbit with plane change

This paper gives a complete analysis of the problem of aeroassisted return from a high Earth orbit to a low Earth orbit with plane change. A discussion of pure propulsive maneuver leads to the necessary change for improvement of the fuel consumption by inserting in the middle of the trajectory an atmospheric phase to obtain all or part of the required plane change. The variational problem is reduced to a parametric optimization problem by using the known results in optimal impulsive transfer and solving the atmospheric turning problem for storage and use in the optimization process. The coupling effect between space maneuver and atmospheric maneuver is discussed. Depending on the values of the plane change i, the ratios of the radii, $\text{n =}\text{r}{}_1\text{r}{}_2$ between the orbits and $\text{a =}\text{r}{}_2\text{R}$ between the low orbit and the atmosphere, and the maximum lift-to-drag ratio E1 of the vehicle, the optimal maneuver can be pure propulsive or aeroassisted. For aeroassisted maneuver, the optimal mode can be parabolic, which requires only drag capability of the vehicle, or elliptic. In the elliptic mode, it can be by one-impulse for deorbit and one or two-impulse in postatmospheric flight, or by two-impulse for deorbit with only one impulse for final circularization. It is shown that whenever an impulse is applied, a plane change is made. The necessary conditions for the optimal split of the plane changes are derived and mechanized in a program routine for obtaining the solution.     

Definition and analysis of properties of optimal three-impulse point-to-orbit transfers under time constraint

This paper deals with energetically optimal multi-impulse transfers of a spacecraft in the central Newtonian gravitational field near a planet. The transfer from a point on initial orbit to the final orbit with the given angular momentum and energy constants is considered. The transfer time is bounded above.With the distance from spacecraft to planet limited and the time free, such parameters of given orbits are chosen that the 3-impulse apsidal transfer T_r is optimal with an intermediate impulse at the maximum distance. On the basis of necessary optimality conditions an algorithm is developed to numerically determine the desired optimal transfer trajectory T_t under time constraint, the apsidal trajectory T_r being taken as initial approach. From the geometry and energy viewpoints, both trajectories T_t and T_r are close to each other. The trajectory T_t is also 3-impulsive, all impulses on it are nonapsidal. The distance from the planet is larger and the sum of impulses is less for this trajectory than for the initial trajectory T_r with the same transfer time.The simplified solution of the problem is constructed producing good approximation to the exact numerical optimization results. The solution asymptotics is found when the transfer time tends to infinity.     

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.     

On Saturn’s rotation relative to a center of mass under the action of the gravitational moments of the Sun and Jupiter

Saturn’s rotation relative to a center of mass is considered within an elliptic restricted three-body problem. It is assumed that Saturn is a solid under the action of gravity of the Sun and Jupiter. The motions of Saturn and Jupiter are considered elliptic with small eccentricities e_S and e_J, respectively; the mean motion of Jupiter n_J is also small. We obtain the averaged Hamiltonian function for a small parameter of ε = n_J and integrals of evolution equations. The main effects of the influence of Jupiter on Saturn’s rotation are described: (α) the evolution of the constant parameters of regular precession for the angular momentum vector I₂; (β) the occurrence of new libration zones of oscillations I₂ near the plane of the celestial equator parallel to the plane of the Jupiter’s orbit; (γ) the occurrence of additional unstable equilibria of vector I₂ at the points of the north and south poles of the celestial sphere and, as a result, the existence of homoclinic trajectories; and (δ) the existence of periodic trajectories with arbitrarily large periods near the homoclinic trajectory. It is shown that the effects of (β), (γ), and (δ) are caused by the eccentricity e of the Jupiter’s orbit and are practically independent of Jupiter’s mass (within satellite approximation).     

Large-size space debris flyby in low earth orbits

the analysis of NORAD catalogue of space objects executed with respect to the overall sizes of upper-stages and last stages of carrier rockets allows the classification of 5 groups of large-size space debris (LSSD). These groups are defined according to the proximity of orbital inclinations of the involved objects. The orbits within a group have various values of deviations in the Right Ascension of the Ascending Node (RAAN). It is proposed to use the RAANs deviations' evolution portrait to clarify the orbital planes’ relative spatial distribution in a group so that the RAAN deviations should be calculated with respect to the concrete precessing orbital plane of the concrete object. In case of the first three groups (inclinations i = 71°, i = 74°, i = 81°) the straight lines of the RAAN relative deviations almost do not intersect each other. So the simple, successive flyby of group’s elements is effective, but the significant value of total ΔV is required to form drift orbits. In case of the fifth group (Sun-synchronous orbits) these straight lines chaotically intersect each other for many times due to the noticeable differences in values of semi-major axes and orbital inclinations. The intersections’ existence makes it possible to create such a flyby sequence for LSSD group when the orbit of one LSSD object simultaneously serves as the drift orbit to attain another LSSD object. This flyby scheme requiring less ΔV was called “diagonal.” The RAANs deviations’ evolution portrait built for the fourth group (to be studied in the paper) contains both types of lines, so the simultaneous combination of diagonal and successive flyby schemes is possible. The value of total ΔV and temporal costs were calculated to cover all the elements of the 4th group. The article is also enriched by the results obtained for the flyby problem solution in case of all the five mentioned LSSD groups. The general recommendations are given concerned with the required reserve of total ΔV and with amount of detachable de-orbiting units onboard the maneuvering platform and onboard the refueling vehicle.     

Results of the early orbit manoeuvres of the ampte-uks

The Active Magnetospheric Particle Tracer Explorers (AMPTE) program consists of three satellites which were launched on 16th August 1984. The scientific aim of the mission is to inject lithium and barium tracer ions inside and outside the Earth's magnetosphere and to detect and monitor these ions as they diffuse through the inner magnetosphere. The first of these three satellites, the U.S. Charge Composition Explorer (CCE) was launched into an elliptical orbit of apogee 8 R_e. The other two satellites are the West German Ion Release Module (IRM) and the U.K. Subsatellite (UKS), both of which were launched on the same vehicle into a highly elliptical orbit of apogee 18 R_e. At discreet intervals during the mission the IRM will release ions into the solar wind, and the movement of these ions will be monitored by the UKS. Depending on the particular scientific requirement, the UKS has to be positioned accurately at a given distance behind the IRM. Initially the UKS has to be located 100 km behind the IRM, and held there for ～9 months. It will then be moved a distance of ～1 R_e behind the IRM. In order to manoeuvre the UKS around its orbit, a cold gas jet system is incorporated on the satellite, allowing impulses to be applied both along and perpendicular to the orbit velocity vector. The orbit control system also has to cater for relative orbit changes due to air drag at perigee, as the IRM and the UKS have different $\text{area}\text{mass}$ ratios. This paper presents an account of the orbit control system implemented on the UKS, together with the mathematical approach adopted, and results from manoeuvres made in the first weeks of the mission.     

A small parameter method in problems of maneuvering space vehicles with aerodynamic efficiency

The motion of a satellite with aerodynamic efficiency along a low near-circular orbit is considered in the paper. The controls of bank angle γ and lift coefficient C_y are used as control functions. The introduction of a small parameter ($\text{? = (}\text{?}{}_0\text{· S · g}{}_0\text{2G}\text{)}$) makes it possible to integrate an adjoint system of equations and to obtain an approximate solution to the complete problem in the class of piecewise-constant control functions. Maximum values for the coordinates of heading angle η and lateral derivation from the plane of a reference orbit ?, which are connected with orbit plane angle by the relation cos i = cos ? · cos h, are used as criteria of maneuvering capability for a satellite with aerodynamic efficiency. Optimal programs for bank angle and incidence variation are derived and the influence of lift-to-drag ratio on the vehicle maneuvering capabilities has been estimated.It is shown that the process of the optimal motion is a special kind of gravitational skipping similar to the Keplerian motion but with continuous descent.     

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 practical approach to the disposal of highly toxic and long-lived spent nuclear fuel waste between Venus and Earth

Extraterrestrial disposal, while not the only alternative (there is at least one very safe terrestrial method), nevertheless assures definite and irreversible removal of the most toxic and long-lived waste from the biosphere. In the foreseeable future, there is far less need to retrieve nuclear waste for later use then to dispose of it in a publicly acceptable manner, beginning in the near future (middle 1980s). It is, therefore, more important to assure safety in the weakest link of this disposal process—transportation into near-Earth orbit—than to engage in a retrievable disposal beyond Earth. The disposal “site” should lie at minimum safe transfer energy level. Primary candidate is the space between Venus and Earth. The number of propulsion phases should be a minimum, preferably only one (out of near-Earth orbit). Lunar gravity assist can be helpful to achieve higher inclination of the heliocentric orbit relative to the ecliptic.Solidified spent fuel isotopes and actinides, sufficient to reduce the residual terrestrial waste to the radiation level of natural uranium deposits after 30–40 yr instead of 1000–1500 yr, is deposited into heliocentric orbits. Transportation systems, requirements, costs and the associated socio-economic benefit potentials of an environmentally more benign and a more vigorous nuclear power generation program are presented.Prior to solidification, an interim storage of 10 yr, following removal from the reactor, may be required. The Shuttle, with one Orbiter modified as Nuclear Waste Carrying Orbiter (NWCO) and an out of near-Earth orbit booster, provides a safe and economic transportation system at (1979$) disposal mission costs from surface to disposal orbit of less than $\text{0.5¢/}\text{kWhe or}\text{? 0.1¢/}\text{kWhe}$ (some 70% of overall exo-disposal cost), depending on level of orbital operations (this at basic Shuttle flight cost of $30M). The orbital operations mode can be selected according to lead time and costs, and can be advanced sequentially, lowering disposal costs while at the same time financing the development of orbital operations techniques needed also for other and larger-scale exo-industrial activities. An average of 10–15 disposal missions of the NWCO is required annually, at the given conditions, to service the spent fuel of 173 reactors 1978 in operation in W. Europe, Japan and North America.     

Families of periodic solutions in three-dimensional restricted three-body problem

Planar orbits of three-dimensional restricted circular three-body problem are considered as a special case of three-dimensional orbits, and the second-order monodromy matrices M (in coordinate z and velocity v _z ) are calculated for them. Semi-trace s of matrix M determines vertical stability of an orbit. If |s| ≤ 1, then transformation of the subspace (z, v _z ) in the neighborhood of solution for the period is reduced to deformation and a rotation through angle φ, cosφ = s. If the angle ? can be rationally expressed through 2π,φ = 2π·p/q, where p and q are integer, then a planar orbit generates the families of three-dimensional periodic solutions that have a period larger by a factor of q (second kind Poincareé periodic solutions). Directions of continuation in the subspace (z, v _z ) are determined by matrix M. If |s| < 1, we have two new families, while only one exists at resonances 1: 1 (s = 1) and 2: 1 (s = ?1). In the course of motion along the family of three-dimensional periodic solutions, a transition is possible from one family of planar solutions to another one, sometimes previously unknown family of planar solutions.     

Applications of multi-body dynamical environments: The ARTEMIS transfer trajectory design

The application of forces in multi-body dynamical environments to permit the transfer of spacecraft from Earth orbit to Sun–Earth weak stability regions and then return to the Earth–Moon libration (L₁ and L₂) orbits has been successfully accomplished for the first time. This demonstrated that transfer is a positive step in the realization of a design process that can be used to transfer spacecraft with minimal Delta-V expenditures. Initialized using gravity assists to overcome fuel constraints; the ARTEMIS trajectory design has successfully placed two spacecrafts into Earth–Moon libration orbits by means of these applications.     

On orbit control using gravity gradient effects

Possibility of orbit control using gravity gradient (GG) effects without any mass expulsion is discussed. For simplicity, a dumb-bell type satellite and circular orbits are mainly considered. It is shown that the GG effects can be applied to convert attitude torques into orbital torques and vice versa. In central gravitational force fields, maximum orbital torques or thrusts are available from the GG force when roll or pitch angle is $\text{±}\text{π}\text{4}$ provided that the attitude angle is null when the dumb-bell axis coincides with the local vertical. Such external torques as geomagnetic or solar wind pressure can be utilized to maintain the $\text{±}\text{π}\text{4}$ attitude, then the orbital torques are available forever. In non-central gravitational fields, without any external torque, the orbital radii of circular orbits can be increased by controlling the satellite attitude using electric energy. The use of the Earth's oblateness effects and the exterior Lunar potential is discussed.     

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.     

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.     

Tides and the evolution of planetary habitability

Tides raised on a planet by the gravity of its host star can reduce the planet's orbital semi-major axis and eccentricity. This effect is only relevant for planets orbiting very close to their host stars. The habitable zones of low-mass stars are also close in, and tides can alter the orbits of planets in these locations. We calculate the tidal evolution of hypothetical terrestrial planets around low-mass stars and show that tides can evolve planets past the inner edge of the habitable zone, sometimes in less than 1 billion years. This migration requires large eccentricities (>0.5) and low-mass stars ( less or similar to 0.35 M(circle)). Such migration may have important implications for the evolution of the atmosphere, internal heating, and the Gaia hypothesis. Similarly, a planet that is detected interior to the habitable zone could have been habitable in the past. We consider the past habitability of the recently discovered, approximately 5 M(circle) planet, Gliese 581 c. We find that it could have been habitable for reasonable choices of orbital and physical properties as recently as 2 Gyr ago. However, when constraints derived from the additional companions are included, most parameter choices that indicate past habitability require the two inner planets of the system to have crossed their mutual 3:1 mean motion resonance. As this crossing would likely have resulted in resonance capture, which is not observed, we conclude that Gl 581 c was probably never habitable.     

Order of magnitude analysis of convective effects in dendritic solidification

The different types of convective phenomena which may occur during the dendritic solidification of metallic alloys are discussed from an order of magnitude analysis. Bulk thermal convection and/or interdendritic solutal convection have to be considered according to the values of the experimental data. Scaling laws for the solute boundary layer resulting from bulk thermal convection have already been derived. It is shown here that the interdendritic flow depends on a solutal Grashof number Gr based on the horizontal density gradient and a characteristic length L_s which is of the order of the liquid channels width. For Gr < 1, which is generally verified in practical cases, the interdendritic flow velocity U_r is proportional to the Grashof number. This a priori law compares favorably with the results of horizontal solidification experiments where the mean interdendritic flow velocity has been estimated from the resulting measured macrosegregation. In these experiments, as well as for most horizontal dendritic solidifications of metallic alloys at 1 g, the ratio $\text{U}{}_{\mathrm{r}}\text{R}$ (R is the growth rate) is of order one. In order to cancel the interdendritic flow effects, this ratio has to be lowered by one order of magnitude. According to our analysis, this can be obtained by performing the experiments either at a slightly reduced g level (～10^?1 g), or at 1 g in a vertical stable configuration with a sufficiently low residual horizontal thermal gradient.