# Analysis of Electric Propulsion System for  Exploration of Saturn

Authors: Carlos Renato Huaura Solórzano; Antonio Fernando Bertachini de Almeida Prado; Alexander Alexandrovich Sukhanov
Journal: Mathematical Problems in Engineering
 (2009)

---

## Abstract

Exploration of the outer planets has experienced new interest with the launch of the Cassini and the New Horizons Missions. At the present time, new technologies are under study for the better use of electric propulsion system in deep space missions. In the present paper, the method of the transporting trajectory is used to study this problem. This approximated method for the flight optimization with power-limited low thrust is based on the linearization of the motion of a spacecraft near a keplerian orbit that is close to the transfer trajectory. With the goal of maximizing the mass to be delivered in Saturn, several transfers were studied using nuclear, radioisotopic and solar electric propulsion systems.


---

## Body

## 1. Introduction

The first mission to Saturn was Pioneer 11 that was launched on April 5, 1973. After making a flyby in Jupiter it determined the mass of the Jupiter's moon Callisto. Looping high above the ecliptic plane and across the Solar System, Pioneer 11 raced toward its appointment with Saturn on September 1, 1979. Voyager 2 was launched on August 20, 1977, and Voyager 1 was launched on September 5, 1977. Voyager 1 reached Saturn on November 12, 1980, followed by Voyager 2 in August 25, 1981. Later, the Cassini-Huygens spacecraft was launched on October 15, 1997. Using the gravity assist technique with the combination Venus-Venus-Earth and Jupiter the spacecraft increased its velocity to a level high enough to reach its final destination. On July 1, 2004, the Cassini-Huygens spacecraft fired its main engine to reduce its speed, allowing the spacecraft to be captured by Saturn's gravity and enter orbit. In May 28, 2008, the Cassini spacecraft passed by Saturn's moon Titan and made its last flyby of the original four-year tour, but Cassini's exploration of Saturn will continue for two more years.With the advent of both the Deep Space 1 mission (Brophy and Noca [1]) and the geosynchronous orbit insertion of a spacecraft, Electric Propulsion (EP) has moved from providing station-keeping capability for spacecrafts to be able to act as its primary propulsion system. More payload mass can be delivered than would have been possible with chemical propulsion due to the propellant reduction achieved. This reduction is due to the specific characteristics of EP systems: lowthrust, high efficiency, and long-burntimes.


## 2. General Information

The performance of a one-stage propulsion system can be roughly characterized by two variables: the maximum ejection velocity, which can be related with the maximum specific impulse and the ratio between the maximum thrust and the engine weight on the ground.Concerning these two parameters (Marec [2]), the propulsion system can be classified into high-thrust system (HT), characterized by a high-thrust acceleration level and a low specific impulse (conventional propulsion), and low-thrust system (LT) characterized by a low-thrust acceleration level and a high specific impulse (electric propulsion).Concerning the operating domain, the propulsion system can be essentially classified into constant ejection velocity (CEV) limited thrust systems, and limited power (LP) variable ejection velocity systems.For the case of the idealized LP system, the only constraint concerns to the power (W≤Wmax). More details are available in Marec [2].In this paper we studied the electric propulsion system of low-thrust type with limited power. Besides, the low-thrust transfers that are studied here have no constraint on the thrust direction. In general, these constraints can be caused by peculiarities of the attitude control system and the mode of the stabilization of the spacecraft.


## 3. Mathematical Model for the Optimization of the Low Thrust Transfers

A description of the mathematical models used to study the low-thrust transfers is now made, in order to explain the procedures used in the present paper.


### 3.1. Optimization of the Low Thrust

