Next Article in Journal
Structural Dynamic Model Updating with Automatic Mode Identification Using Particle Swarm Optimization
Next Article in Special Issue
Adaptive Gaussian Mixture Model for Uncertainty Propagation Using Virtual Sample Generation
Previous Article in Journal
Recommendation Algorithm for Multi-Task Learning with Directed Graph Convolutional Networks
Previous Article in Special Issue
Integrated Design of Moon-to-Earth Transfer Trajectory Considering Re-Entry Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem

1
School of Astronomy and Space Science, Nanjing University, Nanjing 210023, China
2
Institute of Space Environment and Astrodynamics, Nanjing University, Nanjing 210093, China
3
Key Laboratory of Modern Astronomy and Astrophysics, Ministry of Education, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(18), 8952; https://doi.org/10.3390/app12188952
Submission received: 10 July 2022 / Revised: 29 August 2022 / Accepted: 30 August 2022 / Published: 6 September 2022
(This article belongs to the Special Issue Astrodynamics and Celestial Mechanics)

Abstract

:
Resonance transition periodic orbits exist in the chaotic regions where the 1:1 resonance overlaps with nearby interior or exterior resonances in the circular restricted three-body problem (CRTBP). The resonance transition periodic orbits have important applications for tour missions between the interior and the exterior regions of the system. In this work, following the increase of the mass parameter μ in the CRTBP model, we investigate the breakup of the first-order resonant periodic families and their recombination with the resonance transition periodic families. In this process, we can describe in detail how the 1:1 resonance gradually overlaps with nearby first-order resonances with increasing strength of the secondary’s perturbation. Utilizing the continuation method, features of the resonance transition periodic families are discussed and characterized. Finally, an efficient approach to finding these orbits is proposed and some example resonance transition periodic orbits in the Sun–Jupiter system are presented.

1. Introduction

Mean motion resonance (MMR) or orbital resonance is a common phenomenon both in our Solar system and the exoplanetary system that occurs when there is a ratio of integers between the orbital periods of two objects. For example, a resonance occurs between Neptune and Pluto (3:2 resonance) and Jupiter’s moons (1:2:4 resonance) in our Solar system. For MMRs in the exoplanetary system, numerous suspected or confirmed examples exist [1,2].
The resonant motions in the restricted three-body problem have been studied many times in the past years. The resonant periodic families and their stability corresponding to first-, second-, and other higher-order resonances in the CRTBP and the ERTBP model are reviewed by [3,4]. The transition between different resonances is also an interesting problem and has been investigated by many researchers. For example, Oterma, which is Jupiter’s comet, can make a rapid transition from the interior, near 3:2 resonance, to the exterior, near 2:3 resonance, with the secondary. This phenomenon of Oterma is called resonance hopping [5]. The dynamical mechanism of the resonance transition and capture of Jupiter’s comets is explained by [6,7,8]. Resonance transitions associated with the structure of the weak stability boundary are discussed in [9]. Furthermore, investigations about the resonance hopping effect are also conducted in the Neptune–Planet Nine system [10].
Recently, studies [11,12] examined the natural divergence of the first-order resonances for nearly circular orbits at low eccentricities. They indicate the existence of ‘bridges’ between adjacent first-order resonances at low eccentricities that smoothly connect the apocentric libration zone of a first-order resonance and the pericentric libration zone of the neighboring first-order resonance. This phenomenon is also present with the periodic orbits method in this work. Investigations of these bridges and gaps at low-eccentricity first-order resonance under dissipative effects are presented in [13]. In this paper, we conduct an investigation about the breakup of the first-order resonant periodic families by increasing the mass parameter μ . In the chaotic region where the 1:1 resonance overlaps with nearby first-order resonances, specific periodic orbits with part of their trajectories resembling an interior/exterior resonant orbit and part of their trajectories resembling another interior/exterior resonant orbit are found, which are called the resonance transition periodic orbits (RTPOs). The breakup of the first-order resonance with the increasing mass parameter μ and recombination with the resonance transition periodic families (RTPFs) are introduced in this paper.
In recent years, the application of the resonance phenomenon has attracted much interest in engineering and mission design [14,15,16,17,18,19,20,21]. Practical missions using the resonant periodic orbits already exist. For example, in the extension phase of the IBEX (Interstellar Boundary Explorer) mission, the spacecraft was maneuvered into an orbit that is in the 3:1 resonance with the Moon [22,23]. Another mission is the TESS mission (Transiting Exoplanet Survey Satellite), which is in an orbit in the 2:1 resonance with the Moon [24,25]. By choosing a proper phase with the Moon, this resonant periodic orbit can minimize the radiation dose and avoid long eclipses to improve the efficiency of scientific work and sustain long-term stability. The resonant periodic orbits along with the RTPOs investigated in this work may have a wide application for tour missions in the Earth–Moon and other planetary systems in the future.
The remainder of this paper is structured as follows. In Section 2, the model framework is introduced. The structures of the first-order resonant periodic families are introduced in Section 3. The breakup and recombination process of the RTPFs with the first-order resonant periodic families are introduced in Section 4. Starting from unstable orbits in the 1:1 resonance region, an efficient approach to generating the symmetric RTPFs is proposed in Section 5. Section 6 concludes the study.

2. Model Framework

2.1. The CRTBP Model

The circular restricted three-body problem (CRTBP) model has been widely used to approximate the motions of small particles in the Solar system. In the CRTBP model, it is assumed that the two primaries ( P 1 and P 2 ) have circular orbits about their common center of mass due to their mutual gravitational attraction. A massless particle P moves in the gravitational field generated by the two massive primaries: P 1 with mass m 1 and P 2 with mass m 2 [26]. The massless particle P has no influence on the movement of the two primaries. Usually, the dynamics of the CRTBP model are formulated in a rotating frame (synodic frame) and made non-dimensional by the following units:
[ M ] = m 1 + m 2 , [ L ] = a 12 , [ T ] = [ L ] 3 / G [ M ]
where a 12 is the mean distance of the two primaries. The unit of time [ T ] is chosen such that the orbital period of m 1 and m 2 about their center of mass is 2 π . G is the gravitational constant. Equations of motion for P can be expressed as:
r ¨ + 2 y ˙ x ˙ 0 = ( Ω r ) T
Ω = 1 2 ( x 2 + y 2 ) + μ ( 1 μ ) + 1 μ r 1 + μ r 2
where r is the position vector [ x , y , z ] T , and the dot indicates a derivative concerning the non-dimensional time of the rotating frame. r 1 and r 2 indicate the distance of P with respect to the barycenter of P 1 and P 2 , respectively. The mass parameter μ is defined as μ = m 2 / ( m 1 + m 2 ) . In our studies, three systems, the unperturbed two-body ( μ 0 ), the Sun–Earth ( μ S E ), and the Sun = -Jupiter ( μ S J ) systems, are taken as specific examples to display the results. Their mass parameter values are shown in Table 1. For the CRTBP model, there is an integral of motion in the form of
2 Ω ( x ˙ 2 + y ˙ 2 + z ˙ 2 ) = C
where Ω is the gravitational potential and C is the Jacobi constant.

