High-order solutions of invariant manifolds associated with libration point orbits in the elliptic restricted three-body system

High-order analytical solutions of invariant manifolds, associated with Lissajous and halo orbits in the elliptic restricted three-body problem (ERTBP), are constructed in this paper. The equations of motion of ERTBP in the pulsating synodic coordinate system have five equilibrium points, and the three collinear libration points as well as the associated center manifolds are unstable. In our calculation, the general solutions of the invariant manifolds associated with Lissajous and halo orbits around collinear libration points are expressed as power series of five parameters: the orbital eccentricity, two amplitudes corresponding to the hyperbolic manifolds, and two amplitudes corresponding to the center manifolds. The analytical solutions up to arbitrary order are constructed by means of Lindstedt–Poincaré method, and then the center and invariant manifolds, transit and non-transit trajectories in ERTBP are all parameterized. Since the circular restricted three-body problem (CRTBP) is a particular case of ERTBP when the eccentricity is zero, the general solutions constructed in this paper can be reduced to describe the dynamics around the collinear libration points in CRTBP naturally. In order to check the validity of the series expansions constructed, the practical convergence of the series expansions up to different orders is studied.     

Analytical invariant manifolds near unstable points and the structure of chaos

It is known that the asymptotic invariant manifolds around an unstable periodic orbit in conservative systems can be represented by convergent series (Cherry, Proc Lond Math Soc ser 2, 27:151–170, 1926; Moser, Commun Pure Appl Math 9:673, 1956 and 11:257, 1958; Moser, Giorgilli, Discret Contin Dyn Syst 7:855, 2001). The unstable and stable manifolds intersect at an infinity of homoclinic points, generating a complicated homoclinic tangle. In the case of simple mappings it was found (Da Silva Ritter et al., Phys D 29:181, 1987) that the domain of convergence of the formal series extends to infinity along the invariant manifolds. This allows in practice the study of the homoclinic tangle using only series. However in the case of Hamiltonian systems, or mappings with a finite analyticity domain, the convergence of the series along the asymptotic manifolds is also finite. Here, we provide numerical indications that the convergence does not reach any homoclinic points. We discuss in detail the convergence problem in various cases and we find the degree of approximation of the analytical invariant manifolds to the real (numerical) manifolds as (i) the order of truncation of the series increases, and (ii) we use higher numerical precision in computing the coefficients of the series. Then we introduce a new method of series composition, by using action-angle variables, that allows the calculation of the asymptotic manifolds up to an a arbitrarily large extent. This is the first case of an analytic development that allows the computation of the invariant manifolds and their intersections in a Hamiltonian system for an extent long enough to allow the study of homoclinic chaos by analytical means.     

Relative equilibrium solutions in the four body problem

Beyond the casen=3 little was known about relative equilibrium solutions of then-body problem up to recent years. Palmore's work provides in the general case much useful information. In the casen=4 he gives the totality of solutions when the four masses are equal and studies some degeneracies. We present here a survey of solutions for arbitrary masses, discussing the manifolds of degeneracy. The ordering of restricted potentials allows a counting of the number of bifurcation sets and different invariant manifolds. An analysis of linear stability is done in the restricted and general cases. As a result, values of the masses ensuring linear stability are given.     

Combined low-thrust propulsion and invariant manifold trajectories to capture NEOs in the Sun–Earth circular restricted three-body problem

In this paper, a method to capture near-Earth objects (NEOs) incorporating low-thrust propulsion into the invariant manifolds technique is investigated. Assuming that a tugboat-spacecraft is in a rendez-vous condition with the candidate asteroid, the aim is to take the joint spacecraft-asteroid system to a selected periodic orbit of the Sun–Earth restricted three-body system: the orbit can be either a libration point periodic orbit (LPO) or a distant prograde periodic orbit (DPO) around the Earth. In detail, low-thrust propulsion is used to bring the joint spacecraft-asteroid system from the initial condition to a point belonging to the stable manifold associated to the final periodic orbit: from here onward, thanks to the intrinsic dynamics of the physical model adopted, the flight is purely ballistic. Dedicated guided and capture sets are introduced to exploit the combined use of low-thrust propulsion with stable manifolds trajectories, aiming at defining feasible first guess solutions. Then, an optimal control problem is formulated to refine and improve them. This approach enables a new class of missions, whose solutions are not obtainable neither through the patched-conics method nor through the classic invariant manifolds technique.     

