Next Article in Journal
Failure Initiation Analysis of a PRSEUS BWB Wing Subjected to Structural Damage
Next Article in Special Issue
Feasibility Study of the Bare-Photovoltaic-Tether Concept: Prototypes and Experimental Performance Evaluation of the Photovoltaic Tether Segment
Previous Article in Journal
Effect of Working Current on C12A7 Hollow Cathode
Previous Article in Special Issue
E-Sail Optimal Trajectories to Heliostationary Points
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

E-Sail Option for Plunging a Spacecraft into the Sun’s Atmosphere

by
Giovanni Mengali
* and
Alessandro A. Quarta
Department of Civil and Industrial Engineering, University of Pisa, I-56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(4), 340; https://doi.org/10.3390/aerospace10040340
Submission received: 28 January 2023 / Revised: 3 March 2023 / Accepted: 31 March 2023 / Published: 1 April 2023
(This article belongs to the Special Issue Advances in CubeSat Sails and Tethers)

Abstract

:
A close observation of the near-Sun region, with in situ measurements, requires that a scientific probe be placed in a heliocentric orbit with a perihelion distance on the order of a few solar radii only. This is the solution adopted by the Parker Solar Probe (PSP), whose mission design uses a very complex transfer trajectory with seven Venus gravity assists to reach a perihelion radius of roughly 9.9 solar radii in about seven years. This paper aims to discuss the capability of an Electric Solar-Wind Sail (E-sail), i.e., a propellantless propulsion system that exploits the solar wind as a deep-space thrust source using a grid of long and artificially charged tethers, to drive a scientific probe toward a heliocentric orbit with characteristics similar to that considered during the initial design of the PSP mission. The two-dimensional trajectory analysis of an E-sail-based spacecraft is performed in an optimal framework, by considering the physical constraints induced by the thermal loads acting on the propellantless propulsion system when the spacecraft approaches the inner Sun regions. This means that, during the transfer trajectory, the E-sail-based spacecraft must avoid a spherical region around the Sun whose radius depends on the mechanical characteristics of the charged tethers. The paper shows that feasible solutions, in terms of optimal transfer trajectories, are possible even when a medium-performance E-sail is considered in the spacecraft design. In that context, the obtained trajectory can drive a scientific probe on the target (high elliptic) orbit in less than two years, without the use of any intermediate flyby maneuver.

1. Introduction

The study of the Sun and its interactions with the terrestrial planets is crucial for understanding the physical phenomena and the local properties of the near-Earth space environment. For example, a close observation of the near-Sun regions, with in situ measurements, is a fundamental step in the collection of necessary data and the elaboration of a reliable model of the solar magnetic field with an accurate solar-wind description. In this regard, the recent success of NASA’s historic Parker Solar Probe (PSP) mission, which is a milestone in the more ambitious NASA Living with a Star program [1,2], demonstrates the importance of the local observation of the near-Sun space. The PSP mission was conceived to investigate the Sun’s corona from inside and obtain accurate information on the upper solar atmosphere, the occurrence of solar flares and coronal mass ejections, the composition, structure and properties of the solar wind [3], and the influence of solar activity on space weather and the terrestrial environment [4,5].
To meet these scientific objectives, the probe must fly close to the Sun, i.e., its heliocentric orbit must have a very small perihelion distance, in the order of a few solar radii. This requirement makes such a type of scientific mission really a challenging problem due to the velocity change and the launch energy it requires. Not surprisingly, the PSP scenario called for several mission proposals and several changes and new designs [6]. The final trajectory, which is currently flown by the probe since 2018 (when it was launched), has an involved shape [7] and uses seven Venus flybys and approximately seven years to gradually decrease the probe perihelion distance to about 9.9 R or 0.046 au , where R is the reference value of the solar radius. The complexity of the PSP trajectory is the result of several design constraints [7], including the total mission costs (less than $750 M), the lack of any source of nuclear power, the requirement of three solar passes with perihelion distance smaller than 10 R , and a total mission flight time less than 10 years . In particular, the exclusion of radioisotope thermoelectric generators imposed a full mission reappraisal, discarding the idea of using a gravity assist from Jupiter, as planned in early mission design [7]. The renunciation of the Jupiter swingby called for multiple flybys with Venus to reduce the probe velocity enough and allow the spacecraft to approach the Sun at the desired (very close) distance.
This paper aims to analyze, with a simplified two-dimensional approach, a different solution to the transfer problem of a scientific probe to a high elliptic orbit with a very low perihelion radius, by looking for a new trajectory that, while satisfying the main constraints of the original concept behind the PSP mission scenario, can avoid the need for planetary flybys. The proposed approach uses a spacecraft propelled by an Electric Solar-Wind Sail (E-sail), i.e., an innovative propellantless propulsion concept able to produce a deep-space thrust by exploiting the electrostatic interaction between solar-wind ions and a charged tether grid [8,9]. Since an E-sail operates without any propellant stored onboard, it represents a very interesting option for propellantless missions requiring a continuous (and steerable) propulsive acceleration, and, indeed, it has been proposed for different advanced mission scenarios, such as, for example, the reaching of displaced non-Keplerian orbits [10], Solar System escape [11,12], or the maintenance of artificial equilibrium points in the Sun–Earth gravitational field [13]. Another potential application for an E-sail is a mission aimed at reaching high heliolatitudes with a perihelion near the Sun. From that position, it would be possible, through in situ measurements, to analyze electric current structures [14], whose sources are not yet sufficiently studied. Finally, an E-sail propulsion system may also represent an interesting means to insert a scientific probe into a high elliptic orbit with a perihelion radius inside the Sun’s upper atmosphere, as is undertaken by the PSP mission, but using a less involved interplanetary trajectory. Moreover, the E-sail can reach high heliolatitudes with perihelion near the Sun, in order to analyze (through in situ measurements) electric current structures [14], whose sources are not yet sufficiently studied. In the context of a two-dimensional scenario, the transfer trajectory analysis is performed in an optimal framework, by considering the physical constraints induced by the thermal loads acting on the E-sail when it approaches the Sun [15]. In practice, during the transfer, the spacecraft distance from the Sun cannot decrease below a minimum admissible value. Such a mission (path) constraint makes the interplanetary trajectory design a difficult problem to solve, because the E-sail minimum allowable distance from the Sun is one order of magnitude greater than the perihelion radius of the target, high elliptic, orbit.

