Next Article in Journal
Modeling and Simulation of Flight Profile and Power Spectrum for Near-Space Solar-Powered UAV
Next Article in Special Issue
Trajectory Approximation of a Coulomb Drag-Based Deorbiting
Previous Article in Journal
Receding Horizon Trajectory Generation of Stratospheric Airship in Low-Altitude Return Phase
Previous Article in Special Issue
Rapid Evaluation of the Decay Time of a Plasma Brake-Based CubeSat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Circle-to-Ellipse Orbit Transfer for Sun-Facing E-Sail

Department of Civil and Industrial Engineering, University of Pisa, I-56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(11), 671; https://doi.org/10.3390/aerospace9110671
Submission received: 11 October 2022 / Revised: 26 October 2022 / Accepted: 27 October 2022 / Published: 29 October 2022
(This article belongs to the Special Issue Advances in CubeSat Sails and Tethers)

Abstract

:
The transfer between two coplanar Keplerian orbits of a spacecraft with a continuous-thrust propulsion system is a classical problem of astrodynamics, in which a numerical procedure is usually employed to find the transfer trajectory that optimizes (i.e., maximizes or minimizes) a given performance index such as, for example, the delivered payload mass, the propellant mass, the total flight time, or a suitable combination of them. In the last decade, this class of problem has been thoroughly analyzed in the context of heliocentric mission scenarios of a spacecraft equipped with an Electric Solar Wind Sail as primary propulsion system. The aim of this paper is to further extend the existing related literature by analyzing the optimal transfer of an Electric Solar Wind Sail-based spacecraft with a Sun-facing attitude, a particular configuration in which the sail nominal plane is perpendicular to the Sun-spacecraft (i.e., radial) direction, so that the propulsion system is able to produce its maximum propulsive acceleration magnitude. The problem consists in transferring the spacecraft, which initially traces a heliocentric circular orbit, into an elliptic coplanar orbit of given eccentricity with a minimum-time trajectory. Using a classical indirect approach for trajectory optimization, the paper shows that a simplified version of the optimal control problem can be obtained by enforcing the typical transfer constraints. The numerical simulations show that the proposed approach is able to quantify the transfer performance in a parametric and general form, with a simple and efficient algorithm.

1. Introduction

An Electric Solar Wind Sail (or E-sail) [1] is a propellantless device, conceived by Dr. Pekka Janhunen in 2004 [2], whose aim is to navigate the interplanetary space by extracting momentum from the stream of solar wind charged particles. When compared to more conventional propulsion systems, such as chemical or electric thrusters, an E-sail presents interesting peculiarities, including the capability of guaranteing long-lasting space missions or the possibility of conceiving new scenarios and innovative mission concepts.
A very promising class of potential applications for an E-sail consists in the exploration of the outer regions of the Solar System. In this regard, Janhunen et al. [3] analyzed an E-sail-based transfer towards Uranus, while a more ambitious objective is represented by an escape from the Solar System to reach the heliopause [4,5,6,7], the Sun’s gravitational focus [8], or even other stellar systems [9]. Notably, such long transfers can be made faster if the E-sail approaches the Sun in the first phase of its transfer trajectory. Such a strategy, referred to as solar wind assist [5], allows the E-sail to exploit the increased thrust due to a higher plasma density to reach very high orbital speeds and reduce the total mission time.
Other innovative mission scenarios that could be performed by an E-sail include the maintenance of displaced non-Keplerian orbits [10], the generation of artificial equilibrium points in the restricted three-body problem [11], and the asteroid deflection by means of a kinetic impactor [12,13]. Finally, an E-sail could also be exploited for deep space transfers towards planets [14,15], comets [16], asteroids [17], or other targets in the Solar System [18].
Since the E-sail propulsive acceleration depends on the orientation of its nominal plane [19] (that is, the mean plane containing the sail conducting tethers), the attitude control problem of an E-sail is closely related to that of its orbital control [20,21,22,23]. According to some preliminary studies [24,25,26,27], the particular orientation in which the E-sail nominal plane is perpendicular to the Sun-spacecraft line, referred to as Sun-facing configuration, may be passively maintained when the E-sail has an axisymmetric shape. Because of this relative simplicity from the standpoint of attitude control, much effort has been dedicated to the mission analysis of a Sun-facing E-sail configuration [28,29]. The latter is useful, for example, for generating spiral trajectories [30] or to perform phasing maneuvers [31]. A Sun-facing attitude has also been investigated for devices based on photonic propulsion [32,33,34,35], such as solar sails [36] or Smart Dusts [37], to simulate time-optimal transfer trajectories [38], to obtain closed-form solutions to the equations of motion [39], and to calculate the performance requirements for the maintenance of an artificial equilibrium point in the restricted three-body problem [40,41].
This paper focuses on two-dimensional transfers between heliocentric Keplerian orbits using a Sun-facing E-sail. It is assumed that the spacecraft departs from a circular parking orbit with the aim of reaching a coplanar target orbit of given eccentricity. Orbit-to-orbit transfers are analyzed in an optimal framework, where the total flight time is minimized with an indirect approach [42,43,44]. The novelty of the paper is to apply the optimal control theory to the special but important case when the direction of the E-sail propulsive acceleration vector is aligned with the Sun-spacecraft line. The resulting simplification of the control problem due to the reduced form of the thrust vector model, allows some useful analytical results to be obtained in a compact and elegant form. In particular, a set of dimensionless state variables is used in the dynamical model to make the results independent of the radius of the initial orbit. A parametric study is conducted, in which the transfer time is obtained as a function of the final eccentricity and the E-sail performance level. The optimal control law is then applied to simulate some potential heliocentric mission scenarios, including mission towards planets close to the Earth (Mars and Venus), or more ambitious targets, such as Jupiter.
The remainder of the manuscript is structured as follows. Section 2 introduces the dynamical model of a spacecraft subject to the Sun’s gravity and propelled by a Sun-facing E-sail. It also describes the procedure for trajectory optimization. Section 3 shows the results of the parametric analysis and analyzes some potential mission applications. The concluding section summarizes the main outcomes of the paper.