The electric propulsion low-thrust system uses the ionization of a propellant and its subsequent acceleration in an electrostatic field or electromagnetic to generate thrust. In systems with chemical propulsion, the ignition of the engine can run for several minutes, but in the case of the systems that use the electric propulsion, the ignition of the engine needs to run for longer times, up to several months in some cases.The equations of motion are(3.1)r⃗̇=V,⃗V⃗̇=f⃗v+α⃗.
In this way, the general equations of motion for trajectories obtained by the use of the low-thrust propulsion system are the equations of motion of the problem of two bodies with the inclusion of an additional term that represents the acceleration due to the propulsive force. The cost function (i.e., the functional to be minimized) is(3.2)J=x0(t1)=∫t0tff0(r⃗,f⃗v,α⃗)dt.
The Hamiltonian is(3.3)H=p0f0+p⃗rTv⃗+p⃗vTf⃗v+p⃗vTα⃗.
If f⃗v={f0,v⃗,f⃗v,α⃗} does not depend on x0=x0(t) then(3.4)H=-f0+p⃗rTv⃗+p⃗vTf⃗v+p⃗vTα⃗.
The effective power is(3.5)We=ηW=ṁpu22.
Usually We,u are given by(3.6)ṁp=-ṁ=2Weu2,FT=ṁpu=2Weu,α=FTm=2Wemu.
The limited power problem (LP) problem) has only an upper limit of the power given, that is 0≤We≤Wem. the exhaust velocity can be varied inside some limits, that is, umin<u<umax. The acceleration can be arbitrarily varied in the interval αmin≤α<αmax. Equations (3.5)–(8) show that(3.7)ṁp=m2α22We≥m2α22Wem.
In this way, the maximum power provides a minimum propellant consumption. For LP problems (umin<u<umax), due to (3.7) and the relation ṁ=-ṁp(3.8)α22We=ṁm2=ddt(1m)→12∫t0tfα2Wedt=1m1-1m0,
the performance index can be taken as(3.9)J=12∫t0tfα2Wedt→min.
Recalling (3.4), from (3.9) the Hamiltonian is(3.10)H=-α2We+p⃗rTv⃗+p⃗vTf⃗v+p⃗vTα⃗∂H∂α⃗=-α⃗T2We+p⃗vT=0⃗T,
where the optimal thrust is(3.11)α⃗=Wep⃗v.
For the case of constant power, it can be taken as(3.12)J=12∫t0tfα2dt→min,α⃗=p⃗v.
For the solar power, due to the variation with the inverse square of the distance from the Sun, we have(3.13)J=12∫t0tf(rα)2dt→min,α⃗=p⃗vr2.
For the analysis of the low-thrust trajectories, the transporting trajectory method is used.


### 3.2. The Method of the Transporting Trajectory

This approximated method of optimization of flight with ideally controlled small thrust is based on the linearization of the motion of a spacecraft near some reference Keplerian orbit (transporting trajectory). The equation of motion of the spacecraft is(3.14)x⃗̇=f⃗(x⃗)+g⃗,(3.15)f⃗(x⃗)={v⃗,μr3r⃗},g⃗={0⃗,α⃗}
with boundary values(3.16)x⃗(t0)=x⃗0,x⃗(t1)=x⃗1.
Now, considering the associated Keplerian motion described by equation(3.17)x⃗̇=f⃗(x⃗),
and let y⃗=y⃗(t) be a solution of (3.17), with given boundary values(3.18)y⃗(t0)=y⃗0,y⃗(t1)=y⃗1.
Note that the positions in the state vectors y⃗0,y⃗1 are given and the corresponding spacecraft velocities can be found by solving a Lambert’s problem. The solution of (3.14) with boundary values (3.16) is in the form(3.19)x⃗=ξ⃗+y⃗,
due to the low-thrust problem, it is reasonable to assume that(3.20)∥ξ⃗∥≪∥y⃗∥.
A linearization of (3.14) for ξ⃗ yields(3.21)ξ⃗̇=Fξ⃗+g⃗,(3.22)F=∂f⃗∂x⃗=[0IG0],G=μr3(3r⃗r⃗Tr2-I).
Boundary values for ξ⃗ are(3.23)ξ⃗(t0)=x⃗0-y⃗0=ξ⃗0,ξ⃗(t1)=x⃗1-y⃗1=ξ⃗1.
The Keplerian orbit given by y⃗=y⃗(t) is called a reference orbit or a transporting trajectory. Matrices (3.22) are calculated in this orbit.The Hamiltonian for the linearized problem is(3.24)H=-α22We+p⃗TF+p⃗vTα⃗,
where(3.25)p⃗={p⃗r,p⃗v}
is the adjoint (costate) variable.More details of the analysis of the trajectories using low thrust and the transporting trajectory method can be found in Beletsky and Egorov [3], Sukhanov [4, 5].


