Trajectory design strategies applied to temporary comet capture including Poincaré maps and invariant manifolds

Temporary satellite capture (TSC) of Jupiter-family comets has been a focus of investigation within the astronomy community for decades. More recently, TSC has been approached from the perspective of dynamical systems theory, within the context of the circular restricted three-body problem (CR3BP). Thus, this problem serves as a testbed for exploring techniques that support trajectory design in similar dynamical regimes. In particular, an association between the invariant manifolds of libration point orbits and the paths of comets that experience TSC has been explored. In this investigation, TSC is further examined from the perspective of transit, that is, transition through the gateways associated with the collinear libration points, in the three-body problem. Periapsis Poincaré maps, previously employed for trajectory design in several investigations, are used to deliver insight into the nature of transit trajectories for energy levels near those associated with several Jupiter-family comets. The evolution of transit trajectories with increasing energy is explored, and the existence of solutions with similar characteristics to the paths of comets P/1996 R2, 82P/Gehrels 3, and 147P/Kushida–Muramatsu is demonstrated within the context of the planar CR3BP using planar periapsis maps. During TSC, the path of comet 111P/Helin–Roman–Crockett is highly inclined with respect to Jupiter; the motion of this comet is examined relative to invariant manifolds in the spatial CR3BP. A method to display the information contained in higher-dimensional Poincaré maps is also demonstrated, and is employed to locate a trajectory possessing the same qualitative characteristics as the path of 111P/Helin–Roman–Crockett.     

Boundary-value problem formulations for computing invariant manifolds and connecting orbits in the circular restricted three body problem

We demonstrate the remarkable effectiveness of boundary value formulations coupled to numerical continuation for the computation of stable and unstable manifolds in systems of ordinary differential equations. Specifically, we consider the circular restricted three-body problem (CR3BP), which models the motion of a satellite in an Earth–Moon-like system. The CR3BP has many well-known families of periodic orbits, such as the planar Lyapunov orbits and the non-planar vertical and halo orbits. We compute the unstable manifolds of selected vertical and halo orbits, which in several cases leads to the detection of heteroclinic connections from such a periodic orbit to invariant tori. Subsequent continuation of these connecting orbits with a suitable end point condition and allowing the energy level to vary leads to the further detection of apparent homoclinic connections from the base periodic orbit to itself, or the detection of heteroclinic connections from the base periodic orbit to other periodic orbits. Some of these connecting orbits are of potential interest in space mission design.     

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.     

A novel algorithm for generating libration point orbits about the collinear points

This paper presents a numerical algorithm that can generate long-term libration points orbits (LPOs) and the transfer orbits from the parking orbits to the LPOs in the circular-restricted three-body problem (CR3BP) and the full solar system model without initial guesses. The families of the quasi-periodic LPOs in the CR3BP can also be constructed with this algorithm. By using the dynamical behavior of LPO, the transfer orbit from the parking orbit to the LPO is generated using a bisection method. At the same time, a short segment of the target LPO connected with the transfer orbit is obtained, then the short segment of LPO is extended by correcting the state towards its adjacent point on the stable manifold of the target LPO with differential evolution algorithm. By implementing the correction strategy repeatedly, the LPOs can be extended to any length as needed. Moreover, combining with the continuation procedure, this algorithm can be used to generate the families of the quasi-periodic LPOs in the CR3BP.     

Study of lunar gravity assist orbits in the restricted four-body problem