2. Problem Description

Consider a spacecraft that initially traces a heliocentric circular (parking) orbit of assigned radius r 0 . The spacecraft primary propulsion system is an E-sail that, after the deployment at time t 0 0 , is maintained in a Sun-facing condition, so that its propulsive acceleration vector a p is given by the following compact equation
a p = s β μ r 0 2 r 0 r i ^ r s β μ r 0 r i ^ r
where μ is the Sun’s gravitational parameter, r is the Sun-spacecraft distance, i ^ r is the radial (or Sun-spacecraft) unit vector, and s { 0 , 1 } is a dimensionless switching parameter that models the E-sail electron gun on ( s = 1 ) or off ( s = 0 ) condition. Here, s is introduced to account for possible coasting arcs (obtained with s = 0 ) in the spacecraft heliocentric trajectory. In Equation (1), β is a dimensionless positive parameter, defined as the ratio of the maximum propulsive acceleration magnitude a p to the Sun’s gravitational acceleration along the parking orbit, when r = r 0 . Note that β is a sort of reference propulsive acceleration, whose value depends on both the spacecraft mass and the E-sail design characteristics, such as the number and the length of the conducting tethers [45].
Using the propulsive acceleration of Equation (1) and assuming a simplified mission scenario where the external forces on the spacecraft are the E-sail thrust and the Sun’s gravitational pull only, it is easily concluded that the spacecraft propelled trajectory belongs to the parking orbit plane, and the semilatus rectum p of the osculating orbit is a constant of motion, that is
p r 0 for t t 0
Accordingly, the Sun-facing E-sail can potentially reach an elliptic Keplerian (target) orbit of eccentricity e f < 1 and semilatus rectum p f = r 0 , which is coplanar to the circular parking orbit. The aim of this paper is to analyze the optimal transfer trajectory that minimizes the flight time Δ t t f t 0 t f necessary for the E-sail to obtain a rendez-vous with a target coplanar orbit of given eccentricity e f . In particular, the orientation of the target orbit apse line (measured with respect to the Sun-spacecraft line at time t 0 ) is left free and is an output of the optimization process described at the end of this section.

2.1. Spacecraft Dynamics