2.2. Periodic Orbits in the CRTBP Model

The CRTBP model is an autonomous dynamical system. In general, equations of motion (see Equation (2)) describing an autonomous system can be written as
X ˙ = F ( X )
where X denotes the state vector of the small particle P. The solution to this equation is
X ( t ) = Φ ( X 0 , t )
where X 0 is the initial state vector at the initial time t 0 = 0 . The condition that a periodic orbit needs to satisfy is
X ( X 0 , T ) X 0 = 0
where T is the orbital period. To find a periodic orbit is to determine the initial state vector X 0 and the period T that satisfy Equation (7). In total, Equation (7) has seven variables ( X and T), which can be solved by the Newton–Raphson iteration method. Because of the existence of the Jacobi integral (see Equation (4)), only five out of the six equations in Equation (7) are independent. Therefore, we need to fix two variables out of the seven in Equation (7) when solving this equation. Once a solution, i.e., a periodic orbit of Equation (7), is obtained, the whole family of periodic orbits can be generated by the well-known predictor–corrector algorithm [27].
In this work, we focus on two types of periodic orbits in the CRTBP model. The first type of periodic orbit is the planar Lyapunov periodic orbit around the collinear libration points L 1 and L 2 . Families of planar Lyapunov periodic orbits around L 1 and L 2 in the Sun–Earth system are shown as examples in Figure 1a,b, which have been extensively used by previous researchers. Details to compute these periodic orbits can be found in reference [26]. The second type of periodic orbit is the resonant periodic orbit. Orbital resonance is defined by the ratio p:q between the orbital periods, where p and q are incommensurable positive integers. Neglect the secondary’s gravity, i.e., treat the massless particle’s orbit around the primary as an unperturbed Keplerian orbit. We denote the orbital period of the massless particle P as T p and the secondary’s orbital period as T s . A p:q resonant periodic orbit in the synodic frame should strictly satisfy [26,28]
T s T p = p q
When p > q ( p < q ), the resonance is interior (exterior) resonance. Denote k as
k = | p q |
in which k 0 , and k is the resonance order. In this work, k = 1 . One remark is that, in the CRTBP model, the orbital period of the massless particle does not exactly satisfy Equation (8) due to the secondary’s perturbation. Example orbits of the first-order interior resonant periodic orbits (3:2, 4:3, and 5:4 resonances), as well as the exterior resonant periodic orbits (2:3, 3:4, and 4:5 resonances) in the Sun–Earth system are presented in Figure 2.

3. The First-Order Resonant Periodic Families

3.1. The Unperturbed Two-Body Model: μ = 0

In the unperturbed two-body model, two different kinds of symmetric periodic orbits exist in the synodic frame [4,13,29,30]. The periodic orbits related to the circular orbits of the small body are called periodic orbits of the first kind, which belong to the circular family. The periodic orbits related to the elliptic orbits of the small body are called periodic orbits of the second kind, which belong to the resonant periodic family. In Figure 3, the relationships between the circular family and the interior/exterior first-order resonant periodic families are displayed in the T C plane; these are called characteristic curves and are used in previous literature [31,32,33]. For each first-order resonant periodic family, the characteristic curve in Figure 3 represents two branches: one corresponding to the perigee on the positive x-axis and the other to the apogee on the positive x-axis; i.e., they differ only in phase.

3.2. The CRTBP Model: μ > 0

With the increasing of the mass parameter μ , research [31,32,34] indicates that the continuation of circular periodic families of the small body in the CRTBP model is possible in all cases except at the first-order interior resonance ( p + 1 ) :p and exterior resonance p: ( p + 1 ) , p = 1 , 2 , 3 . In the CRTBP model, the two branches of each first-order resonant periodic family split in the T C plane. Example orbits in “Branch-I” (Type-I) and “Branch-II” (Type-II) of the interior 3:2, 4:3, and 5:4 resonances and the exterior 2:3, 3:4, and 4:5 resonances in the Sun–Earth system are presented in Figure 2. Recently, Pan and Hou [4] presented the genealogy of the first-order interior as well as exterior resonant periodic families in the T C plane for the Sun–Earth, Sun–Jupiter, and Earth–Moon systems. Studies indicate that the two branches of each first-order resonant periodic family split and combine with the near-circular branch to form a new family in the T C plane (see Figure 4a,b). This phenomenon is also verified in [11] using the Hamiltonian approach. More specifically, for the interior first-order resonances, the studies [4,11] indicate that starting from p = 1 , “Branch-II” of the ( p + 1 ) :p resonant periodic family and “Branch-I” of the ( p + 2 ) : ( p + 1 ) resonant periodic family combine with part of the near-circular family to form a new periodic family. In the following, we call these new interior periodic families I B ( ( p + 1 ) ( p + 2 ) ) , indicating that they are connecting the ( p + 1 ) :p resonant periodic family with the ( p + 2 ) : ( p + 1 ) resonant periodic family. Similarly, starting from p = 2 , “Branch-I” of the exterior ( p 1 ) :p resonant periodic family is recombined with “Branch-II” of the exterior p: ( p + 1 ) resonant periodic family and part of the near-circular family to form a new periodic family. We denote this new family as family O B ( p ( p + 1 ) ) in this study. Details of the family I B ( 3 4 ) and O B ( 3 4 ) in the Sun–Earth system are presented in Figure 5a,b, respectively.

4. Resonance Transition Periodic Families in the CRTBP Model

In this work, we focus on the RTPFs in the CRTBP model. The RTPOs are the periodic orbits that connect the interior resonance and the exterior resonance by the 1:1 resonance. This kind of periodic orbit has similar orbital behavior to resonance hopping [5,7]. Here, we focus on how the 1:1 resonance overlaps nearby first-order resonances with increasing μ . The breakup of the first-order resonant periodic families and their recombination with the RTPFs in the process of increasing the mass parameter μ are described as follows.

4.1. The RTPFs Related with the 3:2 Resonance

