Escape dynamics and fractal basin boundaries in the planar Earth–Moon system

The escape of trajectories of a spacecraft, or comet or asteroid in the presence of the Earth–Moon system is investigated in detail in the context of the planar circular restricted three-body problem, in a scattering region around the Moon. The escape through the necks around the collinear points \(L_1\) and \(L_2\) as well as the leaking produced by considering collisions with the Moon surface, taking the lunar mean radius into account, were considered. Given that different transport channels are available as a function of the Jacobi constant, four distinct escape regimes are analyzed. Besides the calculation of exit basins and of the spatial distribution of escape time, the qualitative dynamical investigation through Poincaré sections is performed in order to elucidate the escape process. Our analyses reveal the dependence of the properties of the considered escape basins with the energy, with a remarkable presence of fractal basin boundaries along all the escape regimes. Finally, we observe the plentiful presence of stickiness motion near stability islands which plays a remarkable role in the longest escape time behavior. The application of this analysis is important both in space mission design and study of natural systems, given that fractal boundaries are related with high sensitivity to initial conditions, implying in uncertainty between safe and unsafe solutions, as well as between escaping solutions that evolve to different phase space regions.     

Study of the transfer between libration point orbits and lunar orbits in Earth–Moon system

This paper is devoted to the study of the transfer problem from a libration point orbit of the Earth–Moon system to an orbit around the Moon. The transfer procedure analysed has two legs: the first one is an orbit of the unstable manifold of the libration orbit and the second one is a transfer orbit between a certain point on the manifold and the final lunar orbit. There are only two manoeuvres involved in the method and they are applied at the beginning and at the end of the second leg. Although the numerical results given in this paper correspond to transfers between halo orbits around the \(L_1\) point (of several amplitudes) and lunar polar orbits with altitudes varying between 100 and 500 km, the procedure we develop can be applied to any kind of lunar orbits, libration orbits around the \(L_1\) or \(L_2\) points of the Earth–Moon system, or to other similar cases with different values of the mass ratio.     

Effect of Moon perturbation on the energy curves and equilibrium points in the Sun–Earth–Moon system

In this paper, we have considered that the Moon motion around the Earth is a source of a perturbation for the infinitesimal body motion in the Sun–Earth system. The perturbation effect is analyzed by using the Sun–Earth–Moon bi–circular model (BCM). We have determined the effect of this perturbation on the Lagrangian points and zero velocity curves. We have obtained the motion of infinitesimal body in the neighborhood of the equivalent equilibria of the triangular equilibrium points. Moreover, to know the nature of the trajectory, we have estimated the first order Lyapunov characteristic exponents of the trajectory emanating from the vicinity of the triangular equilibrium point in the proposed system. It is noticed that due to the generated perturbation by the Moon motion, the results are affected significantly, and the Jacobian constant is fluctuated periodically as the Moon is moving around the Earth. Finally, we emphasize that this model could be applicable to send either satellite or telescope for deep space exploration.     

A note on the dynamics around the Lagrange collinear points of the Earth–Moon system in a complete Solar System model

In this paper we study the dynamics of a massless particle around the L _1,2 libration points of the Earth–Moon system in a full Solar System gravitational model. The study is based on the analysis of the quasi-periodic solutions around the two collinear equilibrium points. For the analysis and computation of the quasi-periodic orbits, a new iterative algorithm is introduced which is a combination of a multiple shooting method with a refined Fourier analysis of the orbits computed with the multiple shooting. Using as initial seeds for the algorithm the libration point orbits of Circular Restricted Three Body Problem, determined by Lindstedt-Poincaré methods, the procedure is able to refine them in the Solar System force-field model for large time-spans, that cover most of the relevant Sun–Earth–Moon periods.     

Direct and indirect capture of near-Earth asteroids in the Earth–Moon system

Near-Earth asteroids have attracted attention for both scientific and commercial mission applications. Due to the fact that the Earth–Moon \(\hbox {L}_{1}\) and \(\hbox {L}_{2}\) points are candidates for gateway stations for lunar exploration, and an ideal location for space science, capturing asteroids and inserting them into periodic orbits around these points is of significant interest for the future. In this paper, we define a new type of lunar asteroid capture, termed direct capture. In this capture strategy, the candidate asteroid leaves its heliocentric orbit after an initial impulse, with its dynamics modeled using the Sun–Earth–Moon restricted four-body problem until its insertion, with a second impulse, onto the \(\hbox {L}_{2}\) stable manifold in the Earth–Moon circular restricted three-body problem. A Lambert arc in the Sun-asteroid two-body problem is used as an initial guess and a differential corrector used to generate the transfer trajectory from the asteroid’s initial obit to the stable manifold associated with Earth–Moon \(\hbox {L}_{2}\) point. Results show that the direct asteroid capture strategy needs a shorter flight time compared to an indirect asteroid capture, which couples capture in the Sun–Earth circular restricted three-body problem and subsequent transfer to the Earth–Moon circular restricted three-body problem. Finally, the direct and indirect asteroid capture strategies are also applied to consider capture of asteroids at the triangular libration points in the Earth–Moon system.     