## 3.1. Optimization of the Low Thrust

The electric propulsion low-thrust system uses the ionization of a propellant and its subsequent acceleration in an electrostatic field or electromagnetic to generate thrust. In systems with chemical propulsion, the ignition of the engine can run for several minutes, but in the case of the systems that use the electric propulsion, the ignition of the engine needs to run for longer times, up to several months in some cases.The equations of motion are(3.1)r⃗̇=V,⃗V⃗̇=f⃗v+α⃗.
In this way, the general equations of motion for trajectories obtained by the use of the low-thrust propulsion system are the equations of motion of the problem of two bodies with the inclusion of an additional term that represents the acceleration due to the propulsive force. The cost function (i.e., the functional to be minimized) is(3.2)J=x0(t1)=∫t0tff0(r⃗,f⃗v,α⃗)dt.
The Hamiltonian is(3.3)H=p0f0+p⃗rTv⃗+p⃗vTf⃗v+p⃗vTα⃗.
If f⃗v={f0,v⃗,f⃗v,α⃗} does not depend on x0=x0(t) then(3.4)H=-f0+p⃗rTv⃗+p⃗vTf⃗v+p⃗vTα⃗.
The effective power is(3.5)We=ηW=ṁpu22.
Usually We,u are given by(3.6)ṁp=-ṁ=2Weu2,FT=ṁpu=2Weu,α=FTm=2Wemu.
The limited power problem (LP) problem) has only an upper limit of the power given, that is 0≤We≤Wem. the exhaust velocity can be varied inside some limits, that is, umin<u<umax. The acceleration can be arbitrarily varied in the interval αmin≤α<αmax. Equations (3.5)–(8) show that(3.7)ṁp=m2α22We≥m2α22Wem.
In this way, the maximum power provides a minimum propellant consumption. For LP problems (umin<u<umax), due to (3.7) and the relation ṁ=-ṁp(3.8)α22We=ṁm2=ddt(1m)→12∫t0tfα2Wedt=1m1-1m0,
the performance index can be taken as(3.9)J=12∫t0tfα2Wedt→min.
Recalling (3.4), from (3.9) the Hamiltonian is(3.10)H=-α2We+p⃗rTv⃗+p⃗vTf⃗v+p⃗vTα⃗∂H∂α⃗=-α⃗T2We+p⃗vT=0⃗T,
where the optimal thrust is(3.11)α⃗=Wep⃗v.
For the case of constant power, it can be taken as(3.12)J=12∫t0tfα2dt→min,α⃗=p⃗v.
For the solar power, due to the variation with the inverse square of the distance from the Sun, we have(3.13)J=12∫t0tf(rα)2dt→min,α⃗=p⃗vr2.
For the analysis of the low-thrust trajectories, the transporting trajectory method is used.


## 3.2. The Method of the Transporting Trajectory

This approximated method of optimization of flight with ideally controlled small thrust is based on the linearization of the motion of a spacecraft near some reference Keplerian orbit (transporting trajectory). The equation of motion of the spacecraft is(3.14)x⃗̇=f⃗(x⃗)+g⃗,(3.15)f⃗(x⃗)={v⃗,μr3r⃗},g⃗={0⃗,α⃗}
with boundary values(3.16)x⃗(t0)=x⃗0,x⃗(t1)=x⃗1.
Now, considering the associated Keplerian motion described by equation(3.17)x⃗̇=f⃗(x⃗),
and let y⃗=y⃗(t) be a solution of (3.17), with given boundary values(3.18)y⃗(t0)=y⃗0,y⃗(t1)=y⃗1.
Note that the positions in the state vectors y⃗0,y⃗1 are given and the corresponding spacecraft velocities can be found by solving a Lambert’s problem. The solution of (3.14) with boundary values (3.16) is in the form(3.19)x⃗=ξ⃗+y⃗,
due to the low-thrust problem, it is reasonable to assume that(3.20)∥ξ⃗∥≪∥y⃗∥.
A linearization of (3.14) for ξ⃗ yields(3.21)ξ⃗̇=Fξ⃗+g⃗,(3.22)F=∂f⃗∂x⃗=[0IG0],G=μr3(3r⃗r⃗Tr2-I).
Boundary values for ξ⃗ are(3.23)ξ⃗(t0)=x⃗0-y⃗0=ξ⃗0,ξ⃗(t1)=x⃗1-y⃗1=ξ⃗1.
The Keplerian orbit given by y⃗=y⃗(t) is called a reference orbit or a transporting trajectory. Matrices (3.22) are calculated in this orbit.The Hamiltonian for the linearized problem is(3.24)H=-α22We+p⃗TF+p⃗vTα⃗,
where(3.25)p⃗={p⃗r,p⃗v}
is the adjoint (costate) variable.More details of the analysis of the trajectories using low thrust and the transporting trajectory method can be found in Beletsky and Egorov [3], Sukhanov [4, 5].