In this paper, the lunar gravity assist (LGA) orbits starting from the Earth are investigated in the Sun–Earth–Moon–spacecraft restricted four-body problem (RFBP). First of all, the sphere of influence of the Earth–Moon system (SOIEM) is derived. Numerical calculation displays that inside the SOIEM, the effect of the Sun on the LGA orbits is quite small, but outside the SOIEM, the Sun perturbation can remarkably influence the trend of the LGA orbit. To analyze the effect of the Sun, the RFBP outside the SOIEM is approximately replaced by a planar circular restricted three-body problem, where, in the latter case, the Sun and the Earth–Moon barycenter act as primaries. The stable manifolds associated with the libration point orbit and their Poincaré sections on the SOIEM are applied to investigating the LGA orbit. According to our research, the patched LGA orbits on the Poincaré sections can efficiently distinguish the transit LGA orbits from the non-transit LGA orbits under the RFBP. The former orbits can pass through the region around libration point away from the SOIEM, but the latter orbits will bounce back to the SOIEM. Besides, the stable transit probability is defined and analyzed. According to the variant requirement of the space mission, the results obtained can help us select the LGA orbit and the launch window.     

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.     

Capture orbits and melnikov integrals in the planar three-body problem

The separatrix between bounded and unbounded orbits in the three-body problem is formed by the manifolds of forward and backward parabolic orbits. In an ideal problem these manifolds coincide and form the boundary between the sets of bounded and unbounded orbits. As a mass parameter increases the movement of the parabolic manifolds is approximated numerically and by Melnikov's method. The evidence indicates that for positive values of the mass parameter these manifolds no longer coincide, and that capture and oscillatory orbits exist.     

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.     

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.     

Automatic Method for Identifying Photospheric Bright Points and Granules Observed by Sunrise

In this study, we propose methods for the automatic detection of photospheric features (bright points and granules) from ultra-violet (UV) radiation, using a feature-based classifier. The methods use quiet-Sun observations at 214 nm and 525 nm images taken by Sunrise on 9 June 2009. The function of region growing and mean shift procedure are applied to segment the bright points (BPs) and granules, respectively. Zernike moments of each region are computed. The Zernike moments of BPs, granules, and other features are distinctive enough to be separated using a support vector machine (SVM) classifier. The size distribution of BPs can be fitted with a power-law slope ?1.5. The peak value of granule sizes is found to be about 0.5 arcsec². The mean value of the filling factor of BPs is 0.01, and for granules it is 0.51. There is a critical scale for granules so that small granules with sizes smaller than 2.5 arcsec² cover a wide range of brightness, while the brightness of large granules approaches unity. The mean value of BP brightness fluctuations is estimated to be 1.2, while for granules it is 0.22. Mean values of the horizontal velocities of an individual BP and an individual BP within the network were found to be 1.6 km?s^?1 and 0.9 km?s^?1, respectively. We conclude that the effect of individual BPs in releasing energy to the photosphere and maybe the upper layers is stronger than what the individual BPs release into the network.     

Natural and sail-displaced doubly-symmetric Lagrange point orbits for polar coverage

This paper proposes the use of doubly-symmetric, eight-shaped orbits in the circular restricted three-body problem for continuous coverage of the high-latitude regions of the Earth. These orbits, for a range of amplitudes, spend a large fraction of their period above either pole of the Earth. It is shown that they complement Sun-synchronous polar and highly eccentric Molniya orbits, and present a possible alternative to low thrust pole-sitter orbits. Both natural and solar-sail displaced orbits are considered. Continuation methods are described and used to generate families of these orbits. Starting from ballistic orbits, other families are created either by increasing the sail lightness number, varying the period or changing the sail attitude. Some representative orbits are then chosen to demonstrate the visibility of high-latitude regions throughout the year. A stability analysis is also performed, revealing that the orbits are unstable: it is found that for particular orbits, a solar sail can reduce their instability. A preliminary design of a linear quadratic regulator is presented as a solution to stabilize the system by using the solar sail only. Finally, invariant manifolds are exploited to identify orbits that present the opportunity of a ballistic transfer directly from low Earth orbit.     

Nine Years Of Euv Bright Points