Even for very small μ 0 , the RTPFs already exist, although they are not shown in the above characteristic curves of periodic families. In this section, we take the family I B ( 3 4 ) as an example to show the process of family breakup and recombination with the RTPFs. This process can be viewed as the first-order resonances (here the 3:2 and the 4:3 ones) gradually being overlapped by the 1:1 resonance. In Figure 6, the family I B ( 3 4 ) is continued from the Sun–Earth system to the Sun–Jupiter system, and characteristic curves of some example mass parameters are given. Comparing the characteristic curves in Figure 6a with Figure 5a, we can see that the genealogy of family I B ( 3 4 ) becomes more complex and distorted with the increasing μ .
Taking a closer look at the characteristic curves in Figure 6, we can see that there is no significant change of the family I B ( 3 4 ) in the continuation process starting from the Sun–Earth system ( μ S E ) to the system with μ = 10 4 . However, as the mass parameter further increases to μ = 5 × 10 4 , noticeable changes in the characteristic curve occur in the T C plane. The changes are mainly reflected in two aspects: (1) for the branch near the 3:2 resonance, the characteristic curve bends as the Jacobi constant decreases, accompanied by the process of increasing, decreasing, and then increasing orbital period; (2) for the branch near the 4:3 resonance, the orbital period increases first, then decreases when it reaches a certain value, and finally continues to increase. The detailed characteristic curve of the family I B ( 3 4 ) in the system with μ = 5 × 10 4 is separately presented by the orange line in Figure 7a. As the μ value further increases to μ = μ S J , a closed characteristic curve of family I B ( 3 4 ) forms in the T C plane, which is represented by the red line in Figure 6 and Figure 7b.
In the process of increasing μ from μ = 5 × 10 4 to μ S J , most family members of family I B ( 3 4 ) in the system with μ = 5 × 10 4 still belong to the new family I B ( 3 4 ) in the Sun–Jupiter system ( μ S J ). At the same time, small parts of the family I B ( 3 4 ) in the system with μ = 5 × 10 4 , especially for those close to the 3:2 resonance and the 4:3 resonance, combine with the RTPFs to form the new RTPFs in the Sun–Jupiter system. The bifurcation and recombination details are as follows. In the system with μ = 5 × 10 4 , two RTPFs exist (see the green and purple line in Figure 7a). The green family is denoted as RTPF-3:2(a), and the purple family is denoted as RTPF-3:2(b). Example orbits in these two families are shown in Figure 11a,b, respectively. The RTPF-3:2(a) combines with some family members in “Branch-II” of the 3:2 resonance to form a new RTPF-3:2(a) in the Sun–Jupiter system (see the green line in Figure 7b). Similarly, the RTPF-3:2(b) combines with some family members in “Branch-I” of the 4:3 resonance to form a new RTPF-3:2(b) in the Sun–Jupiter system (see the purple line in Figure 7b). The shapes of the orbits in RTPF-3:2(a) and RTPF-3:2(b) are similar. Both of them have four apogees in the synodic frame, two of which are near the collinear libration point L 1 . From the shape of the example RTPOs in Figure 11a,b, we know that part of the orbits around the collinear L 1 point resembles the planar Lyapunov periodic orbits (see Figure 1a), and part of the orbits far away from the collinear point resembles the 3:2 resonant orbit shown in Figure 2a (type-II). In addition, the characteristic curves of RTPF-3:2(a) and RTPF-3:2(b) (see the green and purple line in Figure 7) indicate that the difference in Jacobi constants between these two RTPFs is not significant. Orbits in the two RTPFs are different due to the difference in the number of revolutions around the libration point L 1 (RTPF-3:2(a): two revolutions; RTPF-3:2(b): three revolutions), which can be seen from Figure 11a,b.

4.2. The RTPFs Related with the 4:3 Resonance

Similar to the family I B ( 3 4 ) which connects the 3:2 resonance to the 4:3 resonance, the family I B ( 4 5 ) also has three components: (1) “Branch-II” of the 4:3 resonance, (2) “Branch-I” of the 5:4 resonance, and (3) the branch of the near-circular family. We continue the family I B ( 4 5 ) from the Sun–Earth system to the Sun–Jupiter system. A similar breakup and recombination process happens during the process of increasing μ .
Characteristic curves of some mass parameters during the process of increasing μ from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ) are displayed in Figure 8. In the T C as well as the x 0 C plane, the characteristic curves of the family I B ( 4 5 ) have no significant change from the Sun–Earth system to the system with μ = 5 × 10 5 . However, when the mass parameter further increases to μ = 10 4 , noticeable changes of the characteristic curve occur in the T C plane in two ways: (1) for the branch near the 4:3 resonance, the characteristic curve bends as the Jacobi constant decreases, accompanied with the process of increasing, decreasing, and then increasing orbital period; (2) for the branch near the 5:4 resonance, the orbital period increases first, then decreases when it reaches a certain value, and finally continues to increase. The characteristic curve of the family I B ( 4 5 ) for the system with μ = 10 4 is separately represented by the yellow line in Figure 9a. Similarly, the branches near the 4:3 resonance and the 5:4 resonance come closer to each other with the increasing μ . Different from the case of the family I B ( 3 4 ) , a closed characteristic curve of the family I B ( 4 5 ) already forms in the T C plane as the μ value increases to μ = 5 × 10 4 . This means the breakup and recombination of the family I B ( 4 5 ) happens for smaller μ values when compared with the family I B ( 3 4 ) . When the mass parameter further increases to μ S J , the size of the closed curve becomes smaller, which is represented by the red line in Figure 9b.
In the process of mass parameter continuation from μ = 10 4 to μ S J , most family members of the family I B ( 4 5 ) in the system with μ = 10 4 still belong to the new family I B ( 4 5 ) in the Sun–Jupiter system ( μ S J ). At the same time, small parts of the family I B ( 4 5 ) in the system with μ = 5 × 10 4 , especially for those close to the 4:3 resonance and the 5:4 resonance, combine with the RTPFs to form the new RTPFs in the Sun–Jupiter system. Two RTPFs work in the bifurcation and recombination process, which are denoted as RTPF-4:3(a) and RTPF-4:3(b) (see the green and purple line in Figure 9a). The RTPF-4:3(a) combines with some family members in “Branch-II” of the 4:3 resonance to form a new RTPF-4:3(a) in the Sun–Jupiter system (see the green line in Figure 9b). Similarly, the RTPF-4:3(b) combines with some family members in “Branch-I” of the 5:4 resonance to form a new RTPF-4:3(b)) in the Sun–Jupiter system (see the purple line in Figure 9b). Example orbits from the RTPF-4:3(a) and the RTPF-4:3(b) are presented in Figure 11c,d, respectively. The shapes of the orbits in RTPF-4:3(a) and the RTPF-4:3(b) are similar. Both of them have five apogees in the synodic frame, two of which are near the collinear libration point L 1 . Part of the orbits close to the collinear libration point resembles the planar Lyapunov periodic orbits around the L 1 point and the remaining part of the orbits resembles the 4:3 resonant periodic orbits (see the Type-II in Figure 2b). In addition, characteristic curves (see the green and purple line in Figure 9) of the RTPF-4:3(a) and RTPF-4:3(b) indicate that the difference in Jacobi constants between these two RTPFs is not significant. Orbits in the RTPF-4:3(a) and the RTPF-4:3(b) are different due to the difference in the number of revolutions near the libration point L 1 (RTPF-4:3(a): two revolutions; RTPF-4:3(b): three revolutions), which can be seen from the example orbits in Figure 11c,d.

4.3. The RTPFs Related with the 2:3 Resonance