The two-dimensional spacecraft dynamics can be studied by introducing a heliocentric polar reference frame T ( O ; i ^ r , i ^ θ ) , whose origin coincides with the Sun’s center of mass O, where i ^ θ is the transverse (or circumferential) unit vector in the direction of the spacecraft inertial velocity; see Figure 1.
Bearing in mind that the semilatus rectum of the osculating orbit is a constant of motion and using Equation (1) for the propulsive acceleration vector, the spacecraft equations of motion in the polar reference frame are
d r d t = u
d θ d t = h r 2 μ r 0 r 2
d u d t = μ r 2 + h 2 r 3 + s β μ r 0 2 r 0 r μ r 2 1 + r 0 r + s β r r 0
where u is the radial component of the spacecraft inertial velocity vector, θ is the polar angle measured counterclockwise from the initial Sun-spacecraft line (see Figure 1), and
h = μ p μ r 0
is the (constant) specific angular momentum of the spacecraft osculating orbit. The nonlinear differential Equations (3)–(5) are completed by the following initial conditions
r ( t 0 ) = r 0 , θ ( t 0 ) = 0 , u ( t 0 ) = 0
where, without loss of generality, it is assumed that the initial spacecraft polar angle is zero. The spacecraft dynamics can be more conveniently described in a dimensionless form by introducing the auxiliary (dimensionless) variables
τ t μ r 0 3 , ρ r r 0 , υ u r 0 μ
so that the equations of motion (3)–(5) can be rewritten as
d ρ d τ = υ
d θ d τ = 1 ρ 2
d υ d τ = 1 ρ 2 + 1 ρ 3 + s β ρ
while the initial conditions (7) become
ρ ( τ 0 ) = 1 , θ ( τ 0 ) = 0 , υ ( τ 0 ) = 0
where τ 0 0 is the initial dimensionless time.
The τ -variation of the E-sail switching parameter s is chosen to minimize the flight time Δ τ = τ f τ 0 τ f necessary for the spacecraft to reach an osculating elliptic orbit of given eccentricity e f < 1 . Note that this amounts to reaching a Keplerian orbit with a semimajor axis
a f = r 0 1 e f 2
or a specific mechanical energy
E f = μ 2 a f μ e f 2 1 2 r 0
The latter relation, in a dimensionless form, becomes
E f E f μ / r 0 = e f 2 1 2
The optimal function s = s ( τ ) , with τ [ 0 , τ f ] , and the corresponding value of the minimum flight time τ f are the solution of the minimum-time problem described in the next section.

2.2. Trajectory Optimization

The optimal transfer trajectory is obtained by maximizing the dimensionless performance index
J τ f
Using an indirect approach, introduce the Hamiltonian function
H λ ρ υ + λ θ ρ 2 + λ υ 1 ρ 2 + 1 ρ 3 + s β ρ
where { λ ρ , λ θ , λ υ } are the costates of { ρ , θ , υ } , respectively. The τ -variation of the costates is described by the Euler-Lagrange equations
d λ ρ d τ H ρ = 3 λ υ ρ 4 2 λ υ ρ 3 + s β λ υ ρ 2 + 2 λ θ ρ 3
d λ θ d τ H θ = 0
d λ υ d τ H υ = λ ρ
In particular, Equation (19) states that λ θ is a constant of motion, whose value is found by enforcing the transversality condition. In fact, assuming that the spacecraft final angular position θ f is left free, the transversality condition on λ θ is
λ θ = 0
and Equation (18) becomes
d λ ρ d τ = 3 λ υ ρ 4 2 λ υ ρ 3 + s β λ υ ρ 2
Instead, the initial value of the costate λ υ is derived by exploiting the transversality condition on the final value of the Hamiltonian function, that is
H ( τ f ) = 1
Indeed, observing that H is a constant of motion (as the Hamiltonian (17) does not explicitly depend on τ ) and using Equations (12) and (21), the condition (23) with s = 1 (the thruster is on to start the orbital transfer) gives
H ( τ 0 ) = λ υ ( τ 0 ) β = 1
and the initial value λ υ is
λ υ ( τ 0 ) = 1 β
Finally, the transversality condition gives an additional equation related to the terminal constraint (recall that the transverse component of the dimensionless spacecraft velocity is 1 / ρ ), that is
E ( τ f ) = υ f 2 2 + 1 2 ρ f 2 1 ρ f = E f
where υ f = υ ( τ f ) , ρ f = ρ ( τ f ) , and E f is the final value of the orbital energy given by Equation (15). Note that Equation (26) can be rewritten in standard form as
χ υ f 2 2 + 1 2 ρ f 2 1 ρ f E f = 0
so that the transversality condition gives
ρ f 3 υ f λ ρ ( τ f ) = ρ f 1 λ υ ( τ f )
Using the Pontryagin’s maximum principle, the control function s = s ( τ ) is obtained by maximizing, at any time τ , the portion of the Hamiltonian that explicitly depends on the control s, that is
H λ υ s β ρ
Since ρ and β are both positive parameters, the maximization of H gives a bang-bang control law
s = sign λ υ + 1 2
where sign is the signum function. In other terms, the costate λ υ plays the role of a switching function in determining the value of s. Note that the characteristics of the target orbit appear in Equation (27) through the orbital energy E f defined in Equation (15).
The initial value of the costate λ ρ and the optimal flight time τ f are solutions of a two-point boundary value problem (TPBVP) in which the two final constraints are given by Equations (27) and (28). The TPBVP has been solved, with an absolute error less than 10 8 , through a hybrid numerical technique that uses gradient-based and stochastic methods. Finally, the spacecraft nonlinear equations of motion have been numerically integrated with a variable order Adams-Bashforth-Moulton PECE solver, with absolute and relative errors equal to 10 10 . The numerical results are presented in the next section.

3. Numerical Simulations