2. Mission Description and Mathematical Preliminaries

Consider a two-dimensional heliocentric mission scenario, in which a spacecraft initially traces an Ecliptic circular orbit of radius r 1 au . This situation is consistent with a simplified case where the spacecraft uses a parabolic escape trajectory relative to the Earth, and the eccentricity of the planet’s orbit is neglected. The interplanetary mission requires the scientific probe to be transferred to a coplanar elliptic orbit, with a perihelion distance r p t close to the star to make in situ observations of the Sun’s upper atmosphere. To that end, the scientific probe must be designed to correctly operate under extreme conditions, such as those expected in the near-Sun environment, i.e., in presence of large temperature fluctuations and strong magnetic fields. In addition, the aphelion radius r a t of the target orbit should not be too far from the Sun [7] to allow the solar panels to supply the necessary power to onboard systems without the need for resorting to radioisotope thermoelectric generators. Recalling the main design parameters of NASA’s PSP and using the data reported by Guo [7], it is assumed that the target orbit of the scientific probe has the following characteristics
r p t = 10 R 0.0465 au , r a t = 1 au
that is, the semimajor axis and eccentricity of the heliocentric target orbit are, respectively,
a t = 0.523 au , e t = 0.911
These data correspond to a highly eccentric orbit with an orbital period of about 138 days and a very low perihelion radius. Please note that the nominal final orbit of the PSP has a perihelion radius equal to 9.9 R , an aphelion radius of about 0.73 au (the aphelion point is near the heliocentric orbit of Venus), and an orbital period of about 88 days .
The spacecraft propulsion system is an E-sail, which is used to transfer the scientific probe from the parking circular orbit to the target orbit with elements { a t , e t } defined in Equation (2). Since the E-sail charged tethers cannot operate near the Sun, due to the thermal loads coming from the star, the spacecraft can be ideally thought of as divided into two main macro-components, i.e., the scientific probe and the E-sail-based module, which is jettisoned once the scientific probe is inserted into the target, high elliptic, orbit. In practice, the conceptual design of the E-sail propulsion system [9] prevents the propellantless thruster from operating below a minimum solar distance r m , which depends on the charged tether characteristics. For example, in case of copper tethers, the minimum solar distance is r m = 0.33 au , while aluminum tethers impose the thruster to operate at a higher distance, i.e., r m = 0.5 au from the Sun [15]. In the rest of the paper, a copper tether case is assumed, so that the constraint on the Sun–spacecraft distance during the transfer is
r r m 0.33 au
The latter inequality corresponds to a path constraint to be enforced along the transfer trajectory. Note that, in a two-dimensional mission scenario, Equation (3) describes a forbidden circular region centered at the Sun, which partially contains the target highly elliptic orbit, as sketched in Figure 1. This particular mixture of orbit shape and path constraint makes the transfer trajectory design a difficult task to solve, as discussed below.
Figure 1 also shows the target orbit apse line and a fixed direction in the Ecliptic, which is useful for detecting the spacecraft’s azimuthal position during the transfer. Let ω f [ 0 , 360 ] deg be the argument of periapsis of the target orbit, calculated for the fixed direction, as shown in Figure 2. The apse line orientation is left free, i.e., ω f is unconstrained and is an output of the optimization process used to design the spacecraft transfer trajectory, as discussed later in the paper.