Transport orbits in an equilateral restricted four-body problem

In this paper we consider a restricted equilateral four-body problem where a particle of negligible mass is moving under the Newtonian gravitational attraction of three masses (called primaries) which move on circular orbits around their center of masses such that their configuration is always an equilateral triangle (Lagrangian configuration). We consider the case of two bodies of equal masses, which in adimensional units is the parameter of the problem. We study numerically the existence of families of unstable periodic orbits, whose invariant stable and unstable manifolds are responsible for the existence of homoclinic and heteroclinic connections, as well as of transit orbits traveling from and to different regions. We explore, for three different values of the mass parameter, what kind of transits and energy levels exist for which there are orbits with prescribed itineraries visiting the neighborhood of different primaries.     

Motions close to escapes in the rhomboidal four body problem

In this work we study escape and capture orbits in the planar rhomboidal 4-body problem in a level of constant negative energy. There are only two different values of the masses here. By using numerical analysis, we show certain transversal intersections of the invariant manifolds of parabolic orbits. We then introduce Symbolic Dynamics when the mass ratio is small, and when it is close to one. In the first case the escapes or captures predominate in the direction of one of the diagonals of the rhombus, while in the second case we find solutions escaping or being captured in the direction of both possible diagonals.     

On the Centre Manifold of Collinear Points in the Planar Three-Body Problem

We study the existence of invariant tori in a neighbourhood of the collinear equilibrium points of the planar three-body problem. To this end some properties of the normal form of the Hamiltonian reduced to the 4D central manifold are proved. Using this normal form, we show that the nondegeneracy conditions of KAM theorem are satisfied for all positive masses, including the 2:1 resonance case. The evaluation of the conditions is done numerically.     

Asymptotic orbits in the (<Emphasis Type="Italic">N</Emphasis>+1)-body ring problem

In this paper we study the asymptotic solutions of the (N+1)-body ring planar problem, N of which are finite and ν=N−1 are moving in circular orbits around their center of masses, while the Nth+1 body is infinitesimal. ν of the primaries have equal masses m and the Nth most-massive primary, with m ₀=β m, is located at the origin of the system. We found the invariant unstable and stable manifolds around hyperbolic Lyapunov periodic orbits, which emanate from the collinear equilibrium points L ₁ and L ₂. We construct numerically, from the intersection points of the appropriate Poincaré cuts, homoclinic symmetric asymptotic orbits around these Lyapunov periodic orbits. There are families of symmetric simple-periodic orbits which contain as terminal points asymptotic orbits which intersect the x-axis perpendicularly and tend asymptotically to equilibrium points of the problem spiraling into (and out of) these points. All these families, for a fixed value of the mass parameter β=2, are found and presented. The eighteen (more geometrically simple) families and the corresponding eighteen terminating homo- and heteroclinic symmetric asymptotic orbits are illustrated. The stability of these families is computed and also presented.     

A greedy global search algorithm for connecting unstable periodic orbits with low energy cost.

A method for space mission trajectory design is presented in the form of a greedy global search algorithm. It uses invariant manifolds of unstable periodic orbits and its main advantage is that it performs a global search for the suitable legs of the invariant manifolds to be connected for a preliminary transfer design, as well as the appropriate points of the legs for maneuver application. The designed indirect algorithm bases the greedy choice on the optimality conditions that are assumed for the theoretical minimum transfer cost of a spacecraft when using invariant manifolds. The method is applied to a test case space mission design project in the Earth–Moon system and is found to compare favorably with previous techniques applied to the same project.     

Asymptotic Orbits at the Triangular Equilibria in the Photogravitational Restricted Three-Body Problem

