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.     

Implications of the Galilean satellites ice envelope explosions III. The origin of the Trojans and of some comets

Secondary explosions of the primary ice fragments ejected in the explosion of the electrolyzed massive ice envelopes of the Galilean satellites are capable of imparting velocities of up to ~5km s^–1 to the secondary fragments. As a result, the secondary fragments can enter the orbits of the irregular satellites (Agafonova and Drobyshevski, 1984b) and the Trojan libration orbits. In the latter case a perturbation velocity of V 0.3–2 km s^–1 is sufficient.The primary fragments ejected by the gravitational perturbations due to the Galilean satellites sunward from Jupiter's sphere of action move faster relative to the Sun than Jupiter does and therefore reach their first aphelion ahead of Jupiter in the neighborhood of L ₄. At the same time the fragments propelled from Jupiter's sphere of action beyond the planet's orbit approach it again in their perihelia behind Jupiter in the region of L ₅. The concentration of the fragments and, hence, the frequency of their collisions and explosions at L ₄ turn out to be much greater than those at L ₅. As a result, the number of the secondary fragments of diameter 15 km captured into libration orbits ahead of Jupiter can be as high as many hundreds and should exceed by more than a factor 3.5 that captured behind Jupiter.Since the icy mix of the fragments contains hydrocarbons and particulate material (silicates and the like), after ice sublimation from the surface layers the Trojans should reveal type C and RD spectra typical for Jupiter's irregular satellites, comet nuclei and other distant ice bodies of similar origin. Among the Trojans there cannot be rocky or metallic objects which are known to exist in the main asteroid belt.It is shown that a velocity perturbation of 150–200 m s^–1 resulting from a purely mechanical impact of two bodies may be sufficient to move collision fragments from the orbits of the Trojans to horseshoe-shaped trajectories with a subsequent transfer to the cometary orbits of Jupiter's family.     

Invisible comets on evolutionary track of short-period comets