2.1. Spacecraft Dynamics

During the interplanetary transfer, the spacecraft dynamics may be conveniently described with the aid of a polar reference frame T ( O ; r , θ ) , the origin of which (point O) coincides with the Sun’s center-of-mass, and θ is the polar angle measured counterclockwise from the fixed direction. The latter is assumed to coincide with the Sun–spacecraft line at the initial time t = t 0 0 ; see Figure 3.
The E-sail-induced propulsive acceleration is described by the thrust model proposed by Huo et al. [16], according to which the radial ( a r ) and transverse ( a θ ) components of the propulsive acceleration vector are
a r = τ a c 2 r r 1 + cos α 2
a θ = τ a c 2 r r cos α sin α
where τ { 0 , 1 } is a switching dimensionless parameter, which models the E-sail electron gun operating mode [17], either on (when τ = 1 ) or off (i.e., τ = 0 ), α [ 90 , 90 ] deg is the pitch angle, defined as the angle between the radial line and the normal to the sail nominal plane in the direction opposite to the Sun, and a c is the characteristic acceleration, i.e., the maximum propulsive acceleration magnitude at a solar distance equal to one astronomical unit. Here, the nominal plane is defined as the plane that ideally contains the E-sail charged tethers. The variations of the two components { a r , a θ } and that of the acceleration magnitude a r 2 + a θ 2 with the sail pitch angle are reported (in dimensionless form) in Figure 4 when τ = 1 and r = r .
Using the expressions taken from Equations (4) and (5), the spacecraft equations of motion in the polar reference frame T are
r ˙ = u
θ ˙ = v r
u ˙ = μ r 2 + v 2 r + τ a c 2 r r 1 + cos α 2
v ˙ = u v r + τ a c 2 r r cos α sin α
where u (or v) is the radial (or transverse) component of the spacecraft inertial velocity, and μ is the Sun’s gravitational parameter. The differential system (6)–(9) is completed by four initial conditions. Since the spacecraft’s initial trajectory is a circular orbit of radius r , the initial conditions are
r ( t 0 ) = r , θ ( t 0 ) = 0 , u ( t 0 ) = 0 , v ( t 0 ) = μ r
where the initial polar angle is set equal to zero as it is measured from the fixed direction; see Figure 3.

2.2. Remarks on Trajectory Optimization