Comments on ‘the origin of the Earth—Moon system’

The main points are presented of a new hypothesis of the origin of the Earth—Moon system, developed on the basis of Savi's (1961) theory of the origin of rotation of celestial bodies. The cooling off and contraction due to gravitational attraction on vast particle systems, with the pushing out of electrons from atom shells result in a continually increasing density. Depending on the amount of mass, this pushing out can lead to the expulsion of electrons and the creation of a magnetic field by which a rotational motion is brought about. These conditions are satisfied for the Earth's mass and all larger masses. If the Earth and the Moon formed a unique body, the protoplanet, then once rotational motion had begun, the primeval spherical body must have taken the shape of a large Jacobi ellipsoid. New condensation followed, however no longer solely around the centre of the protoplanet, but also along the edge of the ellipsoid, the process leading to the creation of the dual Earth—Moon system.     

On quasi-periodic motions around the collinear libration points in the real Earth–Moon system

Due to various perturbations, the collinear libration points of the real Earth–Moon system are not equilibrium points anymore. Under the assumption that the Moon’s motion is quasi-periodic, special quasi-periodic orbits called dynamical substitutes exist. These dynamical substitutes replace the geometrical collinear libration points as time-varying equilibrium points. In the paper, the dynamical substitutes of the three collinear libration points in the real Earth–Moon system are computed. For the points L ₁ and L ₂, linearized motions around the dynamical substitutes are described, and the variational equations of the dynamical substitutes are reduced to a form with a near constant coefficient matrix. Then higher order analytical formulae of the central manifolds are constructed. Using these analytical solutions as initial seeds, Lissajous orbits and halo orbits are computed with numerical algorithms.     

Earth–Moon Weak Stability Boundaries in the restricted three and four body problem

This work deals with the structure of the lunar Weak Stability Boundaries (WSB) in the framework of the restricted three and four body problem. Geometry and properties of the escape trajectories have been studied by changing the spacecraft orbital parameters around the Moon. Results obtained using the algorithm definition of the WSB have been compared with an analytical approximation based on the value of the Jacobi constant. Planar and three-dimensional cases have been studied in both three and four body models and the effects on the WSB structure, due to the presence of the gravitational force of the Sun and the Moon orbital eccentricity, have been investigated. The study of the dynamical evolution of the spacecraft after lunar capture allowed us to find regions of the WSB corresponding to stable and safe orbits, that is orbits that will not impact onto lunar surface after capture. By using a bicircular four body model, then, it has been possible to study low-energy transfer trajectories and results are given in terms of eccentricity, pericenter altitude and inclination of the capture orbit. Equatorial and polar capture orbits have been compared and differences in terms of energy between these two kinds of orbits are shown. Finally, the knowledge of the WSB geometry permitted us to modify the design of the low-energy capture trajectories in order to reach stable capture, which allows orbit circularization using low-thrust propulsion systems.     

Secondary craters and ejecta across the solar system: Populations and effects on impact-crater–based chronologies

We review the secondary-crater research over the past decade, and provide new analyses and simulations that are the first to model an accumulation of a combined primary-plus-secondary crater population as discrete cratering events. We develop the secondary populations by using scaling laws to generate ejecta fragments, integrating the trajectories of individual ejecta fragments, noting the location and velocity at impact, and using scaling laws to estimate secondary-crater diameters given the impact conditions. We also explore the relationship between the impactor size–frequency distribution (SFD) and the resulting secondary-crater SFD. Our results from these analyses indicate that the “secondary effect” varies from surface to surface and that no single conclusion applies across the solar system nor at any given moment in time—rather, there is a spectrum of outcomes both spatially and temporally, dependent upon target parameters and the impacting population. Surface gravity and escape speed define the spatial distribution of secondaries. A shallow-sloped impactor SFD will cause proportionally more secondaries than a steeper-sloped SFD. Accounting for the driving factors that define the magnitude and spatial distribution of secondaries is essential to determine the relative population of secondary craters, and their effect on derived surface ages.     