We systematically surveyed the orbits of short-period (SP) comets that show a large change of perihelion distance (q) between 1–2 AU (visible comets) and 4–5 AU (invisible comets) during 4400 years. The data are taken from Cosmo-DICE (Nakamura and Yoshikawa 1991a), which is a long-term orbital evolution project for SP comets. Recognizing that q is the most critical element for observability of comets, an invisibility factor (f), defined as the ratio of unobservable time span to observable span during 4400 years, is calculated for each of the large-q-change comets. A detection limit for each comet is obtained from the heliocentric distance at discovery and/or the absolute magnitude at recent apparitions. A mean f value for 35 SP comets with 2.9 J (J is the Tisserand's invariant) is found to be 19.8. This implies that for each visible SP comet of this J-range, at every epoch of time, there exist about 20 invisible comets near the capture orbits by Jupiter, under the assumptions of steady-state flux and ergodicity for the SP-comet population.     

Homoclinic and heteroclinic orbits in the photogravitational restricted three-body problem

We study numerically the asymptotic homoclinic and heteroclinic orbits around the hyperbolic Lyapunov periodic orbits which emanate from Euler's critical points L ₁ and L ₂, in the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these Lyapunov orbits, are also presented. Poincaré surface of sections of these manifolds on appropriate planes and several homoclinic and heteroclinic orbits for the gravitational case as well as for varying radiation factor q ₁, are displayed. Homoclinic-homoclinic and homoclinic-heteroclinic-homoclinic chains which link the interior with the exterior Hill's regions, are illustrated. We adopt the Sun-Jupiter system and assume that only the larger primary radiates. It is found that for small deviations of its value from the gravitational case (q ₁ = 1), the radiation pressure exerts a significant impact on the Hill's regions and on these asymptotic orbits.     

Temporary satellite capture of comets by Jupiter

This paper studies the dynamical evolution of 97 Jupiter-family comets over an 800-year time period. More than two hundred encounters with Jupiter are investigated, with the observed comets moving during a certain period of time in an elliptic jovicentric orbit. In most cases this is an ordinary temporary satellite capture of a comet in Everhart??s sense, not associated with a transition of the small body into Jupiter??s family of satellites. The phenomenon occurs outside the Hill sphere with comets with a high Tisserand constant relative to Jupiter; the comets?? orbits have a small inclination to the ecliptic plane. An analysis of 236 encounters has allowed the determination within the planar pair two-body problem of a region of orbits in the plane (a, e) whose semimajor axes and eccentricities contribute to the phenomenon under study. Comets with orbits belonging to this region experience a temporary satellite capture during some of their encounters; the jovicentric distance function has several minima; and the encounters are characterized by reversions of the line of apsides and some others features of their combination that are intrinsic to comets in this region. Therefore, this region is called a region of comets with specific features in their encounters with Jupiter. Twenty encounters (out of 236), whereby the comet enters an elliptic jovicentric orbit in the Hill sphere, are identified and investigated. The size and shape of the elliptic heliocentric orbits enabling this transition are determined. It is found that in 11 encounters the motion of small bodies in the Hill sphere has features the most important of which is multiple minima of the jovicentric distance function. The study of these 20 encounters has allowed the introduction of the concept of temporary gravitational capture of a small body into the Hill sphere. An analysis of variations in the Tisserand constant in these (20) encounters of the observable comets shows that their motion is unstable in Hill??s sense.     

Orbital resonance in a dissipative medium

Orbital resonances tend to force bodies into noncircular orbits. If a body is also under the influence of an eccentricity-reducing medium, it will experience a secular change in semimajor axis which may be positive or negative depending on whether its orbit is exterior or interior to that of the perturbing body. Thus a dissipative medium can promote either a loss or a gain in orbital energy. This process may explain the resonant structure of the asteroid belt and of Saturn's rings. For reasonable early solar system parameters, it would clear a gap near the 2:1 resonance with Jupiter on a time scale of a few thousand years; the gap width would be comparable to the Kirkwood gap presently at the location in the asteroid belt. Similarly, a gap comparable in width to Cassini's division would be cleared in Saturn's rings at the 2:1 resonance with Mimas in ～10⁶ yr. Most of the material from the gap would be deposited at the outer edge of ring B. The process would also affect the radial distribution of preplanetary material. Moreover, it provides an explanation for the large amplitude of the Titan-Hyperion libration. Consideration of the effects of dissipation on orbits near the stable L₄ and L₅ points of the restricted three-body problem indicates that energy loss causes particles to move away from these points. This results explains the large amplitude of Trojan asteroids about these points and the possible capture of Trojan into orbit about Jupiter.     

Capture of dust grains in exterior resonances with planets

The temporary capture of the dust grains in the exterior resonances with planets is studied in the frames of the planar circular three-body problem with Poynting-Robertson (PR) drag. For the Earth and particles ~ 10 m the resonances 4/5, 5/6, 6/7, 7/8 are shown to be most effective. The capture is only temporary (of order 10⁵ years) and the position of resonance may be calculated from semi-analytical model using averaged disturbing function. These semi-analytical results are confirmed by numerical integration. For various planet this picture changes as with increasing planetary mass the more exterior resonances become more important. We showed that for Jupiter (at least in the space between Jupiter and Saturn) the resonance 1/2 plays the dominant role. The capture time is here several myr but again eccentricity is evolving to eccentricity e ₀ ~ 0.48 of libration point for this resonance.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation in the present configuration of the two planets

The dynamics of the two Jupiter triangular libration points perturbed by Saturn is studied in this paper. Unlike some previous works that studied the same problem via the pure numerical approach, this study is done in a semianalytic way. Using a literal solution, we are able to explain the asymmetry of two orbits around the two libration points with symmetric initial conditions. The literal solution consists of many frequencies. The amplitudes of each frequency are the same for both libration points, but the initial phase angles are different. This difference causes a temporary spatial asymmetry in the motions around the two points, but this asymmetry gradually disappears when the time goes to infinity. The results show that the two Jupiter triangular libration points should have symmetric spatial stable regions in the present status of Jupiter and Saturn. As a test of the literal solution, we study the resonances that have been extensively studied in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006). The resonance structures predicted by our analytic theory agree well with those found in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006) via a numerical approach. Two kinds of chaotic orbits are discussed. They have different behaviors in the frequency map. The first kind of chaotic orbits (inner chaotic orbits) is of small to moderate amplitudes, while the second kind of chaotic orbits (outer chaotic orbits) is of relatively larger amplitudes. Using analytical theory, we qualitatively explain the transition process from the inner chaotic orbits to the outer chaotic orbits with increasing amplitudes. A critical value of the diffusion rate is given to separate them in the frequency map. In a forthcoming paper, we will study the same problem but keep the planets in migration. The time asymmetry, which is unimportant in this paper, may cause an observable difference in the two Jupiter Trojan groups during a very fast planet migration process.     