## 4. Results of the Mission Analysis

A power-limited low-thrust transfer to Saturn is considered in this paper. For all cases considered below it is assumed that the spacecraft leaves the Earth's sphere of influence with variable velocity and approaches Saturn with zero velocity. Several analyses were considered as a function of the ratio between the final and initial mass, time of flight (TOF), velocity at infinity, and effective specific power. The unit of power used in this paper is the effective specific power that is [efficiency* (total power/total mass)]. Besides, for all cases, the efficiency considered is 100%. In this paper, the solar arrays degradation and decay of the nuclear/radioisotope source are not considered. But the nuclear/radioisotope power is considered constant. On the other way, the variation of the solar power is considered, taking into account the distance from the Sun. The solar power of the solar arrays follows the rule of the inverse of the square of the spacecraft heliocentric distance. The planetary state vector obtained from the planetary ephemerides (J2000) is used, together with Chebyshev’s interpolation. The tolerance used is 10-8. After analyzing several dates of launch between 2020 and 2030, we choose the option in May 06, 2021.


### 4.1. Nuclear Electric Propulsion

Nuclear Electric Propulsion (NEP) uses a reactor power system to provide the electricity for thrusters that ionize and accelerate propellant to produce thrust. The application of nuclear systems in space has advantage in situations where the distances from the Sun are large and the solar power density is too low≳1.5AU) and locations where solar power is not readily or continuously available. One of the key performances for NEP is the power system specific mass, measured in kg/kW. Low values are desired in order to provide the maximum mass allocation for payload or propellant. High values result in minimal delivered payloads.In this section, the study of the NEP system for a trip to Saturn is performed. Figure1 shows the trajectory of the flight Earth-Saturn projected on the ecliptic plane. Figure 2 shows that the use of NEP allows the delivery of larger masses to the target planet, due to the high power of this system.Figure 1
Trajectory of the flight Earth-Saturn with effective specific power (NEP) of 50 W/kg andVinf=10 km/s on Earth’s sphere influence (TOF=12years).Curves for NEP as a function of the ratiomf/m0 for several TOFs. The effective specific powers are (a) 100 W/kg, (b) 80 W/kg, (c) 60 W/kg, (d) 30 W/kg.
(a)(b)(c)(d)When the TOF is larger, the propellant consumption is smaller. These cases showed that the NEP is a good way to transport large payloads to Saturn. Considering the cases withVinf=5 km/s and time of flight of 5 years, the mass relation is between 0.86≤mf/m0≤0.955. Note also that the delivered mass increases significantly with trip time. Increasing the velocity at infinity, the mass ratio also increases. However, the propellant consumption is high when the time of flight is short.The angle of the thrust with the spacecraft velocity is shown in Figure3(a) (continuous lines), however the angles of the thrust with the orbital plane take small values due to the position of the Earth and Saturn near the ecliptic plane (discontinuous lines).(a) Angle to the spacecraft velocity (continuous line) and angle to the orbital plane (discontinuous line). (b) Thrust vector (continuous line) andα/g and specific impulse (discontinuous line). Both for effective specific power of 50 W/kg (NEP), TOF=10 years and Vinf=5 km/s.
(a)(b)The behavior of the thrust is presented in Figure3(b) (continuous lines). It is visible that the thrust has high values at the beginning of the transfers, the specific impulse (discontinuous lines) shows an opposite relative behavior, that is, when the thrust has high values the specific impulse has low values. The propulsion system offers values of variable thrust.


