Next Article in Journal
Pressure Sensors: Working Principles of Static and Dynamic Calibration
Previous Article in Journal
Assessment of Directional–Hemispherical Reflectance of Tablets with Cefuroxime during Storage under Elevated Temperature and Ultraviolet Radiation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deployment of Constellation with Different Inclinations Using the Nodal Precession and Thrust

Graduate School, Space Engineering University, Beijing 101400, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(2), 631; https://doi.org/10.3390/s24020631
Submission received: 12 December 2023 / Revised: 6 January 2024 / Accepted: 16 January 2024 / Published: 19 January 2024
(This article belongs to the Section Remote Sensors)

Abstract

:
Strategy selection is critical for constellation deployment missions, both in terms of energy consumption and time cost. The different effects of impulse thrust and continuous thrust on orbit elements lead to a different choice of strategy. With impulse thrust, constellation types are differentiated according to high and medium-low inclinations. Constellations with high inclination are deployed using a strategy that controls the inclination. Constellations with medium-low inclination are deployed using a strategy that controls the semi-long axis. With continuous thrust, constellations are classified according to high, medium, and low inclination. High inclination constellations are deployed with a strategy of controlling inclination. Medium inclination constellations are deployed with a strategy that controls the semi-long axis. Low inclination constellations are deployed with a strategy of directly applying continuous thrust.

1. Introduction

Satellite constellations are satellite systems consisting of multiple satellites with stable geometrical configurations of space orbits and fixed spatial and temporal relationships between the satellites to fulfil specific space missions, which can give greater play to the role of satellites and expand the forms of satellite applications [1]. It has advantages that are incomparable to those of a single satellite in terms of global coverage, multiplicity, timeliness, and continuity, but at the same time, its construction costs are also incomparable to those of a single satellite, especially at the deployment phase. At present, with the miniaturization of satellites and the enhancement of space maneuver capabilities, the deployment of satellite constellations has gradually become more flexible and diversified, and the existing methods can be divided into two categories: the direct method and the indirect method.
The Direct Method means that satellites are deployed using a direct launch into the target orbit. The satellite is delivered directly into the target orbital plane by a launch vehicle, and the minimum number of launches required is greater than or equal to the number of orbital planes of the constellation, that is, at least one launch is executed per orbital plane [2,3]; this method has been used in many traditional constellations, such as Iridium and Galileo [4]. The huge cost is a key factor determining how successful the plan can be [5]; frequent rocket launches require large capital and time costs, and the traditional constellation deployment method is difficult to apply to large constellations. To address the problem of deploying large or even mega constellations, there is a need to explore new approaches that save both time and energy.
The Indirect Method means that the satellite separates the plane indirectly and then enters the target orbit. A single launch vehicle is used to launch a group of satellites in an initial orbital plane, and subsequently orbital plane separation is achieved as the satellites maneuver to the target plane. The orbital plane separation here mainly refers to the separation of RAAN (Right Ascension of the Ascending Node), and there are two commonly used means at present. One is to apply thrust to change the RAAN, which is to directly change the RAAN by out-of-plane maneuvering to push the satellite into a different orbital plane [6,7,8]. The second is to change the RAAN by controlling the semi-major axis, a process that takes advantage of the different rates of RAAN progression at different semi-major axis due to the non-spherical uptake of the J2 term to achieve an indirect separation of the planes.
The traditional constellation deployment method is difficult to apply to large constellations because of the large capital and time costs associated with frequent rocket launches. To address the problem of deploying large or even mega constellations, it is necessary to explore new approaches that save both time and energy. Therefore, the indirect method is the research focus of this paper. The development of the indirect method has been facilitated by advances in space propulsion technology, such as efficient miniaturized electric thrusters [9,10] and the exploitation of gravitational perturbation. Innovations in methodology and theory reflected in engineering will provide a significant reduction in launch costs. In exceptional cases, a single launch was sufficient to deploy a complete constellation, such as the FORMOSAT-3/COSMIC mission [11]. The mission deployed six satellites in different planes and took about 20 months in total. In 1993, King and Beidleman [12] first patented the use of the Earth’s gravitational perturbation for the deployment of multi-plane satellite constellations. All satellites are first put into a parking orbit. Satellites belonging to the same target orbital plane are placed in a “pallet” and delivered to their target orbits by a carrier. After a certain period, when the RAAN difference between the two orbital planes reaches the desired value due to the Earth’s gravitational perturbation, the second “pallet” is released from the parking orbit, and the process is repeated until all the “pallets” have been deployed. This method was analyzed for a constellation of small satellites in low orbit in [13], where Crisp also considered the effect of atmospheric drag on total mission failure. This method of deploying a multi-plane satellite constellation using RAAN progression is also used in [14]. One satellite changes its semi-major axis while the other satellite maintains its semi-major axis using atmospheric drag compensation with electric thrusters. Eventually a 180° RAAN separation was achieved between the two satellites. Leppinen [15] proposed a similar indirect approach to deploy a constellation of nanosatellites launched on a single carrier rocket to multiple different orbit planes.
Yang et al. [16] proposed a method to deploy a constellation of satellites using a combination of thrust and gravitational perturbation. McGrath et al. used the same strategy of orbit plane separation to solve the deployment problem of circular orbit constellations using a resolved approach [17]. Both Gomez and Yang proposed a combination of thrust and gravitational perturbation for the indirect method of constellation deployment, but the tangential thrust they use is for performing orbital semi-major axis holding of satellites in parked orbits and does not have a direct effect on RAAN separation. Giuseppe [18] simultaneously changes the semi-major axis and inclination to increase the effect of RAAN progression on the target orbit plane, improving the existing deployment method in terms of energy optimality. Zhang Yulin [19], in performing RAAN reconstruction of different orbit planes, similarly noticed that the different inclinations of the orbit planes would lead to different rates of their RAAN progression. Based on the impulse maneuver theory, he found the inclination boundaries for RAAN reconfiguration of different orbital planes: for high inclination orbits, RAAN separation is achieved by using inclination separation; for medium and low inclination orbits, RAAN separation is achieved by using orbital semi-major axis separation. It is necessary to note that, with the development of space electric propulsion, the thrust of electric propulsion can be up to 3 mN–1 N, and the specific impulse can be up to 800–3000 s [20,21]. using continuous thrust to directly change the satellite orbit plane of the RAAN can save a lot of time, but the corresponding energy consumption is pending to be explored.
This paper expands upon the theory of impulse thrust to include continuous thrust and derives an inclination boundary for the efficient separation of RAAN in both thrust modes. It is shown that the velocity increment required to apply continuous thrust to directly change the RAAN is related to inclination. Then, inclination boundary is derived to distinguish between using the method of applying thrust or controlling other orbit elements to achieve RAAN separation. Section 2 derives the relevant mathematical formulas, including the effect of J2 perturbation on orbit elements, the equations of motion under continuous thrust, and the velocity increments required to change different orbit elements for both thrust modes. Section 3 analyzes the deployment strategies under impulse thrust and continuous thrust conditions. Experimental simulations are performed in Section 4 to validate the proposed strategy theory.