Evidence for Lévy Random Walks in the Evolution of Comets from the Oort Cloud

With a simple map model, derived within the framework of the planar circular restricted three-body problem (SunuuJupiteruucomet), we study the dynamical evolution of near-parabolic comets under the perturbation of Jupiter. The commonly adopted random walk assumption about the energy evolution of the comets is examined. Numerical results show that for the comets on Jupiter-crossing orbits, due to the large energy changes with Jupiter per passage, the statistical evolution of the cometary energy follows a Lévy random walk, thus statistically the final energy parameter that a comet reaches is linked to the number of passages by a power law K _f n _f. The mechanism that generates the Lévy random walk is explained in this model.     

On the Stability Regions of the Trojan Asteroids

The orbits of fictitious bodies around Jupiter’s stable equilibrium points L ₄ and L ₅ were integrated for a fine grid of initial conditions up to 100 million years. We checked the validity of three different dynamical models, namely the spatial, restricted three body problem, a model with Sun, Jupiter and Saturn and also the dynamical model with the Outer Solar System (Jupiter to Neptune). We determined the chaoticity of an orbit with the aid of the Lyapunov Characteristic Exponents (=LCE) and used also a method where the maximum eccentricity of an orbit achieved during the dynamical evolution was examined. The goal of this investigation was to determine the size of the regions of motion around the equilibrium points of Jupiter and to find out the dependance on the inclination of the Trojan’s orbit. Whereas for small inclinations (up to i=20°) the stable regions are almost equally large, for moderate inclinations the size shrinks quite rapidly and disappears completely for i>60^°. Additionally, we found a difference in the dynamics of orbits around L ₄ which – according to the LCE – seem to be more stable than the ones around L ₅.     

Relationship of comets and planets

We analyze our earlier data on the numerical integration of the equations of motion for 274 short-period comets (with the period P<200 yr) on a time interval of 6000 yr. As many as 54 comets had no close approaches to planets, 13 comets passed through the Saturnian sphere of action, and one comet passed through the Uranian sphere of action. The orbital elements of these 68 comets changed by no more than ±3 percent in a space of 6000 yr. As many as 206 comets passed close to Jupiter. We confirm Everhart’s conclusion that Jupiter can capture long-period comets with q = 4–6 AU and i < 9° into short-period orbits. We show that nearly parabolic comets cross the solar system mainly in the zone of terrestrial planets. No relationship of nearly parabolic comets and terrestrial planets was found for the epoch of the latest apparition of comets. Guliev’s conjecture about two trans-Plutonian planets is based on the illusory excess of cometary nodes at large heliocentric distances. The existence of cometary nodes at the solar system periphery turns out to be a solely geometrical effect.     

Study of the transfer from the Earth to a halo orbit around the equilibrium pointL 1

