Next Article in Journal
Large-Eddy Simulation of Low-Pressure Turbine Cascade with Unsteady Wakes
Next Article in Special Issue
Three-Axes Attitude Control of Solar Sail Based on Shape Variation of Booms
Previous Article in Journal
The Lunar Radiation Environment: Comparisons between PHITS, HETC-HEDS, and the CRaTER Instrument
Previous Article in Special Issue
\({\mathcal{L}_1}\) Adaptive Loss Fault Tolerance Control of Unmanned Hypersonic Aircraft with Elasticity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Special Types of Orbits around Jupiter

1
State Key Laboratory of Astronautic Dynamics, Xi’an Satellite Control Center, Xi’an 710043, China
2
School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Authors to whom correspondence should be addressed.
Submission received: 19 April 2021 / Revised: 2 July 2021 / Accepted: 4 July 2021 / Published: 8 July 2021
(This article belongs to the Special Issue Spacecraft Dynamics and Control)

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 J 2 and J 4 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 J 2 and J 4 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.

1. Introduction

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 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 magnetosphere and radiation belts are threats to all human spacecraft trying to visit Jupiter. Since the 1970s, several space missions have been launched by NASA, such as Pioneer X, 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 years, for example, the Great Red Spot, Jupiter’s rings, and atmospheric jet streams. As a gaseous planet, Jupiter cannot be explored by landing like a lithospheric planet. We can only use probes to orbit and enter Jupiter’s atmosphere. Therefore, it is necessary to 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 Jovian dust ejected from the Galilean satellites. Research about this will surely contribute to the orbit design of space missions.
Investigating the dynamical environments around planets has been the focus for space missions in the past few decades. For this purpose, much research concerning various special artificial satellite orbits around planets have been conducted. Usually, these special types of orbits include stationary orbits, frozen orbits, sun-synchronous orbits, repeating ground-track orbits, and orbits at the critical inclination. The original idea about 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,14,15,16,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 J 3 and J 2 . 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 J 9 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 J 2 terms of them are dominant among the harmonic coefficients. However, the J 2 term is not as dominant as Earth’s J 2 . The other first few harmonic coefficients are also strong for Mars: about 1–2 orders of magnitude smaller than J 2 ; for the Earth, the other first few harmonic coefficients are about 3–4 orders of magnitude smaller than J 2 .
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 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 Table 1.
In this paper, we investigate some special orbits around Jupiter, considering mainly the effect of the non-spherical perturbation of the gravity field.
From Table 1, one can see that the J2 term is still dominating, J2 is about 25 times the value of J 4 , but is 10 5 times bigger than J 3 . The terms, such as J 3 , J 21 , J 22 , J 5 , J 6 (the values and their uncertainty can be seen in [26]), can be neglected compared with the terms of J 2 and J 4 . Therefore, a good approximation of the gravity model of Jupiter is given by
U = μ r [ 1 J 2 2 ( R r ) 2 ( 3   sin 2 ϕ 1 ) J 4 8 ( R r ) 4 ( 35 sin 4 ϕ 30 sin 2 ϕ + 3 ) ] ,
where μ = G M J is the gravitational constant of Jupiter, G = 6.67428 × 10 11   m 3 k g 1 s 2 , M J 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. Therefore, 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
T = 1 2 ( r ˙ 2 + r 2 cos 2 ϕ λ ˙ 2 + r 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 d t ( T q ˙ ) T q = U q .
More precisely, the equations of motion can be presented as follows:
{ r ¨ r cos 2 ϕ λ ˙ 2 r ϕ ˙ 2 = μ r 2 + 3 μ J 2 R 2 2 r 4 ( 3 sin 2 ϕ 1 ) + 5 μ J 4 R 4 8 r 6 ( 35 sin 4 ϕ 30 sin 2 ϕ + 3 ) , d d t ( r 2 cos 2 ϕ λ ˙ ) = 0 , d d t ( r 2 ϕ ˙ ) + 1 2 r 2 sin ( 2 ϕ ) λ ˙ 2 = [ 3 μ J 2 R 2 2 r 3 + μ J 4 R 4 8 r 5 ( 70 sin 2 ϕ 30 ) ] sin 2 ϕ .
From the second equation of (4), we see that the quantity r 2 cos 2 ϕ λ ˙ 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 perturbations. The presence of these terms increases the radius of the stationary orbit with respect 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 r ˙ = r ¨ = 0 , λ ˙ = n J , λ ¨ = 0 , ϕ ˙ = ϕ ¨ = 0 in these equations; we get
{ n J 2 = μ r 3 + 3 μ J 2 R 2 2 r 5 15 μ J 4 R 4 8 r 7 , sin 2 ϕ = 0 ,
or
{ μ r 2 + 3 μ J 2 R 2 2 r 4 ( 3 sin 2 ϕ 1 ) + 5 μ J 4 R 4 8 r 6 ( 35 sin 4 ϕ 30 sin 2 ϕ + 3 ) + r cos 2 ϕ n J 2 = 0 , 1 2 r 2 n J 2 + 3 μ J 2 R 2 2 r 3 + μ J 4 R 4 8 r 5 ( 70 sin 2 ϕ 30 ) = 0 ,
where n J is the rotational angular velocity of Jupiter. One can verify that the left-hand side of the second equation in (6) is always positive when r is larger than R . Therefore, we only need to analyze solutions of Equation (5). From the second equation of (5), we see that only one meaningful latitude of the stationary orbit exists, i.e., ϕ = 0 ( ϕ = π 2 , π are also roots of sin ( 2 ϕ ) = 0 , but ϕ = π 2 makes no sense for stationary orbits, and ϕ = π corresponds to a stationary orbit which coincides with that of ϕ = 0 ). Introducing the functions f 0 ( r ) = μ r 3 + 3 μ J 2 R 2 2 r 5 15 μ J 4 R 4 8 r 7 n J 2 and f 1 ( r ) = μ r 3 n J 2 , the variations of f 0 ( r ) , f 1 ( r ) , with respect to r , can be seen in the following Figure 1.
It can be proved that f 0 ( r ) vanishes at two real positive values of r , but only one of them is bigger than R . This solution, which is given by r 0 = 2 . 2414 R = 1.6024 × 10 8   m , is denoted by E0 in Figure 1. 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.2381 R , which corresponds to the point E1. Therefore, the existence of the J 2 and J 4 terms increase the altitude of the stationary orbit.
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 r 2 cos 2 ϕ λ ˙ for some given initial conditions. Assuming that the deviations from the stationary orbits are small and writing r = r 0 ( 1 + ε ) , Equation (4) can be linearized as follows (Murray & Dermott [29], Section 11 of Chapter 6):
{ ε ¨ + ( 3 α 1 + α 2 ) ε = 0 , ϕ ¨ + ( α 1 + β 2 ) ϕ + β 1 = 0 , λ h r 0 2 ( 1 + ε ) 2 cos 2 ϕ = 0
where α 1 = 1 r 0 U r ( r 0 , 0 ) , α 2 = 2 U r 2 ( r 0 , 0 ) , β 1 = 1 r 0 2 U ϕ ( r 0 , 0 ) , and β 2 = 1 r 0 2 2 U 2 ϕ ( r 0 , 0 ) . By solving Equation (7), the analytical approximate expressions of ε , ϕ , λ can be formulated as [29]
{ ε = e cos k 1 t , ϕ = β 1 α 1 + β 2 + i cos k 2 t , λ = α 1 t 2 α 1 k 1 e sin k 1 t ,
where k 1 2 = 3 α 1 + α 2 , k 2 2 = α 1 + β 2 , i is the inclination, and e is the eccentricity of the orbit. Note that the average of λ ˙ is α 1 . Let us set k 3 = α 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 k 1 , k 2 , k 3 [29]:
{ k 1 2 = μ r 0 3 [ 1 3 2 J 2 ( R r 0 ) 2 + 45 8 J 4 ( R r 0 ) 4 ] , k 2 2 = μ r 0 3 [ 1 + 9 2 J 2 ( R r 0 ) 2 75 8 J 4 ( R r 0 ) 4 ] , k 3 2 = μ r 0 3 [ 1 + 3 2 J 2 ( R r 0 ) 2 15 8 J 4 ( R r 0 ) 4 ] .
Therefore, the satellites will oscillate in the radial and North-South directions with frequencies k 1 and k 2 , respectively. Moreover, the mean motion frequency, k 3 , 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. Figure 2 shows a stationary orbit and four orbits obtained from r = 0.99 r 0 , ϕ = ± 0.1 deg and r = 1.01 r 0 , ϕ = ± 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 J 2 and J 4 , the three frequencies are usually not equal or even commensurable. Figure 3 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.
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 ( a d ) and central gravity acceleration ( a m ) 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 d a m ) max = 2 m M J ( r 0 ρ ) 3 ,
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 ( a , e , i , Ω , ω , M ) of Io that we adopt here are taken from the 10th China Trajectory Optimization Competition, i.e., ( 422 , 029.687 km , 0.004308 , 0.04 deg , 79.64 deg , 37.991 deg , 4.818 deg ) . Calculation 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.
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 Jupiter is about 122,500 km, and the mass, m r , is about 1.0×1013 kg [32].
For satellites on stationary orbits with position r s , the gravitational acceleration due to the Jovian main ring can be calculated as
F r = G r i n g ( r ( Ρ ) r s ) ρ ( Ρ ) d V | r ( Ρ ) r s | 3 ,
where r ( Ρ ) , ρ ( Ρ ) denote the position vector and the density, respectively, which depend on the position of a point, Ρ , belonging to the ring. We can see that
| r ( Ρ ) r s | r 0 d w 160 , 240 122 , 500 6440 0.4378 R , | F r | G m r ( 0.4378 R ) 2 .
Note that G m r ( 0.4378 R ) 2 / ( G M J r 0 2 ) = m r M J ( r 0 0.4378 R ) 2 1.0 × 10 13 1.9 × 10 27 ( 2.2414 0.4378 ) 2 1.3795 × 10 13 , which is far smaller than J 2 and J 4 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. The magnetic field of Jupiter near stationary orbits when exterior terms are neglected can be written as
B ( r , ϕ , λ ) = V ( r , ϕ , λ ) ,
where
V = R l = 1 m = 1 l ( R r ) l + 1 [ g l , m cos ( m λ ) + h l , m sin ( m λ ) ] P l m ( cos ϕ ) .
Here, g l , m and h l , m are the geomagnetic Gauss coefficients (their values can be found in [33,34]), and P l m ( cos ϕ ) is the normalized Legendre function. The acceleration due to the Lorentz force is
F L = q m ( r ˙ ω × r ) × B ( r , ϕ , λ ) ,
where q = 4 π ε 0 Φ s , ε 0 = 8.854187817 × 10 12   F m 1 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 q = 5.5633 × 10 7 C and | F L | 1.14 × 10 8 N . The Lorentz force compared with the lowest order of gravity provides the ratio
| F L | | G 0 | = ( n n J ) B q r 3 G M m 2.3 × 10 11
This ratio is also far smaller than J 2 and J4. Therefore, it can be concluded that the effects of Lorentz forces on a satellite can be neglected.

3. Sun-Synchronous Orbits

The oblate nature of the primary body can lead to a secular variation of the ascending node of the orbit. However, we can use the orbit perturbations to keep the orientation of the Sun line direction fixed with respect to the orbital plane (for example, perpendicular to it) during one revolution of Jupiter around the Sun.
Based on the mean element theory, the mean nodal precession rate coming from the secular perturbations of the first and second order [36,37] can be described as
Ω . = 3 n J 2 R 2 2 a 2 ( 1 e 2 ) 2 cos i 9 n J 2 2 R 4 4 p 4 cos i { 3 2 + e 2 6 + 1 e 2 sin 2 i ( 5 3 5 e 2 24 + 3 2 1 e 2 ) 35 J 4 18 J 2 2 [ 6 7 + 9 e 2 7 sin 2 i ( 3 2 + 9 e 2 4 ) ] } .
For sun-synchronous orbits, the mean nodal precession rate is equal to the mean motion of Jupiter orbiting around the Sun ( n s ). Thus, we get the following equation:
g ( cos i ) = Γ 1 cos 3 i + Γ 2 cos i + Γ 3 = 0 ,
where
Γ 1 = 9 n J 2 2 R 4 4 p 4 [ 5 3 5 24 e 2 + 3 2 1 e 2 35 J 4 6 J 2 2 ( 1 2 + 3 4 e 2 ) ] ,
Γ 2 = 3 n J 2 R 2 2 p 2 { 1 + 3 J 2 R 2 2 p 2 [ 1 6 + 3 e 2 8 1 2 1 e 2 + 5 J 4 2 J 2 2 ( 1 2 + 3 e 2 4 ) ] } ,
Γ 3 = n s .
Here, we remark that in Equations (17) and (18), p is equal to a(1−e2).
First, setting x = cos i , we note that the equation g ( x ) = 0 has three real roots if the discriminant satisfies the following inequality:
Δ = ( Γ 3 2 4 Γ 1 2 + Γ 2 3 27 Γ 1 3 ) 0 .
Then, it is also necessary to avoid impact with Jupiter. However, the surface of Jupiter cannot be unambiguously defined since it is a gas giant. The critical perijovian distance, d c , is usually much larger than the radius of Jupiter, R , due to radiation safety issues. However, let us take for convenience d c = R , so that the semi-major axis and eccentricity have to satisfy
d c = R < a ( 1 e ) .
The variation of the discriminant with respect to a and e is shown in Figure 4a. We see that the discriminants are always negative for a [ R , 2 R ] and e [ 0 , 1 R a ] , therefore, it can be concluded that there usually exists one meaningful inclination for sun-synchronous orbits when the semi-major axis and eccentricity are given. The inclinations for different semi-major axis and eccentricities are shown in Figure 4b. It was shown that the inclination monotonically increases with respect to the semi-major axis. For the same semi-major axis, the inclination decreases as the eccentricity increases. Moreover, the inclinations with the same altitude ratio, a / R , are relatively smaller than those of near-Earth sun-synchronous orbits.
For satellites on sun-synchronous orbits with medium altitude, perturbations from the gravity of the third body are usually smaller compared to satellites on stationary orbits. However, the risk coming from the Jupiter rings and the magnetic field may increase greatly.
For sun-synchronous orbits, when a = 1.4373 R = 102 , 755.451 km and e = 0.4 , we have i = 90.183 deg; when a = 1.5308 R = 109 , 439.953 km and e = 0.1 , the corresponding inclination is 90.321 deg. The evolution of Ω and of the difference between Δ Ω and n s t over 25T for these two cases are presented in Figure 5. In the first case, this difference is about 0.018 deg; and in the second case, it is about 0.008 deg.

4. Orbits with Critical Inclination

The variations of the eccentricity and the argument of perijove are mainly caused by the equatorial bulge of Jupiter. These usually produce negative effects on space missions to Jupiter. However, we can choose orbits with critical inclinations to avoid these disadvantages.
According to the mean element theory, the mean variation rate of ω caused by perturbations of the first and second order [36,37] can be formulated as follows:
ω ˙ = 3 n J 2 R 2 2 a 2 ( 1 e 2 ) ( 5 2 sin 2 i 2 ) + 9 n J 2 2 R 4 p 4 { 4 + 7 e 2 12 + 2 1 e 2 sin 2 i ( 103 12 + 3 e 2 8 + 11 2 1 e 2 ) + sin 4 i ( 215 48 15 32 e 2 + 15 4 1 e 2 ) 35 J 4 18 J 2 2 [ 12 7 + 27 e 2 14 sin 2 i ( 93 14 + 27 e 2 4 ) + sin 4 i ( 21 4 + 81 e 2 16 ) ] } .
To keep the invariance of the mean argument of perijove, we set ω ˙ = 0 . Therefore, we obtain the following equation
Γ 4 sin 4 i + Γ 5 sin 2 i + Γ 6 = 0 ,
where
Γ 4 = 9 n J 2 2 R 4 4 p 4 [ 215 48 15 e 2 32 + 15 4 1 e 2 35 J 4 18 J 2 2 ( 21 4 + 81 e 2 16 ) ] ,
Γ 5 = 3 n J 2 R 2 4 p 2 [ 5 3 J 2 R 2 p 2 ( 103 12 + 3 e 2 8 + 11 2 1 e 2 ) + 35 J 4 6 J 2 2 ( 31 14 + 9 e 2 4 ) ] ,
Γ 6 = 3 n J 2 R 2 4 p 2 { 4 + 3 J 2 R 2 p 2 [ 4 + 7 e 2 12 + 2 1 e 2 5 J 4 6 J 2 2 ( 4 + 9 e 2 2 ) ] } .
Note that the left hand side of Equation (23) is a second-degree polynomial with respect to sin 2 i . This polynomial in sin 2 i may have two real roots in [0,1], and so we may have up to four values of the inclination in [ 0 , π ] that solve Equation (23). Let us set x = sin 2 i and write Equation (23) as
h ( x ) = Γ 4 x 2 + Γ 5 x + Γ 6 = 0 .
Noting that Γ 4 > 0 , Γ 5 < 0 , Γ 6 > 0 , we find that h ( x ) can have up to two real roots. Moreover, if h ( 1 ) > 0 , they will fall in the interval [0,1]. If there exist four roots in [ 0 , π ] of the critical inclination, then Equation (27) must have two different real roots lying in [0,1]. Observing that Γ 4 > 0 , Γ 5 < 0 , Γ 6 > 0 , we can easily see a necessary condition for Equation (27) to have two different roots lying in [0, 1] is
h ( 0 ) = Γ 6 > 0 , h ( 1 ) = Γ 4 + Γ 5 + Γ 6 > 0
Moreover, the condition (21) should also be satisfied to avoid the impact with Jupiter. The variation of h ( 1 ) for a [ 1.1 R , 9 R ] and e [ 0 , 1 R / a ] are presented in Figure 6a. We note that h ( 1 ) < 0 in this domain of ( a , e ) , therefore, there exist two critical inclinations, and one of them corresponds to a retrograde orbit. The variation of critical inclinations of direct orbits as a increases from R to 2 R are presented in Figure 6b for different values of the eccentricity. We can see that when e is small, for example e < 0.2 , the critical inclination monotonically decreases with respect to the altitude of the orbit. On the other hand, when e > 0.4 , the critical inclination monotonically increases. It can also be concluded that there exist some critical values for e, such that the critical inclination is independent of the altitude of the orbits.
Figure 6b shows that in the case of Jupiter, critical inclinations depend on the semi-major axis and eccentricity, while for the Earth considering only the J 2 , they are fixed (63.435 deg and 116.565 deg). The evolution of ω and e for two different values of the semi-major axis and eccentricity are shown in Figure 7 and Figure 8.
For the first case in Figure 7, the amplitude of ω is about 5 deg over a time of 25 T . The amplitude of e is about 0.013. For the second case, the amplitude of ω is about 1.8 deg over the same time. One can also see from Figure 8 that the amplitude of ω in the first case is much larger. This is reasonable, since the semi-major axis of the second case is larger than that of the first case.

5. Repeating Ground Track Orbits

Repeating ground track orbits are usually used for remote sensing satellites to obtain a good appearance for orbital covering. For satellites on these orbits, the ground tracks can repeat periodically. In this way, the area of Jupiter covered by the satellite can be dynamically monitored to detect the change of the target during this period of time. For example, we can make use of these orbits to detect the change of the Great Red Spot. The interval of adjacent ground track in the equator is
Δ λ = T N ( n J Ω ˙ )
where T N = 2 π M ˙ + ω ˙ is the nodal period of the motion of the spacecraft. Using the mean element theory, we have [36,37]
M ˙ = n + 3 n J 2 R 2 2 p 2 ( 1 3 2 sin 2 i ) 1 e 2 + 9 n J 2 2 R 4 4 p 4 1 e 2 { 1 2 ( 1 3 2 sin 2 i ) 2 1 e 2 + 5 2 + 10 e 2 3 sin 2 i ( 19 3 + 26 e 2 3 ) + sin 4 i ( 233 48 + 103 e 2 12 ) + e 4 1 e 2 ( 35 12 35 4 sin 2 i + 315 32 sin 4 i ) 35 J 4 18 J 2 2 e 2 ( 9 14 45 14 sin 2 i + 45 16 sin 4 i ) } ,
where n is the mean angular velocity. The conditions for repeating ground track orbits can be formulated as
D T N ( n J Ω ˙ ) = D Δ λ = 2 π N
where D and N are positive integers. Equation (31) implies satellites will have a repeating ground track after completing D revolutions in N Jovian days. The ground track repetition parameter, Q = D N , which represents the number of orbital resolutions in a Jovian day, is widely used in engineering practice. Using this parameter combined with the definition of T N , Equation (31) is equivalent to
Q = M ˙ + ω . n J Ω ˙ .
We can see that once Q is fixed, a, e and i will be related by Equation (32). The repeating ground track orbits, which are also sun-synchronous, may play an important role in planetary exploration. For these orbits, Equation (16) was also satisfied. The change of inclinations of sun-synchronous repeating ground track orbits, with respect to the semi-major axis for different values of Q, are presented in Figure 9. One can see that for a given Q, small variations of the semi-major axis produce huge changes in the inclination. Moreover, for orbits with the same semi-major axis, the inclination,i will be larger if Q is larger. For the Earth and Mars, the inclination is not so sensitive to the semi-major axis [9].

6. Conclusions

In this paper, we analyzed some special orbits around Jupiter considering only the J2 and J4 terms of the non-spherical gravitational potential. First, stationary orbits were investigated in spherical coordinates. The radius of stationary orbits was found to be 2.2414 Jupiter radii. The longitude drift due to the non-spherical component of the gravity field did not occur for the lack of J22 terms. Moreover, a small perturbation of r and ϕ would lead to the oscillation phenomenon in the radial and North-South directions. We also showed that the perturbations of the Sun and the Galilean satellites, the main ring and the magnetic field are negligible on short time scales compared to those represented by the terms J2 and J4. However, the obtained solutions may be incorrect on long time scales. Then, sun-synchronous orbits, orbits at critical inclinations, and repeating ground track orbits were discussed based on the mean-element theory. The results showed that only one meaningful inclination exists for sun-synchronous orbits when a and e are fixed in a suitable range. Additionally, only two critical inclinations exist for the critical orbits. Repeating ground track orbits, which are also sun-synchronous, were also calculated here. The traditional frozen orbits were not discussed because the J3 term was not taken into account in our gravity model. However, we remark that there may exist some non-traditional frozen orbits around Jupiter, which will be left for our future investigations.
The dynamics around Jupiter are very complicated and even if some perturbations can be neglected for preliminary orbit design, other effects should be considered. For example, orbits at low-to-moderate inclinations within a few Jovian radii may loiter inside Jupiter’s radiation belts. This can lead to considerable risk for the integrity of a satellite. However, we believe that our work can help to reduce or eliminate the need for station keeping in future space mission around Jupiter.

Author Contributions

Conceptualization, H.Z. and Y.L.; methodology, H.L. and Y.J.; software, Y.J.; validation, H.L., Y.J. and Y.L.; formal analysis, H.L.; investigation, H.Z.; resources, Y.J.; data curation, H.Z. and Y.L.; writing—original draft preparation, Y.L.; writing—review and editing, H.L. and Y.J.; visualization, H.Z. and Y.L.; supervision, H.L. and Y.J.; project administration, H.L. and Y.J.; funding acquisition, Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 11772356).

Acknowledgments

The authors gratefully acknowledge the reviewers for their careful work and thoughtful suggestions that have helped improve this paper substantially.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weibel, W.M.; Kaula, W.M.; Newman, W.I. A computer search for stable orbits between Jupiter and Saturn. Icarus 1990, 83, 382–390. [Google Scholar] [CrossRef]
  2. Jacobson, R.A. The Gravity Field of the Jovian System and the Orbits of the Regular Jovian Satellites. Bull. Am. Astron. Soc. 2001, 33, 11–101. [Google Scholar]
  3. Colwell, J.E.; Horányi, M.; Grün, E. Jupiter’s exogenic dust ring. J. Geophys. Res. Planets 1998, 103, 20023–20030. [Google Scholar] [CrossRef]
  4. Liu, X.; Schmidt, J. Dust in the Jupiter system outside the rings. Astrodynamics 2019, 3, 17–29. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, X.; Sachse, M.; Spahn, F.; Schmidt, J. Dynamics and distribution of Jovian dust ejected from the Galilean satellites. J. Geophys. Res. Planets 2016, 121, 1141–1173. [Google Scholar] [CrossRef] [Green Version]
  6. Clarke, A.C. Extra-Terrestrial Relays: Can Rocket Stations Give World-wide Radio Coverage?—ScienceDirect. Prog. Astronaut. Rocket. 1966, 19, 3–6. [Google Scholar]
  7. Musen, P.; Bailie, A.E. On the motion of a 24-h satellite. J. Geophys. Res. 1962, 67, 1123–1132. [Google Scholar] [CrossRef]
  8. Lara, M.; Elipe, A. Periodic Orbits Around Geostationary Positions. Celest. Mech. Dyn. Astron. 2002, 82, 285–299. [Google Scholar] [CrossRef]
  9. Liu, X.; Baoyin, H.; Ma, X. Five Special Types of Orbits Around Mars. J. Guid. Control Dyn. 2011, 33, 1294–1301. [Google Scholar] [CrossRef]
  10. Liu, X.D.; Baoyin, H.; Ma, X.R. Periodic orbits around areostationary points in the Martian gravity field. Res. Astron. Astrophys. 2012, 12, 551–562. [Google Scholar] [CrossRef] [Green Version]
  11. Orlov, A.A. Obituary: Pochti krugovye periodicheskie dvizheniia materialnoi tochki pod deistviem niutonovskogo pritiazheniia sferoida. Soobshcheniya Gosudarstvennogo Astron. Inst. 1953, 88, 3–38. [Google Scholar]
  12. Brouwer, D. Proceedings of the Celestial Mechanics Conference: Outlines of general theories of the Hill-Brown and Delaunay types for orbits of artificial satellites. Astron. J. 1958, 63, 433. [Google Scholar] [CrossRef]
  13. Jupp, A.H. The problem of the critical inclination revisited. Celest. Mech. 1975, 11, 361–378. [Google Scholar] [CrossRef]
  14. Blitzer, L. Synchronous and Resonant Satellite Orbits Associated With Equatorial Ellipticity. Ars J. 1962, 32, 1016–1019. [Google Scholar] [CrossRef]
  15. Gedeon, G.S. Tesseral resonance effects on satellite orbits. Celest. Mech. Dyn. Astron. 1969, 1, 167–189. [Google Scholar] [CrossRef]
  16. Alessi, E.M.; Buzzoni, A.; Daquin, J.; Carbognani, A.; Tommei, G. Dynamical properties of the Molniya satellite constellation: Long-term evolution of orbital eccentricity—ScienceDirect. Acta Astronaut. 2020, 179, 659–669. [Google Scholar] [CrossRef]
  17. Daquin, J.; Alessi, E.M.; O’Leary, J.; Lemaitre, A.; Buzzoni, A. Dynamical Properties of the Molniya Satellite Constellation: Long-Term Evolution of the Semi-Major Axis. arxiv 2021, arXiv:2103.06251. [Google Scholar]
  18. Cutting, E.; Frautnick, J.C.; Born, G.H. Orbit analysis for Seasat-A. J. Astronaut. Sci. 1978, 26, 315–342. [Google Scholar]
  19. Coffey, S.L.; Deprit, A.; Deprit, E. Frozen orbits for satellites close to an Earth-like planet. Celest. Mech. Dyn. Astron. 1994, 59, 37–72. [Google Scholar] [CrossRef]
  20. Folta, D.; Quinn, D. Lunar Frozen Orbits. In Proceedings of the AIAA/AAS Astrodynamics Specialist Conference and Exhibit, Keystone, CO, USA, 21–24 August 2006. [Google Scholar]
  21. Nie, T.; Gurfil, P. Lunar frozen orbits revisited. Celest. Mech. Dyn. Astron. 2018, 130, 61. [Google Scholar] [CrossRef]
  22. Tresaco, E.; Elipe, A.; Carvalho, J.P. Frozen Orbits for a Solar Sail Around Mercury. J. Guid. Control Dyn. 2016, 39, 1659–1666. [Google Scholar] [CrossRef]
  23. Macdonald, M.; Mckay, R.; Vasile, M.; Frescheville, F.D. Extension of the Sun-Synchronous Orbit. J. Guid. Control Dyn. 2015, 33, 1935–1940. [Google Scholar] [CrossRef] [Green Version]
  24. Lara, M. Repeat Ground Track Orbits of the Earth Tesseral Problem as Bifurcations of the Equatorial Family of Periodic Orbits. Celest. Mech. Dyn. Astron. 2003, 86, 143–162. [Google Scholar] [CrossRef]
  25. Lei, H. Dynamical models for secular evolution of navigation satellites. Astrodynamics 2019, 4, 57–73. [Google Scholar] [CrossRef]
  26. Iess, L.; Folkner, W.M.; Durante, D.; Parisi, M.; Bolton, S.J. Measurement of Jupiter’s asymmetric gravity field. Nature 2018, 555, 220–222. [Google Scholar] [CrossRef] [PubMed]
  27. Lemoine, F.G.; Kenyon, S.C.; Factor, J.K.; Trimmer, R.G.; Pavlis, N.K.; Chinn, D.S.; Cox, C.M.; Klosko, S.M.; Luthcke, S.B.; Torrence, M.H.; et al. The Development of the Joint NASA GSFC and NIMA Geopotential Model EGM96; NASA TM 1998-206861; MIT Press: Cambridge, MA, USA, 1998. [Google Scholar]
  28. Lemoine, F.G.; Smith, D.E.; Rowlands, D.D.; Zuber, M.T.; Neumann, G.A.; Chinn, D.S.; Pavlis, D.E. An improved solution of the gravity field of Mars (GMM-2B) from Mars Global Surveyor. J. Geophys. Res. Planets 2001, 106, 23359–23376. [Google Scholar] [CrossRef]
  29. Murray, C.D.; Dermott, S.F. Solar System Dynamics; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  30. Garmon, J. Moons of the planet Jupiter; College Publishing House: Darya Ganj, IN, USA, 2012. [Google Scholar]
  31. Zhang, H.B. Theories and Methods Of Spacecraft Orbital Mechanics; National Defense Industry Press: Beijing, China, 2015. (In Chinese) [Google Scholar]
  32. Available online: https://solarviews.com/eng/jupiter.htm (accessed on 8 May 2021).
  33. Connerney, J.; Acua, M.H.; Ness, N.F. Voyager 1 assessment of Jupiter’s planetary magnetic field. J. Geophys. Res. Atmos. 1982, 87, 3623–3627. [Google Scholar] [CrossRef] [Green Version]
  34. Connerney, J.; Acuña, M.; Ness, N.F.; Satoh, T. New models of Jupiter’s magnetic field constrained by the Io flux tube footprint. J. Geophys. Res. Space Phys. 1998, 103, 11929–11939. [Google Scholar] [CrossRef]
  35. Garrett, H.B.; Katz, I.; Jun, I.; Kim, W.; Whittlesey, A.C.; Evans, R.W. The Jovian Charging Environment and Its Effects—A Review. Plasma Sci. 2011, 40, 144–154. [Google Scholar] [CrossRef]
  36. Brouwer, D. Solution of the problem of artificial satellite theory without drag. Astron. J. 1959, 64, 378. [Google Scholar] [CrossRef]
  37. Liu, L. Orbit Theory of Spacecraft; National Defense Industry Press: Beijing, China, 2000. (In Chinese) [Google Scholar]
Figure 1. Variations of f 0 ( r ) , f 1 ( r ) as r increases from 2.2 R to 2.3 R .
Figure 1. Variations of f 0 ( r ) , f 1 ( r ) as r increases from 2.2 R to 2.3 R .
Aerospace 08 00183 g001
Figure 2. Orbits around Jupiter: green—stationary orbits; red— r = 0.99 r 0 , ϕ = 0.1 deg; blue— r = 0.99 r 0 and ϕ = 0.1 deg; cyan— r = 1.01 r 0 and ϕ = 0.1 deg; black— r = 1.01 r 0 and ϕ = 0.1 deg.
Figure 2. Orbits around Jupiter: green—stationary orbits; red— r = 0.99 r 0 , ϕ = 0.1 deg; blue— r = 0.99 r 0 and ϕ = 0.1 deg; cyan— r = 1.01 r 0 and ϕ = 0.1 deg; black— r = 1.01 r 0 and ϕ = 0.1 deg.
Aerospace 08 00183 g002
Figure 3. (a,b) Variations of orbital radius ( r ), (c,d) drift of longitude ( Δ λ ), (e,f) evolution of latitude ( ϕ ): green—stationary orbits; red, r = 0.99 r 0 and ϕ = 0.1 deg; blue— r = 0.99 r 0 and ϕ = 0.1 deg; cyan— r = 1.01 r 0 and ϕ = 0.1 deg; black— r = 1.01 r 0 , ϕ = 0.1 deg, and T is a Jovian day.
Figure 3. (a,b) Variations of orbital radius ( r ), (c,d) drift of longitude ( Δ λ ), (e,f) evolution of latitude ( ϕ ): green—stationary orbits; red, r = 0.99 r 0 and ϕ = 0.1 deg; blue— r = 0.99 r 0 and ϕ = 0.1 deg; cyan— r = 1.01 r 0 and ϕ = 0.1 deg; black— r = 1.01 r 0 , ϕ = 0.1 deg, and T is a Jovian day.
Aerospace 08 00183 g003aAerospace 08 00183 g003b
Figure 4. (a) Variation of the discriminant with respect to a and e ; (b) variations of the inclination i as the semi-major axis a increases from R to 2 R .
Figure 4. (a) Variation of the discriminant with respect to a and e ; (b) variations of the inclination i as the semi-major axis a increases from R to 2 R .
Aerospace 08 00183 g004
Figure 5. Evolution of Ω and difference between Δ Ω = Ω Ω 0 and n s t over 25 T , where T is a Jovian day, Ω 0 = 60 deg is the initial condition (the blue curve was indeed obtained by numerical integration of osculating elements). (a) Evolution of Ω for a = 102 , 755.451 km, i = 90.183 deg. (b) Difference between Δ Ω and n s t for a = 102 , 755.451 km, e = 0.4 , i = 90.183 deg. (c) Evolution of for a = 109 , 439.953 km, e = 0.1 , i = 90.321 deg. (d) Difference between Δ Ω and n s t for a = 109 , 439.953 km, e = 0.1 , i = 90.321 deg.
Figure 5. Evolution of Ω and difference between Δ Ω = Ω Ω 0 and n s t over 25 T , where T is a Jovian day, Ω 0 = 60 deg is the initial condition (the blue curve was indeed obtained by numerical integration of osculating elements). (a) Evolution of Ω for a = 102 , 755.451 km, i = 90.183 deg. (b) Difference between Δ Ω and n s t for a = 102 , 755.451 km, e = 0.4 , i = 90.183 deg. (c) Evolution of for a = 109 , 439.953 km, e = 0.1 , i = 90.321 deg. (d) Difference between Δ Ω and n s t for a = 109 , 439.953 km, e = 0.1 , i = 90.321 deg.
Aerospace 08 00183 g005aAerospace 08 00183 g005b
Figure 6. (a) Variation of h ( 1 ) as a function of the semi-major axis a and the eccentricity e . (b) Variation of the critical inclination i as a function of the semi-major axis for different values of the eccentricity.
Figure 6. (a) Variation of h ( 1 ) as a function of the semi-major axis a and the eccentricity e . (b) Variation of the critical inclination i as a function of the semi-major axis for different values of the eccentricity.
Aerospace 08 00183 g006
Figure 7. Evolution of ω and e over 25 T , (a1,a2) a = 1.6832 R = 120 , 335.334 km, i = 63.497 deg; (b1,b2) a = 2.1488 R = 153 , 622.009 km, i = 63.431 deg, where T is a Jovian day.
Figure 7. Evolution of ω and e over 25 T , (a1,a2) a = 1.6832 R = 120 , 335.334 km, i = 63.497 deg; (b1,b2) a = 2.1488 R = 153 , 622.009 km, i = 63.431 deg, where T is a Jovian day.
Aerospace 08 00183 g007
Figure 8. e ~ ω evolution over 25 T , red points correspond to a = 1.6832 R = 120 , 335.334 km and i = 63.497 deg; green points correspond to a = 2.1488 R = 153 , 622.009 km and i = 63.431 deg.
Figure 8. e ~ ω evolution over 25 T , red points correspond to a = 1.6832 R = 120 , 335.334 km and i = 63.497 deg; green points correspond to a = 2.1488 R = 153 , 622.009 km and i = 63.431 deg.
Aerospace 08 00183 g008
Figure 9. Inclinations of sun-synchronous repeating ground track orbits as a function of the semi-major axis, for different values of the repetition parameter Q.
Figure 9. Inclinations of sun-synchronous repeating ground track orbits as a function of the semi-major axis, for different values of the repetition parameter Q.
Aerospace 08 00183 g009
Table 1. Some harmonic coefficients and their ratios of the Earth, Mars, and Jupiter.
Table 1. Some harmonic coefficients and their ratios of the Earth, Mars, and Jupiter.
Planet J 2 ( × 10 3 ) J 3 ( × 10 6 ) J 4 ( × 10 5 ) J 22 ( × 10 6 ) J 3 / J 2 ( × 10 3 ) J 4 / J 2 ( × 10 2 )
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
Table 2. The disturbing acceleration due to the third body.
Table 2. The disturbing acceleration due to the third body.
The Third Body m / M J ( × 10 5 ) ρ / r 0 ( a d / a m ) max ( × 10 6 )
Sun ([32]) 1.0473 × 10 8 4856.80.01823
Io ([30,32])4.70472.6311 5.20000
Europa ([30,32])2.52834.18680.68904
Ganymede ([30,32])7.80566.67750.52379
Callisto ([30,32])5.667311.75110.06989
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, Y.; Jiang, Y.; Li, H.; Zhang, H. Some Special Types of Orbits around Jupiter. Aerospace 2021, 8, 183. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace8070183

AMA Style

Liu Y, Jiang Y, Li H, Zhang H. Some Special Types of Orbits around Jupiter. Aerospace. 2021; 8(7):183. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace8070183

Chicago/Turabian Style

Liu, Yongjie, Yu Jiang, Hengnian Li, and Hui Zhang. 2021. "Some Special Types of Orbits around Jupiter" Aerospace 8, no. 7: 183. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace8070183

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

Article Metrics

Back to TopTop