### 4.2. Radioisotope Electric Propulsion (REP)

The Deep Space 1 mission used a radioisotope thermoelectric generator that was combined with off-the-shelf ion propulsion systems. This combination provides a combined specific mass of almost300 kg/kW. However, advanced radioisotope power system that is under development could achieve the specific mass of 150 kg/kW or lower required for REP. Figure 4 shows the trajectory of the flight Earth-Saturn projected on the ecliptic plane.Figure 4
Trajectory of the flight Earth-Saturn with effective specific power (REP) of10 W/kg and Vinf=10 km/s on Earth’s sphere influence (TOF=12 years).In this section, some simulations for the REP system are shown. Figure5 shows the use of radioisotope power source. When the TOF is larger, the propellant consumption is smaller. These cases showed that the REP is another way to transport payloads to Saturn. For the cases where Vinf=5 km/s and the time of flight is 5 years, the mass ratio are between 0.45≤mf/m0≤0.67.Curves for REP as a function of the ratiomf/m0 for several TOFs. The effective specific powers are (a) 10 W/kg, (b) 8 W/kg, (c) 6 W/kg, (d) 5 W/kg.
(a)(b)(c)(d)The angle of the thrust with the spacecraft velocity is shown in Figure6(a) (continuous lines) as well as the angle that the thrust makes with the orbital plane (discontinuous lines).(a) Angle to the spacraft velocity (continuous line) and angle to the orbital plane (discontinuous line). (b) Thrust vector (continuous line) andα/g and specific impulse (discontinuous line). Both for effective specific power of 10 W/kg (REP), TOF=10 years, and Vinf=5 km/s.
(a)(b)The behavior of the thrust is presented in Figure6(b) (continuous lines) and the specific impulse (discontinuous line).


### 4.3. Solar Electric Propulsion (SEP)

In October 1998, NASA launched the Deep Space 1 that was the first interplanetary mission to be propelled by solar electric propulsion (Rayman and Williams [6], Brophy and Noca [1]). Later, Smart-1 (ESA-2003) and Hayabusa (Japan-2003) were launched. With the successful demonstration made by the Deep Space 1, many studies have been performed to show the applicability and performance of Solar Electric Propulsion (SEP) for interplanetary missions (Brophy and Noca [1], Racca [7]). Previously, we did not consider this option for a trip to Saturn (Solórzano, Prado, and Sukhanov [8]).It is known that the weakness of the SEP technology is the low levels of acceleration that it provides and in the reduced solar irradiance available for photovoltaic power generation at the outer reaches of the solar system. Nevertheless, these drawbacks can be avoided by a suitable design that allows the SEP system to operate efficiently for long periods using a wide range of input powers (Mengali and Quarta [9]).Here, we studied the possible advantages of SEP to reach Saturn. Figure7 shows the trajectory of the flight Earth-Saturn projected on the ecliptic plane for the SEP.Figure 7
Trajectory for the flight Earth-Saturn with effective specific power (SEP) of2 W/kg and Vinf=10 km/s on Earth’s sphere influence (TOF=12 years).Figure8 shows the behavior of the relation final mass/initial mass versus Vinf. It is visible that the relation of mass increases with the effective specific power. For example, when Vinf≅0 km/s and the specific initial power is 2 W/kg, the mass ratio is between 0.027 and 0.064 (5 years ≤TOF≤12 years). For the case when Vinf≅0 km/s and specific initial power is 12 W/kg, the mass relation is between 0.33 and 0.73 (5 years ≤TOF≤12 years).Curves for SEP as a function of the ratiomf/m0 for several TOFs. The effective specific powers are (a) 12 W/kg, (b) 8 W/kg, (c) 4 W/kg, (d) 2 W/kg.
(a)(b)(c)(d)The angle of the thrust with respect of the spacecraft velocity is shown in Figure9(a) (continuous lines), however the angles of the thrust with the orbital plane take small values (discontinuous lines).(a) Angle to the spacraft velocity (continuous line) and angle to the orbital plane (discontinuous line). (b) Thrust vector (continuous line) andα/g and specific impulse (discontinuous line). Both for effective specific power of 2 W/kg (SEP), TOF=10 years, and Vinf=5 km/s.
(a)(b)The behavior of the thrust is presented in Figure9(b) (continuous lines). It is visible that the thrust has high values at the beginning of the transfers. The specific impulse (discontinuous lines) shows an opposite behavior. The maximum peak of the specific impulse happens when the thrust angle and the orbital plane suffer significant changes.


## 4.1. Nuclear Electric Propulsion

Nuclear Electric Propulsion (NEP) uses a reactor power system to provide the electricity for thrusters that ionize and accelerate propellant to produce thrust. The application of nuclear systems in space has advantage in situations where the distances from the Sun are large and the solar power density is too low≳1.5AU) and locations where solar power is not readily or continuously available. One of the key performances for NEP is the power system specific mass, measured in kg/kW. Low values are desired in order to provide the maximum mass allocation for payload or propellant. High values result in minimal delivered payloads.In this section, the study of the NEP system for a trip to Saturn is performed. Figure1 shows the trajectory of the flight Earth-Saturn projected on the ecliptic plane. Figure 2 shows that the use of NEP allows the delivery of larger masses to the target planet, due to the high power of this system.Figure 1
Trajectory of the flight Earth-Saturn with effective specific power (NEP) of 50 W/kg andVinf=10 km/s on Earth’s sphere influence (TOF=12years).Curves for NEP as a function of the ratiomf/m0 for several TOFs. The effective specific powers are (a) 100 W/kg, (b) 80 W/kg, (c) 60 W/kg, (d) 30 W/kg.
(a)(b)(c)(d)When the TOF is larger, the propellant consumption is smaller. These cases showed that the NEP is a good way to transport large payloads to Saturn. Considering the cases withVinf=5 km/s and time of flight of 5 years, the mass relation is between 0.86≤mf/m0≤0.955. Note also that the delivered mass increases significantly with trip time. Increasing the velocity at infinity, the mass ratio also increases. However, the propellant consumption is high when the time of flight is short.The angle of the thrust with the spacecraft velocity is shown in Figure3(a) (continuous lines), however the angles of the thrust with the orbital plane take small values due to the position of the Earth and Saturn near the ecliptic plane (discontinuous lines).(a) Angle to the spacecraft velocity (continuous line) and angle to the orbital plane (discontinuous line). (b) Thrust vector (continuous line) andα/g and specific impulse (discontinuous line). Both for effective specific power of 50 W/kg (NEP), TOF=10 years and Vinf=5 km/s.
(a)(b)The behavior of the thrust is presented in Figure3(b) (continuous lines). It is visible that the thrust has high values at the beginning of the transfers, the specific impulse (discontinuous lines) shows an opposite relative behavior, that is, when the thrust has high values the specific impulse has low values. The propulsion system offers values of variable thrust.


## 4.2. Radioisotope Electric Propulsion (REP)

The Deep Space 1 mission used a radioisotope thermoelectric generator that was combined with off-the-shelf ion propulsion systems. This combination provides a combined specific mass of almost300 kg/kW. However, advanced radioisotope power system that is under development could achieve the specific mass of 150 kg/kW or lower required for REP. Figure 4 shows the trajectory of the flight Earth-Saturn projected on the ecliptic plane.Figure 4
Trajectory of the flight Earth-Saturn with effective specific power (REP) of10 W/kg and Vinf=10 km/s on Earth’s sphere influence (TOF=12 years).In this section, some simulations for the REP system are shown. Figure5 shows the use of radioisotope power source. When the TOF is larger, the propellant consumption is smaller. These cases showed that the REP is another way to transport payloads to Saturn. For the cases where Vinf=5 km/s and the time of flight is 5 years, the mass ratio are between 0.45≤mf/m0≤0.67.Curves for REP as a function of the ratiomf/m0 for several TOFs. The effective specific powers are (a) 10 W/kg, (b) 8 W/kg, (c) 6 W/kg, (d) 5 W/kg.
(a)(b)(c)(d)The angle of the thrust with the spacecraft velocity is shown in Figure6(a) (continuous lines) as well as the angle that the thrust makes with the orbital plane (discontinuous lines).(a) Angle to the spacraft velocity (continuous line) and angle to the orbital plane (discontinuous line). (b) Thrust vector (continuous line) andα/g and specific impulse (discontinuous line). Both for effective specific power of 10 W/kg (REP), TOF=10 years, and Vinf=5 km/s.
(a)(b)The behavior of the thrust is presented in Figure6(b) (continuous lines) and the specific impulse (discontinuous line).