For a given characteristic acceleration a c , the time variation of the sail pitch angle has been obtained by minimizing the flight time Δ t = t f t 0 t f required to transfer the spacecraft from the circular parking orbit of radius r to the (coplanar) elliptic orbit of elements { a t , e t } given by Equation (2). The terminal conditions of Equation (2) can be translated into two equivalent constraints involving the final (i.e., calculated at the unknown time t f ) values of the spacecraft states { r , u , v } , viz.
u ( t f ) 2 + v ( t f ) 2 2 μ r ( t f ) = μ 2 a t
r ( t f ) v ( t f ) = μ a t 1 e t 2
Please note that the right-hand side term in Equation (11) (or in Equation (12)) coincides with the specific mechanical energy (or the specific angular momentum magnitude) of the target orbit. As a result, the problem amounts to looking for the transfer trajectory that minimizes the time necessary to move the dynamical system from the initial state (10) to the state defined by Equations (11) and (12), while enforcing the path constraint on the radial distance given by Equation (3). This problem has been solved with an indirect approach [18], paralleling the procedure described in Ref. [15]. The Hamiltonian function is
H λ r r ˙ + λ θ θ ˙ + λ u u ˙ + λ v v ˙
where { λ r , λ θ , λ u , λ v } are the variables adjoint to { r , θ , u , v } , respectively. The time derivative of a generic adjoint variable is obtained from the Euler-Lagrange equations
λ ˙ i = H i where i = { r , θ , u , v }
the explicit expression of which is here omitted for the sake of brevity. The additional boundary constraints are obtained from the transversality condition [18], i.e.,
H ( t f ) = 1
λ θ ( t f ) = 0
λ u ( t f ) μ r ( t f ) v ( t f ) 2 = r ( t f ) u ( t f ) λ r ( t f ) r ( t f ) λ v ( t f ) v ( t f )
Finally, the optimal control law, i.e., the variation of the sail pitch angle α and that of the switching parameter τ as a function of the adjoint variables is derived from the general results of Ref. [16].
The open-time unconstrained transfer trajectory, i.e., the orbital transfer problem without the path constraint of Equation (3), is a two-point boundary value problem with eight scalar differential equations (the four equations of motion and the four Euler-Lagrange equations) and nine boundary conditions, given by Equations (10)–(12) and (15)–(17). In fact, since the flight time is an output of the optimization process, one of the nine boundary conditions is used to calculate t f .
The approach used to solve the problem is as follows. For a chosen value of the spacecraft characteristic acceleration a c , the open-time unconstrained optimization problem is first considered, i.e., the minimum transfer trajectory with no path constraint is obtained using the previous indirect approach. The corresponding time-history of the Sun–spacecraft distance r = r ( t ) is therefore numerically obtained. If r ( t ) r m at any time, i.e., if both the transfer trajectory and the arrival point sketched in Figure 2 is outside the (circular) forbidden region, the results of the unconstrained transfer problem are consistent with the path constraint on the radial distance, and the obtained flight time coincides with the optimal solution. If, instead, the inequality r ( t ) r m is not met, i.e., there exist points such that r ( t ) < r m , the trajectory analysis is transformed into an optimal control problem with an interior point constraint on one state variable (the distance r in this case), according to the general method detailed in Ref. [18]. In the latter case, the whole transfer trajectory is split into two (or more) arcs, linked at the point where the trajectory is tangent to the circular forbidden zone. The presence of multiple arcs introduces additional (interior point) analytical constraints given, again, by the transversality condition [18]. The two-point boundary value problem has been solved through a hybrid numerical technique that combines genetic algorithms, used to obtain a first guess of the initial adjoint variables, with gradient-based and direct methods to refine the solution. A set of canonical units has been used in the integration of the differential equations to reduce their numerical sensitivity. The differential equations were integrated with double precision with a variable order Adams–Bashforth–Moulton solver scheme with absolute and relative errors of 10 12 .

3. Simulations Results

