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.     

An analytic solution for theJ 2 perturbed equatorial orbit

An analytic solution for theJ ₂ perturbed equatorial orbit is obtained in terms of elliptic functions and integrals. The necessary equations for computing the position and velocity vectors, and the time are given in terms of known functions. The perturbed periapsis and apoapsis distances are determined from the roots of a characteristic cubic.     

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.     

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.     

The use of invariant manifolds for transfers between unstable periodic orbits of different energies

Techniques from dynamical systems theory have been applied to the construction of transfers between unstable periodic orbits that have different energies. Invariant manifolds, trajectories that asymptotically depart or approach unstable periodic orbits, are used to connect the initial and final orbits. The transfer asymptotically departs the initial orbit on a trajectory contained within the initial orbit’s unstable manifold and later asymptotically approaches the final orbit on a trajectory contained within the stable manifold of the final orbit. The manifold trajectories are connected by the execution of impulsive maneuvers. Two-body parameters dictate the selection of the individual manifold trajectories used to construct efficient transfers. A bounding sphere centered on the secondary, with a radius less than the sphere of influence of the secondary, is used to study the manifold trajectories. A two-body parameter, κ, is computed within the bounding sphere, where the gravitational effects of the secondary dominate. The parameter κ is defined as the sum of two quantities: the difference in the normalized angular momentum vectors and eccentricity vectors between a point on the unstable manifold and a point on the stable manifold. It is numerically demonstrated that as the κ parameter decreases, the total cost to complete the transfer decreases. Preliminary results indicate that this method of constructing transfers produces a significant cost savings over methods that do not employ the use of invariant manifolds.     

Analytical solutions for J 2-perturbed unbounded equatorial orbits

While solutions for bounded orbits about oblate spheroidal planets have been presented before, similar solutions for unbounded motion are scarce. This paper develops solutions for unbounded motion in the equatorial plane of an oblate spheroidal planet, while taking into account only the J ₂ harmonic in the gravitational potential. Two cases are distinguished: A pseudo-parabolic motion, obtained for zero total specific energy, and a pseudo-hyperbolic motion, characterized by positive total specific energy. The solutions to the equations of motion are expressed using elliptic integrals. The pseudo-parabolic motion unveils a new orbit, termed herein the fish orbit, which has not been observed thus far in the perturbed two-body problem. The pseudo-hyperbolic solutions show that significant differences exist between the Keplerian flyby and the flyby performed under the the J ₂ zonal harmonic. Numerical simulations are used to quantify these differences.     

Integrable approximation of J 2-perturbed relative orbits

Most existing satellite relative motion theories utilize mean elements, and therefore cannot be used for calculating long-term bounded perturbed relative orbits. The goal of the current paper is to find an integrable approximation for the relative motion problem under the J ₂ perturbation, which is adequate for long-term prediction of bounded relative orbits with arbitrary inclinations. To that end, a radial intermediary Hamiltonian is utilized. The intermediary Hamiltonian retains the original structure of the full J ₂ Hamiltonian, excluding the latitude dependence. This formalism provides integrability via separation, a fact that is utilized for finding periodic relative orbits in a local-vertical local-horizontal frame and determine an initialization scheme that yields long-term boundedness of the relative distance. Numerical experiments show that the intermediary-based computation of orbits provides long-term bounded orbits in the full J ₂ problem for various inclinations. In addition, a test case is shown in which the radial intermediary-based initial conditions of the chief and deputy satellites yield bounded relative distance in a high-precision orbit propagator.     

The transition from elliptic to hyperbolic orbits in the two-body problem by slow loss of mass

Transition from elliptic to hyperbolic orbits in the two-body problem with slowly decreasing mass is investigated by means of asymptotic approximations.Analytical results by Verhulst and Eckhaus are extended to construct approximate solutions for the true anomaly and the eccentricity of the osculating orbit if the initial conditions are nearly-parabolic. It becomes clear that the eccentricity will monotonously increase with time for all mass functions satisfying a Jeans-Eddington relation and even for a larger set of functions. To illustrate these results quantitatively we calculate the eccentricity as a function of time for Jeans-Eddington functionsn=0(1) 5 and 18 nearly-parabolic initial conditions to find that 93 out of 108 elliptic orbits become hyperbolic.     