## 4.3. Solar Electric Propulsion (SEP)

In October 1998, NASA launched the Deep Space 1 that was the first interplanetary mission to be propelled by solar electric propulsion (Rayman and Williams [6], Brophy and Noca [1]). Later, Smart-1 (ESA-2003) and Hayabusa (Japan-2003) were launched. With the successful demonstration made by the Deep Space 1, many studies have been performed to show the applicability and performance of Solar Electric Propulsion (SEP) for interplanetary missions (Brophy and Noca [1], Racca [7]). Previously, we did not consider this option for a trip to Saturn (Solórzano, Prado, and Sukhanov [8]).It is known that the weakness of the SEP technology is the low levels of acceleration that it provides and in the reduced solar irradiance available for photovoltaic power generation at the outer reaches of the solar system. Nevertheless, these drawbacks can be avoided by a suitable design that allows the SEP system to operate efficiently for long periods using a wide range of input powers (Mengali and Quarta [9]).Here, we studied the possible advantages of SEP to reach Saturn. Figure7 shows the trajectory of the flight Earth-Saturn projected on the ecliptic plane for the SEP.Figure 7
Trajectory for the flight Earth-Saturn with effective specific power (SEP) of2 W/kg and Vinf=10 km/s on Earth’s sphere influence (TOF=12 years).Figure8 shows the behavior of the relation final mass/initial mass versus Vinf. It is visible that the relation of mass increases with the effective specific power. For example, when Vinf≅0 km/s and the specific initial power is 2 W/kg, the mass ratio is between 0.027 and 0.064 (5 years ≤TOF≤12 years). For the case when Vinf≅0 km/s and specific initial power is 12 W/kg, the mass relation is between 0.33 and 0.73 (5 years ≤TOF≤12 years).Curves for SEP as a function of the ratiomf/m0 for several TOFs. The effective specific powers are (a) 12 W/kg, (b) 8 W/kg, (c) 4 W/kg, (d) 2 W/kg.
(a)(b)(c)(d)The angle of the thrust with respect of the spacecraft velocity is shown in Figure9(a) (continuous lines), however the angles of the thrust with the orbital plane take small values (discontinuous lines).(a) Angle to the spacraft velocity (continuous line) and angle to the orbital plane (discontinuous line). (b) Thrust vector (continuous line) andα/g and specific impulse (discontinuous line). Both for effective specific power of 2 W/kg (SEP), TOF=10 years, and Vinf=5 km/s.
(a)(b)The behavior of the thrust is presented in Figure9(b) (continuous lines). It is visible that the thrust has high values at the beginning of the transfers. The specific impulse (discontinuous lines) shows an opposite behavior. The maximum peak of the specific impulse happens when the thrust angle and the orbital plane suffer significant changes.


## 5. Conclusions

Low-thrust transfers of the limited power type were considered in this paper. The method of the transporting trajectory was used with the reference orbit composed by a set of short arcs of the keplerian orbits, while the transfer trajectory is subjected to low thrust. Since a maximum power provides a minimum propellant consumption, our goal was to maximize themf/m0 ratio. The NEP system permits to send more payloads, when compared to other options. If larger payloads were required, a nuclear reactor powered system would be needed. NEP is especially applicable for short trip time, large launch masses, and high-energy missions. The REP system appears to be more attractive for higher values of Vinf, due to the fact that, with low values for the specific power it is possible to send larger payloads, when compared to the NEP system. For other values of the Vinf, the REP delivers less science payload with proportional less power available for science instruments. The SEP systems can deliver large payloads for the case where the specific power, time of flight, and Vinf are high. For the SEP systems a better option is a combination with gravity assisted maneuver. In general, maximum power provides a minimum propellant consumption.



---

*Source: 756037-2009-09-15.xml*