The previous optimization procedure has been implemented to simulate the optimal circle-to-ellipse transfer of a Sun-facing E-sail for twenty different mission scenarios, corresponding to all possible combinations of e f { 0.1 , 0.2 , 0.3 , 0.4 , 0.5 } and β { 0.2 , 0.3 , 0.4 , 0.5 } .
The main simulation results are reported in Figure 2, which, in its left-hand side, shows the function τ f = τ f ( e f , β ) obtained with a two-dimensional interpolation procedure of the numerical results (black circles), while the right-hand side of the figure illustrates the level curves of the same function.
The main characteristics of the rendez-vous point, that is, { ρ f , θ f , υ f , ν f } , are summarized in Figure 3 as a function of the design parameters { e f , β } , where ν f [ 0 , 2 π ] rad is the spacecraft true anomaly on the target orbit, calculated with the conic polar equation in dimensionless form, that is
ν f = arccos 1 e f 1 ρ f 1 with ν f [ 0 , π ] rad if υ f 0 ( π , 2 π ] rad if υ f < 0
Finally, Figure 4 reports the semimajor axis a f and apse line rotation angle ω f (see Figure 1) of the elliptic target orbit.
The data of Figure 2, Figure 3 and Figure 4 are useful for obtaining a quick estimate of the overall transfer mission performance in a generic problem, because their values are independent of the radius r 0 of the circular parking.
The twenty analyzed mission scenarios, that is, the optimal transfer trajectories that are representative of the two design parameters ranging in the intervals e f [ 0.1 , 0.5 ] and β [ 0.2 , 0.5 ] , are characterized by the absence of coasting arcs in the transfer trajectory, which is completed in less than one full revolution around the Sun. This aspect is evident in Figure 5 and Figure 6, which show the results of the numerical simulation for the case of e f = 0.3 and β = 0.2 , while Figure 7 shows the optimal transfer trajectory for three sets of design parameters.
A more involved transfer trajectory, with a coasting arc, results when β is sufficiently small. For example, Figure 8 and Figure 9 show the numerical simulations for e f = 0.3 and β = 0.1 . In this case, the spacecraft completes a full revolution around the Sun before reaching the final orbit, and a coating arc appears in the time range τ [ 2 , 6 ] .
The proposed procedure can be used in the study of an interplanetary mission application, as discussed in the next section.

3.1. Mission Application

Consider an interplanetary mission scenario in which the Earth and the target planet orbit are assumed to be coplanar and circular. Starting from a heliocentric orbit of radius r 0 = 1 au , which models a spacecraft deployment along a parabolic escape orbit relative to the Earth, the Sun-facing E-sail is used to insert the vehicle into an elliptic target orbit, whose aphelion (or perihelion) lies on the circular orbit of the target (inner or outer) planet. In this case, the Sun’s gravitational acceleration along the parking orbit is μ / r 0 2 5.93 mm / s 2 , so that the maximum propulsive acceleration can be written as
max a p β × 5.93 mm / s 2 when r = r 0 1 au
while the semimajor axis of the elliptic target orbit, as a function of its eccentricity e f , is
a f = 1 au 1 e f 2
Finally, the perihelion and aphelion radius are, respectively
r p = 1 au 1 + e f
r a = 1 au 1 e f

3.1.1. Mars Case

Consider first a transfer towards Mars, where the planet orbit is approximated as a circular trajectory of radius r = r 1.524 au . In this case, assuming that the aphelion radius r a of the elliptic target orbit is equal to r , Equation (35) gives the value of the (target) eccentricity e f , that is
e f = 1 1 au r 0.3438
Using the eccentricity from the last equation and considering an E-sail with β { 0.2 , 0.3 , 0.4 , 0.5 } , the proposed optimization procedure gives the results summarized in Figure 10. The figure shows that an elliptic orbit with an aphelion radius equal to the heliocentric distance of Mars can be reached with a flight time on the order of one hundred days.
For example, assuming β = 0.3 , which corresponds to a maximum propulsive acceleration of about 1.8 mm / s 2 at 1 au , the optimal transfer trajectory and the transfer orbit characteristics are reported in Figure 11 and Figure 12, respectively.

3.1.2. Venus case

Consider now a circular heliocentric orbit of radius r = r = 0.723 au , that is, a scenario that models a transfer towards Venus. In that case, from Equation (34), the eccentricity of the target orbit is obtained by enforcing the constraint r p = r , and the result is
e f = 1 au r 1 0.3831
Using the usual range of variation of β { 0.2 , 0.3 , 0.4 , 0.5 } , the simulation results are summarized in Figure 13, while Figure 14 and Figure 15 give the transfer characteristics when β = 0.3 .

3.1.3. Jupiter Case