Extension of fast periodic transfer orbits from the Earth–Moon RTBP to the Sun–Earth–Moon Quasi-Bicircular Problem

Starting from 80 families of low-energy fast periodic transfer orbits in the Earth–Moon planar circular Restricted Three Body Problem (RTBP), we obtain by analytical continuation 11 periodic orbits and 25 periodic arcs with similar properties in the Sun–Earth–Moon Quasi-Bicircular Problem (QBCP). A novel and very simple procedure is introduced giving the solar phases at which to attempt continuation. Detailed numerical results for each periodic orbit and arc found are given, including their stability parameters and minimal distances to the Earth and Moon. The periods of these orbits are between 2.5 and 5 synodic months, their energies are among the lowest possible to achieve an Earth–Moon transfer, and they show a diversity of circumlunar trajectories, making them good candidates for missions requiring repeated passages around the Earth and the Moon with close approaches to the last.     

On quasi-periodic motions around the triangular libration points of the real Earth–Moon system

The two triangular libration points of the real Earth–Moon system are not equilibrium points anymore. Under the assumption that the motion of the Moon is quasi-periodic, one special quasi-periodic orbit exists as dynamical substitute for each point. The way to compute the dynamical substitute was discussed before, and a planar approximation was obtained. In this paper, the problem is revisited. The three-dimensional approximation of the dynamical substitute is obtained in a different way. The linearized central flow around it is described.     

Periodic orbits of solar sail equipped with reflectance control device in Earth–Moon system

In this paper, families of Lyapunov and halo orbits are presented with a solar sail equipped with a reflectance control device in the Earth–Moon system. System dynamical model is established considering solar sail acceleration, and four solar sail steering laws and two initial Sun-sail configurations are introduced. The initial natural periodic orbits with suitable periods are firstly identified. Subsequently, families of solar sail Lyapunov and halo orbits around the \(L_{1}\) and \(L_{2}\) points are designed with fixed solar sail characteristic acceleration and varying reflectivity rate and pitching angle by the combination of the modified differential correction method and continuation approach. The linear stabilities of solar sail periodic orbits are investigated, and a nonlinear sliding model controller is designed for station keeping. In addition, orbit transfer between the same family of solar sail orbits is investigated preliminarily to showcase reflectance control device solar sail maneuver capability.     

Cylindrical isomorphic mapping applied to invariant manifold dynamics for Earth–Moon Missions

Several families of periodic orbits exist in the context of the circular restricted three-body problem. This work studies orbital motion of a spacecraft among these periodic orbits in the Earth–Moon system, using the planar circular restricted three-body problem model. A new cylindrical representation of the spacecraft phase space (i.e., position and velocity) is described, and allows representing periodic orbits and the related invariant manifolds. In the proximity of the libration points, the manifolds form a four-fold surface, if the cylindrical coordinates are employed. Orbits departing from the Earth and transiting toward the Moon correspond to the trajectories located inside this four-fold surface. The isomorphic mapping under consideration is also useful for describing the topology of the invariant manifolds, which exhibit a complex geometrical stretch-and-folding behavior as the associated trajectories reach increasing distances from the libration orbit. Moreover, the cylindrical representation reveals extremely useful for detecting periodic orbits around the primaries and the libration points, as well as the possible existence of heteroclinic connections. These are asymptotic trajectories that are ideally traveled at zero-propellant cost. This circumstance implies the possibility of performing concretely a variety of complex Earth–Moon missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining a suitable, convenient end-of-life strategy for spacecraft placed in any of the unstable orbits. The final disposal orbit is an externally confined trajectory, never approaching the Earth or the Moon, and can be entered by means of a single velocity impulse (of modest magnitude) along the right unstable manifold that emanates from the Lyapunov orbit at \(L_2\) .     

Minicomets in the Solar system and in the Earth’s atmosphere and related phenomena

We consider the history of discovery and justify the existence in the Solar system of a new class of bodies—minicomets, i.e., bodies of cometary nature and composition but of low mass. Two classes of minicomets are distinguished: icy ones similar to the Tunguska meteorite, and snow ones, which break up at high altitudes.     

Small impact craters in the lunar regolith — Their morphologies,relative ages,and rates of formation