The purpose of this paper is to study a transfer strategy from the vicinity of the Earth to a halo orbit around the equilibrium pointL ₁ of the Earth-Sun system. The study is done in the real solar system (we use the DE-118 JPL ephemeris in the simulations of motion) although some simplified models, such as the restricted three body problem (RTBP) and the bicircular problem, have been also used in order to clarify the geometrical aspects of the problem. The approach used in the paper makes use of the hyperbolic character of the halo orbits under consideration. The invariant stable manifold of these orbits enables the transfer to be achieved with, theoretically, only one manoeuvre: the one of insertion into the stable manifold. For the total v required, the figures obtained are similar to the ones given by the standard procedures of optimization.     

Chaos in the solar system

We have accumulated thousands of orbits of test particles in the Solar System from the asteroid belt to beyond the orbit of Neptune. We find that the time for an orbit to make a close encounter with a perturbing planet, T _c ,is a function of the Lyapunov time, T _ty .The relation is log (T _c /T _o )= a + b log (T _ly T _o )where T _o is a fiducial period which we have taken as the period of the principal perturber or the period of the asteroid. There are exceptions to this rule interior to the 2/3 resonance with Jupiter. There, at least in the restricted problem, for sufficiently small Jupiter mass, orbits may have a positive Lyapunov exponent and still be blocked from having a close approach to Jupiter by a zero velocity curve. Of more serious concern is whether the relation holds for purely secular resonances, and if it does, how to choose T _o .This is the case of interest for the planets in the solar system.     

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.     

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 the origin of the unusual orbit of Comet 2P/Encke

The orbit of Comet 2P/Encke is difficult to understand because it is decoupled from Jupiter—its aphelion distance is only 4.1 AU. We present a series of orbital integrations designed to determine whether the orbit of Comet 2P/Encke can simply be the result of gravitational interactions between Jupiter-family comets and the terrestrial planets. To accomplish this, we integrated the orbits of a large number of objects from the trans-neptunian region, through the realm of the giant planets, and into the inner Solar System. We find that at any one time, our model predicts that there should be roughly 12 objects in Encke-like orbits. However, it takes roughly 200 times longer to evolve onto an orbit like this than the typical cometary physical lifetime. Thus, we suggest that (i) 2P/Encke became dormant soon after it was kicked inward by Jupiter, (ii) it spent a significant amount of time inactive while rattling around the inner Solar System, and (iii) it only became active again as the ν₆ secular resonance drove down its perihelion distance.     

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.     

Short-Period Comets with a High Tisserand Constant: III. Kinematics of Low-Velocity Encounters

Within the framework of a pair two-body problem (Sun–Jupiter, Sun–comet), the kinematics of the encounter of a minor body with a planet is investigated. The notion of points of low-velocity tangency of the orbits of the comet and Jupiter, as well as the point of Jovicentric velocity and the low-velocity tangent section of a cometary orbit, is introduced. The conditions and definitions of low-velocity and high-velocity encounters are proposed. The systems of inequalities relating the aand eparameters, which make it possible to single out those comets that are likely to be objects with low-velocity encounters, are presented. The regions of orbits that have low-velocity tangent sections, i.e., regions of low-velocity tangency of orbits, are singled out on the (a, e) plane. These regions agree well with the corresponding parameters of the orbits of real comets whose evolution contains low-velocity encounters with Jupiter.     

The Characteristics and Related Problems of the Orbits Around the Earth-Moon Libration Points

Due to the specific dynamics, the probes located at the halo orbits or Lissajous orbits around the Earth-Moon collinear libration point L₁ or L₂ are always studied in the synodic system to understand their trajectories. In fact, they are also orbiting the Earth in a distant Keplerian ellipse. Because of their intrinsic orbital instability, in the orbit prediction the initial errors propagate more prominently than those of the normal orbiting satellites, this requires special attention in the orbit design, maneuver, and control. Despite of all this, they are similar to the normal orbiting satellites in orbit determination and hardly require other special attentions. In this paper, the quantitative results of error propagation under the unstable dynamics, together with the theoretical analysis are presented. The results of precise orbit determination and short-arc orbit predictions are also shown, and compared with the results from the Beijing Aerospace Control Center.