The final mission case refers to a more challenging transfer towards Jupiter, with a maximum propulsive acceleration of 1 mm / s 2 , which, according to Equation (32), corresponds to an E-sail with β 0.1686 . In this case, we assume a circular Jupiter orbit of radius r = r = 5.2 au , so that Equation (35) provides e f = 0.8077 as target eccentricity. This case is characterized by a large value of the target eccentricity and a relatively small value of the E-sail propulsive acceleration. In fact, the optimal transfer trajectory obtained with the proposed approach is more involved and contains a single coasting arc, as shown in Figure 16.
In this case, the transfer trajectory characteristics are summarized in Figure 17.

4. Conclusions

This paper has analyzed the transfer performance of a spacecraft propelled by an Electric Solar Wind Sail with a Sun-facing attitude in a two-dimensional circle-to-ellipse mission scenario. With the introduction of a suitable set of dimensionless spacecraft states, the proposed approach allows the transfer performance of the spacecraft to be evaluated as a function of two parameters only, that is, the final value of the orbital eccentricity and the sail propulsive performance. To that end, the paper presents a set of graphs that can be used to obtain a quick estimate of the minimum flight time and the characteristics of the final orbit without the need of any numerical simulation. The paper also shows a simplified mathematical model for the spacecraft trajectory optimization, which is able to numerically solve the two-point boundary value problem associated with the optimization procedure in a simple and effective way. The potential extensions of this work are related to the analysis of circle-to-ellipse transfers of the Sun-facing Electric Solar Wind Sail in the presence of an intermediate gravity assist maneuver. In that case the flyby could be used, for example, to change the value of the osculating orbit semilatus rectum, thus increasing the degrees of freedom in the design of the final orbit shape.

Author Contributions

Conceptualization, A.A.Q.; methodology, A.A.Q.; software, A.A.Q.; writing—original draft preparation, A.A.Q., M.B., L.N.; writing—review and editing, G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the University of Pisa, Progetti di Ricerca di Ateneo (Grant no. PRA_2022_1).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

aosculating orbit semimajor axis [au]
a p propulsive acceleration vector [mm/s2]
eosculating orbit eccentricity
Especific orbital energy [km2/s2]
E dimensionless specific orbital energy
horbital specific angular momentum magnitude [km2/s]
H dimensionless Hamiltonian function
i ^ r radial unit vector
i ^ θ transverse unit vector
Jdimensionless performance index
OSun’s center of mass
psemilatus rectum [au]
rSun-spacecraft radial distance [au]
sdimensionless switching parameter
ttime [days]
T polar reference frame
uradial component of the spacecraft velocity [km/s]
β dimensionless reference propulsive acceleration magnitude
θ polar angle [rad]
λ ρ dimensionless variable adjoint to ρ
λ υ dimensionless variable adjoint to υ
λ θ dimensionless variable adjoint to θ
μ Sun’s gravitational parameter [km3/s2]
ρ dimensionless radial distance
τ dimensionless time
υ dimensionless radial velocity
χ final dimensionless constraint
ω osculating orbit apse line rotation angle [rad]
Subscripts
0initial, parking orbit
aaphelion
ffinal, target orbit
Jupiter
Mars
pperihelion
Venus