Apparently, there are two types of size-frequency distributions of small lunar craters (1–100 m across): (1) crater production distributions for which the cumulative frequency of craters is an inverse function of diameter to power near 2.8, and (2) steady-state distributions for which the cumulative frequency of craters is inversely proportional to the square of their diameters. According to theory, cumulative frequencies of craters in each morphologic category within the steady-state should also be an inverse function of the square of their diameters. Some data on frequency distribution of craters by morphologic types are approximately consistent with theory, whereas other data are inconsistent with theory.A flux of crater producing objects can be inferred from size-frequency distributions of small craters on the flanks and ejecta of craters of known age. Crater frequency distributions and data on the craters Tycho, North Ray, Cone, and South Ray, when compared with the flux of objects measured by the Apollo Passive Seismometer, suggest that the flux of objects has been relatively constant over the last 100 m.y. (within 1/3 to 3 times of the flux estimated for Tycho).Steady-state frequency distributions for craters in several morphologic categories formed the basis for estimating the relative ages of craters and surfaces in a system used during the Apollo landing site mapping program of the U.S. Geological Survey. The relative ages in this system are converted to model absolute ages that have a rather broad range of values. The range of values of the absolute ages are between about 1/3 to 3 times the assigned model absolute age.     

Sun–Earth L_{1} and L_{2} to Moon transfers exploiting natural dynamics

This paper examines the design of transfers from the Sun–Earth libration orbits, at the \(L_{1}\) and \(L_{2}\) points, towards the Moon using natural dynamics in order to assess the feasibility of future disposal or lifetime extension operations. With an eye to the probably small quantity of propellant left when its operational life has ended, the spacecraft leaves the libration point orbit on an unstable invariant manifold to bring itself closer to the Earth and Moon. The total trajectory is modeled in the coupled circular restricted three-body problem, and some preliminary study of the use of solar radiation pressure is also provided. The concept of survivability and event maps is introduced to obtain suitable conditions that can be targeted such that the spacecraft impacts, or is weakly captured by, the Moon. Weak capture at the Moon is studied by method of these maps. Some results for planar Lyapunov orbits at \(L_{1}\) and \(L_{2}\) are given, as well as some results for the operational orbit of SOHO.     

Indirect transfer to the Earth–Moon L1 libration point

The Earth–Moon L1 libration point is proposed as a human gateway for space transportation system of the future. This paper studies indirect transfer using the perturbed stable manifold and lunar flyby to the Earth–Moon L1 libration point. Although traditional studies indicate that indirect transfer to the Earth–Moon L1 libration point does not save much fuel, this study shows that energy efficient indirect transfer using the perturbed stable manifold and lunar flyby could be constructed in an elegant way. The design process is given to construct indirect transfer to the Earth–Moon L1 libration point. Simulation results show that indirect transfer to the Earth–Moon L1 libration point saves about 420 m/s maneuver velocity compared to direct transfer, although the flight time is about 20 days longer.     

Effect of the Earth’s compression on the physical libration of the Moon

The effect of the Earth??s compression on the physical libration of the Moon is studied using a new vector method. The moment of gravitational forces exerted on the Moon by the oblate Earth is derived considering second order harmonics. The terms in the expression for this moment are arranged according to their order of magnitude. The contribution due to a spherically symmetric Earth proves to be greater by a factor of 1.34 × 10⁶ than a typical term allowing for the oblateness. A linearized Euler system of equations to describe the Moon??s rotation with allowance for external gravitational forces is given. A full solution of the differential equation describing the Moon??s libration in longitude is derived. This solution includes both arbitrary and forced oscillation harmonics that we studied earlier (perturbations due to a spherically symmetric Earth and the Sun) and new harmonics due to the Earth??s compression. We posed and solved the problem of spinorbital motion considering the orientation of the Earth??s rotation axis with regard to the axes of inertia of the Moon when it is at a random point in its orbit. The rotation axes of the Earth and the Moon are shown to become coplanar with each other when the orbiting Moon has an ecliptic longitude of L _? = 90° or L _? = 270°. The famous Cassini??s laws describing the motion of the Moon are supplemented by the rule for coplanarity when proper rotations in the Earth-Moon system are taken into account. When we consider the effect of the Earth??s compression on the Moon??s libration in longitude, a harmonic with an amplitude of 0.03?? and period of T ₈ = 9.300 Julian years appears. This amplitude exceeds the most noticeable harmonic due to the Sun by a factor of nearly 2.7. The effect of the Earth??s compression on the variation in spin angular velocity of the Moon proves to be negligible.