The exterior first-order resonant periodic families have similar breakup and recombination processes as the interior resonances during the process of increasing μ . Taking the family O B ( 3 4 ) as an example, the breakup and recombination process with the RTPFs are as follows. The family O B ( 3 4 ) in the CRTBP model has three components: (1) “Branch-I” of the 2:3 resonance, (2) “Branch-II” of the 3:4 resonance, and (3) the branch of the near-circular family (see Figure 5b). We continue the family O B ( 3 4 ) from the Sun–Earth system to the Sun–Jupiter system. Characteristic curves of some mass parameters during the process of increasing μ from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ) are displayed in the T C plane (see Figure 10a). Similarly, when the mass parameter is small (from the Sun–Earth system to the system with μ = 5 × 10 5 ), there are no significant changes in the characteristic curves. However, when the mass parameter further increases to μ = 10 4 , noticeable changes of the characteristic curve occur in the T C plane in two ways: (1) for the branch near the 2:3 resonance, the characteristic curve bends as the Jacobi constant decreases, accompanied with the process of increasing, decreasing, and then increasing orbital period; (2) for the branch near the 3:4 resonance, the orbital period increases first, then decreases when it reaches a certain value, and finally continues to increase. A closed characteristic curve of family O B ( 3 4 ) forms when μ = 5 × 10 4 (see the orange line in Figure 10a). When the mass parameter further increases to μ S J , the size of the closed curve becomes smaller, which is represented by the red line in Figure 10.
With the increasing mass parameter μ , most family members of the family O B ( 3 4 ) in the system with μ = 10 4 still belong to the new family O B ( 3 4 ) in the Sun–Jupiter system ( μ S J ). At the same time, small parts of the family O B ( 3 4 ) in the system with μ = 10 4 , especially for those close to the 2:3 resonance and the 3:4 resonance, combine with the RTPFs to form the new RTPFs in the Sun–Jupiter system. Two RTPFs work in the bifurcation and recombination process, which are denoted as RTPF-2:3(a) and RTPF-2:3(b). RTPF-2:3(a) combines with some family members in “Branch-I” of the 2:3 resonance to form a new RTPF-2:3(a) in the Sun–Jupiter system (see the green line in Figure 10b). Similarly, RTPF-2:3(b) combines with some family members in “Branch-II” of the 3:4 resonance to form a new RTPF-2:3(b)) in the Sun–Jupiter system (see the purple line in Figure 10b). Example orbits from RTPF-2:3(a) and RTPF-2:3(b) in the Sun–Jupiter system are presented in Figure 11e,f, respectively. The shapes of the orbits in RTPF-2:3(a) and the RTPF-2:3(b) are similar. Both of them have three perigees in the synodic frame, two of which are near the collinear libration point L 2 (see Figure 1b). Part of the orbits close to the collinear libration point resembles the planar Lyapunov periodic orbits around the L 2 point and the remaining part of the orbits resembles the 2:3 resonant periodic orbits (see Figure 2d (Type-I)). In addition, characteristic curves of the two RTPFs in Figure 10b (see the green and purple line) indicate that the difference in Jacobi constants between these two RTPFs is not significant. Orbits in the two RTPFs are different due to the difference in the number of revolutions around the libration point L 2 (RTPF-2:3(a): two revolutions; RTPF-2:3(b): three revolutions), which can be seen from Figure 11e,f.

4.4. Summation

In the CRTBP model, the families I B ( p ( p + 1 ) ) connect the “Branch-II” of the interior p: ( p 1 ) resonance with the “Branch-I” of the interior ( p + 1 ) :p resonance. The families O B ( p ( p + 1 ) ) connect the “Branch-I” of the exterior ( p 1 ) :p resonance with the “Branch-II” of the exterior p: ( p + 1 ) resonance. This genealogy of the periodic families has already been pointed out by some previous studies [4,11,31,32].
With the mass parameter μ increasing, starting from larger values of p, the first-order resonant periodic families I B ( p ( p + 1 ) ) and O B ( p ( p + 1 ) ) gradually disappear because they are overlapped by the 1:1 resonance. The RTPOs have part of the orbit resembling the planar Lyapunov orbits around the L 1 or L 2 point, and the other part resembling the interior or exterior resonant periodic orbits. The collinear libration points along with the orbits around them can be taken as unstable ’saddle points’ of the 1:1 resonance. As a result, the RTPOs can be taken as the connection orbits between the 1:1 resonance region and the interior or the exterior resonance regions. By following the breakup and recombination process of the families I B ( p ( p + 1 ) ) and O B ( p ( p + 1 ) ) with the RTPFs, the details of the resonance overlapping process are revealed by the approach of periodic orbits.

5. Computation of the Symmetric RTPOs

According to the studies above, we know that the strength of the 1:1 resonance increases with the increasing μ . It gradually overlaps with nearby interior or exterior resonances, leading to chaos in these overlapped regions. During this process, the RTPFs play a fundamental role in connecting the 1:1 resonance with the interior or the exterior resonances. The RTPOs presented above have part of their orbits resembling the collinear libration point orbits and the remaining part of the orbits resembling the interior or the exterior resonances. Generally, the regions overlapped by the 1:1 resonance are chaotic. There are also RTPOs connecting both the interior and the exterior resonances, through the orbits around L 1 or the L 2 points. In this section, we propose an efficient way to generate the symmetric RTPOs, and some examples in the Sun–Jupiter system are given.

5.1. Methodology

  • Firstly, we restrict the starting position on the Lyapunov orbit around the collinear L 1 or L 2 point in the Sun–Jupiter system and extend the method in the previous study [35] to the Sun–Jupiter system. The amplitude of the Lyapunov orbit is defined as the distance between its right intersection point (see the yellow dot in Figure 12) with the x-axis and the L 1 or the L 2 point. In our work, it is fixed at 0.05 (dimensionless unit) of the Sun–Jupiter distance. One remark is that a different value can be chosen. Initial states of the planar Lyapunov orbit are denoted as ( x 0 , y 0 , 0 , x 0 ˙ , y 0 ˙ , 0 ) . The Lyapunov orbit is symmetric about the x-axis. The initial position is placed on the positive x-axis and the initial velocity is exactly along the positive y-direction. Therefore, the initial state of the Lyapunov orbit is ( x 0 , 0 , 0 , 0 , y 0 ˙ , 0 ) ;
  • We add a tiny perturbation Δ v along the direction of initial velocity y ˙ 0 on the Lyapunov orbit. The velocity change Δ v is in the range of [−0.01, 0.01] in this study, i.e., the departure velocity y ˙ d at the initial position is restricted to be within [−0.01 + y ˙ 0 , 0.01 + y ˙ 0 ], where 0.01 is a dimensionless velocity in the Sun–Jupiter system. Due to the linear stability of the collinear libration point, the orbit can leave the vicinity of the collinear libration point with this tiny perturbation to the initial velocity. The schematic diagram of generating initials of the RTPOs is displayed in Figure 12;
  • By integration, the trajectory leaves the initial position on the x-axis. If the initial velocity y ˙ d is appropriate, it crosses the x-axis again one or more times after some integration time. Generally, at the new intersection point (see the green point in Figure 12), the velocity is not perpendicular to the x-axis. We record the x component of the velocity at the new intersection point and denote it as x ˙ T / 2 . For a different value of the departure velocity y ˙ d , the value of x ˙ T / 2 is different. By varying the values of y ˙ d [ 0.01 + y ˙ 0 , 0.01 + y ˙ 0 ] and the times of the trajectory intersects with the x-axis after departure, the relationships between x ˙ T / 2 and y ˙ d are obtained;
  • The relationships between x ˙ T / 2 and y ˙ d are presented in Figure 13 (cases for L 1 ) and Figure 14 (cases for L 2 ). According to the well-known symmetry property of the CRTBP model ( x , y , z , x ˙ , y ˙ , z ˙ , t ) ( x , y , z , x ˙ , y ˙ , z ˙ , t ) [26], if x ˙ T / 2 = 0 , the trajectory from t = 0 to t = T / 2 is half of a periodic orbit. From Figure 13 and Figure 14, we can see that there exist some dots satisfying x ˙ T / 2 0 within a certain tolerance, which means that we can use these initial states ( x 0 , 0 , 0 , 0 , y ˙ d , 0 ) where x ˙ T / 2 0 to generate the periodic orbits, especially for the RTPOs. After finding these initial states ( x 0 , 0 , 0 , 0 , y ˙ d , 0 ) , we can use the Newton iteration method to find the exact initial velocity y ˙ d that makes x ˙ T / 2 = 0 . The same algorithm has been used to generate the symmetric horseshoe periodic orbits in the CRTBP model [33,36].