References

  1. Bassetto, M.; Niccolai, L.; Quarta, A.A.; Mengali, G. A comprehensive review of Electric Solar Wind Sail concept and its applications. Prog. Aerosp. Sci. 2022, 128, 100768. [Google Scholar] [CrossRef]
  2. Janhunen, P. Electric sail for spacecraft propulsion. J. Propuls. Power 2004, 20, 763–764. [Google Scholar] [CrossRef]
  3. Janhunen, P.; Lebreton, J.P.; Merikallio, S.; Paton, M.; Mengali, G.; Quarta, A.A. Fast E-sail Uranus entry probe mission. Planet. Space Sci. 2014, 104, 141–146. [Google Scholar] [CrossRef] [Green Version]
  4. Mengali, G.; Quarta, A.A.; Janhunen, P. Considerations of electric sailcraft trajectory design. J. Br. Interplanet. Soc. 2008, 61, 326–329. [Google Scholar]
  5. Quarta, A.A.; Mengali, G. Electric sail mission analysis for outer solar system exploration. J. Guid. Control Dyn. 2010, 33, 740–755. [Google Scholar] [CrossRef]
  6. Sanchez-Torres, A. Propulsive force in electric solar sails for missions in the heliosphere. IEEE Trans. Plasma Sci. 2019, 47, 1657–1662. [Google Scholar] [CrossRef]
  7. Bassetto, M.; Quarta, A.A.; Mengali, G. Locally-optimal electric sail transfer. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 166–179. [Google Scholar] [CrossRef]
  8. Mengali, G.; Quarta, A.A. Trajectory analysis and optimization of Hesperides mission. Universe 2022, 8, 364. [Google Scholar] [CrossRef]
  9. Matloff, G.L. The Solar-Electric Sail: Application to Interstellar Migration and Consequences for SETI. Universe 2022, 8, 252. [Google Scholar] [CrossRef]
  10. Niccolai, L.; Anderlini, A.; Mengali, G.; Quarta, A.A. Electric sail displaced orbit control with solar wind uncertainties. Acta Astronaut. 2019, 162, 563–573. [Google Scholar] [CrossRef]
  11. Niccolai, L.; Caruso, A.; Quarta, A.A.; Mengali, G. Artificial collinear Lagrangian point maintenance with electric solar wind sail. IEEE TRansactions Aerosp. Electron. Syst. 2020, 56, 4467–4477. [Google Scholar] [CrossRef]
  12. Merikallio, S.; Janhunen, P. Moving an asteroid with electric solar wind sail. Astrophys. Space Sci. Trans. 2010, 6, 41–48. [Google Scholar] [CrossRef]
  13. Yamaguchi, K.; Yamakawa, H. Electric solar wind sail kinetic energy impactor for Near Earth Asteroid deflection mission. J. Astronaut. Sci. 2016, 63, 1–22. [Google Scholar] [CrossRef]
  14. Sanchez-Torres, A. Propulsive force in an electric solar sail for outer planet missions. IEEE Trans. Plasma Sci. 2015, 43, 3130–3135. [Google Scholar] [CrossRef]
  15. Huo, M.Y.; Mengali, G.; Quarta, A.A. Optimal planetary rendezvous with an electric sail. Aircr. Eng. Aerosp. Technol. 2016, 88, 515–522. [Google Scholar] [CrossRef]
  16. Quarta, A.A.; Mengali, G.; Janhunen, P. Electric sail option for cometary rendezvous. Acta Astronaut. 2016, 127, 684–692. [Google Scholar] [CrossRef]
  17. Slavinskis, A.; Janhunen, P.; Toivanen, P.; Muinonen, K.; Penttilä, A.; Granvik, M.; Kohout, T.; Gritsevich, M.; Pajusalu, M.; Sunter, I.; et al. Nanospacecraft fleet for multi-asteroid touring with electric solar wind sails. In Proceedings of the IEEE Aerospace Conference Proceedings, Big Sky, MT, USA, 3–10 March 2018; pp. 1–20. [Google Scholar] [CrossRef]
  18. Huo, M.Y.; Zhang, G.; Qi, N.; Liu, Y.; Shi, X. Initial trajectory design of electric solar wind sail based on finite Fourier series shape-based method. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 3674–3683. [Google Scholar] [CrossRef]
  19. Huo, M.Y.; Mengali, G.; Quarta, A.A. Electric sail thrust model from a geometrical perspective. J. Guid. Control Dyn. 2018, 41, 735–741. [Google Scholar] [CrossRef]
  20. Li, G.; Zhu, Z.H.; Du, C.; Meguid, S.A. Characteristics of coupled orbital-attitude dynamics of flexible electric solar wind sail. Acta Astronaut. 2019, 159, 593–608. [Google Scholar] [CrossRef]
  21. Zhao, C.; Huo, M.Y.; Qi, J.; Cao, S.; Zhu, D.; Sun, L.; Sun, H.; Qi, N. Coupled attitude-vibration analysis of an E-sail using absolute nodal coordinate formulation. Astrodynamics 2020, 4, 249–263. [Google Scholar] [CrossRef]
  22. Du, C.; Zhu, Z.H.; Li, G. Analysis of thrust-induced sail plane coning and attitude motion of electric sail. Acta Astronaut. 2021, 178, 129–142. [Google Scholar] [CrossRef]
  23. Du, C.; Zhu, Z.H.; Kang, J. Attitude control and stability analysis of electric sail. IEEE Trans. Aerosp. Electron. Syst. 2022. [Google Scholar] [CrossRef]
  24. Bassetto, M.; Mengali, G.; Quarta, A.A. Thrust and torque vector characteristics of axially-symmetric E-sail. Acta Astronaut. 2018, 146, 134–143. [Google Scholar] [CrossRef]
  25. Bassetto, M.; Mengali, G.; Quarta, A.A. Attitude dynamics of an electric sail model with a realistic shape. Acta Astronaut. 2019, 159, 250–257. [Google Scholar] [CrossRef]
  26. Bassetto, M.; Mengali, G.; Quarta, A.A. Stability and control of spinning E-sail in heliostationary orbit. J. Guid. Control Dyn. 2019, 42, 425–431. [Google Scholar] [CrossRef]
  27. Bassetto, M.; Mengali, G.; Quarta, A.A. E-sail attitude control with tether voltage modulation. Acta Astronaut. 2020, 166, 350–357. [Google Scholar] [CrossRef]
  28. Mengali, G.; Quarta, A.A.; Aliasi, G. A graphical approach to electric sail mission design with radial thrust. Acta Astronaut. 2013, 82, 197–208. [Google Scholar] [CrossRef]
  29. Quarta, A.A.; Mengali, G. Analysis of electric sail heliocentric motion under radial thrust. J. Guid. Control Dyn. 2016, 39, 1431–1435. [Google Scholar] [CrossRef]
  30. Bassetto, M.; Quarta, A.A.; Mengali, G.; Cipolla, V. Spiral trajectories induced by radial thrust with applications to generalized sails. Astrodynamics 2020, 5, 121–137. [Google Scholar] [CrossRef]
  31. Bassetto, M.; Boni, L.; Mengali, G.; Quarta, A.A. Electric sail phasing maneuvers with radial thrust. Acta Astronaut. 2021, 179, 99–104. [Google Scholar] [CrossRef]
  32. Fu, B.; Sperber, E.; Eke, F. Solar sail technology—A state of the art review. Prog. Aerosp. Sci. 2016, 86, 1–19. [Google Scholar] [CrossRef]
  33. Bovesecchi, G.; Corasaniti, S.; Costanza, G.; Tata, M.E. A novel self-deployable solar sail system activated by shape memory alloys. Aerospace 2019, 6, 78. [Google Scholar] [CrossRef] [Green Version]
  34. Zhang, F.; Gong, S.; Baoyin, H. Three-axes attitude control of solar sail based on shape variation of booms. Aerospace 2021, 8, 198. [Google Scholar] [CrossRef]
  35. Zou, J.; Li, D.; Wang, J.; Yu, Y. Experimental study of measuring the wrinkle of solar sails. Aerospace 2022, 9, 289. [Google Scholar] [CrossRef]
  36. McInnes, C.R. Orbits in a generalized two-body problem. J. Guid. Control Dyn. 2003, 26, 743–749. [Google Scholar] [CrossRef] [Green Version]
  37. Niccolai, L.; Bassetto, M.; Quarta, A.A.; Mengali, G. A review of Smart Dust architecture, dynamics, and mission applications. Prog. Aerosp. Sci. 2019, 106, 1–14. [Google Scholar] [CrossRef]
  38. Yamakawa, H. Optimal radially accelerated interplanetary trajectories. J. Spacecr. Rocket. 2006, 43, 116–120. [Google Scholar] [CrossRef]
  39. Bassetto, M.; Quarta, A.A.; Mengali, G.; Cipolla, V. Trajectory analysis of a Sun-facing solar sail with optical degradation. J. Guid. Control Dyn. 2020, 43, 1727–1732. [Google Scholar] [CrossRef]
  40. McInnes, C.R. Solar sail mission applications for non-Keplerian orbits. Acta Astronaut. 1999, 45, 567–575. [Google Scholar] [CrossRef]
  41. Aliasi, G.; Mengali, G.; Quarta, A.A. Artificial equilibrium points for a generalized sail in the elliptic restricted three-body problem. Celest. Mech. Dyn. Astron. 2012, 114, 181–200. [Google Scholar] [CrossRef]
  42. Morante, D.; Sanjurjo Rivo, M.; Soler, M. A Survey on Low-Thrust Trajectory Optimization Approaches. Aerospace 2021, 8, 88. [Google Scholar] [CrossRef]
  43. Vepa, R.; Shaheed, M.H. Optimal Trajectory Synthesis for Spacecraft Asteroid Rendezvous. Symmetry 2021, 13, 1403. [Google Scholar] [CrossRef]
  44. Shi, J.; Wang, J.; Su, L.; Ma, Z.; Chen, H. A Neural Network Warm-Started Indirect Trajectory Optimization Method. Aerospace 2022, 9, 435. [Google Scholar] [CrossRef]
  45. Janhunen, P.; Quarta, A.A.; Mengali, G. Electric solar wind sail mass budget model. Geosci. Instrum. Methods Data Syst. 2013, 2, 85–95. [Google Scholar] [CrossRef]