Reduction and relative equilibria for the two-body problem on spaces of constant curvature

We consider the two-body problem on surfaces of constant nonzero curvature and classify the relative equilibria and their stability. On the hyperbolic plane, for each \(q>0\) we show there are two relative equilibria where the masses are separated by a distance q. One of these is geometrically of elliptic type and the other of hyperbolic type. The hyperbolic ones are always unstable, while the elliptic ones are stable when sufficiently close, but unstable when far apart. On the sphere of positive curvature, if the masses are different, there is a unique relative equilibrium (RE) for every angular separation except \(\pi /2\). When the angle is acute, the RE is elliptic, and when it is obtuse the RE can be either elliptic or linearly unstable. We show using a KAM argument that the acute ones are almost always nonlinearly stable. If the masses are equal, there are two families of relative equilibria: one where the masses are at equal angles with the axis of rotation (‘isosceles RE’) and the other when the two masses subtend a right angle at the centre of the sphere. The isosceles RE are elliptic if the angle subtended by the particles is acute and is unstable if it is obtuse. At \(\pi /2\), the two families meet and a pitchfork bifurcation takes place. Right-angled RE are elliptic away from the bifurcation point. In each of the two geometric settings, we use a global reduction to eliminate the group of symmetries and analyse the resulting reduced equations which live on a five-dimensional phase space and possess one Casimir function.     

Time-varying transformations for Hill–Clohessy–Wiltshire solutions in elliptic orbits

The relative motion of chief and deputy satellites in close proximity with orbits of arbitrary eccentricity can be approximated by linearized time-periodic equations of motion. The linear time-invariant Hill–Clohessy–Wiltshire equations are typically derived from these equations by assuming the chief satellite is in a circular orbit. Two Lyapunov–Floquet transformations and an integral-preserving transformation are here presented which relate the linearized time-varying equations of relative motion to the Hill–Clohessy–Wiltshire equations in a one-to-one manner through time-varying coordinate transformations. These transformations allow the Hill–Clohessy–Wiltshire equations to describe the linearized relative motion for elliptic chief satellites.     

Instability and Diffusion in the Elliptic Restricted Three-body Problem

The importance of the stability characteristics of the planar elliptic restricted three-body problem is that they offer insight about the general dynamical mechanisms causing instability in celestial mechanics. To analyze these concerns, elliptic–elliptic and hyperbolic–elliptic resonance orbits (periodic solutions with lower period) are numerically discovered by use of Newton's differential correction method. We find indications of stability for the elliptic–elliptic resonance orbits because slightly perturbed orbits define a corresponding two-dimensional invariant manifold on the Poincaré surface-section. For the resonance orbit of the hyperbolic–elliptic type, we show numerically that its stable and unstable manifolds intersect transversally in phase-space to induce instability. Then, we find indications that there are orbits which jump from one resonance zone to the next before escaping to infinity. This phenomenon is related to the so-called Arnold diffusion. This revised version was published online in August 2006 with corrections to the Cover Date.     

Analytical short-term orbit predictions withJ 3 andJ 4 in terms of KS elements

Analytical theory for short-term orbit motion of satellite orbits with Earth's zonal harmonicsJ ₃ andJ ₄ is developed in terms of KS elements. Due to symmetry in KS element equations, only two of the nine equations are integrated analytically. The series expansions include terms of third power in the eccentricity. Numerical studies with two test cases reveal that orbital elements obtained from the analytical expressions match quite well with numerically integrated values during a revolution. Typically for an orbit with perigee height, eccentricity and inclination of 421.9 km, 0.17524 and 30 degrees, respectively, maximum differences of 27 and 25 cm in semimajor axis computation are noted withJ ₃ andJ ₄ term during a revolution. For application purposes, the analytical solutions can be used for accurate onboard computation of state vector in navigation and guidance packages.     