5.2. Results

Similar to the computation of the resonant periodic families in the above sections, once a RTPO is found, the whole family can be generated by the predict–correct algorithm. In this section, we limit ourselves to the computation of specific RTPOs, but not the families. Some example orbits are given below.
Examples of the RTPOs related to L 1 as well as L 2 are given in Figure 15 and Figure 16, respectively. In Figure 15, six RTPOs related to L 1 are presented, including (1) a resonance transition between the 1:1 resonance and other resonances (Orbit-1, Orbit-2, and Orbit-4); (2) a resonance transition between the interior and exterior resonances by the 1:1 resonance (Orbit-3); (3) a resonance transition between two interior resonances with the aid of 1:1 resonance (Orbit-5); (4) a resonance transition between two exterior resonances by the 1:1 resonance (Orbit-6). Similarly, six examples of the RTPOs corresponding to L 2 are displayed in Figure 16, which correspond to (1) a resonance transition between the 1:1 and other resonances (Orbit-1, Orbit-2, and Orbit-3); (2) a resonance transition between the interior and exterior resonances by the 1:1 resonance (Orbit-4); (3) a resonance transition between two exterior resonances with the aid of 1:1 resonance (Orbit-5 and Orbit-6).
The RTPOs can stay around L 1 or L 2 for one or more loops before leaving to other resonances. This feature shows that the RTPOs are suitable for missions requiring long-term continuous observations of the second primary. The transition between different resonances by the 1:1 resonance indicates that the RTPOs can be applied to trajectory design for tour missions of the interior and exterior regions of the system. Various RTPOs emanating from the same generating orbit of L 1 or L 2 indicate that the spacecraft can transfer from one RTPO to another by performing a small maneuver in the vicinity of the generating orbit, which may significantly increase the flexibility of the mission design [18]. In our study, the RTPOs are unstable; therefore, station-keeping maneuvers are required if they are applied to practical missions.

6. Conclusions

In summary, details of the breakup and recombination of the interior as well as the exterior resonant periodic families with the RTPFs are revealed in this work with increasing μ . Starting from the first-order resonant periodic families in the Sun–Earth system, we gradually increase the mass parameter μ to the Sun–Jupiter systems. In the process of continuing these periodic families, we observe the gradual changes in the phase space structure of the first-order resonance due to the increasing perturbation strength from the secondary in the CRTBP model. Two phenomena are reported. (1) With the mass parameter μ increasing, the strength of the 1:1 resonance becomes stronger. Starting from large values of p, the interior ( p + 1 ) :p and the exterior p: ( p + 1 ) resonances are gradually overlapped by the 1:1 resonance, in a way that the characteristic curves of these resonance periodic families break up, recombine, and finally disappear. (2) Resonance transition periodic orbits connecting the 1:1 resonance and the interior/exterior resonances appear, forming the bridge connecting the 1:1 resonance region with the interior/exterior resonances. Their appearance indicates that the resonances are already obviously influenced by the 1:1 resonance. Furthermore, starting from the planar Lyapunov orbits in the 1:1 resonance region, we propose an efficient way to compute the symmetric resonance transition periodic orbits, and some examples are presented for the Sun–Jupiter system.
We have to remark that breakup and recombination of interior/exterior resonances with the RTPFs cause difficulties when continuing periodic families from the μ = 0 to large values of μ , because according to our study, the periodic families may undergo many breakups and recombination events during the continuation process, leading to a phase space (for large values of μ ) quite different from that of the two-body problem. Some resonance periodic families may even disappear. We also remark that there are still some limitations of the current work. For example, the regularization method is not used in this work when generating the periodic families, which means that the continuation of the periodic families has to be stopped when the orbit in these families has a close encounter with the primaries. Subsequent work should consider the use of regularization methods to give a more comprehensive understanding of the evolution of the RTPFs. In addition, the symmetry property is now exploited in the calculation of periodic orbits, and it is worthwhile to further investigate some non-symmetric RTPFs, which also have some possible applications for mission designs.

Author Contributions