Figure 1. Reference frame and E-sail-based spacecraft state variables.
Figure 1. Reference frame and E-sail-based spacecraft state variables.
Aerospace 09 00671 g001
Figure 2. Minimum flight time τ f as a function of e f and β .
Figure 2. Minimum flight time τ f as a function of e f and β .
Aerospace 09 00671 g002
Figure 3. Characteristics of the rendez-vous point as a function of e f and β .
Figure 3. Characteristics of the rendez-vous point as a function of e f and β .
Aerospace 09 00671 g003
Figure 4. Geometric characteristics of target elliptic orbit as a function of e f and β .
Figure 4. Geometric characteristics of target elliptic orbit as a function of e f and β .
Aerospace 09 00671 g004
Figure 5. Optimal transfer trajectory when e f = 0.3 and β = 0.2 .
Figure 5. Optimal transfer trajectory when e f = 0.3 and β = 0.2 .
Aerospace 09 00671 g005
Figure 6. Optimal τ -variation of the spacecraft states when e f = 0.3 and β = 0.2 .
Figure 6. Optimal τ -variation of the spacecraft states when e f = 0.3 and β = 0.2 .
Aerospace 09 00671 g006
Figure 7. Optimal transfer trajectory for some set of design parameters.
Figure 7. Optimal transfer trajectory for some set of design parameters.
Aerospace 09 00671 g007
Figure 8. Optimal transfer trajectory when e f = 0.3 and β = 0.1 .
Figure 8. Optimal transfer trajectory when e f = 0.3 and β = 0.1 .
Aerospace 09 00671 g008
Figure 9. Optimal τ -variation of the spacecraft states when e f = 0.3 and β = 0.1 .
Figure 9. Optimal τ -variation of the spacecraft states when e f = 0.3 and β = 0.1 .
Aerospace 09 00671 g009
Figure 10. Optimal transfer performances in a Mars-based mission scenario ( e f = 0.3438 ) as a function of β .
Figure 10. Optimal transfer performances in a Mars-based mission scenario ( e f = 0.3438 ) as a function of β .
Aerospace 09 00671 g010
Figure 11. Optimal transfer trajectory in a Mars-based mission scenario with e f = 0.3438 and β = 0.3 .
Figure 11. Optimal transfer trajectory in a Mars-based mission scenario with e f = 0.3438 and β = 0.3 .
Aerospace 09 00671 g011
Figure 12. Transfer trajectory characteristics in a Mars-based mission scenario with e f = 0.3438 and β = 0.3 .
Figure 12. Transfer trajectory characteristics in a Mars-based mission scenario with e f = 0.3438 and β = 0.3 .
Aerospace 09 00671 g012
Figure 13. Optimal transfer performances in a Venus-based mission scenario ( e f = 0.3831 ) as a function of β .
Figure 13. Optimal transfer performances in a Venus-based mission scenario ( e f = 0.3831 ) as a function of β .
Aerospace 09 00671 g013
Figure 14. Optimal transfer trajectory in a Venus-based mission scenario with e f = 0.3831 and β = 0.3 .
Figure 14. Optimal transfer trajectory in a Venus-based mission scenario with e f = 0.3831 and β = 0.3 .
Aerospace 09 00671 g014
Figure 15. Transfer trajectory characteristics in a Venus-based mission scenario with e f = 0.3831 and β = 0.3 .
Figure 15. Transfer trajectory characteristics in a Venus-based mission scenario with e f = 0.3831 and β = 0.3 .
Aerospace 09 00671 g015
Figure 16. Optimal transfer trajectory in a Jupiter-based mission scenario with e f = 0.8077 and β = 0.1686 .
Figure 16. Optimal transfer trajectory in a Jupiter-based mission scenario with e f = 0.8077 and β = 0.1686 .
Aerospace 09 00671 g016
Figure 17. Transfer trajectory characteristics in a Jupiter-based mission scenario with e f = 0.8077 and β = 0.1686 .
Figure 17. Transfer trajectory characteristics in a Jupiter-based mission scenario with e f = 0.8077 and β = 0.1686 .
Aerospace 09 00671 g017
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Quarta, A.A.; Mengali, G.; Bassetto, M.; Niccolai, L. Optimal Circle-to-Ellipse Orbit Transfer for Sun-Facing E-Sail. Aerospace 2022, 9, 671. https://doi.org/10.3390/aerospace9110671

AMA Style

Quarta AA, Mengali G, Bassetto M, Niccolai L. Optimal Circle-to-Ellipse Orbit Transfer for Sun-Facing E-Sail. Aerospace. 2022; 9(11):671. https://doi.org/10.3390/aerospace9110671

Chicago/Turabian Style

Quarta, Alessandro A., Giovanni Mengali, Marco Bassetto, and Lorenzo Niccolai. 2022. "Optimal Circle-to-Ellipse Orbit Transfer for Sun-Facing E-Sail" Aerospace 9, no. 11: 671. https://doi.org/10.3390/aerospace9110671

APA Style

Quarta, A. A., Mengali, G., Bassetto, M., & Niccolai, L. (2022). Optimal Circle-to-Ellipse Orbit Transfer for Sun-Facing E-Sail. Aerospace, 9(11), 671. https://doi.org/10.3390/aerospace9110671

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