Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Near-optimal capture in the planar circular restricted Pluto-Charon system

Near-optimal capture in the planar circular restricted Pluto-Charon system In this paper, the near-optimal capture problem is numerically investigated to locate the optimal capture points around Charon using the Pluto-Charon planar circular restricted three-body problem (PCRTBP). Capture orbits are divided into the pre- and post-maneuver portions. In the pre-maneuver portion, the gravitational capture conditions are discussed by backward numerical integration. In the post-maneuver portion, the smaller alignment index (SALI) is applied to numerically investigate the orbital characters of long-term capture orbits around Charon. The initial conditions corresponding to three types of motion: (i) bounded (regular or chaotic), (ii) escaping and (iii) collisional are classi ed and studied. Combining results of the pre- and post-maneuver portions, the near-optimal capture method is presented to nd optimal capture points by overlay gures, and then the corresponding capture orbits are constructed in the PCRTBP. The di erent results between the Pluto-Charon system and the Earth-Moon system are compared and analyzed. Our results obtained in this paper, including the optimal capture points and the corresponding capture orbits, could be applied in future space mission design. Keywords: methods: numerical { celestial mechanics { planets and satellites: individual: Charon { planets and satellites: individual: Pluto 1. Introduction the WSB can be extended to interplanetary transfers, such as to Mercury (Jehn et al., 2004) and Mars (Topputo & Capture dynamics is a hot topic in the astronomy and Belbruno, 2015). astronautics. Using Kolmogorov-Arnold-Moser (KAM) tori, In addition, many scholars focus on the capture prob- Astakhov et al. (2003) proposed the chaos-assisted cap- lem of the near-Earth asteroids (NEAs), because the NEAs ture (CAC) mechanism in the circular restricted three- play important roles in the origin of the solar system and body problem (CRTBP) to investigate irregular moons the establishment of early warning of asteroid impact. Baoyin of Jupiter and Saturn, which identi es chaos as the es- et al. (2010) proposed a method to search for low-energy sential feature responsible for initial temporary gravita- NEAs that may be temporarily captured by the Earth. tional trapping within a planet's Hill sphere. Their the- NASA's Asteroid Redirect Robotic Mission (ARRM) was ory can explain why Saturn has more prograde irregular proposed to capture an entire 4  10 m NEA and redi- moons than Jupiter. In the Earth-Moon system, aster- rect it into an orbit around the Moon where astronauts oid 2006 RH was observed to be captured temporarily in the Orion spacecraft could explore it (Strange et al., from September 2006 to June 2007. Granvik et al. (2012) 2013). Yarnoz et al. (2013) found a family of so-called rst proposed a concept of temporary capture (TC) to Easily Retrievable Objects (EROs), which can be trans- describe and investigate this phenomenon in the Earth- ported from accessible heliocentric orbits into the Earth's Moon system. Based on their computation, at any given neighborhood at a ordable costs. Bao et al. (2015) ap- time there should be at least one temporarily-captured ir- plied lunar ybys and Earth ybys to capture NEAs into regular natural Earth satellite (NES) of 1-m diameter or- bounded orbits of the Earth. Based on the CAC mecha- biting the Earth. The well-known gravitational capture or nism in (Astakhov et al., 2003), Verrier & McInnes (2014) weak stability boundary (WSB) theory was rst proposed studied the low-energy capture of asteroids onto KAM tori by Belbruno (1987) to solve the lunar capture problem in in the Sun-Earth system. Qi & de Ruiter (2018) proposed Earth-to-Moon transfers. Then, this theory has been suc- a concept called short-term capture (STC) to investigate cessfully applied to many lunar missions, such as JAXA's temporary capture phenomena of NEAs in the EMS. Space Hiten rescue (Belbruno & Miller, 1990), ESA's SMART-1 and time conditions of STC were de ned and studied. Nu- (Schoenmaekers et al., 2001), and NASA's GRAIL (Ron- merical results indicated that there was a clear association coli & Fujii, 2016). The in-depth research indicates that between the distributions of the time probability of STC and its space conditions. Corresponding author The impulsive maneuver, V , provided by the propul- Email address: evzotos@physics.auth.gr (Euaggelos E. Zotos) Published in Planetary and Space Science April 9, 2019 arXiv:1904.03900v1 [astro-ph.EP] 8 Apr 2019 sion system of the spacecraft is the most practical ma- neuver method in the space mission to the present day. Hence, the optimal capture strategy usually can be evalu- ated by the magnitude of the impulsive maneuver. Qi & Xu (2014) proposed the near-optimal lunar capture theory in the Earth-Moon planar CRTBP to locate the optimal capture point and achieve the permanent lunar capture using the minimum V . This method combined the lu- nar gravitational capture and KAM tori around the Moon. The minimum capture velocity and the minimum capture eccentricity before the maneuver were de ned and anal- ysed. The permanent lunar capture after the maneuver was investigated by the distribution of KAM tori around the Moon. Furthermore, based on the near-optimal lunar capture, Qi & Xu (2016) constructed low-energy transfers from the Earth to permanent lunar orbits. The numerical results indicated that compared with traditional patched- Figure 1: The planar con guration of the circular restricted Pluto- conic transfers, those low-energy transfers can eciently Charon system. save fuel costs, which actually were very close to the the- oretical minimum. In this paper, we apply the near-optimal capture the- ing the pre-maneuver velocity V and increasing the post- ory to investigate the capture problem of Charon in the maneuver velocity V . Under given constraints of the cap- Pluto-Charon planar CRTBP. There are two main di er- ture point and Jacobi constant, there exist a minimum V ences between this paper and the previous research of the and a maximum V . Combining two aspects, the optimal Earth-Moon system in (Qi & Xu, 2014). Firstly, the mass capture point with the minimum V = V V can be 1 2 ratio  of the Pluto-Charon system is signi cantly larger found. than that of the Earth-Moon system, so the gravity of According to the above discussion, the structure of this Charon is more signi cant in its region than the case of the paper is separated into ve parts. Section 2 introduces the Moon. Consequently, we postulate that the near-optimal main properties of the dynamical system. Section 3 focuses capture in the Pluto-Charon planar CRTBP can show dif- on the gravitational capture conditions in pre-maneuver ferent characters and results. Secondly, Qi & Xu (2014) portion and obtain the minimum V . In Section 4, we applied KAM tori to analyze the long-term orbital char- investigate the long-term capture in post-maneuver potion acters after the impulsive maneuver V . Although this and obtain the maximum V . After that, the near-optimal method is simple and visible, the position and coverage of capture of Charon in the Pluto-Charon planar CRTBP is KAM tori are based on the empirical observation rather proposed and analyzed in Section 5. Finally, conclusions than the quantitative analysis. In this paper, the smaller of our investigation are presented in Section 6. alignment index (SALI), a common index for determining the regular or chaotic character of an orbit, is applied to in- 2. Description of the dynamical model vestigate and classify the long-term orbits around Charon more accurately instead of KAM tori. According to the theory of the classical circular re- In Zotos (2015) the orbital dynamics in the planar cir- stricted three-body problem (Szebehely, 1967), the two cular Pluto-Charon system were numerically investigated. primary bodies P (Pluto) and P (Charon) move on cir- 1 2 More speci cally, large sets of initial conditions of orbits cular orbits (Buie et al., 2012), with the same angular were classi ed in an attempt to determine the dynamics velocity, around their common center of mass. The third of a test particle (e.g., a spacecraft, or a comet, or an as- body, also known as the test particle (e.g., a spacecraft, teroid) moving in the vicinity of Charon. It was found or a comet, or an asteroid), moves on the same plane un- that the con guration (x; y) space is covered by a highly der the resultant Newtonian gravitational eld of the two complicated mixture of several types of basins, such as main bodies. Taking into account that the mass of the test basins of escape, basins of collision and basins of bounded particle is extremely small, compared to the masses of the motion. In addition, by varying the value of the Jacobi two primaries, we can reasonably assume that the third constant (or in other words the energy level) he managed body does not perturb the motion of the primaries. The to monitor the parametric evolution of all types of orbits. mass ratio  = m =(m + m ), with m > m de nes the 2 1 2 1 2 In this paper, the quantitative method similar to (Zotos, non-dimensional masses of the two primaries m = 1 2015) will be proposed to analyze the long-term orbital and m = . In this notation, the centers of Pluto and dynamics in the Pluto-Charon system. Based on Qi & Charon are located at (; 0) and (1 ; 0), respectively Xu (2014), the minimum V satisfying the permanent (see Fig. 1). capture can be obtained from two perspectives: decreas- 2 The planar motion of a test particle can be described by the following set of di erential equations @U U (x; y) = = x  2y _; @x @U U (x; y) = = y  + 2x; _ (1) @y where m m  (1 ) 1 1 2 2 2 U (x; y) = + + + x + y ; (2) r r 2 2 1 2 is the time-independent e ective potential function of the system, while r = (x + ) + y ; r = (x +  1) + y ; (3) are the distances to the respective primaries. Similarly, the set of the variational equations is given by (x) = x; _ (y) = y _; Figure 2: Schematic depicting the position of Charon (gray dot) along with the scattering region (cyan annulus) and Charon's sphere (z) = z _; of in uence (red circle). More details, regarding the radii of the 2 2 2 circular regions, are given in the text. The pre-maneuver orbit is @ U @ U @ U shown in dashed purple line, while the post-maneuver orbit is shown (x _ ) = x + y + z + 2y _; @x @x@y @x@z in solid green line. The position of the entrance and capture points 2 2 2 are pinpointed by ve-pointed black stars. (Color gure online). @ U @ U @ U (y _ ) = x + y + z 2x; _ @y@x @y @y@z 2 2 2 @ U @ U @ U velocity V . Similar to Qi & Xu (2014), we also assume (z _) = x + y + z: (4) that capture points are apsis points with respect to Charon, @z@x @z@y @z including periapsides and apoapsides. Since the velocity The set of the equations of motion (1) admits only the direction of the test particle at the capture point is perpen- following integral of motion dicular to the vector from the center of Charon to the test 2 2 J (x; y; x; _ y _ ) = 2U (x; y) x _ + y _ = C; (5) particle, the velocity direction has only two possibilities in the orbital plane: the prograde motion (anticlockwise which is also known as the Jacobi integral. The numeri- motion with respect to Charon) or the retrograde motion cal value C of the Hamiltonian remains constant and it is (clockwise motion with respect to Charon). called the Jacobi constant. The gravitational capture requires that the test parti- According to NASA Space Science Data Coordinated cle starting from the capture point can leave the sphere Archive (Update: 23 December 2016) the value of the mass of in uence (SOI) of Charon under a given ight time by ratio for the Pluto-Charon system is  = 0.10851122058 , backward integration. The radius of the SOI of Charon which is signi cantly larger than that of the Earth-Moon is set as 10000 km (see Fig. 2), about double the dis- planetary system, 0.0121506683 (Qi & Xu, 2014). The tance of L from the center of the smaller primary, which values of the Jacobi constant at the ve Lagrange points is similar to the lunar SOI proposed by Yamakawa et al. (L ; i = 1; :::; 5) are denoted by C , i = 1; :::; 5 and they are i i (1992). The ight time from the boundary of the SOI to equal to: C = 3.717080, C = 3.576079, C = 3.204728, 1 2 3 the capture point is restricted to 15 time units (approx- and C = C = 3. 4 5 imately 15.24 days). For the given position and velocity direction of the capture point, we postulate that a larger 3. Gravitational capture the velocity magnitude can make the escape from Charon (that is the Charon capture in the forward propagation) In this section, the gravitational capture in the pre- easier. Under the limitation of the ight time, we can use maneuver portion is investigated to decease the pre-maneuver the following method to nd the minimum V for the given capture point. It should be noted that in Zotos (2015), the value  = 0:099876695437731 was adopted. However, over the last years the 1. The initial point of backward integration is an apsis ocial value of the mass ratio regarding the Pluto-Charon planetary point with respect to Charon chosen from a scatter- system has been revised. 3 ing region. The scattering region of initial points is in two sides of the Moon; however, in the Pluto-Charon an annulus region around Charon, whose inner and PCRTBP, highly-fractal di used structures of C are max outer radii are equal to the radius of Charon and the more complicated in the Charon realm. distance of L from Charon, i.e., 606.0 and 5848.3 In the calculation of C , the integration time from 1 max km, respectively (see Fig. 2). Inside this annulus the initial point to the boundary of SOI of Charon is de- we de ne a uniform and dense grid of 1024  1024 noted by t . Fig. 4 shows the distributions of t in the esc esc initial conditions. There are two kinds of direction Pluto-Charon rotating frame. Similarly, left and right pan- of velocity: the prograde and retrograde motion, but els correspond to prograde and retrograde capture points, its velocity magnitude in the Pluto-Charon rotating respectively. By looking at the gure, it can be seen that in frame V can be changed. most region, orbits starting from prograde capture points are easier than those from the retrograde points to escape 2. According to Eq. (5), V = 2U (x; y) C , so the the Charon's SOI by backward integration. In addition, minimum V can translate to the maximum C . For we can nd some similarities between structures of t and the given capture point (x; y) and Jacobi constant esc C in Fig. 3. C , its velocity in the Pluto-Charon rotating frame, max According to Eq. (6), the minimum V we seek can be (x; _ y _ ), can be obtained from obtained by C . It should be noted that this minimum max velocity is a magnitude in the Pluto-Charon rotating frame V = 2U (x; y) C; rather than in the Charon-centered inertial frame. We de- (6) x _ = kV y=r ; 1 2 ne the corresponding minimum velocity in the Charon- y _ = kV (x +  1)=r ; 1 2 centered inertial frame as V . Since capture points are min apsides with respect to Charon, we can obtain the follow- where k denotes the direction of the velocity: k = 1 ing relationship, for the prograde motion with respect to Charon; k = 1 for the retrograde motion with respect to Charon. V = 2U (x; y) C + kr : (7) min max 2 r is the distance of the initial point to the center of Fig. 5 shows the distributions of V in the Pluto-Charon min Charon, and can be obtained from Eq. (3). In this dimensionless rotating frame. The left panel corresponds way, the initial state of the capture point (x; y; x; _ y _ ) to prograde capture points, while the right one corresponds in the rotating frame can be derived. to retrograde capture points. 3. We choose the initial C as C . Then, using the From Fig. 5, we observe that, generally, if the capture method in step 2, we can obtain the correspond- point is farther from Charon, the minimum capture veloc- ing initial state (x; y; x; _ y _ ). By backward integra- ity V is smaller. This result is accordant with the fact: min tion, a capture orbit can be obtained from the initial the attraction from Charon at a capture point farther to state. If the orbit is bounded in the SOI within 15 Charon is smaller, which results in a smaller escape veloc- time units or collides with Charon, we postulate that ity. As we can see from the gure, for both prograde and the given C cannot satisfy the gravitational capture. retrograde motions, minimum capture velocity V s are min Then we test whether the backward escape can be mainly distributed in the range of 0.28  2.57. However, achieved for a smaller C . In this way, C is decreas- in Fig. 3 C in the prograde orbits is obviously larger max ing until the backward escape can be satis ed. The than that in the retrograde orbits in most areas. We can corresponding maximum C is de ned as C . max explain this phenomenon by Eq. (7): for the same V at min It should be noted that the above numerical method is a given capture point, C of prograde motion is larger max similar to that in (Qi & Xu, 2014), but in step 3, we use than that of retrograde motion. C is derived from the max C as the variable rather than V in (Qi & Xu, 2014). The state of the rotating frame, so the in uence of the Pluto- advantage of C is that its initial value or upper boundary Charon rotation can in uence C . max is speci c, i.e. C . In the two-body model, we can use the orbital eccen- In Fig. 3 we display the distributions of C in the tricity e to determine whether the orbit is elliptic or hyper- max Pluto-Charon rotating frame, where left and right panels bolic. In the vicinity of a primary body of the PCRTBP, correspond to prograde and retrograde capture points, re- since the gravity of that celestial body is dominant, the spectively. From Fig. 3, we discover that C s of the concept of the eccentricity e can still be valid to study the max prograde motion are larger than those of the retrograde orbital character of the test particle. However, it should be noted that in the PCRTBP, e is time-varying. When e motion in most areas. The lower boundary of C of max is alternated from larger than 1 to smaller than 1, we pos- the retrograde motion is signi cantly smaller than that of tulate that gravitational capture happens. The minimum prograde motion. Compared with distributions of C max capture eccentricity e is de ned as the eccentricity at of Earth-Moon PCRTBP in (Qi & Xu, 2014), structures min the capture point corresponding to the minimum V , of C for retrograde motion are very similar. For the max min and can be obtained from prograde motion, the distribution of C in the Earth- max Moon PCRTBP (Qi & Xu, 2014) is more regular and con- r V min e = 1: (8) min secutive, and there exist two branches with smaller C max 4 Figure 3: Distributions of C in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde max motion. (Color gure online). Figure 4: Distributions of t in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde motion. esc (Color gure online). Figure 5: Distributions of V in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde min motion. (Color gure online). Fig. 6 illustrates distributions of the minimum cap- choose the capture points. As we can see from the gure, ture eccentricity e in the Pluto-Charon rotating frame. in the right gure, all retrograde capture points are avail- min Based on the de nition of the gravitational capture, cap- able, and the structures are similar to those of C . Qi max ture points with e > 1 are infeasible for the gravita- & Xu (2014) found that structures of e of the Earth- min min tional capture. We should avoid these regions when we Moon PCRTBP for prograde motion have two branches 5 Figure 6: Distributions of e in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde motion. min (Color gure online). in two sides of the Moon with e > 1. However, as we is covered either by collision or by regular initial condi- max can see from the left gure, the infeasible prograde capture tions of orbits. As long as C < C the channel around L 1 1 points for gravitational capture are only distributed in a opens and therefore orbits can escape to the Pluto realm. small region on the left edge. The stable areas of regular initial conditions are destroyed abruptly. It should be noted that when C = C two new stability islands emerge at both sides of Charon. We ob- 4. Long-term capture serve that as the value of the Jacobi constant decreases the amount of both regular and collision orbits decrease. In this section, we focus on the maximum post-maneuver On the other hand, the amount of escaping orbits grows velocity V satisfying the long-term capture around Charon. rapidly and for C  C they dominate the entire con gura- According to (Qi & Xu, 2014), a regular orbit in KAM tion plane. In the Earth-Moon PCRTBP (Qi & Xu, 2014), tori around Charon can be bounded in the Charon realm KAM tori rapidly disappear when C decreases to C , but permanently in the Pluto-Charon PCRTBP, otherwise, an apparently, regular orbits in the Pluto-Charon system can orbit located in chaotic areas escapes the Charon realm. still survive when C = C . We postulate that this di er- Therefore, we require that a test particle after V should ence is derived from the larger  of Pluto-Charon system, be inserted into a regular orbit around Charon. In (Qi & because orbits rotating around the smaller primary with a Xu, 2014), KAM tori were applied to nd regular bounded larger  have a greater resistance against the perturbation regions. Although this method is simple and visible, the from the larger primary in the PCRTBP. position and coverage of KAM tori rely on the empiri- When C < C , the abrupt emergence of chaotic regions cal observation rather than the quantitative analysis. The in Fig. 7 looks quite similar to the abrupt bifurcation to SALI is a common index for quantitatively determining fully developed chaotic scattering in (Bleher et al., 1990). the regular or chaotic character of an orbit (Skokos, 2001). However, we postulate that the chaotic phenomenon in The analytical de nition of SALI, as well as the particular Fig. 7 essentially is di erent from the chaotic scattering. threshold values for distinguishing between ordered and Firstly, the fully developed chaotic scattering denotes a sit- chaotic motion are given in Zotos & Jung (2018). The uation in which all periodic orbits are unstable and there classi cations of the long-term orbits around Charon can are no KAM surfaces (i.e., the dynamics is hyperbolic) be obtained by the value of the SALI quantitatively, in- (Bleher et al., 1989), but in Fig. 7, obviously, there ex- cluding (i) bounded (regular or chaotic), (ii) escaping and ist KAM tori in the bounded regular region (denoted by (iii) collision (Zotos, 2015). Essentially, other diagnostic green) when C is smaller than C (Qi & Xu, 2014). In indexes for the orbital chaos, such as the the fast Lyapunov addition, in the CRTBP we consider, not all periodic or- indicator (FLI) (Astakhov & Farrelly, 2004), the correla- bits are unstable, for example the distant retrograde orbits tion dimension (Belbruno et al., 2008), and the frequency around the smaller primary are stable (Henon, 1969; Lam analysis method (Laskar, 1990), can also be applied to & Whi en, 2005). Secondly, Koon et al. (2000) pointed distinguish stable regular orbits from chaotic ones. In this out that the chaotic phenomenon in the region around the section, we use the SALI to investigate bounded basins of smaller primary of the CRTBP is ascribed to the existence the long-term orbits around Charon. of transversal homoclinic and heteroclinic points. We can In Fig. 7 we present the orbital structure of the con g- explain these chaotic phenomena by the symbolic dynam- uration (x; y) plane for the prograde motion. It is seen that ics or the horseshoe-like dynamics (Koon et al., 2000). when C > C all the energetically allowed area of motion It should be noted that trapped chaotic motion is very 6 Figure 7: Basin diagrams on the con guration (x; y) plane for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:68; (d): C = C ; (e): C = C ; (f ): C = C . The color code is as follows: non-escaping regular orbits (green); trapped chaotic orbits (yellow); escaping 2 3 4 orbits (cyan); collision orbits (red); energetically forbidden regions (gray). (Color gure online). limited. This means that the corresponding initial condi- orbital structure of the con guration (x; y) plane for a tions appear only as lonely and isolated points, mainly in xed value of the Jacobi constant (or equivalently, for a the vicinity of the stability islands of regular motion. In- xed orbital energy), for both prograde and retrograde deed, in panel (b) of Fig. 7 one can distinguish several near motion. In an attempt to obtain the same information the boundaries of the bounded basins of regular motion the for a continuous interval of C values, we de ne a plane of presence of several initial conditions which correspond to representation in which the x-coordinate is the abscissa, trapped chaotic motion. while the ordinate is given by the Jacobi constant. This In the same vein, the orbital structure of the con gu- means that all orbits start with initial position (x ; 0), i.e, ration plane for the retrograde motion is depicted in Fig. on the x-axis, with initial velocities (0; y _ ), where the y- 8. It is evident that the orbital structure corresponding to component of the velocity is derived from the third of Eqs. the retrograde motion has many di erences with respect (6). to that of the prograde motion. The most notable ones In Fig. 11 we present the basin diagram on the (x; C ) are: (i) bounded regular motion is possible for much lower plane for C 2 [2:5; 5], for both the prograde and retrograde values of the Jacobi constant, or alternatively, for higher motions. The black continuous line indicates the corre- values of the total orbital energy; (ii) the overall orbital sponding Zero Velocity Curve (ZVC), and it is de ned as structure seems less complicated, without the presence of regions with highly fractal basin boundaries. Those results f (x; C ) = 2U (x; y = 0) = C: (9) are similar to those of the Earth-Moon PCRTBP in(Qi & Xu, 2014). We see that for the prograde motion bounded regular or- The distributions of the corresponding escape and col- bits exist only when C > C , while for the retrograde lision time of the orbits of both the prograde and retro- motion the same type of orbits is still possible even below grade motion are given in Figs. 9 and 10, respectively. As C = C . However, for extremely high and low values of expected, the lowest values of both escape and collision C we encounter the same orbital structure for both types time were measured inside the corresponding basins of es- of motion (prograde and retrograde). In particular, for cape and collision, while the highest times were identi ed extremely low values of the Jacobi constant escaping or- mainly in the fractal basin boundaries, where the degree bits completely dominate, while for relatively high values of the unpredictability of the system maximizes. of C bounded and collision types of orbits share the ener- The color-coded diagrams of Figs. 7 and 8 reveal the getically allowed phase space. In Fig. 12 we present the 7 Figure 8: Basin diagrams on the con guration (x; y) plane for the retrograde motion, where (a): C = 3:9; (b): C = C ; (c): C = C ; (d): 2 4 C = 2:9; (e): C = 2:8; (f ): C = 2:5. The color code is as follows: non-escaping regular orbits (green); trapped chaotic orbits (yellow); escaping orbits (cyan); collision orbits (red); energetically forbidden regions (gray). (Color gure online). Figure 9: Distribution of the escape and collision times of the orbits for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:68; (d): C = C ; (e): C = C ; (f ): C = C . All initial conditions of bounded (regular and chaotic) orbits are shown in white. (Color 2 3 4 gure online). distribution of times employed by the test particle in the and retrograde motions, bounded orbits generally erode escape and collision orbits. or disappear with the decrease of C (see Figs 7 and 8). Based on the previous discussion, for both prograde Hence, for a given capture point (x; y) and its direction 8 Figure 10: Distribution of the escape and collision times of the orbits for the retrograde motion, where (a): C = 3:9; (b): C = C ; (c): C = C ; (d): C = 2:9; (e): C = 2:8; (f ): C = 2:5. All initial conditions of bounded (regular and chaotic) orbits are shown in white. (Color gure online). Figure 11: Basin diagrams on the (x; C)-plane for (a-left): the prograde motion and (b-right): the retrograde motion. The color code is as follows: non-escaping regular orbits (green); trapped chaotic orbits (yellow); escaping orbits (cyan); collision orbits (red); energetically forbidden regions (gray). (Color gure online). of motion, there exists a minimum C satisfying long-term yse the minimum V . According to Qi & Xu (2014), the capture in the Charon region. According to Eq. (5), V key to nd the minimum V is to locate the near-optimal corresponding to the minimum C can be obtained, which capture point. Combining results obtained in previous sec- is the maximum V we seek for. tions, near-optimal capture points can be found by the following method: 1. For the given velocity direction, V corresponding 5. Near-optimal capture of Charon to the capture point (x; y) can be calculated by In the previous two sections, we have investigated how V = jV V j 1 2 to obtain the minimum V and the maximum V , succes- 1 2 p p (10) sively. In this section, we synthesize those results to anal- = j 2U (x; y) C 2U (x; y) Cj max 9 Figure 12: Distribution of the corresponding escape and collision times of the orbits for (a-left): the prograde motion and (b-right): the retrograde motion. All initial conditions of bounded (regular and chaotic) orbits are shown in white. (Color gure online). where C is obtained from Figs. 3 or 4, and C is 14 is similar to that of C for retrograde motion in Fig. max max the Jacobi constant of the post-maneuver portion. It 3 (b). should be noted that if the given capture point (x; y) Figs. 15 and 16 show the overlays of V and the is located in the energetically forbidden regions for bounded basins of di erent C for the prograde and ret- the given C , Eq. (10) is invalid and the correspond- rograde motions, respectively. Gray and cyan regions in ing V does not exist. these gures denote energetically forbidden regions and bounded basins, respectively. Dark blue points are near- 2. Using Eq. (10), the distribution of V in the scat- optimal capture points we chose. Based on the analysis of tering region for given velocity direction and C can Fig. 7, only for a larger C can bounded basins of the pro- be obtained. For example, Figs. 13 and 14 display grade motion exist. Hence, near-optimal prograde capture the distributions of V in the Pluto-Charon rotat- points will disappear with the decrease of C . But near- ing frame corresponding to the cases presented in optimal retrograde capture points can survive in smaller Figs. 7 and 8, respectively. In addition, for the same C s. By looking at these gures, it can be seen that since velocity direction and C , the corresponding basin di- distributions of V and bounded basins of prograde mo- agrams can be obtained using the method in Section tion are both quite di erent from those of retrograde mo- 4 (see Figs. 7 and 8). Then, we put the bounded tion, the distributions of the corresponding near-optimal basins onto the distribution of V , and get an over- capture points have signi cant di erences between two lay gure of them. For example, Figs. 15 and 16 kinds of motions. For example, in Fig. 16, near-optimal display overlays for di erent C s and motions. The retrograde capture points are mostly located in the inner background of the overlay is the distribution of V , boundary of bounded basins near Charon, and similar to and cyan regions denote bounded basins. the distribution of bounded basins, near-optimal capture 3. From the overlay gures, points in bounded basins points gradually move away from Charon with the decrease denote feasible points in terms of the long-term cap- of C . However, situations of near-optimal prograde cap- ture after V . Among these feasible capture points, ture points are more complicated according to Fig. 15, those located in the minimum V regions are cho- where regions with highly fractal basin boundaries exist, sen as near-optimal capture points. For example, especially when C = 3:71. In Fig. 15 (a), near-optimal blue points in Figs. 15 and 16 are near-optimal cap- capture points are distributed in the inner boundary of ture points we chose. The corresponding V is the bounded basin nearest Charon. But when C decreases to minimum V we seek for. 3.71, i.e., necks around L just open, near-optimal capture Figs. 13 and 14 display the distributions of V cor- points are distributed in a highly fractal pattern around responding to Figs. 7 and 8, respectively. Gary regions Charon, which is similar to the structures near Charon in in those gures are the energetically forbidden regions in Figs. 3 and 13. When C decreases to 3.60, corresponding the Charon realm for larger C s. Compared Fig. 13 with bounded basins are distributed in two separated stable is- Fig. 3 (a), the structure of the distribution of V for the lands, we can obtain two near-optimal capture regions in prograde motion actually is quite similar to that of C max two separated bounded basins, respectively. But when C for prograde motion. It can be explained by Eq. (10). decreases to C , near-optimal capture points can only be Similarly, the structure of the distribution of V in Fig. distributed in a small region of the left island. Numerical 10 Figure 13: Distributions of V for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:68; (d): C = C ; (e): C = C ; (f ): 2 3 C = C . (Color gure online). Figure 14: Distributions of V for the retrograde motion, where (a): C = 3:9; (b): C = C ; (c): C = C ; (d): C = 2:9; (e): C = 2:8; (f ): 2 4 C = 2:5. (Color gure online). computation indicates that minimum V s corresponding maneuver portions of capture orbits (or long-term capture to near-optimal capture points in Figs. 15 and 16 are all orbits). Red star points are the optimal capture points (or smaller than 23 m/s. impulsive points) chosen from the near-optimal capture re- Using the near-optimal capture points as the initial gions in Figs. 15 and 16 correspondingly. By looking at state, the corresponding pre- and post-maneuver capture the gures, we observe that both pre- and post-maneuver orbits can be constructed. In Figs. 17 and 18, we illus- portions of capture orbits have signi cant di erences be- trate the capture orbits corresponding to the near-optimal tween prograde and retrograde motions. In Fig. 17, since capture points in Figs. 15 and 16, respectively. The g- the corresponding C s of prograde motion are larger than ures are depicted in the dimensionless Pluto-Charon ro- C , the sizes of post-maneuver orbits are smaller than tating frame. Red dashed lines denote pre-maneuver por- those of retrograde motion in Fig. 18. In addition, all tions of capture orbits, while green solid lines denote post- pre-maneuver portions of prograde near-optimal capture 11 Figure 15: Overlays of V and bounded regular basins, for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:60; (d): C = C . All bounded basins are indicated with cyan color, while the initial conditions of both escape and collision orbits are shown in tones of solar colors, according to the value of V . All near-optimal capture initial conditions are pinpointed by blue dots. (Color gure online). Figure 16: Overlays of V and bounded regular basins, for retrograde motion, where (a): C = C ; (b): C = C ; (c): C = 2:9; (d): C = 2:8. 2 4 All bounded basins are indicated with cyan color, while the initial conditions of both escape and collision orbits are shown in tones of solar colors, according to the value of V . All near-optimal capture initial conditions are pinpointed by blue dots. (Color gure online). 12 Figure 17: Capture obits corresponding to the prograde near-optimal capture points, where the Jacobi constants of the post-maneuver orbit are (a): C = 3:9; (b): C = 3:71; (c): C = 3:60; (d): C = C . The pre-maneuver orbits are shown in dashed, red dashed lines, while the post-maneuver orbits are shown in solid green lines. The black solid lines indicate the boundaries of the scattering region, while the blue solid line delimits the SOI of Charon. The position of the Lagrange points L and L are pinpointed by purple dots. (Color gure online). 1 2 Table 1: Data of prograde near-optimal capture orbits. Table 2: Data of retrograde near-optimal capture orbits. Fig. 17 (x; y) C C V (m/s) Fig. 18 (x; y) C C V (m/s) max max (a) (0:92575062;0:00083565) 3.687 3.90 10.4064 (a) (0:90736622;0:03091922) 3.097 C 22.6293 (b) (0:79538850; 0:00083565) 3.697 3.71 1.2940 (b) (0:89733836; 0:07103064) 2.947 C 3.5683 (c) (0:99093168; 0:15961003) 3.627 3.60 6.0325 (c) (0:88062527;0:09442897) 2.867 2.90 2.5561 (d) (0:82547207;0:00083565) 3.677 C 7.3500 (d) (0:80207374; 0:27158774) 2.777 2.80 3.1574 points pass through the gateway of L . However, in Fig. interplanetary space probe). Although in real space mis- 18, most pre-maneuver portions of retrograde near-optimal sions there are several types of plane perturbations mainly capture points pass through the gateway of L . As we due to the presence of other giant celestial bodies, such as can see Fig. 18, multiple Charon ybys occur during the Jupiter and Saturn, all these perturbations are extremely pre-maneuver portions. Data of prograde and retrograde weak in the vicinity of Charon (which was the scattering near-optimal capture orbits in Figs. 17 and 18 are listed region in our work). Therefore, we argue that the numer- in Tables 1 and 2. From these tables, V s in Fig. 17 are ical results of the Pluto-Charon CRTBP obtained in this smaller than 11 m/s, while V s in Fig. 18 are smaller paper are structurally stable against out of plane pertur- than 23 m/s. Of particular note is that based on the data bations. of Fig. 17 (c) and (d), C s of pre-maneuver portions max are larger than the corresponding C s of post-maneuver 6. Conclusions portions, so the impulsive maneuvers at the near-optimal capture points increase the speed instead of brake. In this paper, the near-optimal capture problem in the The above results we obtained, including optimal cap- Pluto-Charon system was investigated and discussed. In ture points and corresponding capture orbits, could be ap- an attempt to achieve long-term capture around Charon plied in future space mission design (like the New Horizons with the minimum V , optimal capture points were found 13 Figure 18: Capture obits corresponding to the retrograde near-optimal capture points, where the Jacobi constants of the post-maneuver orbit are (a): C = C ; (b): C = C ; (c): C = 2:9; (d): C = 2:8. The pre-maneuver orbits are shown in dashed, red dashed lines, while the 2 4 post-maneuver orbits are shown in solid green lines. The black solid lines indicate the boundaries of the scattering region, while the blue solid line delimits the SOI of Charon. The position of the Lagrange points L and L are pinpointed by purple dots. (Color gure online). 1 2 and located by the near-optimal capture theory. According & Xu, 2014), the structures of C for prograde mo- max to the de nition of V , the capture orbit was divided into tion in the Pluto-Charon PCRTBP are more compli- the pre- and post-maneuver portions. In the pre-maneuver cated in the Charon realm, where highly-fractal dif- portion, the maximum escape Jacobi constant C was fused structures of C exist. Based on the results max max obtained and discussed by backward numerical integra- of t , in most region, orbits starting from prograde esc tion. In the post-maneuver portion, the SALI was applied capture points are easier than those from retrograde to numerically investigate the orbital characters of long- points to escape the Charon's SOI by backward in- term capture orbits around Charon. The initial condi- tegration. tions corresponding to three types of motion: (1) bounded, 2. Based on the de nition of gravitational capture and (2) escaping and (3) collisional were classi ed and distin- distributions of e , we found that all retrograde min guished. Combining the results of pre- and post-maneuver capture points are available for gravitational cap- portions, optimal capture points were located by overlay ture, and its structure of e is similar to that of min gures. Furthermore, corresponding capture orbits were C . Di erent from the two branches in two sides max constructed in the PCRTBP. In addition, the near-optimal of the Moon for prograde motion in the Earth-Moon capture results in the Pluto-Charon system of this paper PCRTBP (Qi & Xu, 2014), the prograde infeasible were compared with those in the Earth-Moon system. We capture points for gravitational capture in the Pluto- postulated that the signi cant di erence of mass ratio Charon PCRTBP are only distributed in a small re- resulted in the di erent results between two planet-moon gion on the left edge. systems. 3. According to the orbital structure of long-term or- The main results of our numerical research can be sum- bits around Charon, we found that the orbital struc- marized as follows: ture corresponding to the retrograde motion has ob- 1. Numerical computation indicated that the C s of vious di erences with respect to that of the prograde max the prograde motion are larger than those of the ret- motion: (1) bounded regular motion is possible for rograde motion in most areas. Compared with distri- much lower values of the Jacobi constant, or alterna- butions of C of the Earth-Moon PCRTBP in (Qi tively, for higher values of the total orbital energy; max 14 (2) the overall orbital structure seems less compli- Astakhov, S.A., Farrelly, D., 2004. Capture and escape in the ellip- tic restricted three-body problem. Monthly Notices of the Royal cated, without the presence of regions with highly Astronomical Society 354, 971-979. fractal basin boundaries. In addition, the lowest val- Bao, C., Yang, H., Barsbold, B., Baoyin, H., 2015. Capturing near- ues of both escape and collision time were measured Earth asteroids into bounded Earth orbits using gravity assist. Astrophysics and Space Science 360, 61. inside the corresponding basins of escape and colli- Baoyin, H-.X., Chen, Y., Li, J.-F., 2010. Capturing near earth ob- sion, while the highest times were identi ed mainly jects. Research in Astronomy and Astrophysics 10, 587-598. in the fractal basin boundaries, where the degree of Belbruno, E.A., AIAA, DGLR, and JSASS, International Electric the unpredictability of the system maximizes. Propulsion Conference, 19th, Colorado Springs, CO, May 11-13, 4. The distributions of near-optimal capture points in- Belbruno, E.A., Miller, J., Jet Propulsion Lab., Pasadena, CA, 1990. dicated that near-optimal prograde capture points Belbruno, E., Topputo, F., Gidea, M., 2008. Resonance transitions disappear rapidly with the decrease of C , but near- associated to weak capture in the restricted three-body problem. optimal retrograde capture points can survive in smaller Advances in Space Research 42, 1330-1351. Bleher, S., Grebogi, C., Ott, E., 1990. Bifurcation to chaotic scat- C s. Near-optimal retrograde capture points are mostly tering. Physica D: Nonlinear Phenomena 46(1), 87-121. located in the inner boundary of bounded basins near Bleher, S., Ott, E., Grebogi, C., 1989. Routes to chaotic scattering. Charon, and gradually move away from Charon with Physical review letters 63(9), 919. the decrease of C . However, situations of prograde Buie, M.W., Tholen, D.J., Grundy, W.M., 2012. The Orbit of Charon is Circular. AJ 144, article id. 15. motion are more complicated, where highly fractal Granvik, M., Vaubaillon, J., Jedicke, R., 2012. The population of patterns around Charon and separated stable islands natural Earth satellites. Icarus 218, 262-277. exist. Hnon, M., 1969. Numerical exploration of the restricted problem, V. Astronomy and Astrophysics 1, 223-238. 5. Numerical results indicated that all pre-maneuver Jehn, R., Campagnola, S., Garcia, D., Kemble, S., 18th International portions of prograde near-optimal capture points pass Symposium on Space Flight Dynamics, 2004. through the gateway of L . However, most pre- 1 Koon W S, Lo M W, Marsden J E, et al., Dynamical systems, the maneuver portions of retrograde near-optimal cap- three-body problem and space mission design, World Scienti c, ture points pass through the gateway of L . Data of Lam, T., Whi en G J., AAS/AIAA Space ight Mechanics Confer- both prograde and retrograde near-optimal capture ence, 2005, paper AAS 05-110. orbits showed that, V s are smaller than 23 m/s. Laskar, J., 1990. The chaotic motion of the solar system: A numerical estimate of the size of the chaotic zones. Icarus 88, 266-291. For the numerical integration (both backward and for- NASA Space Science Data Coordinated Archive: http://nssdc. ward in time) of the equations of motion (1), as well as gsfc.nasa.gov/planetary/factsheet/plutofact.html Press H.P., Teukolsky S.A, Vetterling W.T., Flannery B.P., Numer- of the variational equations (4) we used a double precision ical Recipes in FORTRAN 77, 2nd Ed., Cambridge Univ. Press, Bulirsch-Stoer FORTRAN 77 algorithm (Press et al., 1992), Cambridge, USA, 1992. with a xed time step equal to 10 . In all our calcula- Qi, Y., de Ruiter, A., 2018. Short-term capture of the Earth-Moon system. Monthly Notices of the Royal Astronomical Society 476 tions the numerical error, related to the conservation of (4), 5464-5478. the Jacobi integral (5) was generally smaller than 10 , Qi, Y., Xu, S., 2014. Lunar capture in the planar restricted three- while in most of the cases it was smaller than 10 . All body problem. Celestial Mechanics and Dynamical Astronomy the graphical illustration has been created using the lat- 120, 401-422. Qi, Y., Xu, S., 2016. Earth{Moon transfer with near-optimal lunar est version 11.2 of the software Mathematica (Wolfram, capture in the restricted four-body problem. Aerospace Science 2003). and Technology 55, 282-291. The obtained results in this paper, including optimal Roncoli, R.B., Fujii, K.K., AIAA Meeting Papers on Disc, 2010. capture points and corresponding capture orbits, could be Schoenmaekers, J., Horas, D., Pulido, J.A., 16th International Sym- posium on Space Flight Dynamics, Pasadena, California, 2001. applied in future space mission design. Skokos C., 2001. Alignment indices: a new, simple method for deter- mining the ordered or chaotic nature of orbits.Journal of Physics A 34, 10029-10043. Acknowledgments Strange, N., Landau, D., McElrath, T., Lantoine, G., Lam, T., Overview of mission design for NASA asteroid redirect robotic Yi Qi wishes to acknowledge the support of the Natu- mission concept, Pasadena, CA: Jet Propulsion Laboratory, Na- ral Sciences and Engineering Research Council of Canada tional Aeronautics and Space Administration, 2013. through a Discovery Accelerator Supplement under Grant Szebehely V., Theory of Orbits, Academic Press, New York, 1967. Topputo, F., Belbruno, E., 2015. Earth{Mars transfers with ballistic RGPAS-493042-2016. The authors would like to express capture. Celestial Mechanics and Dynamical Astronomy 121, 329- their warmest thanks to the two anonymous reviewers for the careful reading of the manuscript as well as for all the Verrier, P.E., McInnes, C.R., 2014. Low-Energy Capture of Aster- oids onto Kolmogorov{Arnold{Moser Tori. Journal of Guidance, valuable suggestions and comments which allowed us to Control, and Dynamics 38, 330-335. improve both the quality and the clarity of the paper. Wolfram S., The Mathematica Book, Wolfram Media, Champaign, Yamakawa, H., Kawaguchi, J., Ishii, N., Matsuo, H., A numerical References study of gravitational capture orbit in the Earth-Moon system, Space ight mechanics, 1992, pp. 1113{1132. Astakhov, S.A., Burbanks, A.D., Wiggins, S., Farrelly, D., 2003. Yarnoz, D.G., Sanchez, J.P., McInnes C.R., 2013. Easily retrievable Chaos-assisted capture of irregular moons. Nature 423, 264-267. 15 objects among the NEO population. Celestial Mechanics and Dy- namical Astronomy 116, 367-388. Zotos, E.E., 2015. Orbit classi cation in the planar circular Pluto- Charon system. Astrophysics and Space Science 360, article id. Zotos, E.E., Jung, Ch., 2018. Correlating the escape dynamics and the role of the normally hyperbolic invariant manifolds in a binary system of dwarf spheroidal galaxies. International Journal of Non- Linear Mechanics 99, 182-203. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Astrophysics arXiv (Cornell University)

Near-optimal capture in the planar circular restricted Pluto-Charon system

Zotos, Euaggelos E.; Qi, Yi
Astrophysics , Volume 2019 (1904) – Apr 8, 2019
Download PDF
16 pages

Loading next page...
 
/lp/arxiv-cornell-university/near-optimal-capture-in-the-planar-circular-restricted-pluto-charon-V292siHduB

References (29)

ISSN
0032-0633
eISSN
ARCH-3330
DOI
10.1016/j.pss.2018.11.009
Publisher site
See Article on Publisher Site

Abstract

In this paper, the near-optimal capture problem is numerically investigated to locate the optimal capture points around Charon using the Pluto-Charon planar circular restricted three-body problem (PCRTBP). Capture orbits are divided into the pre- and post-maneuver portions. In the pre-maneuver portion, the gravitational capture conditions are discussed by backward numerical integration. In the post-maneuver portion, the smaller alignment index (SALI) is applied to numerically investigate the orbital characters of long-term capture orbits around Charon. The initial conditions corresponding to three types of motion: (i) bounded (regular or chaotic), (ii) escaping and (iii) collisional are classi ed and studied. Combining results of the pre- and post-maneuver portions, the near-optimal capture method is presented to nd optimal capture points by overlay gures, and then the corresponding capture orbits are constructed in the PCRTBP. The di erent results between the Pluto-Charon system and the Earth-Moon system are compared and analyzed. Our results obtained in this paper, including the optimal capture points and the corresponding capture orbits, could be applied in future space mission design. Keywords: methods: numerical { celestial mechanics { planets and satellites: individual: Charon { planets and satellites: individual: Pluto 1. Introduction the WSB can be extended to interplanetary transfers, such as to Mercury (Jehn et al., 2004) and Mars (Topputo & Capture dynamics is a hot topic in the astronomy and Belbruno, 2015). astronautics. Using Kolmogorov-Arnold-Moser (KAM) tori, In addition, many scholars focus on the capture prob- Astakhov et al. (2003) proposed the chaos-assisted cap- lem of the near-Earth asteroids (NEAs), because the NEAs ture (CAC) mechanism in the circular restricted three- play important roles in the origin of the solar system and body problem (CRTBP) to investigate irregular moons the establishment of early warning of asteroid impact. Baoyin of Jupiter and Saturn, which identi es chaos as the es- et al. (2010) proposed a method to search for low-energy sential feature responsible for initial temporary gravita- NEAs that may be temporarily captured by the Earth. tional trapping within a planet's Hill sphere. Their the- NASA's Asteroid Redirect Robotic Mission (ARRM) was ory can explain why Saturn has more prograde irregular proposed to capture an entire 4  10 m NEA and redi- moons than Jupiter. In the Earth-Moon system, aster- rect it into an orbit around the Moon where astronauts oid 2006 RH was observed to be captured temporarily in the Orion spacecraft could explore it (Strange et al., from September 2006 to June 2007. Granvik et al. (2012) 2013). Yarnoz et al. (2013) found a family of so-called rst proposed a concept of temporary capture (TC) to Easily Retrievable Objects (EROs), which can be trans- describe and investigate this phenomenon in the Earth- ported from accessible heliocentric orbits into the Earth's Moon system. Based on their computation, at any given neighborhood at a ordable costs. Bao et al. (2015) ap- time there should be at least one temporarily-captured ir- plied lunar ybys and Earth ybys to capture NEAs into regular natural Earth satellite (NES) of 1-m diameter or- bounded orbits of the Earth. Based on the CAC mecha- biting the Earth. The well-known gravitational capture or nism in (Astakhov et al., 2003), Verrier & McInnes (2014) weak stability boundary (WSB) theory was rst proposed studied the low-energy capture of asteroids onto KAM tori by Belbruno (1987) to solve the lunar capture problem in in the Sun-Earth system. Qi & de Ruiter (2018) proposed Earth-to-Moon transfers. Then, this theory has been suc- a concept called short-term capture (STC) to investigate cessfully applied to many lunar missions, such as JAXA's temporary capture phenomena of NEAs in the EMS. Space Hiten rescue (Belbruno & Miller, 1990), ESA's SMART-1 and time conditions of STC were de ned and studied. Nu- (Schoenmaekers et al., 2001), and NASA's GRAIL (Ron- merical results indicated that there was a clear association coli & Fujii, 2016). The in-depth research indicates that between the distributions of the time probability of STC and its space conditions. Corresponding author The impulsive maneuver, V , provided by the propul- Email address: evzotos@physics.auth.gr (Euaggelos E. Zotos) Published in Planetary and Space Science April 9, 2019 arXiv:1904.03900v1 [astro-ph.EP] 8 Apr 2019 sion system of the spacecraft is the most practical ma- neuver method in the space mission to the present day. Hence, the optimal capture strategy usually can be evalu- ated by the magnitude of the impulsive maneuver. Qi & Xu (2014) proposed the near-optimal lunar capture theory in the Earth-Moon planar CRTBP to locate the optimal capture point and achieve the permanent lunar capture using the minimum V . This method combined the lu- nar gravitational capture and KAM tori around the Moon. The minimum capture velocity and the minimum capture eccentricity before the maneuver were de ned and anal- ysed. The permanent lunar capture after the maneuver was investigated by the distribution of KAM tori around the Moon. Furthermore, based on the near-optimal lunar capture, Qi & Xu (2016) constructed low-energy transfers from the Earth to permanent lunar orbits. The numerical results indicated that compared with traditional patched- Figure 1: The planar con guration of the circular restricted Pluto- conic transfers, those low-energy transfers can eciently Charon system. save fuel costs, which actually were very close to the the- oretical minimum. In this paper, we apply the near-optimal capture the- ing the pre-maneuver velocity V and increasing the post- ory to investigate the capture problem of Charon in the maneuver velocity V . Under given constraints of the cap- Pluto-Charon planar CRTBP. There are two main di er- ture point and Jacobi constant, there exist a minimum V ences between this paper and the previous research of the and a maximum V . Combining two aspects, the optimal Earth-Moon system in (Qi & Xu, 2014). Firstly, the mass capture point with the minimum V = V V can be 1 2 ratio  of the Pluto-Charon system is signi cantly larger found. than that of the Earth-Moon system, so the gravity of According to the above discussion, the structure of this Charon is more signi cant in its region than the case of the paper is separated into ve parts. Section 2 introduces the Moon. Consequently, we postulate that the near-optimal main properties of the dynamical system. Section 3 focuses capture in the Pluto-Charon planar CRTBP can show dif- on the gravitational capture conditions in pre-maneuver ferent characters and results. Secondly, Qi & Xu (2014) portion and obtain the minimum V . In Section 4, we applied KAM tori to analyze the long-term orbital char- investigate the long-term capture in post-maneuver potion acters after the impulsive maneuver V . Although this and obtain the maximum V . After that, the near-optimal method is simple and visible, the position and coverage of capture of Charon in the Pluto-Charon planar CRTBP is KAM tori are based on the empirical observation rather proposed and analyzed in Section 5. Finally, conclusions than the quantitative analysis. In this paper, the smaller of our investigation are presented in Section 6. alignment index (SALI), a common index for determining the regular or chaotic character of an orbit, is applied to in- 2. Description of the dynamical model vestigate and classify the long-term orbits around Charon more accurately instead of KAM tori. According to the theory of the classical circular re- In Zotos (2015) the orbital dynamics in the planar cir- stricted three-body problem (Szebehely, 1967), the two cular Pluto-Charon system were numerically investigated. primary bodies P (Pluto) and P (Charon) move on cir- 1 2 More speci cally, large sets of initial conditions of orbits cular orbits (Buie et al., 2012), with the same angular were classi ed in an attempt to determine the dynamics velocity, around their common center of mass. The third of a test particle (e.g., a spacecraft, or a comet, or an as- body, also known as the test particle (e.g., a spacecraft, teroid) moving in the vicinity of Charon. It was found or a comet, or an asteroid), moves on the same plane un- that the con guration (x; y) space is covered by a highly der the resultant Newtonian gravitational eld of the two complicated mixture of several types of basins, such as main bodies. Taking into account that the mass of the test basins of escape, basins of collision and basins of bounded particle is extremely small, compared to the masses of the motion. In addition, by varying the value of the Jacobi two primaries, we can reasonably assume that the third constant (or in other words the energy level) he managed body does not perturb the motion of the primaries. The to monitor the parametric evolution of all types of orbits. mass ratio  = m =(m + m ), with m > m de nes the 2 1 2 1 2 In this paper, the quantitative method similar to (Zotos, non-dimensional masses of the two primaries m = 1 2015) will be proposed to analyze the long-term orbital and m = . In this notation, the centers of Pluto and dynamics in the Pluto-Charon system. Based on Qi & Charon are located at (; 0) and (1 ; 0), respectively Xu (2014), the minimum V satisfying the permanent (see Fig. 1). capture can be obtained from two perspectives: decreas- 2 The planar motion of a test particle can be described by the following set of di erential equations @U U (x; y) = = x  2y _; @x @U U (x; y) = = y  + 2x; _ (1) @y where m m  (1 ) 1 1 2 2 2 U (x; y) = + + + x + y ; (2) r r 2 2 1 2 is the time-independent e ective potential function of the system, while r = (x + ) + y ; r = (x +  1) + y ; (3) are the distances to the respective primaries. Similarly, the set of the variational equations is given by (x) = x; _ (y) = y _; Figure 2: Schematic depicting the position of Charon (gray dot) along with the scattering region (cyan annulus) and Charon's sphere (z) = z _; of in uence (red circle). More details, regarding the radii of the 2 2 2 circular regions, are given in the text. The pre-maneuver orbit is @ U @ U @ U shown in dashed purple line, while the post-maneuver orbit is shown (x _ ) = x + y + z + 2y _; @x @x@y @x@z in solid green line. The position of the entrance and capture points 2 2 2 are pinpointed by ve-pointed black stars. (Color gure online). @ U @ U @ U (y _ ) = x + y + z 2x; _ @y@x @y @y@z 2 2 2 @ U @ U @ U velocity V . Similar to Qi & Xu (2014), we also assume (z _) = x + y + z: (4) that capture points are apsis points with respect to Charon, @z@x @z@y @z including periapsides and apoapsides. Since the velocity The set of the equations of motion (1) admits only the direction of the test particle at the capture point is perpen- following integral of motion dicular to the vector from the center of Charon to the test 2 2 J (x; y; x; _ y _ ) = 2U (x; y) x _ + y _ = C; (5) particle, the velocity direction has only two possibilities in the orbital plane: the prograde motion (anticlockwise which is also known as the Jacobi integral. The numeri- motion with respect to Charon) or the retrograde motion cal value C of the Hamiltonian remains constant and it is (clockwise motion with respect to Charon). called the Jacobi constant. The gravitational capture requires that the test parti- According to NASA Space Science Data Coordinated cle starting from the capture point can leave the sphere Archive (Update: 23 December 2016) the value of the mass of in uence (SOI) of Charon under a given ight time by ratio for the Pluto-Charon system is  = 0.10851122058 , backward integration. The radius of the SOI of Charon which is signi cantly larger than that of the Earth-Moon is set as 10000 km (see Fig. 2), about double the dis- planetary system, 0.0121506683 (Qi & Xu, 2014). The tance of L from the center of the smaller primary, which values of the Jacobi constant at the ve Lagrange points is similar to the lunar SOI proposed by Yamakawa et al. (L ; i = 1; :::; 5) are denoted by C , i = 1; :::; 5 and they are i i (1992). The ight time from the boundary of the SOI to equal to: C = 3.717080, C = 3.576079, C = 3.204728, 1 2 3 the capture point is restricted to 15 time units (approx- and C = C = 3. 4 5 imately 15.24 days). For the given position and velocity direction of the capture point, we postulate that a larger 3. Gravitational capture the velocity magnitude can make the escape from Charon (that is the Charon capture in the forward propagation) In this section, the gravitational capture in the pre- easier. Under the limitation of the ight time, we can use maneuver portion is investigated to decease the pre-maneuver the following method to nd the minimum V for the given capture point. It should be noted that in Zotos (2015), the value  = 0:099876695437731 was adopted. However, over the last years the 1. The initial point of backward integration is an apsis ocial value of the mass ratio regarding the Pluto-Charon planetary point with respect to Charon chosen from a scatter- system has been revised. 3 ing region. The scattering region of initial points is in two sides of the Moon; however, in the Pluto-Charon an annulus region around Charon, whose inner and PCRTBP, highly-fractal di used structures of C are max outer radii are equal to the radius of Charon and the more complicated in the Charon realm. distance of L from Charon, i.e., 606.0 and 5848.3 In the calculation of C , the integration time from 1 max km, respectively (see Fig. 2). Inside this annulus the initial point to the boundary of SOI of Charon is de- we de ne a uniform and dense grid of 1024  1024 noted by t . Fig. 4 shows the distributions of t in the esc esc initial conditions. There are two kinds of direction Pluto-Charon rotating frame. Similarly, left and right pan- of velocity: the prograde and retrograde motion, but els correspond to prograde and retrograde capture points, its velocity magnitude in the Pluto-Charon rotating respectively. By looking at the gure, it can be seen that in frame V can be changed. most region, orbits starting from prograde capture points are easier than those from the retrograde points to escape 2. According to Eq. (5), V = 2U (x; y) C , so the the Charon's SOI by backward integration. In addition, minimum V can translate to the maximum C . For we can nd some similarities between structures of t and the given capture point (x; y) and Jacobi constant esc C in Fig. 3. C , its velocity in the Pluto-Charon rotating frame, max According to Eq. (6), the minimum V we seek can be (x; _ y _ ), can be obtained from obtained by C . It should be noted that this minimum max velocity is a magnitude in the Pluto-Charon rotating frame V = 2U (x; y) C; rather than in the Charon-centered inertial frame. We de- (6) x _ = kV y=r ; 1 2 ne the corresponding minimum velocity in the Charon- y _ = kV (x +  1)=r ; 1 2 centered inertial frame as V . Since capture points are min apsides with respect to Charon, we can obtain the follow- where k denotes the direction of the velocity: k = 1 ing relationship, for the prograde motion with respect to Charon; k = 1 for the retrograde motion with respect to Charon. V = 2U (x; y) C + kr : (7) min max 2 r is the distance of the initial point to the center of Fig. 5 shows the distributions of V in the Pluto-Charon min Charon, and can be obtained from Eq. (3). In this dimensionless rotating frame. The left panel corresponds way, the initial state of the capture point (x; y; x; _ y _ ) to prograde capture points, while the right one corresponds in the rotating frame can be derived. to retrograde capture points. 3. We choose the initial C as C . Then, using the From Fig. 5, we observe that, generally, if the capture method in step 2, we can obtain the correspond- point is farther from Charon, the minimum capture veloc- ing initial state (x; y; x; _ y _ ). By backward integra- ity V is smaller. This result is accordant with the fact: min tion, a capture orbit can be obtained from the initial the attraction from Charon at a capture point farther to state. If the orbit is bounded in the SOI within 15 Charon is smaller, which results in a smaller escape veloc- time units or collides with Charon, we postulate that ity. As we can see from the gure, for both prograde and the given C cannot satisfy the gravitational capture. retrograde motions, minimum capture velocity V s are min Then we test whether the backward escape can be mainly distributed in the range of 0.28  2.57. However, achieved for a smaller C . In this way, C is decreas- in Fig. 3 C in the prograde orbits is obviously larger max ing until the backward escape can be satis ed. The than that in the retrograde orbits in most areas. We can corresponding maximum C is de ned as C . max explain this phenomenon by Eq. (7): for the same V at min It should be noted that the above numerical method is a given capture point, C of prograde motion is larger max similar to that in (Qi & Xu, 2014), but in step 3, we use than that of retrograde motion. C is derived from the max C as the variable rather than V in (Qi & Xu, 2014). The state of the rotating frame, so the in uence of the Pluto- advantage of C is that its initial value or upper boundary Charon rotation can in uence C . max is speci c, i.e. C . In the two-body model, we can use the orbital eccen- In Fig. 3 we display the distributions of C in the tricity e to determine whether the orbit is elliptic or hyper- max Pluto-Charon rotating frame, where left and right panels bolic. In the vicinity of a primary body of the PCRTBP, correspond to prograde and retrograde capture points, re- since the gravity of that celestial body is dominant, the spectively. From Fig. 3, we discover that C s of the concept of the eccentricity e can still be valid to study the max prograde motion are larger than those of the retrograde orbital character of the test particle. However, it should be noted that in the PCRTBP, e is time-varying. When e motion in most areas. The lower boundary of C of max is alternated from larger than 1 to smaller than 1, we pos- the retrograde motion is signi cantly smaller than that of tulate that gravitational capture happens. The minimum prograde motion. Compared with distributions of C max capture eccentricity e is de ned as the eccentricity at of Earth-Moon PCRTBP in (Qi & Xu, 2014), structures min the capture point corresponding to the minimum V , of C for retrograde motion are very similar. For the max min and can be obtained from prograde motion, the distribution of C in the Earth- max Moon PCRTBP (Qi & Xu, 2014) is more regular and con- r V min e = 1: (8) min secutive, and there exist two branches with smaller C max 4 Figure 3: Distributions of C in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde max motion. (Color gure online). Figure 4: Distributions of t in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde motion. esc (Color gure online). Figure 5: Distributions of V in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde min motion. (Color gure online). Fig. 6 illustrates distributions of the minimum cap- choose the capture points. As we can see from the gure, ture eccentricity e in the Pluto-Charon rotating frame. in the right gure, all retrograde capture points are avail- min Based on the de nition of the gravitational capture, cap- able, and the structures are similar to those of C . Qi max ture points with e > 1 are infeasible for the gravita- & Xu (2014) found that structures of e of the Earth- min min tional capture. We should avoid these regions when we Moon PCRTBP for prograde motion have two branches 5 Figure 6: Distributions of e in the dimensionless Pluto-Charon rotating frame. (a-left): Prograde motion and (b-right): retrograde motion. min (Color gure online). in two sides of the Moon with e > 1. However, as we is covered either by collision or by regular initial condi- max can see from the left gure, the infeasible prograde capture tions of orbits. As long as C < C the channel around L 1 1 points for gravitational capture are only distributed in a opens and therefore orbits can escape to the Pluto realm. small region on the left edge. The stable areas of regular initial conditions are destroyed abruptly. It should be noted that when C = C two new stability islands emerge at both sides of Charon. We ob- 4. Long-term capture serve that as the value of the Jacobi constant decreases the amount of both regular and collision orbits decrease. In this section, we focus on the maximum post-maneuver On the other hand, the amount of escaping orbits grows velocity V satisfying the long-term capture around Charon. rapidly and for C  C they dominate the entire con gura- According to (Qi & Xu, 2014), a regular orbit in KAM tion plane. In the Earth-Moon PCRTBP (Qi & Xu, 2014), tori around Charon can be bounded in the Charon realm KAM tori rapidly disappear when C decreases to C , but permanently in the Pluto-Charon PCRTBP, otherwise, an apparently, regular orbits in the Pluto-Charon system can orbit located in chaotic areas escapes the Charon realm. still survive when C = C . We postulate that this di er- Therefore, we require that a test particle after V should ence is derived from the larger  of Pluto-Charon system, be inserted into a regular orbit around Charon. In (Qi & because orbits rotating around the smaller primary with a Xu, 2014), KAM tori were applied to nd regular bounded larger  have a greater resistance against the perturbation regions. Although this method is simple and visible, the from the larger primary in the PCRTBP. position and coverage of KAM tori rely on the empiri- When C < C , the abrupt emergence of chaotic regions cal observation rather than the quantitative analysis. The in Fig. 7 looks quite similar to the abrupt bifurcation to SALI is a common index for quantitatively determining fully developed chaotic scattering in (Bleher et al., 1990). the regular or chaotic character of an orbit (Skokos, 2001). However, we postulate that the chaotic phenomenon in The analytical de nition of SALI, as well as the particular Fig. 7 essentially is di erent from the chaotic scattering. threshold values for distinguishing between ordered and Firstly, the fully developed chaotic scattering denotes a sit- chaotic motion are given in Zotos & Jung (2018). The uation in which all periodic orbits are unstable and there classi cations of the long-term orbits around Charon can are no KAM surfaces (i.e., the dynamics is hyperbolic) be obtained by the value of the SALI quantitatively, in- (Bleher et al., 1989), but in Fig. 7, obviously, there ex- cluding (i) bounded (regular or chaotic), (ii) escaping and ist KAM tori in the bounded regular region (denoted by (iii) collision (Zotos, 2015). Essentially, other diagnostic green) when C is smaller than C (Qi & Xu, 2014). In indexes for the orbital chaos, such as the the fast Lyapunov addition, in the CRTBP we consider, not all periodic or- indicator (FLI) (Astakhov & Farrelly, 2004), the correla- bits are unstable, for example the distant retrograde orbits tion dimension (Belbruno et al., 2008), and the frequency around the smaller primary are stable (Henon, 1969; Lam analysis method (Laskar, 1990), can also be applied to & Whi en, 2005). Secondly, Koon et al. (2000) pointed distinguish stable regular orbits from chaotic ones. In this out that the chaotic phenomenon in the region around the section, we use the SALI to investigate bounded basins of smaller primary of the CRTBP is ascribed to the existence the long-term orbits around Charon. of transversal homoclinic and heteroclinic points. We can In Fig. 7 we present the orbital structure of the con g- explain these chaotic phenomena by the symbolic dynam- uration (x; y) plane for the prograde motion. It is seen that ics or the horseshoe-like dynamics (Koon et al., 2000). when C > C all the energetically allowed area of motion It should be noted that trapped chaotic motion is very 6 Figure 7: Basin diagrams on the con guration (x; y) plane for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:68; (d): C = C ; (e): C = C ; (f ): C = C . The color code is as follows: non-escaping regular orbits (green); trapped chaotic orbits (yellow); escaping 2 3 4 orbits (cyan); collision orbits (red); energetically forbidden regions (gray). (Color gure online). limited. This means that the corresponding initial condi- orbital structure of the con guration (x; y) plane for a tions appear only as lonely and isolated points, mainly in xed value of the Jacobi constant (or equivalently, for a the vicinity of the stability islands of regular motion. In- xed orbital energy), for both prograde and retrograde deed, in panel (b) of Fig. 7 one can distinguish several near motion. In an attempt to obtain the same information the boundaries of the bounded basins of regular motion the for a continuous interval of C values, we de ne a plane of presence of several initial conditions which correspond to representation in which the x-coordinate is the abscissa, trapped chaotic motion. while the ordinate is given by the Jacobi constant. This In the same vein, the orbital structure of the con gu- means that all orbits start with initial position (x ; 0), i.e, ration plane for the retrograde motion is depicted in Fig. on the x-axis, with initial velocities (0; y _ ), where the y- 8. It is evident that the orbital structure corresponding to component of the velocity is derived from the third of Eqs. the retrograde motion has many di erences with respect (6). to that of the prograde motion. The most notable ones In Fig. 11 we present the basin diagram on the (x; C ) are: (i) bounded regular motion is possible for much lower plane for C 2 [2:5; 5], for both the prograde and retrograde values of the Jacobi constant, or alternatively, for higher motions. The black continuous line indicates the corre- values of the total orbital energy; (ii) the overall orbital sponding Zero Velocity Curve (ZVC), and it is de ned as structure seems less complicated, without the presence of regions with highly fractal basin boundaries. Those results f (x; C ) = 2U (x; y = 0) = C: (9) are similar to those of the Earth-Moon PCRTBP in(Qi & Xu, 2014). We see that for the prograde motion bounded regular or- The distributions of the corresponding escape and col- bits exist only when C > C , while for the retrograde lision time of the orbits of both the prograde and retro- motion the same type of orbits is still possible even below grade motion are given in Figs. 9 and 10, respectively. As C = C . However, for extremely high and low values of expected, the lowest values of both escape and collision C we encounter the same orbital structure for both types time were measured inside the corresponding basins of es- of motion (prograde and retrograde). In particular, for cape and collision, while the highest times were identi ed extremely low values of the Jacobi constant escaping or- mainly in the fractal basin boundaries, where the degree bits completely dominate, while for relatively high values of the unpredictability of the system maximizes. of C bounded and collision types of orbits share the ener- The color-coded diagrams of Figs. 7 and 8 reveal the getically allowed phase space. In Fig. 12 we present the 7 Figure 8: Basin diagrams on the con guration (x; y) plane for the retrograde motion, where (a): C = 3:9; (b): C = C ; (c): C = C ; (d): 2 4 C = 2:9; (e): C = 2:8; (f ): C = 2:5. The color code is as follows: non-escaping regular orbits (green); trapped chaotic orbits (yellow); escaping orbits (cyan); collision orbits (red); energetically forbidden regions (gray). (Color gure online). Figure 9: Distribution of the escape and collision times of the orbits for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:68; (d): C = C ; (e): C = C ; (f ): C = C . All initial conditions of bounded (regular and chaotic) orbits are shown in white. (Color 2 3 4 gure online). distribution of times employed by the test particle in the and retrograde motions, bounded orbits generally erode escape and collision orbits. or disappear with the decrease of C (see Figs 7 and 8). Based on the previous discussion, for both prograde Hence, for a given capture point (x; y) and its direction 8 Figure 10: Distribution of the escape and collision times of the orbits for the retrograde motion, where (a): C = 3:9; (b): C = C ; (c): C = C ; (d): C = 2:9; (e): C = 2:8; (f ): C = 2:5. All initial conditions of bounded (regular and chaotic) orbits are shown in white. (Color gure online). Figure 11: Basin diagrams on the (x; C)-plane for (a-left): the prograde motion and (b-right): the retrograde motion. The color code is as follows: non-escaping regular orbits (green); trapped chaotic orbits (yellow); escaping orbits (cyan); collision orbits (red); energetically forbidden regions (gray). (Color gure online). of motion, there exists a minimum C satisfying long-term yse the minimum V . According to Qi & Xu (2014), the capture in the Charon region. According to Eq. (5), V key to nd the minimum V is to locate the near-optimal corresponding to the minimum C can be obtained, which capture point. Combining results obtained in previous sec- is the maximum V we seek for. tions, near-optimal capture points can be found by the following method: 1. For the given velocity direction, V corresponding 5. Near-optimal capture of Charon to the capture point (x; y) can be calculated by In the previous two sections, we have investigated how V = jV V j 1 2 to obtain the minimum V and the maximum V , succes- 1 2 p p (10) sively. In this section, we synthesize those results to anal- = j 2U (x; y) C 2U (x; y) Cj max 9 Figure 12: Distribution of the corresponding escape and collision times of the orbits for (a-left): the prograde motion and (b-right): the retrograde motion. All initial conditions of bounded (regular and chaotic) orbits are shown in white. (Color gure online). where C is obtained from Figs. 3 or 4, and C is 14 is similar to that of C for retrograde motion in Fig. max max the Jacobi constant of the post-maneuver portion. It 3 (b). should be noted that if the given capture point (x; y) Figs. 15 and 16 show the overlays of V and the is located in the energetically forbidden regions for bounded basins of di erent C for the prograde and ret- the given C , Eq. (10) is invalid and the correspond- rograde motions, respectively. Gray and cyan regions in ing V does not exist. these gures denote energetically forbidden regions and bounded basins, respectively. Dark blue points are near- 2. Using Eq. (10), the distribution of V in the scat- optimal capture points we chose. Based on the analysis of tering region for given velocity direction and C can Fig. 7, only for a larger C can bounded basins of the pro- be obtained. For example, Figs. 13 and 14 display grade motion exist. Hence, near-optimal prograde capture the distributions of V in the Pluto-Charon rotat- points will disappear with the decrease of C . But near- ing frame corresponding to the cases presented in optimal retrograde capture points can survive in smaller Figs. 7 and 8, respectively. In addition, for the same C s. By looking at these gures, it can be seen that since velocity direction and C , the corresponding basin di- distributions of V and bounded basins of prograde mo- agrams can be obtained using the method in Section tion are both quite di erent from those of retrograde mo- 4 (see Figs. 7 and 8). Then, we put the bounded tion, the distributions of the corresponding near-optimal basins onto the distribution of V , and get an over- capture points have signi cant di erences between two lay gure of them. For example, Figs. 15 and 16 kinds of motions. For example, in Fig. 16, near-optimal display overlays for di erent C s and motions. The retrograde capture points are mostly located in the inner background of the overlay is the distribution of V , boundary of bounded basins near Charon, and similar to and cyan regions denote bounded basins. the distribution of bounded basins, near-optimal capture 3. From the overlay gures, points in bounded basins points gradually move away from Charon with the decrease denote feasible points in terms of the long-term cap- of C . However, situations of near-optimal prograde cap- ture after V . Among these feasible capture points, ture points are more complicated according to Fig. 15, those located in the minimum V regions are cho- where regions with highly fractal basin boundaries exist, sen as near-optimal capture points. For example, especially when C = 3:71. In Fig. 15 (a), near-optimal blue points in Figs. 15 and 16 are near-optimal cap- capture points are distributed in the inner boundary of ture points we chose. The corresponding V is the bounded basin nearest Charon. But when C decreases to minimum V we seek for. 3.71, i.e., necks around L just open, near-optimal capture Figs. 13 and 14 display the distributions of V cor- points are distributed in a highly fractal pattern around responding to Figs. 7 and 8, respectively. Gary regions Charon, which is similar to the structures near Charon in in those gures are the energetically forbidden regions in Figs. 3 and 13. When C decreases to 3.60, corresponding the Charon realm for larger C s. Compared Fig. 13 with bounded basins are distributed in two separated stable is- Fig. 3 (a), the structure of the distribution of V for the lands, we can obtain two near-optimal capture regions in prograde motion actually is quite similar to that of C max two separated bounded basins, respectively. But when C for prograde motion. It can be explained by Eq. (10). decreases to C , near-optimal capture points can only be Similarly, the structure of the distribution of V in Fig. distributed in a small region of the left island. Numerical 10 Figure 13: Distributions of V for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:68; (d): C = C ; (e): C = C ; (f ): 2 3 C = C . (Color gure online). Figure 14: Distributions of V for the retrograde motion, where (a): C = 3:9; (b): C = C ; (c): C = C ; (d): C = 2:9; (e): C = 2:8; (f ): 2 4 C = 2:5. (Color gure online). computation indicates that minimum V s corresponding maneuver portions of capture orbits (or long-term capture to near-optimal capture points in Figs. 15 and 16 are all orbits). Red star points are the optimal capture points (or smaller than 23 m/s. impulsive points) chosen from the near-optimal capture re- Using the near-optimal capture points as the initial gions in Figs. 15 and 16 correspondingly. By looking at state, the corresponding pre- and post-maneuver capture the gures, we observe that both pre- and post-maneuver orbits can be constructed. In Figs. 17 and 18, we illus- portions of capture orbits have signi cant di erences be- trate the capture orbits corresponding to the near-optimal tween prograde and retrograde motions. In Fig. 17, since capture points in Figs. 15 and 16, respectively. The g- the corresponding C s of prograde motion are larger than ures are depicted in the dimensionless Pluto-Charon ro- C , the sizes of post-maneuver orbits are smaller than tating frame. Red dashed lines denote pre-maneuver por- those of retrograde motion in Fig. 18. In addition, all tions of capture orbits, while green solid lines denote post- pre-maneuver portions of prograde near-optimal capture 11 Figure 15: Overlays of V and bounded regular basins, for the prograde motion, where (a): C = 3:9; (b): C = 3:71; (c): C = 3:60; (d): C = C . All bounded basins are indicated with cyan color, while the initial conditions of both escape and collision orbits are shown in tones of solar colors, according to the value of V . All near-optimal capture initial conditions are pinpointed by blue dots. (Color gure online). Figure 16: Overlays of V and bounded regular basins, for retrograde motion, where (a): C = C ; (b): C = C ; (c): C = 2:9; (d): C = 2:8. 2 4 All bounded basins are indicated with cyan color, while the initial conditions of both escape and collision orbits are shown in tones of solar colors, according to the value of V . All near-optimal capture initial conditions are pinpointed by blue dots. (Color gure online). 12 Figure 17: Capture obits corresponding to the prograde near-optimal capture points, where the Jacobi constants of the post-maneuver orbit are (a): C = 3:9; (b): C = 3:71; (c): C = 3:60; (d): C = C . The pre-maneuver orbits are shown in dashed, red dashed lines, while the post-maneuver orbits are shown in solid green lines. The black solid lines indicate the boundaries of the scattering region, while the blue solid line delimits the SOI of Charon. The position of the Lagrange points L and L are pinpointed by purple dots. (Color gure online). 1 2 Table 1: Data of prograde near-optimal capture orbits. Table 2: Data of retrograde near-optimal capture orbits. Fig. 17 (x; y) C C V (m/s) Fig. 18 (x; y) C C V (m/s) max max (a) (0:92575062;0:00083565) 3.687 3.90 10.4064 (a) (0:90736622;0:03091922) 3.097 C 22.6293 (b) (0:79538850; 0:00083565) 3.697 3.71 1.2940 (b) (0:89733836; 0:07103064) 2.947 C 3.5683 (c) (0:99093168; 0:15961003) 3.627 3.60 6.0325 (c) (0:88062527;0:09442897) 2.867 2.90 2.5561 (d) (0:82547207;0:00083565) 3.677 C 7.3500 (d) (0:80207374; 0:27158774) 2.777 2.80 3.1574 points pass through the gateway of L . However, in Fig. interplanetary space probe). Although in real space mis- 18, most pre-maneuver portions of retrograde near-optimal sions there are several types of plane perturbations mainly capture points pass through the gateway of L . As we due to the presence of other giant celestial bodies, such as can see Fig. 18, multiple Charon ybys occur during the Jupiter and Saturn, all these perturbations are extremely pre-maneuver portions. Data of prograde and retrograde weak in the vicinity of Charon (which was the scattering near-optimal capture orbits in Figs. 17 and 18 are listed region in our work). Therefore, we argue that the numer- in Tables 1 and 2. From these tables, V s in Fig. 17 are ical results of the Pluto-Charon CRTBP obtained in this smaller than 11 m/s, while V s in Fig. 18 are smaller paper are structurally stable against out of plane pertur- than 23 m/s. Of particular note is that based on the data bations. of Fig. 17 (c) and (d), C s of pre-maneuver portions max are larger than the corresponding C s of post-maneuver 6. Conclusions portions, so the impulsive maneuvers at the near-optimal capture points increase the speed instead of brake. In this paper, the near-optimal capture problem in the The above results we obtained, including optimal cap- Pluto-Charon system was investigated and discussed. In ture points and corresponding capture orbits, could be ap- an attempt to achieve long-term capture around Charon plied in future space mission design (like the New Horizons with the minimum V , optimal capture points were found 13 Figure 18: Capture obits corresponding to the retrograde near-optimal capture points, where the Jacobi constants of the post-maneuver orbit are (a): C = C ; (b): C = C ; (c): C = 2:9; (d): C = 2:8. The pre-maneuver orbits are shown in dashed, red dashed lines, while the 2 4 post-maneuver orbits are shown in solid green lines. The black solid lines indicate the boundaries of the scattering region, while the blue solid line delimits the SOI of Charon. The position of the Lagrange points L and L are pinpointed by purple dots. (Color gure online). 1 2 and located by the near-optimal capture theory. According & Xu, 2014), the structures of C for prograde mo- max to the de nition of V , the capture orbit was divided into tion in the Pluto-Charon PCRTBP are more compli- the pre- and post-maneuver portions. In the pre-maneuver cated in the Charon realm, where highly-fractal dif- portion, the maximum escape Jacobi constant C was fused structures of C exist. Based on the results max max obtained and discussed by backward numerical integra- of t , in most region, orbits starting from prograde esc tion. In the post-maneuver portion, the SALI was applied capture points are easier than those from retrograde to numerically investigate the orbital characters of long- points to escape the Charon's SOI by backward in- term capture orbits around Charon. The initial condi- tegration. tions corresponding to three types of motion: (1) bounded, 2. Based on the de nition of gravitational capture and (2) escaping and (3) collisional were classi ed and distin- distributions of e , we found that all retrograde min guished. Combining the results of pre- and post-maneuver capture points are available for gravitational cap- portions, optimal capture points were located by overlay ture, and its structure of e is similar to that of min gures. Furthermore, corresponding capture orbits were C . Di erent from the two branches in two sides max constructed in the PCRTBP. In addition, the near-optimal of the Moon for prograde motion in the Earth-Moon capture results in the Pluto-Charon system of this paper PCRTBP (Qi & Xu, 2014), the prograde infeasible were compared with those in the Earth-Moon system. We capture points for gravitational capture in the Pluto- postulated that the signi cant di erence of mass ratio Charon PCRTBP are only distributed in a small re- resulted in the di erent results between two planet-moon gion on the left edge. systems. 3. According to the orbital structure of long-term or- The main results of our numerical research can be sum- bits around Charon, we found that the orbital struc- marized as follows: ture corresponding to the retrograde motion has ob- 1. Numerical computation indicated that the C s of vious di erences with respect to that of the prograde max the prograde motion are larger than those of the ret- motion: (1) bounded regular motion is possible for rograde motion in most areas. Compared with distri- much lower values of the Jacobi constant, or alterna- butions of C of the Earth-Moon PCRTBP in (Qi tively, for higher values of the total orbital energy; max 14 (2) the overall orbital structure seems less compli- Astakhov, S.A., Farrelly, D., 2004. Capture and escape in the ellip- tic restricted three-body problem. Monthly Notices of the Royal cated, without the presence of regions with highly Astronomical Society 354, 971-979. fractal basin boundaries. In addition, the lowest val- Bao, C., Yang, H., Barsbold, B., Baoyin, H., 2015. Capturing near- ues of both escape and collision time were measured Earth asteroids into bounded Earth orbits using gravity assist. Astrophysics and Space Science 360, 61. inside the corresponding basins of escape and colli- Baoyin, H-.X., Chen, Y., Li, J.-F., 2010. Capturing near earth ob- sion, while the highest times were identi ed mainly jects. Research in Astronomy and Astrophysics 10, 587-598. in the fractal basin boundaries, where the degree of Belbruno, E.A., AIAA, DGLR, and JSASS, International Electric the unpredictability of the system maximizes. Propulsion Conference, 19th, Colorado Springs, CO, May 11-13, 4. The distributions of near-optimal capture points in- Belbruno, E.A., Miller, J., Jet Propulsion Lab., Pasadena, CA, 1990. dicated that near-optimal prograde capture points Belbruno, E., Topputo, F., Gidea, M., 2008. Resonance transitions disappear rapidly with the decrease of C , but near- associated to weak capture in the restricted three-body problem. optimal retrograde capture points can survive in smaller Advances in Space Research 42, 1330-1351. Bleher, S., Grebogi, C., Ott, E., 1990. Bifurcation to chaotic scat- C s. Near-optimal retrograde capture points are mostly tering. Physica D: Nonlinear Phenomena 46(1), 87-121. located in the inner boundary of bounded basins near Bleher, S., Ott, E., Grebogi, C., 1989. Routes to chaotic scattering. Charon, and gradually move away from Charon with Physical review letters 63(9), 919. the decrease of C . However, situations of prograde Buie, M.W., Tholen, D.J., Grundy, W.M., 2012. The Orbit of Charon is Circular. AJ 144, article id. 15. motion are more complicated, where highly fractal Granvik, M., Vaubaillon, J., Jedicke, R., 2012. The population of patterns around Charon and separated stable islands natural Earth satellites. Icarus 218, 262-277. exist. Hnon, M., 1969. Numerical exploration of the restricted problem, V. Astronomy and Astrophysics 1, 223-238. 5. Numerical results indicated that all pre-maneuver Jehn, R., Campagnola, S., Garcia, D., Kemble, S., 18th International portions of prograde near-optimal capture points pass Symposium on Space Flight Dynamics, 2004. through the gateway of L . However, most pre- 1 Koon W S, Lo M W, Marsden J E, et al., Dynamical systems, the maneuver portions of retrograde near-optimal cap- three-body problem and space mission design, World Scienti c, ture points pass through the gateway of L . Data of Lam, T., Whi en G J., AAS/AIAA Space ight Mechanics Confer- both prograde and retrograde near-optimal capture ence, 2005, paper AAS 05-110. orbits showed that, V s are smaller than 23 m/s. Laskar, J., 1990. The chaotic motion of the solar system: A numerical estimate of the size of the chaotic zones. Icarus 88, 266-291. For the numerical integration (both backward and for- NASA Space Science Data Coordinated Archive: http://nssdc. ward in time) of the equations of motion (1), as well as gsfc.nasa.gov/planetary/factsheet/plutofact.html Press H.P., Teukolsky S.A, Vetterling W.T., Flannery B.P., Numer- of the variational equations (4) we used a double precision ical Recipes in FORTRAN 77, 2nd Ed., Cambridge Univ. Press, Bulirsch-Stoer FORTRAN 77 algorithm (Press et al., 1992), Cambridge, USA, 1992. with a xed time step equal to 10 . In all our calcula- Qi, Y., de Ruiter, A., 2018. Short-term capture of the Earth-Moon system. Monthly Notices of the Royal Astronomical Society 476 tions the numerical error, related to the conservation of (4), 5464-5478. the Jacobi integral (5) was generally smaller than 10 , Qi, Y., Xu, S., 2014. Lunar capture in the planar restricted three- while in most of the cases it was smaller than 10 . All body problem. Celestial Mechanics and Dynamical Astronomy the graphical illustration has been created using the lat- 120, 401-422. Qi, Y., Xu, S., 2016. Earth{Moon transfer with near-optimal lunar est version 11.2 of the software Mathematica (Wolfram, capture in the restricted four-body problem. Aerospace Science 2003). and Technology 55, 282-291. The obtained results in this paper, including optimal Roncoli, R.B., Fujii, K.K., AIAA Meeting Papers on Disc, 2010. capture points and corresponding capture orbits, could be Schoenmaekers, J., Horas, D., Pulido, J.A., 16th International Sym- posium on Space Flight Dynamics, Pasadena, California, 2001. applied in future space mission design. Skokos C., 2001. Alignment indices: a new, simple method for deter- mining the ordered or chaotic nature of orbits.Journal of Physics A 34, 10029-10043. Acknowledgments Strange, N., Landau, D., McElrath, T., Lantoine, G., Lam, T., Overview of mission design for NASA asteroid redirect robotic Yi Qi wishes to acknowledge the support of the Natu- mission concept, Pasadena, CA: Jet Propulsion Laboratory, Na- ral Sciences and Engineering Research Council of Canada tional Aeronautics and Space Administration, 2013. through a Discovery Accelerator Supplement under Grant Szebehely V., Theory of Orbits, Academic Press, New York, 1967. Topputo, F., Belbruno, E., 2015. Earth{Mars transfers with ballistic RGPAS-493042-2016. The authors would like to express capture. Celestial Mechanics and Dynamical Astronomy 121, 329- their warmest thanks to the two anonymous reviewers for the careful reading of the manuscript as well as for all the Verrier, P.E., McInnes, C.R., 2014. Low-Energy Capture of Aster- oids onto Kolmogorov{Arnold{Moser Tori. Journal of Guidance, valuable suggestions and comments which allowed us to Control, and Dynamics 38, 330-335. improve both the quality and the clarity of the paper. Wolfram S., The Mathematica Book, Wolfram Media, Champaign, Yamakawa, H., Kawaguchi, J., Ishii, N., Matsuo, H., A numerical References study of gravitational capture orbit in the Earth-Moon system, Space ight mechanics, 1992, pp. 1113{1132. Astakhov, S.A., Burbanks, A.D., Wiggins, S., Farrelly, D., 2003. Yarnoz, D.G., Sanchez, J.P., McInnes C.R., 2013. Easily retrievable Chaos-assisted capture of irregular moons. Nature 423, 264-267. 15 objects among the NEO population. Celestial Mechanics and Dy- namical Astronomy 116, 367-388. Zotos, E.E., 2015. Orbit classi cation in the planar circular Pluto- Charon system. Astrophysics and Space Science 360, article id. Zotos, E.E., Jung, Ch., 2018. Correlating the escape dynamics and the role of the normally hyperbolic invariant manifolds in a binary system of dwarf spheroidal galaxies. International Journal of Non- Linear Mechanics 99, 182-203.

Journal

Astrophysics – arXiv (Cornell University)

Published: Apr 8, 2019

There are no references for this article.