The capability of an E-sail-based spacecraft to reach the highly elliptic orbit defined by Equation (2) has been checked by simulating the optimal (constrained) transfer trajectories for fixed values of the characteristic acceleration within the set a c { 1.5 , 1.6 , 1.7 , 1.8 , 1.9 , 2 } mm / s 2 . The chosen range of variation of characteristic acceleration, consistent with a medium (i.e., a c 1.5 mm / s 2 ) or a high (i.e., a c 2 mm / s 2 ) performance E-sail, allows the reader to appreciate the topology variation of the optimal transfer trajectory and the corresponding change of total flight time. For example, a medium–low performance E-sail with a payload mass of about 100 kg has a total mass of about 400 kg , with a power budget of about 500 W .
Using the previously described optimization procedure, the results of the numerical simulations have been reassumed in Table 1, which reports the transfer time ( t f ), some characteristics of the optimal transfer trajectory (perihelion r p and aphelion r a distance), the spacecraft conditions at the arrival (radial distance r f r ( t f ) , velocity components u f u ( t f ) and v f v ( t f ) , final polar angle θ f θ ( t f ) , and true anomaly along the target orbit ν f ), the argument of periapsis of the target orbit ω f , and the number N of full revolutions around the Sun during the interplanetary transfer. According to Table 1, the variation of the perihelion distance with the characteristic acceleration indicates that the path constraint of Equation (3) becomes active when a c 1.8 mm / s 2 . Moreover, the variation of N with a c shows that a close passage of the spacecraft to the forbidden region takes place when a c 1.6 mm / s 2 . More precisely, if a c > 1.6 mm / s 2 , the optimal transfer trajectory is made of a single arc connecting the circular parking orbit and the elliptic target orbit. The endpoint of such a trajectory arc is just at the intersection point between the elliptic target orbit and the circular forbidden zone (i.e., the final distance is exactly r f = 0.33 au ) if a c { 1.8 , 1.7 } mm / s 2 . If, instead, a c 1.6 mm / s 2 the optimal transfer trajectory shows two (or, in general, more) arcs, and the trajectory itself is tangent (at one or more points) to the circular forbidden region.
This behavior is clearly illustrated in Figure 5, which shows the optimal transfer trajectory as a function of the characteristic acceleration within the selected set of values. In particular, the graphs reassumed in Figure 5 highlight the differences between the case when a c = 1.7 mm / s 2 and that when a c = 1.6 mm / s 2 in terms of trajectory topology; see Figure 5d,e, respectively. Note that, for the assigned characteristics of the target orbit, a c = 1.6 mm / s 2 represents a sort of threshold value of a c below which the optimal transfer trajectory requires at least a full revolution around the Sun and the path constraint of Equation (3) gives a tangent condition at one (or more) points.
For exemplary purposes, consider a spacecraft with a characteristic acceleration below such a threshold value, i.e., a c = 1.5 mm / s 2 . In this case, the time variation of the spacecraft states are those reported in Figure 6a, where the red dash line in the upper graph of the figure (corresponding to the radial distance r) represents the constraint on the minimum solar radius of Equation (3). According to Figure 5f, the Sun–spacecraft distance is just equal to r m at two different points, one coinciding with the final spacecraft position (i.e., the point occupied at t = t f ), and the other corresponding to the spacecraft position at about t = t f / 2 , when the transfer trajectory becomes tangent to the forbidden region, as sketched in Figure 5f. From Figure 6a, it is clear that the maximum transverse component of the spacecraft velocity v is reached just at that tangent point. This is by no means a surprising result, because the inner part of the trajectory (when r < r ) is used by the spacecraft to increase the osculating orbit aphelion radius r a , so that the succeeding part of the transfer trajectory has a behavior similar to that obtained for high values of a c ; see Figure 5a,b.
Figure 6b shows the time variation of the characteristics of the spacecraft osculating orbit, in term of semimajor axis a, eccentricity e, perihelion radius r p , and aphelion radius r a . In particular, the first two graphs of Figure 6b confirm that the elliptic target orbit is reached at the end of the transfer, while the other two plots indicate that the optimal transfer trajectory has a coasting phase (that is, a phase in which τ = 0 ) in the neighborhood of the tangent point with the circular forbidden zone. The presence of such a coasting phase is confirmed by Figure 7, which shows the time variation of the two control variables α and τ . Please note that the duration of the coasting phase is small (about 33 days ) when compared with the total flight time of 707 days .

4. Conclusions

The use of an advanced propellantless propulsion concept such as the E-sail allows complex interplanetary transfer trajectories and high-energy mission scenarios to become a viable option and a promising alternative to more conventional choices, usually based on chemical systems. In this context, this paper has analyzed the problem of reaching a highly elliptic heliocentric orbit with a perihelion distance of a few solar radii, to complete a campaign of in situ observation of the Sun’s upper atmosphere. This particular orbit, whose characteristics are consistent with the guidelines of the current PSP mission, can be reached with a flight time smaller than 2 years using a medium-performance E-sail. The preliminary design of the (two-dimensional) interplanetary transfer trajectory has been obtained as the solution to a minimum-time control problem that constraints the spacecraft to not enter inside a specified spherical region around the Sun to avoid excessive thermal loads on the E-sail charged tethers. This kind of (interior) path constraint makes the trajectory design a non-trivial problem to solve. In fact, the presence of a forbidden region that includes a non-negligible part of the target orbit reduces the admissible variation of the rendezvous point on the final orbit to a narrow range of (final) true anomalies. Nevertheless, feasible solutions are possible and the simulation results show that an E-sail with a characteristic acceleration of about 1.5 mm / s 2 may send the probe on the required target orbit without the need for several Venus flybys, as in the case of the PSP mission.

Author Contributions

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

Funding

This work is partly 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

The following abbreviations are used in this manuscript:
E-sailElectric Solar-Wind Sail
PSPParker Solar Probe

Notation

