aerospace
Article Some Special Types of Orbits around Jupiter
Yongjie Liu 1, Yu Jiang 1,*, Hengnian Li 1,2,* and Hui Zhang 1
1 State Key Laboratory of Astronautic Dynamics, Xi’an Satellite Control Center, Xi’an 710043, China; [email protected] (Y.L.); [email protected] (H.Z.) 2 School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China * Correspondence: [email protected] (Y.J.); [email protected] (H.L.)
Abstract: This paper intends to show some special types of orbits around Jupiter based on the mean element theory, including stationary orbits, sun-synchronous orbits, orbits at the critical inclination,
and repeating ground track orbits. A gravity model concerning only the perturbations of J2 and J4 terms is used here. Compared with special orbits around the Earth, the orbit dynamics differ greatly: (1) There do not exist longitude drifts on stationary orbits due to non-spherical gravity
since only J2 and J4 terms are taken into account in the gravity model. All points on stationary orbits are degenerate equilibrium points. Moreover, the satellite will oscillate in the radial and North-South directions after a sufficiently small perturbation of stationary orbits. (2) The inclinations of sun-synchronous orbits are always bigger than 90 degrees, but smaller than those for satellites around the Earth. (3) The critical inclinations are no-longer independent of the semi-major axis and eccentricity of the orbits. The results show that if the eccentricity is small, the critical inclinations will decrease as the altitudes of orbits increase; if the eccentricity is larger, the critical inclinations will increase as the altitudes of orbits increase. (4) The inclinations of repeating ground track orbits are monotonically increasing rapidly with respect to the altitudes of orbits.
Keywords: mean element theory; orbit dynamics; longitude drift; critical inclinations
Citation: Liu, Y.; Jiang, Y.; Li, H.; Zhang, H. Some Special Types of Orbits around Jupiter. Aerospace 2021, 8, 183. https://doi.org/10.3390/ 1. Introduction aerospace8070183 Jupiter is the most massive planet in the solar system. It is also a gas giant planet. Its mass is 2.5 times the mass of other planets in the solar system. The large size of Jupiter also Academic Editor: Fanghua Jiang makes it relatively easy to be observed. As a result, it was discovered very early. Jupiter has been one of the major targets for planetary exploration. However, Jupiter’s powerful Received: 19 April 2021 magnetosphere and radiation belts are threats to all human spacecraft trying to visit Jupiter. Accepted: 4 July 2021 Since the 1970s, several space missions have been launched by NASA, such as Pioneer X, Published: 8 July 2021 Voyager 1, Galileo, and Juno. However, there is still a long way to go to explore Jupiter. There exist many interesting phenomena that have attracted people to explore for many Publisher’s Note: MDPI stays neutral years, for example, the Great Red Spot, Jupiter’s rings, and atmospheric jet streams. As with regard to jurisdictional claims in a gaseous planet, Jupiter cannot be explored by landing like a lithospheric planet. We published maps and institutional affil- can only use probes to orbit and enter Jupiter’s atmosphere. Therefore, it is necessary to iations. investigate the orbit dynamics around Jupiter. Weibel et al. [1] researched stable orbits between Jupiter and the Sun. Jacobson [2] investigated the gravity field of Jupiter and the orbits of its Galilean satellites. Colwell et al. [3] studied the exogenic dust ring. Recently, Liu et al. [4,5] discussed the dust in the Jupiter system outside the rings and distribution of Copyright: © 2021 by the authors. Jovian dust ejected from the Galilean satellites. Research about this will surely contribute Licensee MDPI, Basel, Switzerland. to the orbit design of space missions. This article is an open access article Investigating the dynamical environments around planets has been the focus for distributed under the terms and space missions in the past few decades. For this purpose, much research concerning conditions of the Creative Commons various special artificial satellite orbits around planets have been conducted. Usually, these Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ special types of orbits include stationary orbits, frozen orbits, sun-synchronous orbits, 4.0/). repeating ground-track orbits, and orbits at the critical inclination. The original idea about
Aerospace 2021, 8, 183. https://doi.org/10.3390/aerospace8070183 https://www.mdpi.com/journal/aerospace Aerospace 2021, 8, 183 2 of 15
geostationary orbits was first put up by Clarke [6]. He pointed out that satellites with an altitude of 36,000 km above the equator of the Earth would have the same rotation rates as the Earth and stay stationary relative to an observer on the equator. Therefore, geostationary satellites are often used for the sake of communications and navigations. Four equilibrium solutions for geostationary orbits were shown to exist by Musen & Bailie [7]. Moreover, two of them were stable while the other two were unstable. Lara & Elipe [8] calculated periodic orbits around equilibrium points in the Earth second degree and order gravity field. For Mars, the stationary orbits, also known as areostationary orbits, and the equilibrium points were studied by Liu et al. [9]. The periodic orbits around the equilibrium points were also calculated by Liu et al. [10]. For orbits at the critical inclination, the eccentricity and argument of perigee are invariant on average. The concept of the Earth critical inclination was first introduced by Orlov [11]. Brouwer [12] used canonical transformations to eliminate short-period terms. Coffey et al. gave a geometrical interpretation of the critical inclination for satellites by investigating the averaged Hamilton system [8]. Representatives of orbits at the critical inclination are the Russian Molniya satellites. The combined effects of the critical inclination and the 2:1 mean motion resonance of a Molniya orbit have been intensively studied since then, for example, in [13–17]. Similarly, frozen orbits are characterized by the invariance of average eccentricity and argument of perigee. Frozen orbits are not limited to specific inclinations. They may exist at any inclination. Usually, the argument of perigee is equal to 90 or 270 deg, depending on the sign of the ratio of the harmonic coefficients J3 and J2. Frozen orbits were first proposed by Cutting et al. [18] for orbit analysis of the Earth satellite SEASAT-A. Coffey et al. [19] showed that there exist three families of frozen orbits in the averaged zonal problem up to J9 in the gravity field of an Earth-like planet. The frozen orbits around the moon in the full gravity model were considered by Folta & Quinn [20], and Nie & Gurfil [21]. Some researchers also view orbits at the critical inclination as frozen orbits, for example [19,22]. Sun-synchronous orbits are defined with a precession rate of the orbital plane equal to the revolution angular velocity around the sun. Generally, remote sensing satellites are placed into these orbits. Macdonald et al. [23] used an undefined, non-orientation- constrained, low-thrust propulsion system to consider an extension of the sun-synchronous orbits. For repeating ground track orbits, the trajectory ground track repeats after a whole number of revolutions within some days. Orbits of this type are widely used to achieve better coverage properties. Lara showed that orbits repeating their ground track on the surface of the Earth were members of periodic-orbit families of the tesseral problem of the Earth artificial satellite [24]. Lei [25] considered the leading terms of the Earth’s oblateness and the luni-solar gravitational perturbations to describe the secular dynamics of navigation satellites moving in the medium Earth orbit and geosynchronous orbit regions. Liu et al. [9] calculated these five types of special orbits around Mars with analytical formulations and numerical simulations. In fact, the gravity field of Mars shares many similarities with that of the Earth. The J2 terms of them are dominant among the harmonic coefficients. However, the J2 term is not as dominant as Earth’s J2. The other first few harmonic coefficients are also strong for Mars: about 1–2 orders of magnitude smaller than J2; for the Earth, the other first few harmonic coefficients are about 3–4 orders of magnitude smaller than J2. The situation is rather different for Jupiter compared with the Earth and Mars. Due to the difficulty of determining the gravity field of Jupiter, there exist few studies about special types of orbits around Jupiter, as far as we know. However, the situation has greatly improved since Juno’s gravitational measurements were conducted. Iess et al. [26] provided measurements of Jupiter’s gravity harmonics (both even and odd) through precise Doppler tracking of the Juno spacecraft. Moreover, they pointed out a North-South asymmetry, which is a signature of atmospheric and interior flows. Here, we mainly use the results Aerospace 2021, 8, 183 3 of 15
in [26] to build a simplified gravity model of Jupiter. Some harmonic coefficients of the gravity model of the Earth, Mars, and Jupiter can be seen in Table1.
Table 1. Some harmonic coefficients and their ratios of the Earth, Mars, and Jupiter.
−3 −6 −5 −6 −3 −2 Planet J2(×10 ) J3(×10 ) J4(×10 ) J22(×10 ) J3/J2(×10 ) J4/J2(×10 ) Earth ([27]) 1.08263 −2.53266 −1.61962 1.81534 −2.3394 −1.4960 Mars ([28]) 1.95545 31.4498 −1.53774 63.0692 16.083 −0.7864 Jupiter ([26]) 14.6965 0.042 −58.661 0 0.0021378 −3.9923
In this paper, we investigate some special orbits around Jupiter, considering mainly the effect of the non-spherical perturbation of the gravity field. From Table1, one can see that the J2 term is still dominating, J2 is about 25 times the 5 value of J4, but is 10 times bigger than J3. The terms, such as J3, J21, J22, J5, J6 (the values and their uncertainty can be seen in [26]), can be neglected compared with the terms of J2 and J4. Therefore, a good approximation of the gravity model of Jupiter is given by " # µ J R 2 J R 4 U = 1 − 2 3 sin2 φ − 1 − 4 35 sin4 φ − 30 sin2 φ + 3 , (1) r 2 r 8 r
−11 3 −1 −2 where µ = GMJ is the gravitational constant of Jupiter, G = 6.67428 × 10 m · kg · s , MJ is the mass of Jupiter, R is the radius of Jupiter, r is the distance of the satellite relative to the center of mass of Jupiter, φ is the latitude of the satellite. Equation (1) indicates that the gravity model of Jupiter that we use here is symmetrical with respect to the z-axis. This leads to different characteristics of satellites orbiting around Jupiter and the Earth or Mars. In the next sections, we adopt the gravity model represented by Equation (1) and use it to study some mean features of orbits around Jupiter.
2. Stationary Orbits Satellites on stationary orbits are well-known for their stationary ground track. There- fore, stationary orbits are preferred for designing communications and navigation satellites. There exist numerous studies on stationary orbits of the Earth and Mars. However, the gravity field of Jupiter is significantly different. In this subsection, we will calculate the stationary orbit of Jupiter and investigate their stability in a spherical coordinate system. In the spherical coordinates of an inertial frame, O − r, λ, φ, where O is the center of mass of Jupiter, λ is the jovicentric longitude, r and φ are the same as those in Equation (11), and the kinetic energy of the spacecraft can be written as 1 . 2 . 2 . 2 T = r + r2 cos2 φλ + r2φ , (2) 2
. . . where r, λ, φ are the derivatives with respect to time. From the expressions of T and U, one can see that λ is a cyclic variable. Let us introduce q = [r, λ, φ]T, the Lagrangian can be written as L = T − U. By Lagrange equations, we have
d ∂T ∂T ∂U . − = . (3) dt ∂q ∂q ∂q
More precisely, the equations of motion can be presented as follows: .. . 2 . 2 2 4 r − r cos2 φλ − rφ = − µ + 3µJ2R 3 sin2 φ − 1 + 5µJ4R 35 sin4 φ − 30 sin2 φ + 3 , r2 2r4 8r6 . d 2 2 dt r cos φλ = 0, (4) . . 2 h 2 4 i d r2 + 1 r2 ( ) = − 3µJ2R + µJ4R 2 − dt φ 2 sin 2φ λ 2r3 8r5 70 sin φ 30 sin 2φ. Aerospace 2021, 8, x FOR PEER REVIEW 4 of 17
.. μ 35μμJR24 JR r−−=−+−+−+rrcos22φλ φ 224() 3sin 2 φ 1() 35sin 4 φ 30 sin 2 φ 3 , rr2428 r 6 d 22 ()r cos φλ = 0, (4) dt d 1 3μμJR24 JR ()22φφλ+=−+−() 224() 2 φφ rrsin270sin30s. 35 in 2 dt 2 28rr
From the second equation of (4), we see that the quantity r 22cos φλ is invariant, which can also be obtained by conservation of the angular momentum along the z-axis. Moreover, the zonal terms of the gravity field only lead to radial and North-South pertur- bations. The presence of these terms increases the radius of the stationary orbit with re- spect to the case of a spherical planet with the same mass of Jupiter. When the orbital ⋅ plane coincides with the equatorial plane, namely φφ==0, 0 , the vertical perturbation ==λλφφ = = == vanishes. In order to find stationary orbits for Jupiter, let rr0, nJ , 0, 0 in these equations; we get μ 315μμJ RJR24 n2 =+24 − , J rr3528 r 7 (5) sin 2φ = 0, Aerospace 2021, 8, 183 4 of 15 or
24. μ 35μμJR JR2 2 From the second equation−+ of (4),24 we see that24222 theφφφφ −+ quantity r cos φλ is invariant, − which ++ = ()3sin 1() 35sin 30sin 3rncos J 0 , can also be obtained by conservationrr2428 of the angular momentum r along 6 the z-axis. Moreover, the zonal terms of the gravity field only lead to radial and North-South perturbations. The (6) 1 3μμJR24 JR presence of these terms increasesrn22+= the radius24 of+− the stationary()70 sin orbit 2φ with 30 respect 0, to the case of a spherical planet with theJ same mass35 of Jupiter. When the orbital plane coincides with 2 28· rr the equatorial plane, namely φ = 0, φ = 0, the vertical perturbation vanishes. In order to ...... find stationary orbits for Jupiter, let r = r = 0, λ = n , λ = 0, φ = φ = 0 in these equations; where nJ is the rotational angularJ velocity of Jupiter. One can verify that the left-hand we get ( 2 4 side of the second equationn2 = µ + in3µJ 2(6)R − is15 µalwaysJ4R , positive when r is larger than R . Therefore, J r3 2r5 8r7 (5) we only need to analyzesin 2 φsolution= 0, s of Equation (5). From the second equation of (5), we see or π φ = φ π that only one2 meaningful4 latitude of the stationary orbit exists, i.e., 0 ( =, are − µ + 3µJ2R 3 sin2 φ − 1 + 5µJ4R 35 sin4 φ − 30 sin2 φ + 3 + r cos2 φn 2 = 0, r2 2r4 8r6 J 2 2 4 (6) 1 r2n 2 + 3µJ2R + µJ4R 70 sin2 φ − 30 = 0, π also2 rootsJ 2rof3 sin8r5 () 2φ = 0 , but φ= makes no sense for stationary orbits, and φ=π where nJ is the rotational angular velocity of Jupiter. One2 can verify that the left-hand side ofcorresponds the second equation to ina (6)stationary is always positive orbit when whichr is larger coincides than R. Therefore, with wethat only of φ=0 ). Introducing the func- need to analyze solutions of Equation (5). From the second equation of (5), we see that only one meaningful latitudeμ of theμμ stationary24 orbit exists, i.e., φ = 0 (φ = π , π are also rootsμ 315JR24 JR 2 2 2 of sin(s2φ) = 0, but=+φ = π makes no sense − for stationary − orbits, and φ = π correspondsf ()rn=− tion f0 ()rn352 7J and 1 3 J , the variations of to a stationary orbit whichrr coincides28 with that r of φ = 0). Introducing the functionsr 2 4 f (r) = µ + 3µJ2R − 15µJ4R − n2 and f (r) = µ − n2, the variations of f (r), f (r), with 0f01(),rfrr3 ()2r5 , with8r7 respectJ to1 r , rcan3 Jbe seen in the0 following1 Figure 1. respect to r, can be seen in the following Figure1.
Figure 1. Variations of f0(r), f1(r) as r increases from 2.2R to 2.3R. Figure 1. Variations of f01(),rfr () as r increases from 2.2R to 2.3R . It can be proved that f0(r) vanishes at two real positive values of r, but only one of 8 them is bigger than R. This solution, which is given by r0 = 2.2414R = 1.6024 × 10 m, is denoted by E0 in Figure1. If the perturbation due to the non-spherical gravity field is not considered, one can find that the radius for stationary orbits is about 2.2381R, which corresponds to the point E1. Therefore, the existence of the J2 and J4 terms increase the altitude of the stationary orbit. Aerospace 2021, 8, 183 5 of 15
To study the stability of the stationary orbits under small perturbations, we mainly used the epicyclic theory (Murray & Dermott [29]). We can denote by the constant h, the . value taken by r2 cos2 φλ for some given initial conditions. Assuming that the deviations from the stationary orbits are small and writing r = r0(1 + ε), Equation (4) can be linearized as follows (Murray & Dermott [29], Section 11 of Chapter 6): .. ε + (3α1 + α2)ε = 0, .. φ + (α1 + β2)φ + β1 = 0, · (7) − h = λ 2 2 2 0 r0 (1+ε) cos φ
1 ∂U ∂2U 1 ∂U 1 ∂2U where α1 = − (r0, 0), α2 = − 2 (r0, 0), β1 = − 2 (r0, 0), and β2 = − 2 2 (r0, 0). r0 ∂r ∂r r0 ∂φ r0 ∂ φ By solving Equation (7), the analytical approximate expressions of ε, φ, λ can be formulated as [29] ε = e cos k t, 1 β1 φ = − + i cos k2t, α1+β2 √ (8) √ 2 α λ = α t − 1 e sin k t, 1 k1 1 2 2 where k = 3α1 + α2, k = α1 + β2, i is the inclination, and e is the eccentricity of the orbit. 1 2 . √ √ Note that the average of λ is α1. Let us set k3 = α1. The first two equations in (7) or (8) show that the radial and North-South motions are uncoupled. Using the expression (1), we can get a precise form of the three frequencies k1, k2, k3 [29]:
2 4 k 2 = µ − 3 J R + 45 J R 1 r 3 1 2 2 r 8 4 r , 0 0 0 2 4 2 µ 9 R 75 R k2 = 3 1 + 2 J2 r − 8 J4 r , (9) r0 0 0 2 4 2 µ 3 R 15 R k3 = 3 1 + J2 − J4 . r0 2 r0 8 r0
Therefore, the satellites will oscillate in the radial and North-South directions with frequencies k1 and k2, respectively. Moreover, the mean motion frequency, k3, in the West- East direction is larger than in the Keplerian case. Namely, for a given semi-major axis, the satellite moves faster than the rate expected at that location in the Keplerian case. In the following, we give some numerical examples to illustrate the above characteristics. We computed the evolution of r, φ, and λ when we selected initial values of r and φ close to r0 and 0, respectively. Figure2 shows a stationary orbit and four orbits obtained from r = 0.99r0, φ = ±0.1 deg and r = 1.01r0, φ = ±0.1 deg. We can see that these four disturbed orbits are no longer periodic. This can be explained by using Equation (9). Due to the presence of the terms that contain J2 and J4, the three frequencies are usually not equal or even commensurable. Figure3 illustrates that the satellite oscillates in the radial and North-South directions. On the other hand, the satellite on stationary orbits does not drift . in these directions. It only moves with a constant λ along the orbit. Therefore, the drift of longitude would not occur for satellites on stationary orbits since only zonal harmonics are taken into account in the gravity model. Furthermore, there are no significant differences between points on the stationary orbit. As a result of the symmetry of the gravity field, it can be concluded that the points on stationary orbits are degenerate equilibrium points. Here, an equilibrium point is degenerate if, and only if, the matrix of the linearized equation for (4) is degenerate. This is the major difference with respect to stationary orbits of other planets, such as the Earth and Mars. For the Earth and Mars, there exist four equilibrium points on stationary orbits, among which two are stable and the other two are unstable. Aerospace 2021, 8, x FOR PEER REVIEW 6 of 17
Aerospace 2021, 8, x FOR PEER REVIEW 6 of 17
==±φ ==±φ obtained from rr0.990 , 0.1 deg and rr1.010 , 0.1 deg. We can see that these four disturbed orbits are no longer periodic. This can be explained by using Equation (9). ==±φ ==±φ obtained from rr0.990 , 0.1 deg and rr1.010 , 0.1 deg. We can see that these Due to the presence of the terms that contain J 2 and J 4 , the three frequencies are usu- allyfour not disturbed equal or orbits even arecommensurable. no longer periodic. Figure This 3 illustrates can be explained that the bysatellite using oscillatesEquation in(9). theDue radial to the and presence North-South of the directions. terms that On contain the ot herJ 2 hand,and J the4 , the satellite three on frequencies stationary are orbits usu- doesally not not drift equal in orthese even directions. commensurable. It only moves Figure with 3 illustrates a constant that λ the along satellite the orbit. oscillates There- in fore,the radialthe drift and of North-South longitude would directions. not occur On the fo rot satellitesher hand, on the stationary satellite on orbits stationary since onlyorbits zonaldoes harmonics not drift in are these taken directions. into account It only in moves the gravity with a model. constant Furthermore, λ along the there orbit. are There- no significantfore, the driftdifferences of longitude between would points not on occur the st foationaryr satellites orbit. on As stationary a result of orbits the symmetry since only ofzonal the gravity harmonics field, are it can taken be concludedinto account that in the the points gravity on model. stationary Furthermore, orbits are degeneratethere are no equilibriumsignificant points.differences Here, between an equilibrium points on point the st isationary degenerate orbit. if, As and a resultonly if, of the the matrixsymmetry of theof linearizedthe gravity equation field, it can for be(4) concluded is degenerate. that theThis points is the on major stationary difference orbits with are respectdegenerate to stationaryequilibrium orbits points. of other Here, planets, an equilibrium such as thepoint Earth is degenerate and Mars. if, For and the only Earth if, the and matrix Mars, of therethe linearizedexist four equilibriumequation for points (4) is degenerate.on stationary This orbits, is the among major which difference two arewith stable respect and to thestationary other two orbits are unstable. of other planets, such as the Earth and Mars. For the Earth and Mars, Aerospace 2021, 8, 183 6 of 15 there exist four equilibrium points on stationary orbits, among which two are stable and the other two are unstable.
− − = φ =− − Figure 2. Orbits around Jupiter: green stationary orbits; red rr0.99 0 , 0.1 deg; blue = φ = − = φ =− − = φ = rr0.99 0 and 0.1 deg; cyan rr1.01 0 and 0.1 deg; black rr1.01 0 and 0.1 Figure 2. Orbits around Jupiter: green—stationary− orbits; red—− r = 0.99r , φφ ==−−0.1 deg; blue—− deg.Figure 2. Orbits around Jupiter: green stationary orbits; red rr0.99 00, 0.1 deg; blue r ==0.99r0 and φφ==0.1 deg; cyan—r−= 1.01= r0 and φ =φ−=−0.1 deg; black—r −= 1.01= r0 and φ = 0.1φ = deg. rr0.99 0 and 0.1 deg; cyan rr1.01 0 and 0.1 deg; black rr1.01 0 and 0.1 deg.
Aerospace 2021, 8, x FOR PEER REVIEW 7 of 17
(a) (b)
(a) (b)
(c) (d)
(e) (f)
FigureFigure 3. 3.(a,b(a,)b Variations) Variations of of orbital orbital radius radius ( (rr),), ((cc,d,d)) driftdrift ofof longitudelongitude ( ∆( Δλλ), (),e ,(fe,f) evolution) evolution of of latitude lati- φ − = φ =− − = tude(φ): ( green—stationary): green stationary orbits; orbits; red, r =red,0.99 rrr0 and0.99φ0=and−0.1 deg;0.1 blue— deg;r =blue0.99 r0 andrr0.99φ = 00.1 anddeg; φ = = − = = − φ =− = =− = φ = cyan—0.1 deg;r 1.01cyanr0 andrrφ 1.010.10 and deg; black—0.1r deg;1.01 blackr0, φ 0.1rr deg,1.01 and0 , T is a0.1 Joviandeg, day.and T is a Jovian day.
In the following, we show that the effects of third bodies, the Jovian ring and the
magnetic field of Jupiter, are negligible compared to those of the J 2 and J 4 terms. Based on the data of the Sun and moons around Jupiter (see for example, [30]), one
can calculate the maximal ratio of the disturbed acceleration ( ad ) and central gravity ac-
celeration ( am ) for satellites on stationary orbits, which is achieved when the satellite, Ju- piter, and the third body are in a straight line. Therefore, the maximal ratio can be calcu- lated as [31]
3 ard m 0 = 2 , (10) aMρ mJmax where m is the mass of the third body, and ρ is the distance between Jupiter and the third body. Values of the maximum ratio (109) for different bodies are reported in Table
2. Since the J 4 term of Jupiter is 100 times bigger than the ratio of acceleration for Io (which gives the highest value among the Galilean satellites), it is reasonable to ignore these effects when investigating the qualitive character of stationary orbits on short time scales. To verify the validity of these solutions, we use the numerical integration method to see the effect of the Galilean moon Io in 800 Jovian days. The orbital elements ()aei,,,Ω ,ω , M of Io that we adopt here are taken from the 10th China Trajectory Optimization Competition, i.e., ()422029.687km, 0. 004308, 0.04 deg,− 79.64 deg,37. 9 91deg, 4.818 deg . Calculation results show that the drift of longitude ( Δλ ) for the satellite on stationary orbits is less than 0.1
Aerospace 2021, 8, 183 7 of 15
In the following, we show that the effects of third bodies, the Jovian ring and the magnetic field of Jupiter, are negligible compared to those of the J2 and J4 terms. Based on the data of the Sun and moons around Jupiter (see for example, [30]), one can calculate the maximal ratio of the disturbed acceleration (ad) and central gravity acceleration (am) for satellites on stationary orbits, which is achieved when the satellite, Jupiter, and the third body are in a straight line. Therefore, the maximal ratio can be calculated as [31] a m r 3 d = 2 0 , (10) am max MJ ρ where m is the mass of the third body, and ρ is the distance between Jupiter and the third body. Values of the maximum ratio (109) for different bodies are reported in Table2. Since the J4 term of Jupiter is 100 times bigger than the ratio of acceleration for Io (which gives the highest value among the Galilean satellites), it is reasonable to ignore these effects when investigating the qualitive character of stationary orbits on short time scales. To verify the validity of these solutions, we use the numerical integration method to see the effect of the Galilean moon Io in 800 Jovian days. The orbital elements (a, e, i, Ω, ω, M) of Io that we adopt here are taken from the 10th China Trajectory Optimization Competi- tion, i.e., 422, 029.687 km , 0.004308, 0.04 deg, −79.64 deg, 37.991 deg, 4.818 deg. Cal- culation results show that the drift of longitude (∆λ) for the satellite on stationary orbits is less than 0.1 deg in 800 Jovian days. The inclination changes no more than 1.0 × 10−4 deg. The oscillation of the semi-major axis is less than 0.13% of the obital radius. However, it should be noted that the obtained solutions may not be correct on long time scales.
Table 2. The disturbing acceleration due to the third body.
−5 −6 The Third Body m/MJ(×10 ) ρ/r0 (ad/am)max(×10 ) Sun ([32]) 1.0473 × 108 4856.8 0.01823 Io ([30,32]) 4.7047 2.6311 5.20000 Europa ([30,32]) 2.5283 4.1868 0.68904 Ganymede ([30,32]) 7.8056 6.6775 0.52379 Callisto ([30,32]) 5.6673 11.7511 0.06989
The Jovian main ring is about 6440 km wide and probably less than 30 km thick. We denote this width by w. Moreover, the distance, d, of the ring from the center of mass of 13 Jupiter is about 122,500 km, and the mass, mr, is about 1.0×10 kg [32]. For satellites on stationary orbits with position rs, the gravitational acceleration due to the Jovian main ring can be calculated as
ZZZ (r(P) − r )ρ(P)dV = s Fr G 3 , (11) ring |r(P) − rs|
where r(P), ρ(P) denote the position vector and the density, respectively, which depend on the position of a point, P, belonging to the ring. We can see that
Gmr |r(P) − rs|≥ r0 − d − w ≈ 160, 240 − 122, 500 − 6440 ≈ 0.4378R, |Fr| ≤ . (0.4378R)2
13 2 Gmr GMJ mr r0 2 1.0×10 2.2414 −13 Note that 2 / 2 = · ≈ 27 · ≈ 1.3795 × 10 , (0.4378R) r0 MJ 0.4378R 1.9×10 0.4378 which is far smaller than J2 and J4 terms. Therefore, the effects of the Jovian main ring on stationary orbits can be neglected. Another perturbation that may affect stationary orbits is the magnetic field of Jupiter when the satellite is charged. Here, we briefly analyze the effects under some assumptions. Aerospace 2021, 8, 183 8 of 15
The magnetic field of Jupiter near stationary orbits when exterior terms are neglected can be written as B(r, φ, λ) = −∇V(r, φ, λ), (12) where ∞ l l+1 R m V = R∑ ∑ [gl,m cos(mλ) + hl,m sin(mλ)]Pl (cos φ). (13) l=1 m=1 r
Here, gl,m and hl,m are the geomagnetic Gauss coefficients (their values can be found m in [33,34]), and Pl (cos φ) is the normalized Legendre function. The acceleration due to the Lorentz force is q . F = r − ω × r × B(r, φ, λ), (14) L m −12 −1 where q = 4πε0Φs,ε0 = 8.854187817 × 10 F · m is the vacuum permittivity, Φ is the surface potential of satellites, s is the radius of the satellite, ω is the angular velocity vector of Jupiter. For a charged satellite with Φ = 500 V (an in-depth study of the effect of spacecraft charging at Jupiter can be found in [35]), s = 10 m, and m = 100 kg, we have −7 −8 q = 5.5633 × 10 C and FL ≤ 1.14 × 10 N. The Lorentz force compared with the lowest order of gravity provides the ratio