2. Problem Formulation

2.1. Effects of J2 Perturbation

Neglecting the effect of the short-periodic term, the effect of the J2 perturbation [22,23] on the classical orbital elements is as follows.
d a d t = 0 d e d t = 0 d i d t = 0 d Ω d t = 3 J 2 R e 2 2 p 2 n cos i d M d t = 3 J 2 R e 2 2 p 2 n 2 5 2 sin 2 i d ω d t = 3 J 2 R e 2 2 p 2 n 1 3 2 sin 2 i 1 e 2
where a is the semi-major axis, e is the eccentricity, i is the inclination, Ω is the RAAN, M is the mean anomaly, ω is the argument of perigee, n is the angular velocity, p = a(1 − e2) and p is the semilatus rectum of the conic curve. Under the J2 perturbation, the semi-major axis, eccentricity, and inclination do not change, and the change in RAAN is related to the semi-major axis, eccentricity, and inclination.
The control of eccentricity is complicated and energy consuming, which is beyond the study scope of this paper. In this paper, we focus on analyzing the influence of semi-major axis and inclination on the variation law of RAAN. Under the non-spherical J2 perturbation of the Earth, the change rate of RAAN is derived from the semi-major axis and inclination as follows.
Ω ˙ J 2 a = Δ Ω ˙ J 2 Δ a = 21 J 2 R e 2 4 a p 2 n cos i
Ω ˙ J 2 i = Δ Ω ˙ J 2 Δ i = 3 J 2 R e 2 2 p 2 n sin i
For near-circular orbits, there are:
Δ a = 4 a 3 Δ Ω ˙ J 2 21 J 2 R e 2 n cos i
Δ i = 2 a 2 Δ Ω ˙ J 2 3 J 2 R e 2 n sin i

2.2. Equations of Motion for Continuous Thrust

The orbital equation of motion under continuous thrust is the key theory for constellation deployment. The perturbation equation [24] under continuous thrust is as follows.
d a d t = 2 a 2 h e f R sin ν + p r f C d e d t = 1 h p f R sin ν + f C ( p + r ) cos ν + r e d i d t = r cos ϑ h f N d Ω d t = r sin ϑ h sin i f N d ω d t = 1 h e p f R cos ν + p + r f C sin ν r sin ϑ h sin i f N cos i
where ϑ = ν + ω, ν is the true proximity angle, h is the orbital angular momentum, and μ is the Earth’s gravitational parameter. The thrust (fR, fC, fN) applied by the electric thrusters can be expressed in the RCN (Radial Circumferential Normal) reference system. The radial component fR is in the orbit plane, along the zenith; the normal component fN is perpendicular to the orbit plane, in the same direction as R × V ; and the circumferential component fC is in the orbit plane, perpendicular to the radial vector in the orbit plane.
Under the influence of continuous thrust, there is a complex coupling relationship between the orbit elements. As can be seen from Equation (6), the variation of the semi-major axis of the orbit is directly related to the radial component fR and the circumferential component fC; the variation of the inclination and the RAAN is directly related to the normal component fN. The direct correlation with the semi-major axis implies that there is some correlation between inclination and RAAN with both the radial component fR and the circumferential component fC. The eccentricity and argument of perigee are also affected by forces in all three directions. As shown in Figure 1, the thrust acceleration vector f is used to describe the applied thrust. Introducing an azimuth angle α (−π < α < π) and a pitch angle β (−π/2 < β < π/2), the acceleration component is expressed as:
f R   =   f cos β sin α f C   =   f cos β cos α f N   =   f sin β
From this, the equations of motion for the orbit after the introduction of the thrust acceleration vector can be obtained.
d a d t = 2 a 2 h f cos β e sin ν sin α + p r cos α d e d t = 1 h f cos β p sin ν sin α + ( p + r ) cos ν + r e cos α d i d t = f r h cos ϑ sin β d Ω d t = f r h sin ϑ sin i sin β d ω d t = 1 h e f cos β p cos ν sin α + p + r sin ν sin α r sin ϑ cos i h sin i sin β
Equation (8) represents the analytical relationship between the thrust angle (α, β) and the instantaneous change in each orbit element. Observing Equation (8), it can be observed that if the normal thrust does not change direction during an orbital period, the change in inclination and RAAN is almost 0. To maximize the change in inclination or RAAN, the normal thrust needs to change direction at a phase-symmetric position (ϑ = ±π/2 or ϑ = 0, π). For near-circular orbits, Bruce [25] and Zhao [26] gives equations for the control of the continuous thrust that changes direction at ϑ = ±π/2 and ϑ = 0, π.
(9) The   normal   thrust   changes   sign   at   ϑ = ± π / 2 d a d t   =   f 2 a v cos β cos α d i d t   =   f 2 π v sin β d Ω d t = 0 d ϑ d t = v 3 μ
(10) The   normal   thrust   changes   sign   at   ϑ = 0 ,   π . d a d t   =   f 2 a v cos β cos α d i d t = 0 d Ω d t   =   f 2 π v sin i sin β d ϑ d t = v 3 μ

2.3. Velocity Increment Required to Change the Orbit Elements

Orbital parameters can be changed using both impulsive and continuous thrusts. King [12] and A. Ruggiero [24] derived the velocity increments required to change different orbit elements using impulsive and continuous thrusts for near-circular orbital conditions.
Velocity   increment   with impulse   thrust   Δ V impulse (11) Δ V i m p u l s e a = 1 2 μ a 1 a Δ a (12) Δ V i m p u l s e i = μ a Δ i (13) Δ V i m p u l s e R A A N = μ a Δ Ω
Equation (11) represents the velocity increment required to change the semi-major axis by impulse thrust, Equation (12) represents the velocity increment required to change the inclination by impulse thrust, and Equation (13) represents the velocity increment required to change the RAAN deficit by impulse thrust.
Figure 2 shows the velocity increments required to change the semi-major axis and inclination of the spacecraft at different inclinations under impulse thrust conditions. Observing the solid blue line, it can be found that when the inclination is close to 90° (perpendicular to the equatorial plane), a large velocity increment is required to change the semi-major axis, while a small amount of energy is needed to change the inclination; and observing the solid yellow line, it can be found that when the inclination is close to 0 or 180° (parallel to the equatorial plane), a large velocity increment is required to change the inclination of the orbit, while a small amount of energy is needed to change the semi-major axis of the orbit. The velocity increments required to change the semi-major axis and inclination using impulse thrust have distinct inclination boundaries, and a detailed derivation of the process will be given in the next section. Spacecraft are capable of carrying not only impulse thrust loads but also continuous thrust loads, and the velocity increments required to change the orbital elements under continuous thrust conditions are given below.
Velocity   increment   with continuous   thrust   Δ V const (14) Δ V c o n s t a = 1 2 μ a 1 a Δ a (15) Δ V c o n s t i = π 2 μ a Δ i (16) Δ V c o n s t R A A N = π 2 μ a sin i Δ Ω
Equation (14) represents the velocity increment required to change the semi-major axis for continuous thrust, Equation (15) represents the velocity increment required to change the inclination for continuous thrust, and Equation (16) represents the velocity increment required to change the RAAN deficit for continuous thrust.
Figure 3 shows the velocity increments required to change the semi-major axis and inclination for the spacecraft at different inclinations under continuous thrust. Comparison of Equations (13) and (16) shows that the velocity increment required to change the RAAN under continuous thrust is also affected by the inclination. As can be seen in Figure 3, there is an intersection of the green solid line with the blue solid line. In other words, when the inclination is small enough, it is more energy efficient to apply thrust directly to control the RAAN.
For changing the semi-major axis of the orbit, the velocity increment required is the same for impulse thrust and continuous thrust. However, for changing the inclination of the orbit, the velocity increment required to use continuous thrust will be π/2 times the velocity increment required to use impulse thrust. In addition, for changing the RAAN, the situation is more complicated, and the velocity increment required to change the RAAN is related to the inclination.
Δ V c o n s t R A A N < Δ V i m p u l s e R A A N 0 i < 39.5402 ,   140.4598 < i 180 Δ V c o n s t R A A N = Δ V i m p u l s e R A A N i = 39.5402 ,   i = 140.4598 Δ V c o n s t R A A N > Δ V i m p u l s e R A A N 39.5402 < i < 140.4598
Compared to the velocity increment required to change the orbital inclination, applying continuous thrust to change the RAAN requires an additional term sin(i) in the equation, as can also be seen from Equation (4) in Equation (8). Here it is important to note that in terms of the properties of the sinusoidal function, the direct application of thrust to change the RAAN at low inclination will be more fuel and cost efficient than the velocity increments required by the other strategies.

3. Multi-Satellite Multiplane Deployment Strategy

3.1. Strategy of RAAN Separation

The simultaneous deployment of multiple satellites in different orbital planes is the main mode of deployment of mega-constellations, and the capability of launching multiple satellites with a single rocket makes this deployment mode feasible. The essence of multi-satellite heterodyne deployment is to realize RAAN separation between multiple satellites, and the selection of an appropriate RAAN separation strategy is the key to constellation deployment. As shown in Figure 4, the RAAN separation strategy is classified into “Indirect control of RAAN” and “Direct control of RAAN”.
The strategy of “indirect control of RAAN” can be achieved by changing the semi-major axis or inclination as shown in Figure 4a,b. The strategy of “direct control of the RAAN” is to apply a direct normal thrust to the spacecraft, as shown in Figure 4c. Whether impulse or continuous thrust is used, if the method of “indirect RAAN separation by changing orbit elements” is used, the constellation would be deployed using the “group delay” method. The launch vehicle puts all the satellites into a parking orbit, divides all the satellites into three groups, and then conducts the first maneuver, the second maneuver and the third maneuver, with RAAN separations generated during the maneuver time intervals. The schematic of the group delay transfer strategy is shown in Figure 5.

3.2. Inclination Boundary with Impulse Thrust

Substitute Equations (4) and (5) into Equations (11) and (12).
Δ V i m p u l s e a = 1 2 μ a 4 a 2 Δ Ω ˙ J 2 21 J 2 R e 2 n cos i
Δ V i m p u l s e i = μ a 2 a 2 Δ Ω ˙ J 2 3 J 2 R e 2 n sin i
Equations (18) and (19) indicate that an impulse thrust Δ V i m p u l s e a is required to achieve a RAAN drift rate Δ Ω ˙ J 2 by adjusting the semi-major axis, and an impulse thrust Δ V i m p u l s e i is required by adjusting the inclination of the orbit. Comparing the velocity increments required by the two methods, there is
Δ V i m p u l s e a Δ V i m p u l s e i = sin i 7 cos i
When tan i = 7 , Δ V i m p u l s e a / Δ V i m p u l s e i = 1 , so i = 81.8699° or 98.1301°. As shown in Figure 6, when 0 ≤ i < 81.87° or 98.13° < i < 180°, adjusting the semi-major axis to separate the orbit planes requires a smaller velocity increment. However, when 81.8699° ≤ i ≤ 98.1301°, it is more fuel efficient to complete the separation of the orbit planes by adjusting the inclination. In practical constellation configurations, the difference between the RAANs of different orbital planes is generally large, and it can be seen from Equation (13) that directly changing the RAAN by applying a normal thrust requires a large velocity increment, so this strategy is not considered here.

3.3. Inclination Boundary with Continuous Thrust

Substitute Equations (4) and (5) into Equations (14) and (15).
Δ V c o n s t a = 1 2 μ a 1 a 4 a 3 Δ Ω ˙ J 2 21 J 2 R e 2 n cos i
Δ V const i = π 2 μ a 2 a 2 Δ Ω ˙ J 2 3 J 2 R e 2 n sin i
Δ V const a Δ V const i = 2 sin i 7 π cos i
Same analytical approach as in Section 3.1, when Δ V c o n s t a / Δ V c o n s t i = 1 , i = 84.8035° or 95.1965°. Therefore, when 84.8035° ≤ i ≤ 95.1965°, it is more fuel efficient to complete the separation of the orbit planes by adjusting the inclination.
From the fourth equation of Equation (8), it can be seen that the efficiency of changing the RAAN by normal continuous thrust would increase for low inclination. Section 2.3 also analyzes that for low inclination, compared with impulse thrust, the velocity increment required to control the RAAN with continuous thrust will decrease.
According to the characteristics of continuous thrust dynamics [26], the total velocity increment required from a0 to af is Δ V c o n s t a .
Δ V c o n s t a = μ a 0 μ a f
where a0 is the semi-major axis of the parking orbit and af is the semi-major axis of the target orbit. Equation (24) is equivalent to Equation (14). Meanwhile, Equation (16) can be expressed as
Δ V c o n s t R A A N = π 2 μ a f sin i Δ Ω
By comparing Equations (24) and (25), we can analyze how to more efficiently change the RAAN under continuous thrust. Observation of the above equation shows that there exists a critical inclination iC that enables Δ V c o n s t a / Δ V c o n s t R A A N = 1 . As shown in Figure 7, when 0 ≤ i < iC or 180° − iC < i < 180°, a smaller velocity increment is required when thrust is applied directly to achieve RAAN separation. When iCi < 84.8035° or 95.1965° < i < 180° − iC, fuel savings are better achieved by adjusting the semi-major axis. When 84.8035° ≤ i ≤ 95.1965°, adjusting the inclination is more fuel efficient.
The values of iC in different cases are analyzed below.
As can be seen in Figure 8, the smaller the required RAAN separation, the larger the value of iC is taken; the lower the height of the parking orbit, the smaller the value of iC is taken.

4. Simulation Verification

Theoretical validation needs full-model simulation in real scenarios. The following parameters are set to validate the satellite constellation deployment strategy with two thrust modes (impulse thrust and continuous thrust) proposed in this study.

4.1. Multi-Satellite Multiplane Deployment

The space environment simulation parameter settings are shown in Table 1.
The target orbit and satellite thrust parameters are shown in Table 2.
The following is an example of controlling the semi-major axis (parking orbit height of 350 km) to achieve RAAN separation under continuous thrust condition, to illustrate the indirect control RAAN strategy in detail. The specific operation flow is shown in Table 3, and the process schematic is shown in Figure 9.
As shown in Figure 9, between T1 and T2 is the process of achieving RAAN separation for Group 2 satellites, and between T1 and T4 is the process of achieving RAAN separation for Group 3 satellites. The total velocity increment and time required to deploy the mission at different inclinations are analyzed below. Where the total velocity increment is the sum of the energy required for the three maneuvers and the time required is the total duration of the deployment mission, which is T6.

4.2. Low Inclination

Figure 10 shows the velocity increment (left Y-axis) and deployment time (right Y-axis) required for constellation deployment at inclination i = 3.2°. The energy and time required to control the semi-major axis to indirectly achieve RAAN separation under impulse thrust conditions are indicated in Figure 10a. The horizontal coordinate indicates the altitude of the parking orbit, and the target orbit height is 550 km. When the parking orbit height is 550 km, the velocity increment is 0, and the RAAN separation cannot be achieved, so the required time is infinite. The larger the height difference between the parking orbit and the target orbit, the larger the required velocity increment; but the required deployment time is shorter. Figure 10b represents the energy and time required to indirectly achieve RAAN separation by controlling the inclination under impulse thrust conditions, and the horizontal coordinate represents the orbital inclination of the parking orbit. Similarly, Figure 10c,d represent the energy and time required to control the semi-major axis and inclination under continuous thrust conditions, respectively. The time required for deployment using the controlled orbital semi-major axis under small inclination conditions is much less than the time required to control the inclination for both impulse and continuous thrust.
In this case, applying thrust to directly change the RAAN requires a velocity increment of 1317.85 m/s with impulse thrust and only 115.55 m/s with continuous thrust. Comparing the two scenarios of controlling the semi-major axis and directly applying thrust: (1) when the velocity increment is 115.55 m/s, the semi-major axis of the parking orbit is 343.94 km and the deployment time is 37.41 days; (2) With continuous thrust, the same velocity increment of 115.55 m/s is required for deployment, which takes only 26.74 days, saving 28.51% of the deployment time compared with the former.

4.3. Medium Inclination

Figure 11 shows the velocity increment (left Y-axis) and deployment time (right Y-axis) required for constellation deployment with impulse thrust and continuous thrust, respectively, when i = 53.7°. From the right Y-axis of the plots in Figure 11, for the medium inclination constellations, the efficiencies of the RAAN changing for controlling the semi-major axis and the inclination are approximated. In this case, the direct application of thrust to change the RAAN requires a velocity increment of 1668.23 m/s and a deployment time of 193.08 days. The strategy of directly applying thrust is not applicable to the deployment of medium inclination constellations, but it may be feasible to simultaneously control the semi-major axis and inclination.

4.4. High Inclination

Figure 12 shows the velocity increment (left Y-axis) and deployment time (right Y-axis) required for constellation deployment with impulse thrust and continuous thrust, respectively, at i = 87.9°. From the right Y-axis, for both impulse thrust and continuous thrust, the time required for deploying the constellation with controlled inclination at a small inclination is much smaller than that required for controlling the semi-major axis, indicating that in this case, the deployment of the constellation with controlled inclination is more efficient at high orbital inclination. Similarly, applying direct thrust to change the RAAN of the satellite requires a velocity increment of 2068.64 m/s and takes 239.42 days, which is not applicable to high inclination constellations.

5. Conclusions

The study can be applied to the selection of constellation deployment strategies at different inclinations and to the selection of constellation reconfiguration strategies at different inclinations. In the simulations, the proposed method is analyzed and validated and the following conclusions are obtained:
Numbered lists can be added as follows:
  • Due to the effect of integration, the velocity increment required for continuous thrust to change orbit elements is different from the velocity increment required for impulse thrust, and the impact is different for different orbit elements;
  • In the case of impulse thrust, there are orbital inclination boundaries that distinguish between high and medium-low inclinations. For the medium-low inclination constellations, it is more efficient to control the semi-major axis in orbit to achieve RAAN separation; for the high angle constellations, controlling the inclination is more efficient to achieve RAAN separation;
  • In the case of continuous thrust, there exist inclination boundaries distinguishing between high, medium, and low inclinations. For low inclination, it is more efficient to apply thrust directly to change the RAAN. For medium inclination, controlling the semi-major axis to achieve RAAN separation is more efficient. For high inclination, the efficiency of RAAN separation is higher by controlling the inclination;
  • In constellation deployment missions, controlling orbit elements is no longer the only viable option to achieve RAAN separation. The use of continuous thrust to directly change the RAAN at small inclinations saves both time and energy.
Future work may include the optimization of transfer orbits and the derivation of more sophisticated orbital dynamics. The next step to be investigated is the effect of other perturbation forces on the constellation deployment strategy and the use of additional perturbation forces to optimize the deployment strategy.

Author Contributions

Conceptualization, S.Z. and Y.Z.; Methodology, S.Z. and X.T.; Validation, S.Z. and Q.Z.; Writing—original draft, S.Z.; Writing—review & editing, X.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wang, K.; El-Mowafy, A. Proposed Orbital Products for Positioning Using Mega-Constellation LEO Satellites. Sensors 2020, 20, 5806. [Google Scholar] [CrossRef] [PubMed]
  2. Budianto, I.A.; Olds, J.R. Design and deployment of a satellite constellation using collaborative optimization. J. Spacecr. Rockets 2004, 41, 956–963. [Google Scholar] [CrossRef]
  3. Buzzi, P.G.; Selva, D.; Hitomi, N.; Blackwell, W.J. Assessment of constellation designs for earth observation: Application to the TROPICS mission. Acta Astronaut. 2019, 161, 166–182. [Google Scholar] [CrossRef]
  4. Bartolomé, J.P.; Maufroid, X.; Hernández, I.F.; López Salcedo, J.A.; Granados, G.S. Overview of Galileo System. In GALILEO Positioning Technology Signals and Communication Technology; Nurmi, J., Lohan, E.S., Sand, S., Hurskainen, H., Eds.; Springer: Dordrecht, The Netherlands, 2015; pp. 9–33. [Google Scholar] [CrossRef]
  5. Jia, L.; Zhang, Y. Cost Estimation Model for Mega-Constellation Deployment Missions. IEEE Access 2021, 9, 156778–156788. [Google Scholar] [CrossRef]
  6. Bainum, P.M.; Strong, A.; Tan, Z.; Capó-Lugo, P.A. Techniques for deploying elliptically orbiting constellations in along-track formation. Acta Astronaut. 2005, 57, 685–697. [Google Scholar] [CrossRef]
  7. Carlo, M.; Di Ricciardi, L.; Vasile, M. Multi-objective optimisation of constellation deployment using low-thrust propulsion. In AIAA/AAS Astrodynamics Specialist Conference; Aerospace Research Central: Reston, VA, USA, 2016; pp. 1–16. [Google Scholar] [CrossRef]
  8. Lee, H.W.; Jakob, P.C.; Ho, K.; Shimizu, S.; Yoshikawa, S. Optimization of satellite constellation deployment strategy considering uncertain areas of interest. Acta Astronaut. 2018, 153, 213–228. [Google Scholar] [CrossRef]
  9. Wijnen, M.; Correyero Plaza, S.; Pérez Grande, D.; de Saavedra, B.; Villegas Prado, D. ATHENA a modular, customizable electrospray thruster for pico-, nano- and micro-satellites. In Proceedings of the International Electric Propulsion Conference (IEPC), Cambridge, MA, USA, 19–23 June 2022. [Google Scholar]
  10. Grimaud, L.; Krejci, D.; Seifert, B. The IFM Micro FEEP thruster: A modular design for smallsat propulsion. In Proceedings of the International Electric Propulsion Conference (IEPC), Vienna, Austria, 15–20 September 2019; Available online: https://electricrocket.org/2019/675.pdf (accessed on 11 December 2023).
  11. Fong, C.J.; Shiau, W.T.; Lin, C.T.; Kuo, T.C.; Chu, N.H.; Yang, S.K.; Yen, N.L.; Chen, S.S.; Kuo, Y.H.; Liou, Y.A.; et al. Constellation deployment for the FORMOSAT-3/COSMIC mission. IEEE Trans. Geosci. Remote Sens. 2008, 46, 3367–3379. [Google Scholar] [CrossRef]
  12. King, J.A.; Beidleman, N.J. Method and Apparatus for Deploying a Satellite Network. U.S. Patent 5,199,672, 6 April 1993. [Google Scholar]
  13. Crisp, N.H.; Smith, K.; Hollingsworth, P. Launch and deployment of distributed small satellite systems. Acta Astronaut. 2015, 114, 65–78. [Google Scholar] [CrossRef]
  14. Gomez Jenkins, M.; Krejci, D.; Lozano, P. CubeSat constellation management using Ionic Liquid Electrospray Propulsion. Acta Astronaut. 2018, 151, 243–252. [Google Scholar] [CrossRef]
  15. Leppinen, H. Deploying a single-launch nanosatellite constellation to several orbital planes using drag maneuvers. Acta Astronaut. 2016, 121, 23–28. [Google Scholar] [CrossRef]
  16. Yang, H.; Jiang, Y.; Baoyin, H. Fuel efficient control strategy for constellation orbital deployment. Aircr. Eng. Aerosp. Technol. 2016, 88, 159–167. [Google Scholar] [CrossRef]
  17. McGrath, C.N.; Macdonald, M. General perturbation method for satellite constellation deployment using nodal precession. J. Guid. Control. Dyn. 2020, 43, 814–824. [Google Scholar] [CrossRef]
  18. Di Pasquale, G.; Sanjurjo-Rivo, M.; Pérez Grande, D. Optimization of constellation deployment using on-board propulsion and Earth nodal regression. Adv. Space Res. 2022, 70, 3281. [Google Scholar] [CrossRef]
  19. Zhang, Y. Satellite Constellation Theory and Design; Science Press: Beijing, China, 2008. [Google Scholar]
  20. Hang, G.R.; Li, S.N.; Kang, X.L.; Jin, Y.Z.; Sun, W.X. Current Status and Future Prospects of Hall Electric Propulsion Space Applications. Propuls. Technol. 2023, 44, 28–41. [Google Scholar]
  21. Tirila, V.-G.; Demairé, A.; Ryan, C.N. Review of Alternative Propellants in Hall Thrusters. Acta Astronaut. 2023, 212, 284–306. [Google Scholar] [CrossRef]
  22. Vallado, D.A. Fundamentals of Astrodynamics and Applications, 4th ed.; Microcosm Press: Portland, OR, USA, 2013. [Google Scholar]
  23. Bakhtiari, M.; Abbasali, E.; Danesh, K. Minimum Cost Perturbed Multi-impulsive Maneuver Methodology to Accomplish an Optimal Deployment Scheduling for a Satellite Constellation. J. Astronaut. Sci. 2023, 70, 18. [Google Scholar] [CrossRef]
  24. Ruggiero, A.; Pergola, P.; Marcuccio, S.; Andrenucci, M. Low-Thrust Maneuvers for the Efficient Correction of Orbital Elements. In Proceedings of the 32nd International Electric Propulsion Conference, Wiesbaden, Germany, 11–15 September 2011; IEPC-2011-102. Available online: https://api.semanticscholar.org/CorpusID:29040684 (accessed on 11 December 2023).
  25. Conway, B.A. Spacecraft Trajectory Optimization; Cambridge University Press: Cambridge, UK, 2010; pp. 112–138. ISBN 978-0-521-51850-5. [Google Scholar]
  26. Zhao, S.; Tao, X.; Li, Z. Initial Identification of Thrust and Orbit Elements for Continuous Thrust Spacecraft in Circular Orbit. Aerospace 2023, 10, 1012. [Google Scholar] [CrossRef]
Figure 1. RCN reference system and thrust Angle.
Figure 1. RCN reference system and thrust Angle.
Sensors 24 00631 g001
Figure 2. Velocity increment with impulse thrust.
Figure 2. Velocity increment with impulse thrust.
Sensors 24 00631 g002
Figure 3. Velocity increment with continuous thrust.
Figure 3. Velocity increment with continuous thrust.
Sensors 24 00631 g003
Figure 4. Strategy for RAAN separation.
Figure 4. Strategy for RAAN separation.
Sensors 24 00631 g004
Figure 5. Group delay strategy.
Figure 5. Group delay strategy.
Sensors 24 00631 g005
Figure 6. Constellation deployment strategy selection for different inclinations with impulse thrust.
Figure 6. Constellation deployment strategy selection for different inclinations with impulse thrust.
Sensors 24 00631 g006
Figure 7. Constellation deployment strategy selection for different inclinations with continuous thrust.
Figure 7. Constellation deployment strategy selection for different inclinations with continuous thrust.
Sensors 24 00631 g007
Figure 8. Values of iC for different scenarios (target orbital height 1200 km).
Figure 8. Values of iC for different scenarios (target orbital height 1200 km).
Sensors 24 00631 g008
Figure 9. Deployment process.
Figure 9. Deployment process.
Sensors 24 00631 g009
Figure 10. Velocity increments and deployment times required for different thrust modes (i = 3.2°). (a) Change semi-major axis under impulse thrust; (b) Change inclination under impulse thrust; (c) Change semi-major axis under continuous thrust; (d) Change inclination under continuous thrust.
Figure 10. Velocity increments and deployment times required for different thrust modes (i = 3.2°). (a) Change semi-major axis under impulse thrust; (b) Change inclination under impulse thrust; (c) Change semi-major axis under continuous thrust; (d) Change inclination under continuous thrust.
Sensors 24 00631 g010
Figure 11. Velocity increments and deployment times required for different thrust modes (i = 53.7°). (a) Change semi-major axis under impulse thrust; (b) Change inclination under impulse thrust; (c) Change semi-major axis under continuous thrust; (d) Change inclination under continuous thrust.
Figure 11. Velocity increments and deployment times required for different thrust modes (i = 53.7°). (a) Change semi-major axis under impulse thrust; (b) Change inclination under impulse thrust; (c) Change semi-major axis under continuous thrust; (d) Change inclination under continuous thrust.
Sensors 24 00631 g011
Figure 12. Velocity increments and deployment times required for different thrust modes (i = 87.9°). (a) Change semi-major axis under impulse thrust; (b) Change inclination under impulse thrust; (c) Change semi-major axis under continuous thrust; (d) Change inclination under continuous thrust.
Figure 12. Velocity increments and deployment times required for different thrust modes (i = 87.9°). (a) Change semi-major axis under impulse thrust; (b) Change inclination under impulse thrust; (c) Change semi-major axis under continuous thrust; (d) Change inclination under continuous thrust.
Sensors 24 00631 g012
Table 1. Perturbation parameter settings.
Table 1. Perturbation parameter settings.
PerturbationModels and Values
Three-body gravitySun, Moon and Major Planets: JPL DE405
TidesSolid Tides: IERS Conventions 2003
Non-spherical gravitational fieldGravity Modelling: EGM2008
Degree: 21; Order: 21
RelativityIERS Conventions 2003
Solar pressureShadow Models: Dual Cone
Optical Pressure Coefficient:1.00
Surface-to-Mass Ratio: 0.02 m2/kg
Atmospheric dragDensity Model: Jacchia-Roberts
Drag coefficient: 2.20
Surface-to-mass ratio: 0.02 m2/kg
Table 2. Target orbit elements and satellite thrust parameter settings.
Table 2. Target orbit elements and satellite thrust parameter settings.
Simulation ParametersParameter Value
Target orbit height (km)550
Number of target orbit planes3
RAAN separation (°)10
Mass of the satellite (kg)500
Thruster thrust (mN)50
Table 3. Deployment operation flow.
Table 3. Deployment operation flow.
StepOperationTime (Days)
Step 1After entering the parking orbit, all satellites were divided into three groups, labelled Group 1, Group 2, and Group 3;T1 = 0
Step 2Group 1 satellites began maneuvering (towards the target orbit);T1 = 0
Step 3Group 2 satellites began maneuvering (towards the target orbit);T2 = 12.41
Step 4Group 1 satellites reach the target orbit;T3 = 12.97
Step 5Group 3 satellites began maneuvering (towards the target orbit);T4 = 24.82
Step 6Group 2 satellites reach the target orbit;T5 = 25.38
Step 7Group 3 satellites reached the target orbit and deployment was completed.T6 = 37.79
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, S.; Zhu, Q.; Tao, X.; Zhang, Y. Deployment of Constellation with Different Inclinations Using the Nodal Precession and Thrust. Sensors 2024, 24, 631. https://doi.org/10.3390/s24020631

AMA Style

Zhao S, Zhu Q, Tao X, Zhang Y. Deployment of Constellation with Different Inclinations Using the Nodal Precession and Thrust. Sensors. 2024; 24(2):631. https://doi.org/10.3390/s24020631

Chicago/Turabian Style

Zhao, Shuailong, Qinyu Zhu, Xuefeng Tao, and Yasheng Zhang. 2024. "Deployment of Constellation with Different Inclinations Using the Nodal Precession and Thrust" Sensors 24, no. 2: 631. https://doi.org/10.3390/s24020631

APA Style

Zhao, S., Zhu, Q., Tao, X., & Zhang, Y. (2024). Deployment of Constellation with Different Inclinations Using the Nodal Precession and Thrust. Sensors, 24(2), 631. https://doi.org/10.3390/s24020631

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