We discuss early results derived from an algorithm that automates the detection, cataloging, and analysis of extreme-ultraviolet (EUV) “bright points” (BP) from 9 years of data acquired by the Extreme-ultraviolet Imaging Telescope (EIT) on the Solar and Heliospheric Observatory (SOHO). The algorithm relies upon the computation of a map of “intensity significance”; this then contains the location of the EUV BPs. By mapping the location of BPs in each image and linking them through long sequences of EIT images we can describe the temporal and spatial variations of the 1.3× 10⁸ EUV BPs observed by SOHO to date. We suggest that there is a considerable amount of physical information about the solar coronal plasma that can be readily gleamed from the BP detection database. In this paper we discuss only a small portion of the possible correlations, but we point to the possibility of BP lifetime distributions that are well described by modified power-laws; the components of which vary with the different temperature filters and with time over the present solar cycle. Dedicated to the memory of John (Ian) Hamilton (1938–2004).     

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.     

Spacecraft trajectories to the L<Subscript>3</Subscript> point of the Sun–Earth three-body problem

Of the three collinear libration points of the Sun–Earth Circular Restricted Three-Body Problem (CR3BP), L₃ is that located opposite to the Earth with respect to the Sun and approximately at the same heliocentric distance. Whereas several space missions have been launched to the other two collinear equilibrium points, i.e., L₁ and L₂, taking advantage of their dynamical and geometrical characteristics, the region around L₃ is so far unexploited. This is essentially due to the severe communication limitations caused by the distant and permanent opposition to the Earth, and by the gravitational perturbations mainly induced by Jupiter and the close passages of Venus, whose effects are more important than those due to the Earth. However, the adoption of a suitable periodic orbit around L₃ to ensure the necessary communication links with the Earth, or the connection with one or more relay satellites located at L₄ or L₅, and the simultaneous design of an appropriate station keeping-strategy, would make it possible to perform valuable fundamental physics and astrophysics investigations from this location. Such an opportunity leads to the need of studying the ways to transfer a spacecraft (s/c) from the Earth’s vicinity to L₃. In this contribution, we investigate several trajectory design methods to accomplish such a transfer, i.e., various types of two-burn impulsive trajectories in a Sun-s/c two-body model, a patched conics strategy exploiting the gravity assist of the nearby planets, an approach based on traveling on invariant manifolds of periodic orbits in the Sun–Earth CR3BP, and finally a low-thrust transfer. We examine advantages and drawbacks, and we estimate the propellant budget and time of flight requirements of each.     

Study of resonances for the restricted 3‐body problem

Our aim is to identify and classify mean‐motion resonances (MMRs) for the coplanar circular restricted three‐body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the coplanar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the previously established resonances for the solar system, i.e., specifically, in regards to the asteroid belt. Previous work employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to observed or proposed resonances, of extra‐solar planets in binary stellar systems (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

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

Escapes and Recurrence in a Simple Hamiltonian System

Many physical systems can be modeled as scattering problems. For example, the motions of stars escaping from a galaxy can be described using a potential with two or more escape routes. Each escape route is crossed by an unstable Lyapunov orbit. The region between the two Lyapunov orbits is where the particle interacts with the system. We study a simple dynamical system with escapes using a suitably selected surface of section. The surface of section is partitioned in different escape regions which are defined by the intersections of the asymptotic manifolds of the Lyapunov orbits with the surface of section. The asymptotic curves of the other unstable periodic orbits form spirals around various escape regions. These manifolds, together with the manifolds of the Lyapunov orbits, govern the transport between different parts of the phase space. We study in detail the form of the asymptotic manifolds of a central unstable periodic orbit, the form of the escape regions and the infinite spirals of the asymptotic manifolds around the escape regions. We compute the escape rate for different values of the energy. In particular, we give the percentage of orbits that escape after a finite number of iterations. In a system with escapes one cannot define a Poincaré recurrence time, because the available phase space is infinite. However, for certain domains inside the lobes of the asymptotic manifolds there is a finite minimum recurrence time. We find the minimum recurrence time as a function of the energy.This revised version was published online in October 2005 with corrections to the Cover Date.     

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.