asemimajor axis [au]
a c       characteristic acceleration [mm/s 2 ]
a r radial component of propulsive acceleration [mm/s 2 ]
a θ transverse component of propulsive acceleration [mm/s 2 ]
eeccentricity
H Hamiltonian function
OSun’s center of mass
rSun–spacecraft radial distance [au]
R solar radius [au]
ttime [days]
T polar reference frame
uradial component of spacecraft velocity [km/s]
vtransverse component of spacecraft velocity [km/s]
α pitch angle [deg]
θ polar angle [deg]
λ r variable adjoint to r
λ u variable adjoint to u
λ v variable adjoint to v
λ θ variable adjoint to θ
μ Sun’s gravitational parameter [km 3 /s 2 ]
ν true anomaly [deg]
τ switching parameter
ω argument of periapsis [deg]
Subscripts
Earth
0initial, parking orbit
aaphelion
ffinal
mminimum
pperihelion
ttarget

References

  1. Brewer, D.A.; Barth, J.L.; Label, K.A.; Kauffman, W.J.; Giffin, G. NASA’s Living with a Star Program: Science with relevance. Acta Astronaut. 2002, 51, 609–616. [Google Scholar] [CrossRef]
  2. Withbroe, G.; Guhathakurta, M.; Hoeksema, J. Origins of the international living with a star program. Adv. Space Res. 2005, 35, 40–43. [Google Scholar] [CrossRef]
  3. Perrone, D.; Perri, S.; Bruno, R.; Stansby, D.; D’Amicis, R.; Jagarlamudi, V.; Laker, R.; Toledo-Redondo, S.; Stawarz, J.; Telloni, D.; et al. Evolution of coronal hole solar wind in the inner heliosphere: Combined observations by Solar Orbiter and Parker Solar Probe. Astron. Astrophys. 2022, 668, 1–13. [Google Scholar] [CrossRef]
  4. McComas, D.J.; Acton, L.W.; Balat-Pichelin, M.; Bothmer, V.; Dirling, R.B.; Feldman, W.C.; Gloeckler, G.; Habbal, S.R.; Hassler, D.M.; Mann, I.; et al. Solar Probe Plus: Report of the Science and Technology Definition Team; Technical Memorandum NASA/TM–2008-214161; NASA: Washington, DC, USA, 2008. [Google Scholar]
  5. Guo, Y.; Thompson, P.; Wirzburger, J.; Pinkine, N.; Bushman, S.; Goodson, T.; Haw, R.; Hudson, J.; Jones, D.; Kijewski, S.; et al. Execution of Parker Solar Probe’s unprecedented flight to the Sun and early results. Acta Astronaut. 2021, 179, 425–438. [Google Scholar] [CrossRef]
  6. Guo, Y. Solar Probe Plus: Mission design challenges and trades. Acta Astronaut. 2010, 67, 1063–1072. [Google Scholar] [CrossRef]
  7. Guo, Y. Trajectory Design of Solar Probe+ Using Multiple Venus Gravity Assists. In Proceedings of the AIAA/AAS Astrodynamics Specialist Conference and Exhibit, Honolulu, HI, USA, 18–21 August 2008. Paper AIAA 2008-7365. [Google Scholar] [CrossRef]
  8. Janhunen, P. Electric sail for spacecraft propulsion. J. Propuls. Power 2004, 20, 763–764. [Google Scholar] [CrossRef]
  9. 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]
  10. Mengali, G.; Quarta, A.A. Non-Keplerian orbits for electric sails. Celest. Mech. Dyn. Astron. 2009, 105, 179–195. [Google Scholar] [CrossRef]
  11. Matloff, G.L. The Solar-Electric Sail: Application to Interstellar Migration and Consequences for SETI. Universe 2022, 8, 252. [Google Scholar] [CrossRef]
  12. Mengali, G.; Quarta, A.A.; Janhunen, P. Considerations of electric sailcraft trajectory design. J. Br. Interplanet. Soc. 2008, 61, 326–329. [Google Scholar]
  13. Aliasi, G.; Mengali, G.; Quarta, A.A. Artificial Equilibrium Points for an Electric Sail with Constant Attitude. J. Spacecr. Rocket. 2013, 50, 1295–1298. [Google Scholar] [CrossRef]
  14. Khabarova, O.V.; Malova, H.V.; Kislov, R.A.; Zelenyi, L.M.; Obridko, V.N.; Kharshiladze, A.F.; Tokumaru, M.; Sokół, J.M.; Grzedzielski, S.; Fujiki, K. High-latitude Conic Current Sheets in the Solar Wind. Astrophys. J. 2017, 836, 108. [Google Scholar] [CrossRef]
  15. 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]
  16. 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]
  17. Janhunen, P.; Toivanen, P.K.; Polkko, J.; Merikallio, S.; Salminen, P.; Haeggström, E.; Seppänen, H.; Kurppa, R.; Ukkonen, J.; Kiprich, S.; et al. Electric solar wind sail: Toward test missions. Rev. Sci. Instrum. 2010, 81, 111301. [Google Scholar] [CrossRef] [PubMed]
  18. Bryson, A.E., Jr.; Ho, Y.C. Applied Optimal Control; Hemisphere Publishing Corporation: New York, NY, USA, 1975. [Google Scholar]