Optimal transfers between unstable periodic orbits using invariant manifolds

This paper presents a method to construct optimal transfers between unstable periodic orbits of differing energies using invariant manifolds. The transfers constructed in this method asymptotically depart the initial orbit on a trajectory contained within the unstable manifold of the initial orbit and later, asymptotically arrive at the final orbit on a trajectory contained within the stable manifold of the final orbit. Primer vector theory is applied to a transfer to determine the optimal maneuvers required to create the bridging trajectory that connects the unstable and stable manifold trajectories. Transfers are constructed between unstable periodic orbits in the Sun–Earth, Earth–Moon, and Jupiter-Europa three-body systems. Multiple solutions are found between the same initial and final orbits, where certain solutions retrace interior portions of the trajectory. All transfers created satisfy the conditions for optimality. The costs of transfers constructed using manifolds are compared to the costs of transfers constructed without the use of manifolds. In all cases, the total cost of the transfer is significantly lower when invariant manifolds are used in the transfer construction. In many cases, the transfers that employ invariant manifolds are three times more efficient, in terms of fuel expenditure, than the transfer that do not. The decrease in transfer cost is accompanied by an increase in transfer time of flight.     

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

On topological stability in the general three-body problem

The confining curves in the general three-body problem are studied; the role of the integralc ² h (angular momentum squared times energy) as bifurcation parameter is established in a very simple way by using symmetries and changes of scale. It is well known (Birkhoff, 1927) that the bifurcations of the level manifolds of the classical integrals occur at the Euler-Lagrange relative equilibrium configurations. For small values of the mass ratio ε=m ₃/m ₂ both the positions of the collinear equilibrium points and thec ² h integral are expanded in power series of ε. In this way the relationship is found between the confining curves resulting from thec ² h integral in the general problem, and the zero velocity curves given by the Jacobi integral in the corresponding restricted problem. For small values of ε the singular confining curves in the general and in the restricted problem are very similar, but they do not correspond to each other: the offset of the two bifurcation values is, in the usual, system of units of the restricted problem, about one half of the eccentricity squared of the orbits of the two larger bodies. This allows the definition of an approximate stability criterion, that applies to the systems with small ε, and quantifies the qualitatively well known destabilizing effect of the eccentricity of the binary on the third body. Because of this destabilizing effect the third body cannot be bounded by any topological criterion based on the classical integrals unless its mass is larger than a minimum value. As an example, the three-body systems formed by the Sun, Jupiter and one of the small planets Mercury, Mars, Pluto or anyone of the asteroids are found to be ‘unstable’, i.e. there is no way of proving, with the classical integrals, that they cannot cross the orbit of Jupiter. This can be reliably checked with the approximate stability criterion, that given for the most important three-body subsystems of the Solar System essentially the same information on ‘stability’ as the full computation of thec ² h integral and of the bifurcation values.     

Mars interplanetary trajectory design via Lagrangian points

With the increase in complexities of interplanetary missions, the main focus has shifted to reducing the total delta-V for the entire mission and hence increasing the payload capacity of the spacecraft. This paper develops a trajectory to Mars using the Lagrangian points of the Sun-Earth system and the Sun-Mars system. The whole trajectory can be broadly divided into three stages: (1) Trajectory from a near-Earth circular parking orbit to a halo orbit around Sun-Earth Lagrangian point L₂. (2) Trajectory from Sun-Earth L₂ halo orbit to Sun-Mars L₁ halo orbit. (3) Sun-Mars L₁ halo orbit to a circular orbit around Mars. The stable and unstable manifolds of the halo orbits are used for halo orbit insertion. The intermediate transfer arcs are designed using two-body Lambert’s problem. The total delta-V for the whole trajectory is computed and found to be lesser than that for the conventional trajectories. For a 480 km Earth parking orbit, the total delta-V is found to be 4.6203 km/s. Another advantage in the present approach is that delta-V does not depend upon the synodic period of Earth with respect to Mars.     