Conceptualization, S.P. and X.H.; methodology, S.P. and X.H.; formal analysis, S.P. and X.H.; investigation, S.P.; data curation, S.P.; writing—original draft preparation, S.P.; writing—review and editing, X.H.; supervision, X.H.; funding acquisition, S.P. and X.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 11773017) and the China Scholarship Council (No. 202006190153).

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the reviewers and editors for their comments and suggestions, which have greatly improved the presentation of our work. This work was supported by the National Natural Science Foundation of China (No. 11773017) and the China Scholarship Council (No. 202006190153).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rivera, E.J.; Laughlin, G.; Butler, R.P.; Vogt, S.S.; Haghighipour, N.; Meschiari, S. The Lick-Carnegie Exoplanet Survey: A Uranus-Mass Fourth Planet for GJ 876 in an Extrasolar Laplace Configuration. Astrophys. J. 2010, 719, 890–899. [Google Scholar] [CrossRef]
  2. MacDonald, M.G.; Ragozzine, D.; Fabrycky, D.C.; Ford, E.B.; Holman, M.J.; Isaacson, H.T.; Lissauer, J.J.; Lopez, E.D.; Mazeh, T.; Rogers, L.; et al. A Dynamical Analysis of the Kepler-80 System of Five Transiting Planets. Astron. J. 2016, 152, 105. [Google Scholar] [CrossRef]
  3. Hadjidemetriou, J.D. Periodic Orbits and Stability; IAU Colloq. 96: The Few Body Problem; Valtonen, M.J., Ed.; Springer: Berlin/Heidelberg, Germany, 1988; p. 31. [Google Scholar] [CrossRef]
  4. Pan, S.; Hou, X. Review Article: Resonant families of periodic orbits in the restricted three-body problem. Res. Astron. Astrophys. 2022, 22, 072002. [Google Scholar] [CrossRef]
  5. Belbruno, E.; Marsden, B.G. Resonance Hopping in Comets. Astron. J. 1997, 113, 1433. [Google Scholar] [CrossRef]
  6. Howell, K.; Marchand, B.; Lo, M. Temporary Satellite Capture of Short-Period Jupiter Family Comets from the Perspective of Dynamical Systems. J. Astronaut. Sci. 2000, 49. [Google Scholar] [CrossRef]
  7. Koon, W.S.; Lo, M.W.; Marsden, J.E.; Ross, S.D. Heteroclinic connections between periodic orbits and resonance transitions in celestial mechanics. Chaos 2000, 10, 427–469. [Google Scholar] [CrossRef]
  8. Koon, W.S.; Lo, M.W.; Marsden, J.E.; Ross, S.D. Resonance and Capture of Jupiter Comets. Celest. Mech. Dyn. Astron. 2001, 81, 27–38. [Google Scholar] [CrossRef]
  9. Belbruno, E.; Topputo, F.; Gidea, M. Resonance transitions associated to weak capture in the restricted three-body problem. Adv. Space Res. 2008, 42, 1330–1351. [Google Scholar] [CrossRef]
  10. Khain, T.; Becker, J.; Adams, F. The Resonance Hopping Effect in the Neptune-planet Nine System. Publ. Astron. Soc. Pac. 2020, 132, 124401. [Google Scholar] [CrossRef]
  11. Malhotra, R.; Zhang, N. On the divergence of first-order resonance widths at low eccentricities. Mon. Not. R. Astron. Soc. 2020, 496, 3152–3160. [Google Scholar] [CrossRef]
  12. Malhotra, R. New results on orbital resonances. IAU Symp. 2022, 364, 85–101. [Google Scholar] [CrossRef] [PubMed]
  13. Antoniadou, K.I.; Libert, A.S. Bridges and gaps at low-eccentricity first-order resonances. Mon. Not. R. Astron. Soc. 2021, 506, 3010–3017. [Google Scholar] [CrossRef]
  14. Lo, M.; Parker, J. Unstable Resonant Orbits near Earth and Their Applications in Planetary Missions. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Providence, RI, USA, 16–19 August 2004. [Google Scholar] [CrossRef]
  15. Anderson, R.L. Tour Design Using Resonant-Orbit Invariant Manifolds in Patched Circular Restricted Three-Body Problems. J. Guid. Control Dyn. 2021, 44, 106–119. [Google Scholar] [CrossRef]
  16. Vaquero, M.; Howell, K.C. Leveraging Resonant-Orbit Manifolds to Design Transfers Between Libration-Point Orbits. J. Guid. Control Dyn. 2014, 37, 1143–1157. [Google Scholar] [CrossRef]
  17. Canales, D.; Gupta, M.; Park, B.; Howell, K.C. A transfer trajectory framework for the exploration of Phobos and Deimos leveraging resonant orbits. Acta Astronaut. 2022, 194, 263–276. [Google Scholar] [CrossRef]
  18. Lei, H.; Xu, B. Resonance transition periodic orbits in the circular restricted three-body problem. Astrophys. Space Sci. 2018, 363, 70. [Google Scholar] [CrossRef]
  19. Lei, H. Dynamical models for secular evolution of navigation satellites. Astrodynamics 2019, 4, 57–73. [Google Scholar] [CrossRef]
  20. McCarthy, B.; Howell, K. Leveraging quasi-periodic orbits for trajectory design in cislunar space. Astrodynamics 2021, 5, 139–165. [Google Scholar] [CrossRef]
  21. Lian, Y.; Huo, Z.; Cheng, Y. On the dynamics and control of the Sun-Earth L2 tetrahedral formation. Astrodynamics 2021, 5, 331–346. [Google Scholar] [CrossRef]
  22. Mccomas, D.; Carrico, J.; Hautamaki, B.; Intelisano, M.; Lebois, R.; Loucks, M.; Policastri, L.; Reno, M.; Scherrer, J.; Schwadron, N.; et al. A new class of long-term stable lunar resonance orbits: Space weather applications and the Interstellar Boundary Explorer. Space Weather 2011, 9, 11002. [Google Scholar] [CrossRef]
  23. Dichmann, D.; Lebois, R.; Carrico, J. Dynamics of Orbits Near 3:1 Resonance in the Earth-Moon System. J. Astronaut. Sci. 2014, 60, 51–86. [Google Scholar] [CrossRef]
  24. Gangestad, J.W.; Henning, G.A.; Persinger, R.R.; Ricker, G.R. A High Earth, Lunar Resonant Orbit for Lower Cost Space Science Missions. arXiv 2013, arXiv:1306.5333. [Google Scholar]
  25. Ricker, G.R. The Transiting Exoplanet Survey Satellite (TESS): Discovering New Earths and Super-Earths in the Solar Neighborhood. In AAS/Division for Extreme Solar Systems Abstracts; American Astronomical Society: Washington, DC, USA, 2015. [Google Scholar]
  26. Szebehely, V. Theory of Orbits: The Restricted Problem of Three Bodies. Am. J. Phys. 1968, 36, 375. [Google Scholar] [CrossRef]
  27. Lara, M.; Peláez, J. On the numerical continuation of periodic orbits. An intrinsic, 3-dimensional, differential, predictor-corrector algorithm. Astron. Astrophys. 2002, 389, 692–701. [Google Scholar] [CrossRef]
  28. Murray, C.D.; Dermott, S.F. Solar System Dynamics; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar] [CrossRef]
  29. Hadjidemetriou, J.D. Resonant Motion in the Restricted Three Body Problem. Celest. Mech. Dyn. Astron. 1993, 56, 201–219. [Google Scholar] [CrossRef]
  30. Hou, X.; Xin, X.; Feng, J. Genealogy and stability of periodic orbit families around uniformly rotating asteroids. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 93–114. [Google Scholar] [CrossRef]
  31. Hadjidemetriou, J.D.; Ichtiaroglou, S. A qualitative study of the Kirkwood gaps in the asteroids. Astron. Astrophys. 1984, 131, 20–32. [Google Scholar]
  32. Hadjidemetriou, J.D. Periodic Orbits. Celest. Mech. 1984, 34, 379–393. [Google Scholar] [CrossRef]
  33. Hou, X.H.; Liu, L. The Symmetric Horseshoe Periodic Families and the Lyapunov Planar Family Around L 3. Astron. J. 2008, 136, 67. [Google Scholar] [CrossRef]
  34. Colombo, G.; Franklin, F.A.; Munford, C.M. Saturn’s Rings. Astron. J. 1968, 73, 111. [Google Scholar] [CrossRef]
  35. Zhang, H.; Li, Y.; Zhang, K. A novel method of periodic orbit computation in circular restricted three-body problem. Sci. China E Technol. Sci. 2011, 54, 2197–2203. [Google Scholar] [CrossRef]
  36. Schanzle, A.F. Horseshoe-shaped orbits in the Jupiter-sun restricted problem. Astron. J. 1967, 72, 149. [Google Scholar] [CrossRef]