We study numerically the asymptotic homoclinic and heteroclinic orbits associated with the triangular equilibrium points L ₄ and L ₅, in the gravitational and the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these critical points, are also presented. Hundreds of asymptotic orbits for equal mass of the primaries and for various values of the radiation pressure are computed and the most interesting of them are illustrated. In the Copenhagen case, which the problem is symmetric with respect to the x- and y-axis, we found and present non-symmetric heteroclinic asymptotic orbits. So pairs of heteroclinic connections (from L ₄ to L ₅ and vice versa) form non-symmetric heteroclinic cycles. The termination orbits (a combination of two asymptotic orbits) of all the simple families of symmetric periodic orbits, in the Copenhagen case, are illustrated.     

The Evolution of the Stable and Unstable Manifold of an Equilibrium Point

We consider the evolution of the stable and unstable manifolds of an equilibrium point of a Hamiltonian system of two degrees of freedom which depends on a parameter,ν. The eigenvalues of the linearized system are complex for ν < 0 and purely imaginary for ν > 0. Thus for ν < 0 the equilibrium has a two‐dimensional stable manifold and a two‐dimensional unstable manifold, but for ν > 0 these stable and unstable manifolds are gone. We study the system defined by the truncated generic normal form in this situation. One of two things happens depending on the sign of a certain quantity in the normal form expansion. In one case the two families detach as a single invariant manifold and recedes from the equilibrium as ν tends away from 0 through positive values. In the other case the stable and unstable manifold are globally connected for ν < 0 and the whole structure of these manifolds shrinks to the equilibrium as ν → 0 and disappears. These considerations have interesting implications about Strömgren's conjecture in celestial mechanics and the blue sky catastrophe of Devaney.     

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\) .     

The coalescence of invariant manifolds and the spiral structure of barred galaxies

In a previous paper (Voglis et al., Paper I), we demonstrated that, in a rotating galaxy with a strong bar, the unstable asymptotic manifolds of the short-period family of unstable periodic orbits around the Lagrangian points L ₁ or L ₂ create correlations among the apocentric positions of many chaotic orbits, thus supporting a spiral structure beyond the bar. In this paper, we present evidence that the unstable manifolds of all the families of unstable periodic orbits near and beyond corotation contribute to the same phenomenon. Our results refer to a N -body simulation, a number of drawbacks of which, as well as the reasons why these do not significantly affect the main results, are discussed. We explain the dynamical importance of the invariant manifolds as due to the fact that they produce a phenomenon of 'stickiness' slowing down the rate of chaotic escape in an otherwise non-compact region of the phase space. We find a stickiness time of the order of 100 dynamical periods, which is sufficient to support a long-living spiral structure. Manifolds of different families become important at different ranges of values of the Jacobi constant. The projections of the manifolds of all the different families in the configuration space produce a pattern due to the 'coalescence' of the invariant manifolds. This follows closely the maxima of the observed   m = 2  component near and beyond corotation. Thus, the manifolds support both the outer edge of the bar and the spiral arms.     

Lagrangian coherent structures in the planar elliptic restricted three-body problem

This study investigates Lagrangian coherent structures (LCS) in the planar elliptic restricted three-body problem (ER3BP), a generalization of the circular restricted three-body problem (CR3BP) that asks for the motion of a test particle in the presence of two elliptically orbiting point masses. Previous studies demonstrate that an understanding of transport phenomena in the CR3BP, an autonomous dynamical system (when viewed in a rotating frame), can be obtained through analysis of the stable and unstable manifolds of certain periodic solutions to the CR3BP equations of motion. These invariant manifolds form cylindrical tubes within surfaces of constant energy that act as separatrices between orbits with qualitatively different behaviors. The computation of LCS, a technique typically applied to fluid flows to identify transport barriers in the domains of time-dependent velocity fields, provides a convenient means of determining the time-dependent analogues of these invariant manifolds for the ER3BP, whose equations of motion contain an explicit dependency on the independent variable. As a direct application, this study uncovers the contribution of the planet Mercury to the Interplanetary Transport Network, a network of tubes through the solar system that can be exploited for the construction of low-fuel spacecraft mission trajectories. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

On final evolutions in the restricted planar parabolic three-body problem

In this paper, we prove the existence of special type of motions in the restricted planar parabolic three-body problem, of the type exchange, emission–capture, and emission–escape with close passages to collinear and equilateral triangle configuration, among others. The proof is based on a gradient-like property of the Jacobian function when equations of motion are written in a rotating–pulsating reference frame, and the extended phase space is compactified in the time direction. Thus a phase space diffeomorphic to -coordinates (θ, ζ, ζ′) is obtained with the boundary manifolds θ = ± π/2 corresponding to escapes of the binaries when time tends to ± ∞. It is shown there exists exactly five critical points on each boundary, corresponding to classic homographic solutions. The connections of the invariant manifolds associated to the collinear configurations, and stable/unstable sets associated to binary collision on the boundary manifolds, are obtained for arbitrary masses of the primaries. For equal masses extra connections are obtained, which include equilateral configurations. Based on the gradient-like property, a geometric criterion for capture is proposed and is compared with a criterion introduced by Merman (1953b) in the fifties, and an example studied numerically by Kocina (1954).     

Orbits and manifolds near the equilibrium points around a rotating asteroid

We study the orbits and manifolds near the equilibrium points of a rotating asteroid. The linearised equations of motion relative to the equilibrium points in the gravitational field of a rotating asteroid, the characteristic equation and the stable conditions of the equilibrium points are derived and discussed. First, a new metric is presented to link the orbit and the geodesic of the smooth manifold. Then, using the eigenvalues of the characteristic equation, the equilibrium points are classified into 8 cases. A theorem is presented and proved to describe the structure of the submanifold as well as the stable and unstable behaviours of a massless test particle near the equilibrium points. The linearly stable, the non-resonant unstable, and the resonant equilibrium points are discussed. There are three families of periodic orbits and four families of quasi-periodic orbits near the linearly stable equilibrium point. For the non-resonant unstable equilibrium points, there are four relevant cases; for the periodic orbit and the quasi-periodic orbit, the structures of the submanifold and the subspace near the equilibrium points are studied for each case. For the resonant equilibrium points, the dimension of the resonant manifold is greater than 4, and we find at least one family of periodic orbits near the resonant equilibrium points. As an application of the theory developed here, we study relevant orbits for the asteroids 216 Kleopatra, 1620 Geographos, 4769 Castalia and 6489 Golevka.     

A two-level perturbation method for connecting unstable periodic orbits with low fuel cost and short time of flight: application to a lunar observation mission

The proposed method connects two unstable periodic orbits by employing trajectories of their associated invariant manifolds that are perturbed in two levels. A first level of velocity perturbations is applied on the trajectories of the discretized manifolds at the points where they approach the nominal unstable periodic orbit in order to accelerate them. A second level of structured velocity perturbations is applied to trajectories that have already been subjected to first level perturbations in order to approximately meet the necessary conditions for a low \(\varDelta \text {V}\) transfer. Due to this two-level perturbation approach, the number of the trajectories obtained is significantly larger compared with approaches that employ traditional invariant manifolds. For this reason, the problem of connecting two unstable periodic orbits through perturbed trajectories of their manifolds is transformed into an equivalent discrete optimization problem that is solved with a very low computational complexity algorithm that is proposed in this paper. Finally, the method is applied to a lunar observation mission of practical interest and is found to perform considerably better in terms of \(\varDelta \text {V}\) cost and time of flight when compared with previous techniques applied to the same project.     

On magnetostatic equilibrium in a stratified atmosphere

This is a study of the relationship between a magnetic field and its embedding plasma in static equilibrium in a uniform gravity. The ideal gas law is assumed. A system invariant in a given direction is treated first. We show that an exact integral of the equation for force balance across field lines can be derived in a closed form. Using this integral, exact solutions can be generated freely by integrating directly for the distributions of pressure, density and temperature necessary to keep a given magnetic field in equilibrium. Particular solutions are presented for illustration with the solar atmosphere in mind. Extending the treatment to the general system depending on all three spatial coordinates, we arrive at the general form of a theorem of Parker that a magnetic field in static equilibrium must possess certain symmetries. We derive an equation involving the Euler potentials of the magnetic field stipulating these necessary symmetries. Only those magnetic fields satisfying this equation can be in static equilibrium and for these fields, the endowed symmetries make the construction of exact solutions an essentially two dimensional problem as exemplified by the special case of invariance in a given direction.