Some properties of the dumbbell satellite attitude dynamics

The dumbbell satellite is a simple structure consisting of two point masses connected by a massless rod. We assume that it moves around the planet whose gravity field is approximated by the field of the attracting center. The distance between the point masses is assumed to be much smaller than the distance between the satellite’s center of mass and the attracting center, so that we can neglect the influence of the attitude dynamics on the motion of the center of mass and treat it as an unperturbed Keplerian one. Our aim is to study the satellite’s attitude dynamics. When the center of mass moves on a circular orbit, one can find a stable relative equilibrium in which the satellite is permanently elongated along the line joining the center of mass with the attracting center (the so called local vertical). In case of elliptic orbits, there are no stable equilibrium positions even for small values of the eccentricity. However, planar periodic motions are determined, where the satellite oscillates around the local vertical in such a way that the point masses do not leave the orbital plane. We prove analytically that these planar periodic motions are unstable with respect to out-of-plane perturbations (a result known from numerical investigations cf. Beletsky and Levin Adv Astronaut Sci 83, 1993). We provide also both analytical and numerical evidences of the existence of stable spatial periodic motions.     

Distant quasi-periodic orbits around Mercury

In this paper, distant quasi-periodic orbits around Mercury are studied for future Mercury missions. All of these orbits have relatively large sizes, with their altitudes near or above the Mercury sphere of influence. The research is carried out in the framework of the elliptic restricted three-body problem (ER3BP) to account for the planet’s non-negligible orbital eccentricity. Retrograde and prograde quasi-periodic trajectories in the planar ER3BP are generalized from periodic orbits in the CR3BP by the homotopy algorithm, and the shape evolution of such quasi-periodic trajectories around Mercury is investigated. Numerical simulations are performed to evaluate the stability of these distant orbits in the long term. These two classes of orbits present different characteristics: retrograde orbits can maintain shape stability with a large size, although the trajectories in some regions may oscillate with larger amplitudes; for prograde orbits, the range of existence is much smaller, and their trajectories easily move away from the vicinity of Mercury when the orbits become larger. Distant orbits can be used to explore the space environment in the vicinity of Mercury, and some orbits can be taken as transfer orbits for low-cost Mercury return missions or other programs for their high maneuverability.     

Long-term orbit computations with KS uniformly regular canonical elements with oblateness

This paper concerns with the study of KS uniformly regular canonical elements with Earth's oblateness. These elements, ten in number, are all constant in the unperturbed motion and even in the perturbed motion, the substitution is straightforward and elementary due to the transformation laws being explicit and closed expression. By utilizing the recursion formulas of Legendre's polynomials, we are able to include any number of Earth's zonal harmonics J _n in the package and also economize the computations. A fixed step-size fourth-order Runge-Kutta-Gill method is employed for numerical integration of the canonical equations.Utilizing 5 test cases covering a large range of semimajor axis and eccentricity, we have carried out computations to study the effects of Earth's zonal harmonics (up to J ₃₆) and integration step-size variation. Bilinear relations and energy equation are used for checking the accuracies of numerical integration. From the application point of view, the package is utilized to study the behaviour of 900 km height near-circular sun-synchronous satellite orbit over a longer duration of 220 days time (nearly 3078 revolutions) and the necessity of including more number of Earth's zonal harmonic terms is noticed. The package is also used to study the effect of higher zonal harmonics on three 900 km height near-circular orbits with inclinations of 60, 63.2, and 65 degrees, by including Earth's zonal harmonics up to J ₂₄. The mean eccentricity (e _m) is found to have long-periods of 459.6, 6925.1 and 1077.6 days, respectively. Sharp changes in the variation of _m near the minima to e_m are noticed. The values of _m are found to be very near to +-90 degrees at the extrema of e_m. The same orbit is employed to study the effect of variation of inclination from 0 to 180 degrees on long-period (T) of eccentricity with J ₂ to J ₂₄ terms. T is found to increase rapidly as we proceed towards the critical inclinations.