Figure 1. The planar Lyapunov families in the Sun–Earth system: (a) the Lyapunov family around L 1 ; (b) the Lyapunov family around L 2 .
Figure 1. The planar Lyapunov families in the Sun–Earth system: (a) the Lyapunov family around L 1 ; (b) the Lyapunov family around L 2 .
Applsci 12 08952 g001
Figure 2. Example resonant periodic orbits in the Sun–Earth system ( μ S E ): (a) 3:2 resonance; (b) 4:3 resonance; (c) 5:4 resonance; (d) 2:3 resonance; (e) 3:4 resonance; (f) 4:5 resonance.
Figure 2. Example resonant periodic orbits in the Sun–Earth system ( μ S E ): (a) 3:2 resonance; (b) 4:3 resonance; (c) 5:4 resonance; (d) 2:3 resonance; (e) 3:4 resonance; (f) 4:5 resonance.
Applsci 12 08952 g002
Figure 3. Characteristic curves of the circular families and the first-order resonant periodic families in the unperturbed model. The solid lines indicate the circular families and the dashed lines indicate the resonant periodic families. (a): Interior circular family and the interior first-order resonant periodic families (3:2, 4:3, and 5:4 resonances). (b): Exterior circular family and the exterior first-order resonant periodic families (2:3, 3:4, and 4:5 resonances). The abscissa is half the orbital period ( T / 2 ) in the synodic frame and the ordinate is the Jacobi constant C.
Figure 3. Characteristic curves of the circular families and the first-order resonant periodic families in the unperturbed model. The solid lines indicate the circular families and the dashed lines indicate the resonant periodic families. (a): Interior circular family and the interior first-order resonant periodic families (3:2, 4:3, and 5:4 resonances). (b): Exterior circular family and the exterior first-order resonant periodic families (2:3, 3:4, and 4:5 resonances). The abscissa is half the orbital period ( T / 2 ) in the synodic frame and the ordinate is the Jacobi constant C.
Applsci 12 08952 g003
Figure 4. Characteristic curves of the circular families and the first-order resonant periodic families in the CRTBP model of the Sun–Earth system ( μ S E ). (a): Interior circular family and the interior first-order resonant periodic families (3:2, 4:3, and 5:4 resonances). (b): Exterior circular family and the exterior first-order resonant periodic families (2:3, 3:4, and 4:5 resonances). The abscissa is half the orbital period ( T / 2 ) in the synodic frame and the ordinate is the Jacobi constant C.
Figure 4. Characteristic curves of the circular families and the first-order resonant periodic families in the CRTBP model of the Sun–Earth system ( μ S E ). (a): Interior circular family and the interior first-order resonant periodic families (3:2, 4:3, and 5:4 resonances). (b): Exterior circular family and the exterior first-order resonant periodic families (2:3, 3:4, and 4:5 resonances). The abscissa is half the orbital period ( T / 2 ) in the synodic frame and the ordinate is the Jacobi constant C.
Applsci 12 08952 g004
Figure 5. The family I B ( 3 4 ) (a) and the family O B ( 3 4 ) (b) in the Sun–Earth system ( μ S E ), shown in the T C plane. The abscissa is half the orbital period ( T / 2 ) in the synodic frame and the ordinate is the Jacobi constant C.
Figure 5. The family I B ( 3 4 ) (a) and the family O B ( 3 4 ) (b) in the Sun–Earth system ( μ S E ), shown in the T C plane. The abscissa is half the orbital period ( T / 2 ) in the synodic frame and the ordinate is the Jacobi constant C.
Applsci 12 08952 g005
Figure 6. Characteristic curves of the family I B ( 3 4 ) continued from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ). (a) Family I B ( 3 4 ) shown in the T C plane; (b) family I B ( 3 4 ) shown in the x 0 C plane. T / 2 is half the orbital period in the synodic frame. C is the Jacobi constant. x 0 is the initial state of the periodic orbits on the x-axis. Different colors indicate different mass parameters.
Figure 6. Characteristic curves of the family I B ( 3 4 ) continued from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ). (a) Family I B ( 3 4 ) shown in the T C plane; (b) family I B ( 3 4 ) shown in the x 0 C plane. T / 2 is half the orbital period in the synodic frame. C is the Jacobi constant. x 0 is the initial state of the periodic orbits on the x-axis. Different colors indicate different mass parameters.
Applsci 12 08952 g006
Figure 7. Characteristic curves of the family I B ( 3 4 ) and two RTPFs in the system for μ = 5 × 10 4 (a) and μ = μ S J (b).
Figure 7. Characteristic curves of the family I B ( 3 4 ) and two RTPFs in the system for μ = 5 × 10 4 (a) and μ = μ S J (b).
Applsci 12 08952 g007
Figure 8. Characteristic curves of the family I B ( 4 5 ) continued from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ). (a) Family I B ( 4 5 ) shown in the T C plane; (b) family I B ( 4 5 ) shown in the x 0 C plane. T / 2 is half the orbital period in the synodic frame. C is the Jacobi constant. x 0 is the initial state of the periodic orbits on the x-axis. Different colors indicate different mass parameters.
Figure 8. Characteristic curves of the family I B ( 4 5 ) continued from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ). (a) Family I B ( 4 5 ) shown in the T C plane; (b) family I B ( 4 5 ) shown in the x 0 C plane. T / 2 is half the orbital period in the synodic frame. C is the Jacobi constant. x 0 is the initial state of the periodic orbits on the x-axis. Different colors indicate different mass parameters.
Applsci 12 08952 g008
Figure 9. Characteristic curves of the family I B ( 4 5 ) and two RTPFs for the system with μ = 10 4 (a) and μ = μ S J (b).
Figure 9. Characteristic curves of the family I B ( 4 5 ) and two RTPFs for the system with μ = 10 4 (a) and μ = μ S J (b).
Applsci 12 08952 g009
Figure 10. (a) Characteristic curves of the family O B ( 3 4 ) continued from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ). Different colors indicate different mass parameters. (b) Characteristic curves of the family O B ( 3 4 ) and two RTPFs in the system with μ S J . T / 2 is half the orbital period in the synodic frame. C is the Jacobi constant. x 0 is the initial state of the periodic orbits on the x-axis.
Figure 10. (a) Characteristic curves of the family O B ( 3 4 ) continued from the Sun–Earth system ( μ S E ) to the Sun–Jupiter system ( μ S J ). Different colors indicate different mass parameters. (b) Characteristic curves of the family O B ( 3 4 ) and two RTPFs in the system with μ S J . T / 2 is half the orbital period in the synodic frame. C is the Jacobi constant. x 0 is the initial state of the periodic orbits on the x-axis.
Applsci 12 08952 g010
Figure 11. (a) Example orbit in RTPF-3:2(a); (b) Example orbit in RTPF-3:2(b); (c) Example orbit in RTPF-4:3(a); (d) Example orbit in RTPF-4:3(b); (e) Example orbit in RTPF-2:3(a); (f) Example orbit in RTPF-2:3(b). The initials of these example orbits are given in Table 2.
Figure 11. (a) Example orbit in RTPF-3:2(a); (b) Example orbit in RTPF-3:2(b); (c) Example orbit in RTPF-4:3(a); (d) Example orbit in RTPF-4:3(b); (e) Example orbit in RTPF-2:3(a); (f) Example orbit in RTPF-2:3(b). The initials of these example orbits are given in Table 2.
Applsci 12 08952 g011
Figure 12. Schematic diagram of generation of the initials of RTPOs in the Sun–Jupiter system.
Figure 12. Schematic diagram of generation of the initials of RTPOs in the Sun–Jupiter system.
Applsci 12 08952 g012
Figure 13. Curves of the velocity along the x-direction when the trajectory intersects with the x-axis for the first (a), second (b), third (c), and fourth (d) time. The abscissa is the initial velocity along the y-axis. The reference orbit is a planar Lyapunov orbit around the L 1 point with a Jacobi constant of C = 3.00391283019118 . The purple dots indicate a positive Δ v and the red points indicate a negative Δ v .
Figure 13. Curves of the velocity along the x-direction when the trajectory intersects with the x-axis for the first (a), second (b), third (c), and fourth (d) time. The abscissa is the initial velocity along the y-axis. The reference orbit is a planar Lyapunov orbit around the L 1 point with a Jacobi constant of C = 3.00391283019118 . The purple dots indicate a positive Δ v and the red points indicate a negative Δ v .
Applsci 12 08952 g013
Figure 14. Curves of the velocity along the x-direction when the trajectory intersects with the x-axis for the first (a), second (b), third (c), and fourth (d) time. The abscissa is the initial velocity along the y-axis. The reference orbit is a planar Lyapunov orbit around the L 2 point with a Jacobi constant of C = 2.99747681824366 . The purple dots indicate a positive Δ v and the red points indicate a negative Δ v .
Figure 14. Curves of the velocity along the x-direction when the trajectory intersects with the x-axis for the first (a), second (b), third (c), and fourth (d) time. The abscissa is the initial velocity along the y-axis. The reference orbit is a planar Lyapunov orbit around the L 2 point with a Jacobi constant of C = 2.99747681824366 . The purple dots indicate a positive Δ v and the red points indicate a negative Δ v .
Applsci 12 08952 g014
Figure 15. Example RTPOs passing through the L 1 region in the Sun–Jupiter system. (a) Orbit-1: initials generated by the second intersection of the trajectory with the x-axis (see Figure 13b). (b) Orbit-2 and (c) Orbit-3: initials generated by the third intersection of the trajectory with the x-axis (see Figure 13c). (d) Orbit-4, (e) Orbit-5 and (f) Orbit-6: initials generated by the fourth intersection of the trajectory with the x-axis (see Figure 13d). The larger and the smaller blue dots represent the positions of P 1 and P 2 , respectively.
Figure 15. Example RTPOs passing through the L 1 region in the Sun–Jupiter system. (a) Orbit-1: initials generated by the second intersection of the trajectory with the x-axis (see Figure 13b). (b) Orbit-2 and (c) Orbit-3: initials generated by the third intersection of the trajectory with the x-axis (see Figure 13c). (d) Orbit-4, (e) Orbit-5 and (f) Orbit-6: initials generated by the fourth intersection of the trajectory with the x-axis (see Figure 13d). The larger and the smaller blue dots represent the positions of P 1 and P 2 , respectively.
Applsci 12 08952 g015
Figure 16. Example RTPOs passing through the L 2 region in the Sun–Jupiter system. (a) Orbit-1: initials generated by the second intersection of the trajectory with the x-axis (see Figure 14b). (b) Orbit-2 and (c) Orbit-3: initials generated by the third intersection of the trajectory with the x-axis (see Figure 14c). (d) Orbit-4, (e) Orbit-5 and (f) Orbit-6: initials generated by the fourth intersection of the trajectory with the x-axis (see Figure 14d). The larger and the smaller blue dots represent the positions of P 1 and P 2 , respectively.
Figure 16. Example RTPOs passing through the L 2 region in the Sun–Jupiter system. (a) Orbit-1: initials generated by the second intersection of the trajectory with the x-axis (see Figure 14b). (b) Orbit-2 and (c) Orbit-3: initials generated by the third intersection of the trajectory with the x-axis (see Figure 14c). (d) Orbit-4, (e) Orbit-5 and (f) Orbit-6: initials generated by the fourth intersection of the trajectory with the x-axis (see Figure 14d). The larger and the smaller blue dots represent the positions of P 1 and P 2 , respectively.
Applsci 12 08952 g016aApplsci 12 08952 g016b
Table 1. The mass parameter values of the unperturbed two-body system ( μ 0 ), the Sun–Earth system ( μ S E ), and the Sun–Jupiter ( μ S J ) system.
Table 1. The mass parameter values of the unperturbed two-body system ( μ 0 ), the Sun–Earth system ( μ S E ), and the Sun–Jupiter ( μ S J ) system.
System/NameValues
Two-body system ( μ 0 )0
Sun–Earth system ( μ S E )3.003480575402412 × 10 6
Sun–Jupiter system ( μ S J )9.538811803631013 × 10 4
Table 2. Initial states of the example RTPOs shown in Figure 11.
Table 2. Initial states of the example RTPOs shown in Figure 11.
Orbit/Name x 0 y ˙ 0 T / 2 C
RTPF-3:2(a)-orbit ( μ = 5 × 10 4 )0.996693105698827−0.60672168269537013.5726320536319882.986678114083724
RTPF-3:2(b)-orbit ( μ = 5 × 10 4 )0.8398073560072940.32469898590257417.3717726920481652.985500211312612
RTPF-4:3(a)-orbit ( μ = 10 4 )0.999071125547079−0.49806795270147317.1474731118844692.992921442199951
RTPF-4:3(b)-orbit ( μ = 10 4 )0.8723996282744390.25329322122785021.2741748490527912.990627628838839
RTPF-2:3(a)-orbit ( μ = μ S J )1.0035482073430150.65698811218420816.1248462542868622.989316084875661
RTPF-2:3(b)-orbit ( μ = μ S J )1.240062798333267−0.43260417313864421.2237833594118042.969522450671411
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pan, S.; Hou, X. Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem. Appl. Sci. 2022, 12, 8952. https://doi.org/10.3390/app12188952

AMA Style

Pan S, Hou X. Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem. Applied Sciences. 2022; 12(18):8952. https://doi.org/10.3390/app12188952

Chicago/Turabian Style

Pan, Shanshan, and Xiyun Hou. 2022. "Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem" Applied Sciences 12, no. 18: 8952. https://doi.org/10.3390/app12188952

APA Style

Pan, S., & Hou, X. (2022). Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem. Applied Sciences, 12(18), 8952. https://doi.org/10.3390/app12188952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop