Next Article in Journal
Model-Based Assessment of Elastic Material Parameters in Rheumatic Heart Disease Patients and Healthy Subjects
Next Article in Special Issue
Comparative Study of Metaheuristic Optimization of Convolutional Neural Networks Applied to Face Mask Classification
Previous Article in Journal / Special Issue
Observer-Based State Estimation for Recurrent Neural Networks: An Output-Predicting and LPV-Based Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelization of Low-Cost Maneuvers for an Areostationary Preliminary Mission Design

by
Marta M. Sánchez-García
1,
Gonzalo Barderas
1,2,* and
Pilar Romero
1,2
1
U.D. Astronomía y Geodesia, Facultad de Matemáticas, Universidad Complutense de Madrid, Plaza de las Ciencias, 3, E-28040 Madrid, Spain
2
Instituto de Matemática Interdisciplinar, Facultad de Matemáticas, Universidad Complutense de Madrid, Plaza de las Ciencias, 3, E-28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2023, 28(6), 105; https://doi.org/10.3390/mca28060105
Submission received: 20 September 2023 / Revised: 15 October 2023 / Accepted: 25 October 2023 / Published: 27 October 2023

Abstract

:
The aim of this paper is to analyze the determination of interplanetary trajectories from Earth to Mars to evaluate the cost of the required impulse magnitudes for an areostationary orbiter mission design. Such analysis is first conducted by solving the Lambert orbital boundary value problem and studying the launch and arrival conditions for various date combinations. Then, genetic algorithms are applied to investigate the minimum-energy transfer orbit. Afterwards, an iterative procedure is used to determine the heliocentric elliptic transfer orbit that matches at the entry point of Mars’s sphere of influence with an areocentric hyperbolic orbit imposing specific conditions on inclination and periapsis radius. Finally, the maneuvers needed to obtain an areostationary orbit are numerically computed for different objective condition values at the Mars entry point to evaluate an areostationary preliminary mission cost for further study and characterization. Results show that, for the dates of the minimum-energy Earth–Mars transfer trajectory, a low value for the maneuvers to achieve an areostationary orbit is obtained for an arrival hyperbola with the minimum possible inclination and a capture into an elliptical trajectory with a low periapsis radius and an apoapsis at the stationary orbit. For a 2026 mission with a TOF of 304 for the minimum-energy Earth–Mars transfer trajectory, for a capture with a periapsis of 300 km above the Mars surface the value achieved will be 2.083 km/s.

1. Introduction

Small relay satellites in areostationary orbit are considered the most efficient candidates to support the telecommunication needs in the 2020s [1,2,3,4,5,6,7]. Areostationary orbiters, like geostationary satellites for Earth [8,9], can provide continuous access at very high data rates to remotely supervise a significant population of probes and robotic missions on the Martian surface. The determination of transfer trajectories from Earth to Mars aimed at lowering costs in terms of impulses has become a key factor in mission planning, allowing for larger payloads to be transported at a minimum energy cost.
In this work, we analyze the design of an interplanetary Earth–Mars transfer to reach the areostationary orbit with the minimum impulsive maneuvers cost. Several authors have studied the optimization of interplanetary trajectories: in [10,11] transfer trajectories to the Moon and Jupiter, respectively, passing close to a Lagrangian point, are considered; in [12], a method is developed to obtain approximate near-optimal low-thrust interplanetary transfers using solar electric propulsion spacecraft; in [13], the optimization is performed with a cost function with variable coefficients; in [14], launch constraints are imposed for the optimization. We derive the heliocentric elliptic transfer characterizing the launch windows using an heuristic optimization method for determining an optimal time of flight (TOF) that minimizes the characteristic energies [15,16]. We will analyze the sensibility of this parameter in the optimization of impulsive maneuvers.
The first step consists in solving the Lambert problem [17] for various combinations of departure and arrival dates. Departure characteristic energy and hyperbolic arrival velocity plots are usually examined to investigate possible transfer windows [18]. We use genetic algorithms [19] to simultaneously minimize these two key parameters within these launch windows, comparing their performance.
Then, we match this interplanetary transfer with an entry hyperbola around Mars. The classic patched conic problem has been used to achieve a continuous trajectory composed of the trajectory between two planets and the planetocentric trajectory [20,21,22,23,24]. We use the iterative procedure [25] with imposed conditions on the periapsis distance, the arrival hyperbolic inclination, and a fixed radius for the Mars sphere of influence (SOI), and analyze the changes in the B-plane [26] due to the variations in the arrival asymptote direction. This iterative procedure enables the evaluation of these selected parameters in order to minimize fuel consumption for planning an areostationary mission. Once the fully matched trajectory to arrive at Mars is obtained, the maneuvers necessary to capture the orbiter and to place it in the areostationary orbit are analyzed.
The paper is organized as follows. In Section 2, we describe the dynamical model of the minimum-energy launch window problem. The determination of the Earth–Mars transfer trajectory with imposed hyperbolic arrival trajectory conditions is presented in Section 3. Section 4 describes the maneuvers performed to capture the spacecraft into an areostationary orbit and the numerical simulations to evaluate these maneuvers for different conditions. Finally, in Section 5, we briefly summarize the main conclusions.

2. Minimum-Energy Launch Window for Earth–Mars Transfer Trajectories

We first analyze the determination of interplanetary trajectories from Earth to Mars by minimizing the required energy at Earth departure and Mars arrival. We assume point mass gravitational forces for Earth and Mars within their respective spheres of influence and an unperturbed Keplerian orbit around the Sun.
The key parameters commonly used [27,28] to analyze the Earth–Mars mission launch opportunities are the characteristic energy at departure from Earth, C 3 = V E 2 , the hyperbolic excess velocity to escape from Earth, V E , and Mars arrival hyperbolic excess velocity, V M . In order to obtain these two parameters, it is necessary to first solve the Lambert orbital boundary value problem for the heliocentric spacecraft position, r s S ,
r ¨ s S = μ S r s S r s S 3 ,
constrained by two points, P 1 and P 2 , and an elapsed TOF, t 2 = t 1 + TOF , as illustrated in Figure 1,
r s S ( t 1 ) = r E ,
r s S ( t 2 ) = r M ,
where r E and r M are the heliocentric position vectors for the Earth at t 1 and Mars at t 2 , respectively.
The solution to the Lambert problem results in an elliptic conic section connecting P 1 and P 2 . We consider the shortcut solution that satisfies the boundary conditions, based on the iterative procedure, choosing the time transfer function introduced by Lancaster [29] as parameter for the iteration. The solution results in an elliptic conic section connecting P 1 and P 2 , with departure and arrival velocities, V s S 1 and V s S 2 , at t 1 and t 2 , respectively.
For each Earth departure and Mars arrival combination of dates, V E and V M change as V s S 1 and V s S 2 change according to
V E = V s S 1 V E S 1 ,
V M = V s S 2 V M S 2 ,
where V E S 1 and V M S 2 are the heliocentric velocity vectors for the Earth at t 1 and Mars at t 2 , respectively.
To analyze launch and arrival window opportunities, we first focus on data visualization [15,18] of the departure characteristic energy C 3 and the hyperbolic arrival velocity V M for various combinations of departure and arrival dates. We use Matlab software available from [30] to solve the Lambert problem, first obtaining a reduced launch window for the minimum-energy solution. The porkchop plots [31,32] shown in Figure 2 depict the contour lines of constant C 3 in km 2 / s 2 and V M in km / s for the 2019–2029 departure and the 2020–2030 arrival time frames. It is possible to observe that the launch and arrival windows that give the minimum values approximately repeat every Mars synodic period of about 780 days. In more detail, Figure 3a shows the departure characteristic energy and hyperbolic arrival velocity contour plots for the Earth–Mars transfer covering 27 months, from 1 July 2019 to 1 November 2021, and Figure 3b presents these plots for the departure window on July 2020 and the arrival window on February 2021.
Now, we search for the solution minimizing the equation:
C = C 3 + V M .
The minimum C in Equation (6) tends to give lower values of the impulsive maneuvers required, first to obtain an Earth escape velocity, and after, at the Mars arriving hyperbolic orbit, to reduce the hyperbolic excess velocity to capture the probe.
To minimize Equation (6), applied to the reduced windows previously estimated from porkchop plots, we now use a Matlab optimizer [33] that includes a library dedicated to genetic algorithms with different implementations of selection and crossover functions (see [19]). In order to analyze the accuracy of the genetic algorithms when applied to this problem, we compare the performance of the Remainder and the Stochastic Uniform functions as selection functions to select the individuals that contribute to the population at the next generation, and the Heuristic, the Scattered, and the Single point rules as crossover functions to combine two individuals to form the next generation for populations of 100, 500, and 2000 individuals.
Table 1 summarizes the key results to compare the genetic algorithms performances. It can be concluded that the considered selection and crossover functions do not significantly change the results. CPU time depends on the population size, leading to equivalent results.
Moreover, the results in Table 1 are in agreement with the trajectories defined by the different missions launched in the year 2020. The Mars 2020 mission (EEUU) [34] was launched on 30 July 2020, and its rover, the Perseverance rover, landed on 18 February 2021, with TOF = 203 days; The Tianwen-1 mission (China) [35] was launched on 23 July 2020 and arrived in Mars on 10 February 2021, with a TOF of 202 days. The Emirates Mars Mission (UAE) [36] was launched on 19 July 2020, and arrived at the orbit around Mars on 9 February 2021, with a TOF of 205 days.
In order to analyze how a lower TOF impacts V E and V M , we compare the heliocentric Earth–Mars optimal transfer orbit with the resulting orbit when a TOF 31 days less than the optimal minimum-energy trajectory is imposed. We choose a population of 500 individuals, the Stochastic Uniform function as the selection function, and the Heuristic as the crossover function. The resulting heliocentric elliptical orbits are illustrated in Figure 4, and their orbital elements are listed in Table 2. For the optimal minimum-energy orbit with a TOF of 197 days, departure on 20 July 2020 and arrival on 1 February 2021, we obtain values of V E = 3.6361 km/s and V M = 2.7682 km/s. When reducing the TOF to 166 days, departure and arrival dates change to 22 July 2020 and 4 January 2021, respectively, resulting in a more eccentric orbit, with V M = 3.5888 km/s significantly increased.
For the analysis henceforth, we focus on the optimal orbit for a transference in the launch window of 2026. We consider the window with the earliest launch date on 1 March 2026 and latest arrival date on 1 November 2027. We consider the same parameters for the genetic algorithms as considered in the case of the launch window of 2020: a population of 500 individuals, the Stochastic Uniform function as the selection function, and the Heuristic as crossover function. The obtained trajectory has its launch date on 31 October 2026 and arrival date on 31 August 2027. Table 3 presents the heliocentric elliptic orbital elements for this transference. The results, particularly the higher TOF value with respect to 2020, are in good agreement with the mission analysis referenced for year 2026 in [28].

3. Determination of Earth–Mars Trajectories with Hyperbolic Orbital Objective Values

Once the launch and arrival dates for the optimal minimum-energy solution have been determined, we move on to deal with the determination of the orbit for the entry point at the Mars SOI. To this end, we implement a procedure for matching the Earth–Mars elliptic transfer orbit and the Mars arrival hyperbolic orbit, fixing the periapsis distance, r p s M , the arrival hyperbolic inclination, i s M , and the radius for the SOI. We first present the determination of the arrival hyperbolic orbit and subsequent computation of the position and velocity at the matching point, which will be iterated with the Earth–Mars elliptic transfer afterwards. The goal is to analyze the changes in the B-plane [26] due to the variations in the arrival asymptote direction.
The iterative procedure used to compute the entry point at the SOI, in which both the heliocentric elliptic transfer orbit for the obtained TOF and the areocentric hyperbolic arrival orbit match, is formulated as follows:
r s M ( i + 1 ) = g f r s M ( i ) .
The function f provides the areocentric ecliptic transfer velocity, V M , at the Mars SOI entry point, r s M :
f ( r s M ( i ) ) = V M ( i ) ,
obtaining first V s S 2 , by solving Lambert’s problem (1) with a modified condition (3),
r s S ( t 2 ) = r M + r s M ,
and then using Equation (5).
The function g in (7) provides the position vector, r s M , in the areocentric ecliptic reference frame:
g ( V M ( i ) ) = r s M ( i + 1 ) .
This function gathers a set of expressions to obtain r s M at the SOI for the objective values of i s M and r p s M , based on [22,24,25], fixing the radius of the SOI instead of the time that the spacecraft is inside the SOI, as in [22].
To this end, we first determine the areoequatorial coordinates α , δ of the arrival asymptote, given by V s M (obtained by transforming V M to the areocentric areoequatorial reference frame), to obtain the parameter σ as
σ = arcsin tan δ tan i s M ,
defining the minimum inclination as the value of δ .
For a direct orbit, the right ascension of the ascending node, Ω s M , can be computed in two different ways:
Ω s M = α σ ( V s Mz > 0 ) ,
Ω s M = α + σ + π ( V s Mz < 0 ) .
The value of V s M determines the semimajor axis, a s M . Then, the eccentricity of the hyperbola, e s M , for a periapsis radius, r p s M , is fixed as
e s M = 1 + r p s M a s M .
The true anomaly, v s M , of the spacecraft at the SOI is determined using the standard procedure (i.e., [37,38]).
According to Figure 5, the unitary vectors defining the local reference system for the arrival asymptote, u T , u B , u H , are obtained as
u T = V s M | | V s M | | ,
u H = sin i s M sin Ω s M sin i s M cos Ω s M cos i s M ,
u B = u H × u T .
The plane containing u B , the planet center, and perpendicular to the arrival orbit is known as the B-plane, which is considered a fundamental tool in analyzing planetary arrivals [26]. With this iterative procedure, we update the B-plane according to the objective values imposed (see Figure 6).
The vectors (15)–(17) are then transformed to a local reference system in the periapsis, u R p , u T p , u H , according to
u R p = sin η u T + cos η u N ,
u T p = cos η u T sin η u N ,
where
η = arcsin 1 e s M .
Then, the areocentric areoequatorial position R s M = ( X s M , Y s M , Z s M ) vector components at the SOI are determined in terms of the base vectors u R p , u T p , u H , in a classical way, as
R s M = r s M cos v s M u R p + sin v s M u T p .
Finally, the position vector R s M is transformed to obtain r s M in the areocentric ecliptic reference frame, as a result of the function g in (10).
We consider that the iterative method defined by (7) converges, for a given tolerance τ , when the following condition is achieved:
r s M ( i + 1 ) r s M ( i ) < τ .
For comparison with Table 3, in Table 4, we present the results for the elliptic transfer orbit in the heliocentric ecliptic reference frame for the minimum-energy TOF ( 304.4618 days), the minimum arrival hyperbolic inclination ( i s M m i n = 16 . 1167 ), r p s M = 20,428 km, and a SOI of r I M = 577,239 km. We note that V M changes by 0.015 km/s approximately. Also, the resulting areocentric hyperbolic arrival trajectory orbital elements in the areoequatorial reference frame are shown in Table 4.

4. Mars Arrival Maneuvers Evaluation for an Areostationary Mission

We now conduct a preliminary evaluation of the total impulsive maneuvers Δ V M needed to capture the probe and to place it in an areostationary orbit:
Δ V M = Δ V c + Δ V P + Δ V A + Δ V i ,
where Δ V c is the capture maneuver to avoid the probe leaving the SOI on a flyby trajectory; Δ V P and Δ V A are the two Hohmann transfer maneuvers, at the perigee and apogee, respectively, designed to obtain a transfer orbit from the capture orbit to the target areostationary orbit; and Δ V i is the inclination correction maneuver to reach the final zero desired inclination.
If we consider a Δ V c maneuver at the periapsis of the hyperbolic orbit in order to obtain a circular capture orbit, its magnitude would be
Δ V c = V M 2 + 2 μ M r p s M μ M r p s M .
The hyperbolic periapsis distance, r p s M , is calculated as follows:
r p s M = a s M ( e s M 1 ) ,
with a s M as the hyperbolic semimajor axis and e s M as the hyperbolic eccentricity.
The two Hohmann transfer maneuvers required to achieve an orbit at the areostationary semimajor axis, r A = 20,428 km, in the case that r p s M < r A , are calculated as
Δ V P = μ M r p s M 2 r A r p s M + r A 1 ,
Δ V A = μ M r A 1 2 r p s M r p s M + r A .
Finally, the inclination maneuver is performed at the stationary distance in order to minimize its magnitude. For a maneuver at the node, this is obtained as
Δ V i = 2 μ M r A 2 μ M r A cos i s M ,
where i s M is the hyperbolic inclination.
In the case that r p s M > r A , the inclination maneuver at the node would be the first to be performed at the r p s M distance:
Δ V i = 2 μ M r p s M 2 μ M r p s M cos i s M .
Then, the two Hohmann maneuvers would be performed according to the following equations:
Δ V P = μ M r p s M 1 2 r A r p s M + r A ,
Δ V A = μ M r A 2 r p s M r p s M + r A 1 .
As can be seen from Equations (24) to (31), the periapsis distance, r p s M , and the orbital inclination of the arrival trajectory with respect to Mars, i s M , are the key design parameters in order to minimize the required impulses. Corrections to change the approach asymptotic plane are not considered, and neither are the combinations of the capture maneuver with the first Hohmann maneuver or the inclination maneuver with the second Hohmann maneuver.
Several numerical simulations are carried out to analyze the magnitude of these impulses due to the variations in the B-plane when different imposed values for TOF, i s M , and r p s M are considered.
We fix the launch window for 2026 according to the parameters given in Table 3 for the launch and arrival dates connecting Earth and Mars that minimize V M . In the first set of simulations, we consider different imposed objective values for the inclination, i s M , and the periapsis radius, r p s M , to determine the entry point at the SOI, using the iterative procedure defined in (7) to evaluate the total impulse Δ V M in Equation (23). We consider different values for i s M , ranging in the interval i s M m i n , 90 , and for r p s M , in the interval 15 , 000 km , 25 , 000 km , with the arrival velocity, V M , recalculated for the different entry points.
Figure 7a,b represent the results obtained for V M with different values of r p s M and i s M , respectively. For the objective conditions of i s M = 16 . 1167 , which corresponds to the i s M m i n for r p s M = r A = 20,428 km, a value of V M = 2.5763 km/s is obtained. As can be observed in Figure 7a for i s M m i n , V M decreases linearly when increasing r p s M from 2.5768 km/s to 2.5759 km/s. Figure 7b shows variations of about 2 × 10 3 km/s in V M , for the areostationary distance, as the objective inclination increases from i s M m i n to 90 . Table 5 summarizes the values for the C 3 , V M , and total cost C of Equation (6) for the minimum and maximum values of r p s M and i s M considered.
Figure 8a,b represent the values of Δ V c , Δ V i , Δ V A , Δ V P , and the total Δ V M for different values of r p s M and i s M to achieve the circular areostationary orbit at r A with zero areoequatorial inclination. In the case that r p s M < r A , the maneuvers are computed using (24) and (26)–(28). In the other case, the maneuvers are computed by means of (24) and (29)–(31).
As can be seen in Figure 8a, with the above-described strategy, the values for the capture maneuver Δ V c change, for i s M m i n , from 1.8246 km/s for r p s M = 15,000 km to 1.8631 km/s for r p s M = 25,000 km, depending on the different objective periapsis radius considered. Obviously, the Hohmann maneuver vanishes for a capture at r A , and its value varies from 0.2404 km/s for r p s M = 15,000 km to 0.1387 km/s for r p s M = 25,000. The inclination maneuver, Δ V i , to achieve zero inclination, from the i s M m i n for each r p s M , has a slight variation from 0.4062 km/s to 0.3673 km/s. The minimum for the total amount, Δ V M , has a value of 2.2493 km/s, which corresponds to r p s M = r A . This magnitude rises to 2.4712 km/s for r p s M = 15,000 km and to 2.3691 km/s for r p s M = 25,000 km.
Figure 8b shows that, for r p s M = r A , with Δ V H = 0 km/s, Δ V M increases when different values for the objective hyperbolic inclination are considered. From the same minimum of Δ V M = 2.2493 km/s for i s M m i n = 16 . 1167 , the total quantity of impulses rises to 3.8922 km/s, as illustrated in Figure 8a.
In a second set of simulations, we consider a different strategy to evaluate the Oberth effect that, due to the potential energy for a capture near the Mars surface, allows one to reduce the required Δ V M . To this end, we consider a capture into an elliptical trajectory with a low periapsis radius and an apoapsis at the stationary orbit. Then, a circularization process is performed at the apoapsis. For a capture at a periapsis of r p s M = r M + 300 km into an ellipse with an apoapsis of r A , the value of the impulse required to obtain a circular orbit at the areostationary distance is decreased to 1.6764 km/s, with Δ V c = 1.0294 km/s and Δ V A = 0.6470 km/s. The total Δ V M = 2.0834 km/s reduces 0.1659 km/s with respect to the previous strategy.
Finally, we analyze the influence of the TOF in Δ V M for a possible launch date delay. We consider up to two weeks delay from 31 October 2026 for the launch date and up to two months after 31 August 2027 for the arrival date for the minimum i s M . In the first case, a capture into a circular orbit of r p s M = r A is considered, as shown in Figure 9a. In the second case, a capture into an ellipse with r p s M = r M + 300 and an apoapsis of r A is considered, as shown in Figure 9b. It can be noticed, by comparing Figure 9a,b, that the second case is less sensitive to TOF variations. For case (a), a delay of 14 days at departure and 60 at arrival, increases the value of 2.25 km/s by 0.75 km/s, corresponding to the minimum-energy transfer TOF of 304 days, instead of an increase of 0.52 km/s from the value 2.08 km/s in case (b). In both cases, a launch delay of 14 days maintaining a TOF of 304 days would increase Δ V M by 0.25 km/s.

5. Conclusions

In this paper, we conducted a preliminary analysis of the impulsive maneuvers cost to transfer a spacecraft from Earth to Mars and to position it in an areostationary orbit. We first obtained the minimum-energy interplanetary transfer from Earth to Mars by applying genetic algorithms to select launch and arrival dates. Several simulations were carried out to analyze the performance of the genetic algorithms. Results show that differences in the final energy cost and the time of flight are negligible, and the only significant change is the CPU time needed to converge, which is dependent on the population size.
With the optimized launch and arrival dates, an iterative procedure was used to match the interplanetary trajectory obtained with the genetic algorithms, with an entry hyperbola defined by imposing objective conditions for the inclination and the periapsis radius of the orbit. Two different strategies were computed to evaluate the cost of the mission, Δ V M : The first includes a capture maneuver to a circular orbit at different periapsis radii, as well as Hohmann maneuvers and an inclination maneuver; The second includes a capture maneuver to an elliptic orbit with a low periapsis and an apoapsis at the areostationary orbit. Simulations with different imposed conditions on the entry hyperbola were conducted depending on the two key parameters: the hyperbolic inclination, i s M , and the periapsis radius, r p s M . For a circular capture at the stationary radius, results show that for a 2026 mission with a TOF of 304 days for the minimum-energy Earth–Mars transfer trajectory, the values achieved are Δ V c = 1.84 , Δ V H = 0 and Δ V i = 0.41 , being the total impulse Δ V M = 2.25 km/s for the minimum possible inclination i s M = 16 . 1167 and r p s M = 20,428 km corresponding to an areostationary radius. For an elliptical capture with a periapsis of r p s M = r M + 300 km and an apoapsis of 20,428, the values achieved are Δ V c = 1.03 , Δ V A = 0.65 , and Δ V i = 0.41 , the total impulse amounting to Δ V M = 2.08 km/s. A launch delay of two weeks would increase this minimum value of Δ V M by 0.25 km/s in both cases. The capture into an ellipse with a delay of 14 days at departure and 60 at arrival is less sensitive to TOF variations, increasing the total Δ V M by 0.52 km/s, instead of 0.75 for the direct circular capture.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CPUCentral Processing Unit
SOISphere Of Influence
TOFTime Of Flight

References

  1. Edwards, C.; Arnold, B.; DePaula, R.; Kazz, G.; Lee, C.; Noreen, G. Relay communications strategies for Mars exploration through 2020. Acta Astronaut. 2006, 59, 310–318. [Google Scholar] [CrossRef]
  2. Edwards, C.; DePaula, R. Key telecommunications technologies for increasing data return for future Mars exploration. Acta Astronaut. 2007, 61, 131–138. [Google Scholar] [CrossRef]
  3. Jentsch, C.; Rathke, A.; Wallner, O. Interplanetary communication: A review of future missions. In Proceedings of the 2009 International Workshop on Satellite and Space Communications, Siena, Italy, 9–11 September 2009; pp. 291–294. [Google Scholar] [CrossRef]
  4. Podnar, G.; Dolan, J.; Elfes, A. Telesupervised robotic systems and the human exploration of Mars. J. Cosmol. 2010, 12, 4058–4067. [Google Scholar]
  5. Romero, P.; Pablos, B.; Barderas, G. Analysis of orbit determination from Earth-based tracking for relay satellites in a perturbed areostationary orbit. Acta Astronaut. 2017, 136, 434–442. [Google Scholar] [CrossRef]
  6. Montabone, L.; Heavens, N.; Babuscia, A.; Barba, N.; Battalio, J.; Bertrand, T.; Edwards, C.; Guzewich, S.; Kahre, M.; Kass, D.; et al. Observing Mars from Areostationary orbit: Benefits and applications. In Proceedings of the Mars Exploration Program Analysis Group (MEPAG) #38, Virtual, 15–17 April 2020. [Google Scholar]
  7. Montabone, L.; Heavens, N.; Alvarellos, J.L.; Lillis, R.; Aye, M.; Liuzzi, G.; Babuscia, A.; Mischna, M.A.; Barba, N.; Newman, C.E.; et al. Observing Mars from Areostationary Orbit: Benefits and Applications. Available online: https://doi.org/10.13140/RG.2.2.21498.72643 (accessed on 20 September 2023).
  8. Romero, P.; Gambi, J. Optimal control in the east/west station-keeping manoeuvres for geostationary satellites. Aerosp. Sci. Technol. 2004, 8, 729–734. [Google Scholar] [CrossRef]
  9. Romero, P.; Gambi, J.; Patiño, E. Stationkeeping manoeuvres for geostationary satellites using feedback control techniques. Aerosp. Sci. Technol. 2007, 11, 229–237. [Google Scholar] [CrossRef]
  10. Prado, A.; Broucke, R. Transfer orbits in the Earth-Moon system using a regularized model. J. Guid. Control. Dyn. 1996, 19, 929–933. [Google Scholar] [CrossRef]
  11. Broucke, R.; Prado, A. Jupiter swing-by trajectories passing near the Earth. Adv. Astronaut. Sci. 1993, 82, 1159–1176. [Google Scholar]
  12. Kluever, C. Efficient Computation of Optimal Interplanetary Trajectories Using Solar Electric Propulsion. J. Guid. Control. Dyn. 2014, 38, 5. [Google Scholar] [CrossRef]
  13. Li, X.; Qiao, D.; Chen, H. Interplanetary transfer optimization using cost function with variable coefficients. Astrodynamics 2019, 3, 173–188. [Google Scholar] [CrossRef]
  14. Chen, L.; Li, J. Optimization of Earth-Mars transfer trajectories with launch constraints. Astrophys. Space Sci. 2022, 367, 12. [Google Scholar] [CrossRef]
  15. Woolley, R.; Whetsel, C. On the nature of Earth-Mars Porkchop plots. Adv. Astronaut. Sci. 2013, 148, 413–426. [Google Scholar]
  16. Sanchez-Garcia, M.M.; Barderas, G.; Romero, P. Analysis of the optimization for an Earth to Mars areostationary mission. In Proceedings of the European Planetary Science Congress, Virtual, 21 September–9 October 2020; Volume 14, p. EPSC2020-134. [Google Scholar] [CrossRef]
  17. Gooding, R.H. A procedure for the solution of Lambert’s orbital boundary-value problem. Celest. Mech. Dyn. Astron. 1990, 48, 145–165. [Google Scholar] [CrossRef]
  18. Conte, D. Survey of Earth-Mars trajectories using Lambert’s Problem and Applications. Ph.D. Thesis, The Pennsylvania State University, State College, PA, USA, 2014. [Google Scholar]
  19. Matlab. Genetic Algorithms Options. 2022. Available online: https://es.mathworks.com/help/gads/genetic-algorithm-options.html (accessed on 20 September 2023).
  20. Clarke, V.C., Jr.; Bollman, W.E.; Feitis, P.H.; Roth, R.Y. Design Parameters for Ballistic Interplanetary Trajectories, Part II: One-Way Transfers to Mercury and Jupiter; Technical Report 32–77, JPL; US Gov.: Washington, DC, USA, 1966.
  21. Cornelisse, J.W. Trajectory analysis for interplanetary missions. ESA J. 1978, 2, 131–144. [Google Scholar]
  22. Parvathi, S.P.; Ramanan, R.V. Direct Transfer Trajectory Design Options for Interplanetary Orbiter Missions using an Iterative Patched Conic Method. Adv. Space Res. 2016, 59, 1763–1774. [Google Scholar] [CrossRef]
  23. Parvathi, S.P.; Ramanan, R.V. Direct interplanetary trajectory design with a precise V-infinity targeting technique. In Proceedings of the 2017 First International Conference on Recent Advances in Aerospace Engineering (ICRAAE), Coimbatore, India, 3–4 March 2017; pp. 1–6. [Google Scholar] [CrossRef]
  24. Iwabuchi, M.; Satoh, S.; Yamada, K. Smooth and continuous interplanetary trajectory design of spacecraft using iterative patched-conic method. Acta Astronaut. 2021, 185, 58–69. [Google Scholar] [CrossRef]
  25. Bond, V.R. Matched-Conic Solutions to Round-Trip Interplanetary Trajectory Problems That Insure State-Vector Continuity at All Boundaries; Technical Note D-4942; NASA: Washington, DC, USA, 1969. [Google Scholar]
  26. Farnocchia, D.; Eggl, S.; Chodas, P.; Giorgini, J.; Chesley, S. Planetary encounter analysis on the B-plane: A comprehensive formulation. Celest. Mech. Dyn. Astron. 2019, 131, 36. [Google Scholar] [CrossRef]
  27. George, L.E.; Kos, L.D. Interplanetary Mission Design Handbook: Earth-to-Mars Mission Opportunities and Mars-to-Earth Return Opportunities, 2009–2024; Technical Report TM-1998-208533; NASA: Washington, DC, USA, 1998. [Google Scholar]
  28. Burke, L.M.; Falck, R.D.; McGuire, M.L. Interplanetary Mission Design Handbook: Earth-to-Mars Mission Opportunities 2026 to 2045; Technical Report TM-2010-216764; NASA: Washington, DC, USA, 2010. [Google Scholar]
  29. Lancaster, E.R.; Blanchard, R.C. A Unified Form of Lambert’s Theorem; Technical Note D-5368; NASA: Washington, DC, USA, 1969. [Google Scholar]
  30. Oldenhuis, R. Robust Solver for Lambert’s Orbital-Boundary Value Problem. 2017. Available online: http://es.mathworks.com/matlabcentral/fileexchange/26348-robust-solver-for-lambert-s-orbital-boundary-value-problem?? (accessed on 20 September 2023).
  31. Conte, D.; Di Carlo, M.; Ho, K.; Spencer, D.; Vasile, M. Earth-Mars transfers through Moon Distant Retrograde Orbits. Acta Astronaut. 2017, 143, 372–379. [Google Scholar] [CrossRef]
  32. Conte, D.; Spencer, D. Mission Analysis for Earth to Mars-Phobos Distant Retrograde Orbits. Acta Astronaut. 2018, 151, 761–771. [Google Scholar] [CrossRef]
  33. Trajectory Optimization Tool, v2.1.1. 2011. Available online: http://www.orbithangar.com/searchid.php?ID=5418 (accessed on 20 September 2023).
  34. Farley, K.A.; Williford, K.H.; Stack, K.M.; Bhartia, R.; Chen, A.; de la Torre, M.; Hand, K.P.; Goreva, Y.; Herd, C.D.K.; Hueso, R.; et al. Mars 2020 Mission Overview. Space Sci. Rev. 2020, 216, 142. [Google Scholar] [CrossRef]
  35. Jiang, X.; Yang, B.; Li, S. Overview of China’s 2020 Mars mission design and navigation. Astrodynamics 2018, 2, 1–11. [Google Scholar] [CrossRef]
  36. Sharaf, O.; Amiri, S.; AlDhafri, S.; Withnell, P.; Brain, D. Sending Hope to Mars. Nat. Astron. 2020, 4, 722. [Google Scholar] [CrossRef]
  37. Capderou, M. Satellites. Orbits and Missions; Springer: Paris, France, 2005. [Google Scholar]
  38. Montenbruck, O.; Gill, E. Satellite Orbits: Models, Methods, and Applications; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
Figure 1. Lambert orbital boundary value problem and the heliocentric elliptic trajectory.
Figure 1. Lambert orbital boundary value problem and the heliocentric elliptic trajectory.
Mca 28 00105 g001
Figure 2. Departure characteristic energy and hyperbolic arrival velocity contour plots for the Earth–Mars transfer spanning 10 years.
Figure 2. Departure characteristic energy and hyperbolic arrival velocity contour plots for the Earth–Mars transfer spanning 10 years.
Mca 28 00105 g002
Figure 3. Departure characteristic energy and hyperbolic arrival velocity contour plots for the Earth–Mars transfer (a) from July 2019 to November 2021 and (b) from July 2020 to February 2021.
Figure 3. Departure characteristic energy and hyperbolic arrival velocity contour plots for the Earth–Mars transfer (a) from July 2019 to November 2021 and (b) from July 2020 to February 2021.
Mca 28 00105 g003
Figure 4. Minimum-energy optimal trajectory for the Earth–Mars transfer for departure date on 1 July 2019 and arrival date on 1 November 2021 (top) in comparison with a trajectory with 31 fewer days of flight (bottom).
Figure 4. Minimum-energy optimal trajectory for the Earth–Mars transfer for departure date on 1 July 2019 and arrival date on 1 November 2021 (top) in comparison with a trajectory with 31 fewer days of flight (bottom).
Mca 28 00105 g004
Figure 5. Local reference system vectors u T , u B and u R p , u T p in the arrival hyperbolic orbit plane defined by the normal vector u H .
Figure 5. Local reference system vectors u T , u B and u R p , u T p in the arrival hyperbolic orbit plane defined by the normal vector u H .
Mca 28 00105 g005
Figure 6. Changes in the B-plane orientation due to the arrival asymptotic velocity variations.
Figure 6. Changes in the B-plane orientation due to the arrival asymptotic velocity variations.
Mca 28 00105 g006
Figure 7. V M values for (a) the i s M m i n for different periapsis radii, r p s M , ranging in the interval [15,000 km, 25,000 km] and (b) the objective values of r p s M = r A = 20,428 km for different inclination values, i s M , ranging in the interval i s M m i n = 16 . 1167 , 90 .
Figure 7. V M values for (a) the i s M m i n for different periapsis radii, r p s M , ranging in the interval [15,000 km, 25,000 km] and (b) the objective values of r p s M = r A = 20,428 km for different inclination values, i s M , ranging in the interval i s M m i n = 16 . 1167 , 90 .
Mca 28 00105 g007
Figure 8. Δ V c , Δ V i , Δ V A , Δ V P , and the total Δ V M values for (a) the i s M m i n for different periapsis radii, r p s M , ranging in the interval [15,000 km, 25,000 km] and (b) the objective values of r p s M = r A = 20,428 km for different inclinations, i s M , ranging in the interval i s M m i n = 16 . 1167 , 90 .
Figure 8. Δ V c , Δ V i , Δ V A , Δ V P , and the total Δ V M values for (a) the i s M m i n for different periapsis radii, r p s M , ranging in the interval [15,000 km, 25,000 km] and (b) the objective values of r p s M = r A = 20,428 km for different inclinations, i s M , ranging in the interval i s M m i n = 16 . 1167 , 90 .
Mca 28 00105 g008
Figure 9. Δ V values for a combination of dates ranging from 31 October 2026 to 13 November 2026 for the launch date and from 31 August 2027 to 30 October 2027 for the arrival date, for the following objective conditions: (a) minimum i s M and r p s M = r A = 20,428 km; (b) minimum i s M and a capture into an ellipse with a periapsis of r p s M = r M + 300 km and an apoapsis of r A .
Figure 9. Δ V values for a combination of dates ranging from 31 October 2026 to 13 November 2026 for the launch date and from 31 August 2027 to 30 October 2027 for the arrival date, for the following objective conditions: (a) minimum i s M and r p s M = r A = 20,428 km; (b) minimum i s M and a capture into an ellipse with a periapsis of r p s M = r M + 300 km and an apoapsis of r A .
Mca 28 00105 g009
Table 1. Simulation scenarios in a launch window in 2020 to compare the genetic algorithm performances minimizing the cost function, C = C 3 + V M , of the required C 3 energy and the arrival velocity, V M .
Table 1. Simulation scenarios in a launch window in 2020 to compare the genetic algorithm performances minimizing the cost function, C = C 3 + V M , of the required C 3 energy and the arrival velocity, V M .
Pop.CrossoverSelectionCPUDepartureTOF C 3 V M C
Time (s) Time (Days) ( km 2 / s 2 ) ( km / s )
Departure Date 20 July 2020
100Heuristicremainder10.8601:19:08196.939713.22162.768115.9897
stoch. unif.10.3601:05:27196.925313.22122.768515.9897
Scatteredremainder10.3001:04:01196.950113.22132.768415.9897
stoch. unif.10.4901:01:17196.928113.22122.768515.9897
Single ptremainder10.1401:06:02196.965713.22182.767915.9897
stoch. unif.10.3501:14:40196.956613.22172.768015.9897
500Heuristicremainder48.3001:06:27196.928813.22122.768515.9897
stoch. unif.47.8601:13:05196.942013.22152.768215.9897
Scatteredremainder48.3101:12:22196.942013.22152.768215.9897
stoch. unif.47.6001:21:17196.953713.22182.767915.9897
Single ptremainder47.3201:24:10196.940713.22172.768115.9898
stoch. unif.47.1901:16:41196.945013.22162.768115.9897
2000Heuristicremainder188.0601:04:09196.939313.22142.768315.9897
stoch. unif.186.5701:07:45196.933213.22132.768415.9897
Scatteredremainder222.9801:16:41196.945013.22162.768115.9897
stoch. unif.229.0301:10:38196.936113.22142.768315.9897
Single ptremainder238.3801:08:54196.934313.22142.768315.9897
stoch. unif.237.4901:26:54196.948013.22182.767915.9897
Table 2. Comparison of Earth–Mars optimal transfer heliocentric orbital parameters with respect to the mean ecliptic and equinox of J2000 for (a) the optimal launch and arrival dates in the window from 1 July 2019 to 1 November 2021 and (b) the launch and arrival dates in the window from 1 July 2019 to 1 November 2021 with the constraint of having 31 days of flight less than the optimal.
Table 2. Comparison of Earth–Mars optimal transfer heliocentric orbital parameters with respect to the mean ecliptic and equinox of J2000 for (a) the optimal launch and arrival dates in the window from 1 July 2019 to 1 November 2021 and (b) the launch and arrival dates in the window from 1 July 2019 to 1 November 2021 with the constraint of having 31 days of flight less than the optimal.
Parameter(a)(b)
Departure Date, t 1 20 July 2020 01:13:0522 July 2020 15:20:14
Arrival Date, t 2 1 February 20214 January 2021
Arrival time23:49:34 h15:21:40 h
Semimajor axis, a s S (km)198,312,598.97202,972,264.04
Eccentricity, e s S (unitless)0.233460.25130
Inclination, i s S (deg)1.736260.72898
Ascending node long., Ω s S (deg)297.4399299.8067
Arg. of the perihelium, ω s S (deg)359.6131358.1402
True anomaly, v s S (deg)0.46522.0420
V E (km/s)3.63613.7803
V M (km/s)2.76823.5888
Time of flight, TOF (days)196.9420166.0010
Total cost, C15.989717.8795
Table 3. Elliptic heliocentric orbital elements for the optimal launch and arrival dates in the window from 1 March 2026 to 1 November 2027.
Table 3. Elliptic heliocentric orbital elements for the optimal launch and arrival dates in the window from 1 March 2026 to 1 November 2027.
ParameterValue
Departure Date, t 1 31 October 2026 05:42:13 h
Arrival Date, t 2 31 August 2027 16:47:12 h
Semimajor axis, a s S (km)189,961,652.992134
Eccentricity, e s S (unitless)0.218496
Inclination, i s S (deg)0.8695
Periapsis argument, ω s S (deg)4.1768
Right ascension node longitude, Ω s S (deg)37.5815
True anomaly, v s S (deg)197.9132
V E (km/s)3.0311
V M (km/s)2.5913
Time of flight, TOF (days)304.4618
Table 4. Elliptic heliocentric orbital elements for (a) the Earth–Mars transfer orbit at the Mars sphere of influence entry point and final hyperbolic areocentric orbital elements and (b) launch date on 31 October 2026 and arrival date on 31 August 2027, for the objective values i s M = i s M m i n = 16 . 1167 and r p s M = 20,428 km.
Table 4. Elliptic heliocentric orbital elements for (a) the Earth–Mars transfer orbit at the Mars sphere of influence entry point and final hyperbolic areocentric orbital elements and (b) launch date on 31 October 2026 and arrival date on 31 August 2027, for the objective values i s M = i s M m i n = 16 . 1167 and r p s M = 20,428 km.
OrbitParameterUnitsValue
(a)Departure date, t 1 (UT)31 October 2026 05:42:13 h
Arrival date, t 2 (UT)31 August 2027 16:47:12 h
Semimajor axis, a s S (km)189,905,238.422086
Eccentricity, e s S (unitless)0.218286
Inclination, i s S (deg)0.9311
Periapsis argument, ω s S (deg)4.3094
Right ascension node longitude, Ω s S (deg)37.5723
True anomaly, v s S (deg)197.9302
V E (km/s)3.0333
V M (km/s)2.5763
Time of flight, TOF(days)304.4618
(b)Semimajor axis, a s M (km)6600.229103
Eccentricity, e s M (unitless)4.0950441
Inclination, i s M (deg)16.1167
Periapsis argument, ω s M (deg)194.1344
Right ascension node longitude, Ω s M (deg)162.7116
True anomaly at t 2 , v s M (deg)258.4533
Periapsis radius, r p s M (km)20,467.9232
Arrival at periapsis date, t p (UT)31 August 2027 20:03:55
Table 5. Interplanetary transfer energy variations for the considered combinations of i s M and the maximum and minimum r p s M and of r p s M = r A with the maximum and minimum i s M .
Table 5. Interplanetary transfer energy variations for the considered combinations of i s M and the maximum and minimum r p s M and of r p s M = r A with the maximum and minimum i s M .
i s M (deg) r ps M (km) C 3 ( km 2 / s 2 ) V M (km/s)C
i s M m i n = 16.1158 15,0009.20142.576811.7782
i s M m i n = 16.1175 25,0009.20082.575911.7767
i s M m i n = 16.1167 r A = 20 , 428 9.20112.576311.7774
90 r A = 20 , 428 9.21232.578111.7904
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

Sánchez-García, M.M.; Barderas, G.; Romero, P. Modelization of Low-Cost Maneuvers for an Areostationary Preliminary Mission Design. Math. Comput. Appl. 2023, 28, 105. https://doi.org/10.3390/mca28060105

AMA Style

Sánchez-García MM, Barderas G, Romero P. Modelization of Low-Cost Maneuvers for an Areostationary Preliminary Mission Design. Mathematical and Computational Applications. 2023; 28(6):105. https://doi.org/10.3390/mca28060105

Chicago/Turabian Style

Sánchez-García, Marta M., Gonzalo Barderas, and Pilar Romero. 2023. "Modelization of Low-Cost Maneuvers for an Areostationary Preliminary Mission Design" Mathematical and Computational Applications 28, no. 6: 105. https://doi.org/10.3390/mca28060105

APA Style

Sánchez-García, M. M., Barderas, G., & Romero, P. (2023). Modelization of Low-Cost Maneuvers for an Areostationary Preliminary Mission Design. Mathematical and Computational Applications, 28(6), 105. https://doi.org/10.3390/mca28060105

Article Metrics

Back to TopTop