Figure 1. E-sail forbidden region and target orbit perihelion distance.
Figure 1. E-sail forbidden region and target orbit perihelion distance.
Aerospace 10 00340 g001
Figure 2. Target orbit Ecliptic position.
Figure 2. Target orbit Ecliptic position.
Aerospace 10 00340 g002
Figure 3. Reference frame and orbit characteristics.
Figure 3. Reference frame and orbit characteristics.
Aerospace 10 00340 g003
Figure 4. Variation of the propulsive acceleration (dimensionless) components and magnitude with the sail pitch angle when τ = 1 and r = r .
Figure 4. Variation of the propulsive acceleration (dimensionless) components and magnitude with the sail pitch angle when τ = 1 and r = r .
Aerospace 10 00340 g004
Figure 5. Optimal transfer trajectories as a function of the characteristic acceleration. Blue line → parking orbit; green line → target orbit; black line → transfer trajectory; black circle → starting point; black square → arrival point; filled pink circle → E-sail forbidden region.
Figure 5. Optimal transfer trajectories as a function of the characteristic acceleration. Blue line → parking orbit; green line → target orbit; black line → transfer trajectory; black circle → starting point; black square → arrival point; filled pink circle → E-sail forbidden region.
Aerospace 10 00340 g005
Figure 6. Time histories of spacecraft states and osculating orbit characteristics when a c = 1.5 mm / s 2 . Black circle → starting point; black square → arrival point.
Figure 6. Time histories of spacecraft states and osculating orbit characteristics when a c = 1.5 mm / s 2 . Black circle → starting point; black square → arrival point.
Aerospace 10 00340 g006
Figure 7. Time variation of the sail pitch angle α and switching parameter τ when a c = 1.5 mm / s 2 . Black circle → starting point; black square → arrival point.
Figure 7. Time variation of the sail pitch angle α and switching parameter τ when a c = 1.5 mm / s 2 . Black circle → starting point; black square → arrival point.
Aerospace 10 00340 g007
Table 1. Optimal transfer trajectory characteristics as a function of the E-sail propulsive performance.
Table 1. Optimal transfer trajectory characteristics as a function of the E-sail propulsive performance.
a c [mm/s 2 ] t f [Days] r p [au] r a [au] θ f [deg] r f [au] u f [km/s] v f [km/s] ν f [deg] ω f [deg]N
2.05100.41851.78101930.4185−45.7521.21210.17354.280
1.95460.37251.86222010.3725−49.9923.83213.31345.270
1.85870.33001.86222090.3300−54.3726.90216.68357.480
1.76390.33001.86222140.3300−54.3726.90216.68352.160
1.66660.33001.59855620.3300−54.3726.90216.68347.481
1.57070.33001.70175710.3300−54.3726.90216.68343.271
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

Mengali, G.; Quarta, A.A. E-Sail Option for Plunging a Spacecraft into the Sun’s Atmosphere. Aerospace 2023, 10, 340. https://doi.org/10.3390/aerospace10040340

AMA Style

Mengali G, Quarta AA. E-Sail Option for Plunging a Spacecraft into the Sun’s Atmosphere. Aerospace. 2023; 10(4):340. https://doi.org/10.3390/aerospace10040340

Chicago/Turabian Style

Mengali, Giovanni, and Alessandro A. Quarta. 2023. "E-Sail Option for Plunging a Spacecraft into the Sun’s Atmosphere" Aerospace 10, no. 4: 340. https://doi.org/10.3390/aerospace10040340

APA Style

Mengali, G., & Quarta, A. A. (2023). E-Sail Option for Plunging a Spacecraft into the Sun’s Atmosphere. Aerospace, 10(4), 340. https://doi.org/10.3390/aerospace10040340

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