SciELO - Scientific Electronic Library Online

 
vol.58 número1Precise Measurement of Orbital and Physical Parameters of Bright Detached Solar Analog Eclipsing BinariesOff-Axis Exact Ray Tracing Algorithm for Zero Coma Point Determination in Classical and Non-Classical Reflective Telescopes índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de astronomía y astrofísica

versión impresa ISSN 0185-1101

Rev. mex. astron. astrofis vol.58 no.1 Ciudad de México abr. 2022  Epub 23-Ene-2023

https://doi.org/10.22201/ia.01851101p.2022.58.01.02 

Articles

A Study on Different Attitude Strategies and Mission Parameters Based on Lightsail-2

L. G. Meireles1 

A. F. B. A. Prado1 

C. F. de Melo2 

M. C. Pereira2 

1Space Mechanics and Control Division, National Institute for Space Research, Brazil.

2Mechanical Engineering Department, Federal University of Minas Gerais, Brazil.


ABSTRACT

The Planetary Society’s LightSail-2 mission successfully validated the orbital maneuvering capability of a solar radiation pressure (SRP) propelled spacecraft. This paper presents a study on two alternative attitude strategies for the orientation of a solar sail. The goal is to increase the effect of the SRP acceleration over the spacecraft’s orbital trajectory, with the intention of maintaining or even gaining altitude over time. Furthermore, one of these strategies was employed while varying a few of the mission’s parameters to determine if it would be viable to maintain the spacecraft’s average altitude. Results show that it is possible to increase the average altitude of the spacecraft over time while still reducing the number of maneuvers necessary to change the spacecraft’s attitude. With that result in hand, it is also possible to change some of the mission parameters without compromising the solar sailing performance.

Key Words: celestial mechanics; methods; numerical; minor planets; asteroids; general; space vehicles

RESUMEN

La misión LightSail-2 de la Planetary Society validó con éxito la capacidad de maniobra orbital de una nave impulsada por la presión de la radiación solar (PRS). Presentamos un estudio sobre dos estrategias distintas para la orientación de una vela solar. El objetivo es aumentar el efecto de la aceleración por la PRS a lo largo de la trayectoria orbital de la nave, con la intención de mantener o incrementar con el tiempo la altitud. Además, una de las estrategias se usó para determinar la posibilidad de mantener la altitud promedio de la nave al variar algunos parámetros de la misión. Los resultados muestran que es posible aumentar con el tiempo la altitud promedio de la nave, mientras que se reduce el número de maniobras necesarias para cambiar la orientación de la nave. Con este resultado se ve que también es posible modificar algunos parámetros de la misión sin afectar el desempeño de la vela solar.

1. INTRODUCTION

Since the beginning of the space exploration age, space agencies and, most recently, private companies around the world, have faced the grand problem of vehicle autonomy and the limitations it imposes on a mission’s duration. So far, the most traditional and reliable method of space vehicle’s propulsion is chemical combustion. It requires that part of the overall space and total mass of the vehicle be reserved for the storage of fuel, which in turn will be used for required mission maneuvers. In this manner, it serves as a finite resource and is the main limiting factor of a mission’s length over time. Alternative propulsion systems have long been of interest for agencies and organizations with the purpose of overcoming these problems. A method of propulsion which uses an abundant energy available throughout the whole space in the inner Solar System was conceived in the early 1920s by Soviet scientists Konstantin Tsiolkovsky and Friedrich Zander (McInnes 2004). This energy is called the solar radiation pressure (SRP) and is simply the momentum of the photons emitted by the Sun. In spite of its simple concept, this technology requires a tremendous effort in materials engineering, in order to build structures large enough to absorb sufficient SRP while still having a limited total mass which facilitates its maneuverability (Vulpetti et al. 2015). These structures are denominated solar sails and have increasingly become target of studies (Forward 1990; Angrilli & Bortolami 1990; McInnes 1999; Wie 2007; Guerman et al. 2008; Zeng et al. 2015, 2016; Zhang & Zhou 2017; D’Ambrosio et al. 2019; Meireles 2019) and space missions (Johnson et al. 2011; Fernandez et al. 2014; Palla et al. 2017; Betts et al. 2017a; Viquerat et al. 2015; McNutt et al. 2014; Heiligers et al. 2019; Mori et al. 2020) from the 1990s until today, in order to validate and further explore their propelling concept and promising mission possibilities (McInnes 2003a,b). A special mention must be given to the Japan Aerospace Exploration Agency (JAXA), which has been responsible for significant contributions to the advances and validation of solar sailing and the first ever mission to employ a solar sail: IKAROS (Mori et al. 2010; Tsuda et al. 2011; Mori et al. 2012; Funase et al. 2012; Tsuda et al. 2013). Lastly, the most recent successful mission which employs a solar sail is LightSail-2, from The Planetary Society (Betts et al. 2017b, 2019; Mansell et al. 2020; Spencer et al. 2021). Despite the mission’s successful demonstration of the solar sailing concept, its attitude configuration requires abrupt changes in the spacecraft’s orientation. Considering the difficulty in maneuvering the large surface of the sail with a control system with components sufficiently small to be embarked in the mission, these attitude changes require a greater settling time than ideal for the employed strategy. This causes a SRP acceleration in unwanted directions and jeopardizes the solar sail’s full potential (Plante et al. 2017). Based on this difficulty, this study proposes some alternative approaches, as presented in McInnes (2004), for the use of the solar sail by smoothing the attitude’s rate of change during a maneuver, while still maintaining a similar solar sailing performance.

2. SYSTEM DYNAMICS

The LightSail-2 spacecraft has a total mass m of approximately 5kg and its sail fully deployed has a surface area A of 32m2. Its initial orbit is circular with an altitude of 720km and an inclination of 24 (The Planetary Society 2020). Therefore, considering the Earth’s mean radius as 6378km, the initial osculating orbital elements are: semi-major axis of 7098km, eccentricity of 0.0 and inclination of 24. The three remaining orbital elements are considered to be zero, due to a lack of further information and their reduced importance for our study.

The dynamical model assumes a 4-body problem, including the Sun, the Earth and the Moon, besides the spacecraft. This means that some of the forces acting on the spacecraft are of gravitational nature, from the Sun, the Earth and the Moon. Nevertheless, a few other perturbation forces were considered in the dynamics. The gravitational potential of a non-spherical Earth was included by using up to the 4th order zonal harmonics perturbation acceleration (Bate et al. 1971). Atmospheric effects were also included to determine the drag forces acting on the spacecraft by taking into account the 1976 U.S. Standard Atmosphere atmospheric model (National Oceanic and Atmospheric Administration (NOAA);National Aeronautics and Space Administration (NASA); United States Air Force (USAF) 1976). Furthermore, a solar radiation pressure (SRP) net force was also considered. It is the propelling force, whose model is presented in details in the next section (Vulpetti et al. 2015). Another widely used model for the SRP force can be found in McInnes (2004). This choice was made based on the existence of legacy codes from previous works developed with the format presented in Vulpetti et al. (2015). Finally, the Earth’s shadow was considered as having a cylinder shaped umbra and no penumbra region to determine a light exposure coefficient for the SRP net force. The orbit was numerically integrated with the use of Cowell’s method (Bate et al. 1971) and a RADAU integrator (Everhart 1985).

LightSail-2 mission’s Attitude Determination and Control System (ADCS) is composed of 5 sun sensors, 4 magnetometers, 2 gyro sensors, 3 torque rods and 1 momentum wheel (Plante et al. 2017). However, this study does not take into consideration the delays from the actuators nor the errors from the sensors. The attitude values used in the orbital simulations are equal to the values proposed by the strategy, in order to compare the difference in their influence in the solar sailing performance.

2.1. Solar Sail Dynamics

To derive the dynamics of a solar sail propelled spacecraft it is convenient to define a Spacecraft Oriented Frame (SOF), as illustrated in Figure 1. The SOF’s X axis points in the radial direction of a heliocentric inertial frame. It has the same direction of the incoming sunlight, which is indicated by the unit vector u. The SOF’s Z axis is defined as having the same direction of the spacecraft’s orbital angular momentum, indicated by the unit vector h. The Y axis is defined in agreement with a coordinate system in dextrorotation. The sail’s orientation is represented by the unit vector n orthogonal to the solar sail surface. Two important angles are consequently defined, the azimuth α and the elevation δ, necessary to resolve n in the SOF. An important consideration is made: the n unit vector exists only in the opposite semi-space of the sunlight beam or, in other words, opposite to the sunlit layer of the solar sail. This limits both the azimuth and elevation angles to the interval [−90,90].

Fig. 1 Spacecraft oriented frame (SOF) referenced at the spacecraft’s barycenter, adapted from Vulpetti et al. (2015)

To determine the thrust T resulting from SRP, as a function of the sail’s orientation, it is furthermore useful to define the lightness vector L, conceived in the SOF as the impulsive acceleration normalized by the Sun’s local gravitational acceleration and defined as

L=(12σcσ)n×2rspecn×+frdiff+kαn+α+rdiffu (1)

where σ c is a constant referred to as critical loading (equation 2), σ is the sail loading (equation 3), n x is the SOF’s X axis n versor’s component, r spec is the specular reflectance coefficient, r diff is the diffuse reflectance coefficient, χ is the emission/diffusion coefficient (the subscript f refers to the front side of the solar sail), κ is a net thrust dimensionless factor that results from the absorbed and re-emitted power on both sides of the sail (equation 4), and a is the absorptivity coefficient.

The previously mentioned critical loading σ c (equation 2), sail loading σ (equation 3) and net thrust dimensionless factor κ (equation 4), are defined as

σc2I1AUc g1AU1.5368 g/m2, (2)

σmA, (3)

kXfϵfT-Xbϵb(T)ϵfT+ϵb(T) (4)

where I 1AU = 1366W/m2 is the energy flux emitted by the Sun at 1 Astronomic Unit, the average Sun-Earth distance, equal to 149597870700m, c ≈ 2,9979 × 108 m/s is the speed of light in a vacuum, g 1UA ≈ 5,930 × 10−3 m/s2 is the gravitational acceleration of the Sun at 1 Astronomical Unit, m is the spacecraft’s total mass, A is the solar sail surface area, T is the solar sail’s temperature and is the emittance coefficient as a function of T (again, the subscript f refers to the front side of the sail, while the subscript b refers to the back side). From equation 3 and LightSail-2 mission’s specifications, the sail loading is σ = 156.25g/m2.

Further considerations are made in regard of the optical coefficients absorptivity α, diffuse reflectance r diff and specular reflectance r spec. They are all considered to be a function of the angle of the incident light beam, as portrayed in Figure 2.

Fig. 2 Absorptivity a, diffuse reflectance r diff and specular reflectance r spec coefficients as a function of the incident light angle for an aluminium sail with a root mean square roughness equal to 20nm (Vulpetti & Scaglione 1999). The color figure can be viewed online. 

An incidence angle of 0 represents the case where the unit vectors u and n are parallel. Since both the reflectance coefficients are also a function of the sail’s reflective layer roughness, the worst case scenario for the graphs of a value of 20nm root mean square roughness was considered for the simulations.

The Jet Propulsion Laboratory (JPL) published some of its solar sail models coefficients obtained from experimental studies developed in the early 2000’s, which are presented in Table 1 (McInnes 2004). The values of and χ of a square solar sail were used for this study.

TABLE 1 JPL’S MODEL SOLAR SAIL OPTICAL COEFFICIENTS 

r rspec/r ϵf ϵb χf χb
Ideal 1 1 0 0 2/3 2/3
Square 0,88 0,94 0,05 0,55 0,79 0,55

3. ATTITUDE STRATEGIES

This section serves the purpose of explaining in details each of the three different strategies, as presented in McInnes (2004) and used to simulate the solar sail’s orientation and their consequences to the spacecraft’s orbital trajectory. The results obtained from the Planetary Society strategy are used as a reference for comparison to the other two strategies analyzed by this study.

3.1. Strategy 1: The Planetary Society

The Planetary Society initial strategy consists of operating the solar sail at an “on-off” regime (Mansell et al. 2020). It consists of maintaining the solar sail at full exposure to the sunlight whenever the spacecraft is traveling away from the Sun. In other words, the unit vector orthogonal to the sail n is parallel to u and, consequently, both the azimuth and elevation angles α and δ are equal to zero. On the other hand, the “off” regime consists of maintaining the solar sail without any exposure to the sunlight whenever the spacecraft is traveling in the direction of the Sun. This means that the unit vectors n and u are orthogonal between themselves and α and/or δ are equal to 90. The schematics of this strategy is illustrated in Figure 3.

Fig. 3 LightSail-2 orientation schematics, adapted from Betts et al. (2017b). The color figure can be viewed online. 

In order to maximize the thrust component in the SOF XY-plane, δ is always equal to zero. Consequently, only the value of α changes. The first day of simulation (with the purpose of zooming in the X axis) for this strategy is illustrated in Figure 4.

Fig. 4 Azimuth angle α as a function of time for Strategy 1. The color figure can be viewed online. 

It is possible to observe that two attitude maneuvers are necessary at every orbital revolution.

The first 200 days of the orbital dynamics of this system were simulated using Strategy 1, as well as the spacecraft’s altitude evolution, are presented in Figure 5.

Fig. 5 Spacecraft’s altitude as function of time for Strategy 1. The color figure can be viewed online. 

The blue region of the plot is in fact the altitude over time. It is seen as a blur due to the large scale of time covered in the X axis. As expected, it behaves as a sinusoidal wave enveloped by two other curves of the maximum and minimum altitude reached at every revolution, which in turn are seen as the dashed black lines. The dashed red line is a simple average of the envelope curves and is considered to be the average altitude as a function of the simulation instant in days h av(t[days]). It starts at a value of 711.6867km and decreases 5.5204km to 706.1663km after 200 days. The dashed orange line indicates h av(t) for a simulation with disabled solar sailing, meaning a null SRP net force, to demonstrate the importance of solar sailing for maintaining the spacecraft’s average altitude over time. After 200 days, h av(200) = 671.6718km, having lost an extra 34.4945km, or 7.25 times more, of average altitude in comparison with the situation where the solar sailing is active.

Figure 6 presents the actual highest and lowest points of the orbit’s altitude values for the LightSail-2 mission, publicly made available and directly taken from The Planetary Society’s link for the LightSail-2 Mission Control (The Planetary Society 2020). A similar chart can be seen in the same link. There are clear differences from the results seen in Figure 5, despite the fact that the general behavior of these figures are the same. This happens due to the simplified orbital model used in the simulation and to limited access to initial conditions information for the spacecraft. In addition, the model considered overestimates the solar radiation pressure net thrust in comparison to the drag force. Therefore, in order to approximate the simulation results to the LightSail-2’s mission published data, the SRP net thrust had to be scaled down by a factor of 7×10−2.

Fig. 6 LightSail-2’s apogee, perigee and average altitude curves as a function of time from LightSail-2 Mission Control (https://www.planetary.org/explore/projects/lightsail-solar-sailing/lightsail-mission-control.html). The color figure can be viewed online. 

The initial average altitude of LightSail-2 mission’s at sail’s deployment, on 23 July 2019, is 717.4105km. On 3 February 2020, 200 days after sail deployment, the average altitude is 711.6255km, in contrast to the simulated value of 706.1663km. Although a divergence from the simulated and actual trajectory data is not ideal, it is not a vital requirement for this study that they match in value. As is the case, the main result is a comparison from all the attitude strategies, and since they are all simulated using the same orbital models, the results from Strategy 1 serve mainly as a comparison reference for other strategies.

3.2. Strategy 2: Maximum Ɛ˙

Solar radiation pressure propulsion is a specific type of low-thrust. Therefore, the thrust is a function of the spacecraft’s orientation, but for solar sails, its magnitude is also a function of the spacecraft’s position. Therefore, an analysis of its trajectory as a function of its orientation is possible from a traditional approach (Keaton 1986) given this extra constraint |T| = f(n , r), where r is the spacecraft’s position in the heliocentric inertial frame (HIF).

From equation 5, it is known that the specific orbital energy’s rate of change over time is maximized when the net thrust component in the direction of the spacecraft’s velocity vector has its greatest magnitude. This does not necessarily mean (and often it is the case) that the sail and/or the net thrust vector should be pointed in the direction of travel. In equation 5 is the spacecraft’s velocity vector in the HIF

ε˙=r˙T1m. (5)

To determine this optimal strategy, it is necessary to search for the orientation which gives the maximum component of T at the direction. Considering one point of the simulation at a time, that is to say, having a fixed position r and velocity , this is a well bound one dimensional convex problem, with α constrained in the interval [−90,90]. A simple hill climbing search algorithm is sufficient, and was used in this study to find the solutions. The first day of simulation (again, with the purpose of zooming in the X axis) for this strategy is illustrated in Figure 7.

Fig. 7 Azimuth angle α as a function of time for Strategy 2. The color figure can be viewed online. 

In contrast to the limited number of attitude maneuvers per orbital revolution from Strategy 1, this approach requires a constant control over the sail’s orientation. Nevertheless, it is possible to observe an interesting characteristic of this strategy’s behavior: the sail completes half a rotation (the azimuth angle α rises from −90 to 90) for every orbital revolution. This is a crucial information considered in the development of Strategy 3. It is also important to note that, from the definition of unit vector n (it exists in the semi-space opposite to the incoming sunbeam and, put differently, it always points away from the sail’s side not exposed to the sunlight), the transition of value from α = 90 to α = −90 does not mean a change of 180 in the sail’s orientation. In this situation, the unit vector n merely changes its side back to the one in the shadow. Furthermore, for this consideration to work with the simulations in this study, it is necessary for the solar sail to have equally reflective sides.

Figure 8 displays the same quantities presented in Figure 5, but this time, with the employment of Strategy 2.

Fig. 8 Spacecraft’s altitude as function of time for Strategy 2. The color figure can be viewed online. 

There is a noticeable quantitative difference between the results of Strategies 1 and 2. The latter presents an average altitude after 200 days of 719.9000km, which is 13.7337km, or 1.94%, higher. But there is also a remarkable contrast between both simulations. While Strategy 1 decreases its average altitude over time, Strategy 2 actually increases its value. Starting from 711.7190km, it gains 8.1810km over the course of 200 days. This result is of great value given the objective of LightSail-2 mission in demonstrating the hidden potential of SRP in maneuvering a spacecraft. The more the solar sail is capable of increasing the average altitude of the spacecraft, the more the concept of solar sailing gains relevance.

3.3. Strategy 3: Constant α˙

From the desire of maintaining a performance in altitude gain close to the one obtained from Strategy 2 while trying to reduce the number of attitude maneuvers, in order to maintain a similarity to LightSail-2 mission’s approach, a third strategy was analyzed. It consists of keeping the sail at a constant rate of change of its azimuth angle α˙, with a value that guarantees it completes half a rotation for every orbital revolution, as observed in Strategy 2. But there is a minor problem with this consideration. Since the goal is to increase the spacecraft’s altitude and, consequently, its specific orbital energy, then the period of revolution will also increase. This would mean that the α˙ implemented is no longer able to maintain the half rotation per orbital revolution desired. Therefore, it is necessary to perform attitude maneuvers once in a while to correct α˙ and α. More specifically, an attitude maneuver will only be performed after a complete revolution, which means after an integer number of revolutions N rev. As a reference, the maneuvers are performed at the point where the spacecraft is traveling directly in the direction of the Sun. In other words, when its velocity vector points in the opposite direction of its heliocentric position vector. At this instant, the desired value of α is always −90.

Figure 9 illustrates some examples of the proposed strategy. Every instant of an attitude maneuver is indicated with a red cross in the graph. Once again, only the first day of simulation is presented with the intention of zooming on the X axis values and having a better view for further analysis. The top graph corresponds to an attitude maneuver once every orbital revolution (N rev = 1). The bottom graph corresponds to an attitude maneuver once every five orbital revolutions (N rev = 5).

Fig. 9 Azimuth angle α as a function of time for Strategy 3. N rev = 1 for the upper graph. N rev = 5 for the lower graph. The color figure can be viewed online. 

Figure 10 presents the azimuth angle rate of change α˙ from the examples in Figure 9, in the same time interval. Once again, the instant of the maneuver is indicated by a red cross.

Fig. 10 Azimuth angle rate of change α˙ as a function of time for Strategy 3. N rev = 1 for the upper graph. N rev = 5 for the lower graph. The color figure can be viewed online. 

From the data presented in Figure 10, it is interesting to observe the general behavior of α˙ throughout the simulation (Figure 11).

Fig. 11 Azimuth angle rate of change α˙ as a function of time for Strategy 3. N rev = 1 for the upper graph. N rev = 20 for the lower graph. The color figure can be viewed online. 

As ε increases, so does the orbital period. Therefore, α˙ is supposed to decrease in this situation. The inverse is also true. As ε oscillates, so does α˙. In the case of N rev = 20, is it possible to observe an overall increase in the average value of α˙, which is an indication of an overall loss of ε and, consequently, average altitude. An example of N rev = 20 was considered instead of N rev = 5 to make this even clearer in the graphical representation of the discrete range of values of α˙. As N rev increases, the longer the sail maintains its azimuth angle rate of change.

It became interesting to simulate a range of N rev to evaluate for how long the spacecraft can be left without any attitude maneuver and still have a better performance than Strategy 1. The average altitude after 200 days of simulation h av(200) and its gain compared to the initial average altitude h gain (equation 6) are presented in Table 2, as a function of N rev

hgaint=hav(t)hav(0)-1×100. (6)

TABLE 2 EVALUATION OF HAV(200) AS A FUNCTION OF NREV FOR STRATEGY 3 

Nrev hav(200) (km) hgain(200) (%)
1 711.9803 +0.0367
2 711.9994 +0.0394
3 712.0116 +0.0411
4 712.0119 +0.0412
5 712.0039 +0.0400
10 711.8236 +0.0147
15 711.3687 -0.0492
20 710.6479 -0.1505
25 709.6213 -0.2947
30 708.3483 -0.4736
35 706.7824 -0.6936
40 704.9903 -0.9454
45 702.9148 -1.2370
50 700.5899 -1.5637

Remarkably, only values of N rev greater than 35 resulted in h av(200) smaller than 706.1663km, obtained by Strategy 1. This means that over 70 times fewer maneuvers could be performed without a loss in the solar sail’s performance in keeping the spacecraft’s average altitude. An unexpected result is having a maximum value for N rev = 4. The expected maximum was found for N rev = 1. Despite this fact, the whole range of N rev from 1 up to 10 presents similar results for h av(200), all of which, as is the case of Strategy 2, show an increase of h av over time. This indicates that other conditions are maybe more significant for this range of values. Results from Table 2 are presented graphically in Figure 12.

Fig. 12 Average altitude after 200 days as a function of N rev for Strategy 3. The color figure can be viewed online. 

The blue dashed line at the top is the result obtained from Strategy 2. It can be considered as a maximum, yet unattainable, goal for Strategy 3. On the other hand, the red dashed line at the bottom is the result obtained from Strategy 1. Since it was already implemented in the LightSail-2 mission, it can be considered as a minimum to overcome. For values of N rev smaller than or equal to 35, Strategy 3 proved to be a successful attempt to improve the performance of a solar sail to gain altitude.

4. MISSION PARAMETERS ANALYSIS

In this section, a set of mission parameter values were modified to investigate their general behavior. The use of Strategy 2 is justified as being a theoretical limit to the maneuvering potential of solar sailing. Namely, it is the best case scenario.

An important consideration is made for the average altitude’s change of value over time. It is assumed to have a linear variation along the 200 days simulated. Therefore, it was possible to implement a linear least-squares regression for h av(t), as indicated in Figure 13 by the yellow dashed line. The slope of these regression lines is henceforth referred to as hav(200)t. It is noted that this yellow line follows very closely the red line.

Fig. 13 The spacecraft altitude as a function of time for Strategy 2. The color figure can be viewed online. 

As final considerations, after each sweep of the mission parameters, the regression line’s slope values are fitted to a given function, and are presented in the corresponding section.

4.1. Initial Altitude Sweep

An interesting analysis comes from the sweep of the mission’s initial altitude h(0). As indicated by Figure 5, the drag forces are greatly responsible for the spacecraft’s altitude decay. An also known fact is that, as the altitude decreases, the drag forces increase in magnitude. In turn, this would decrease the slope of the regression line. The inverse is also expected to be true. Therefore, a sweep of values in the [690,750]km interval was made (Figure 14). In addition, given the exponential nature of the atmospheric model considered, the data is fitted to the function presented in equation 7

f(x) = c1e-c2(x-c3) + c4, (7)

where x is the independent variable from known data, y is the dependent variable for the fitted curve and c n are the function’s parameters.

Fig. 14 Average altitude slope after 200 days as a function of the initial altitude h(0) for Strategy 2. The parameters of equation 7 are also displayed. The color figure can be viewed online. 

Figure 14 shows that Strategy 2 can avoid the decay of the spacecraft for initial altitudes close to 700km. Since the atmospheric model only considers that the atmosphere’s density is not null up to an altitude of approximately 865km, and therefore drag forces only exist until this altitude, it is interesting to extend the interval of analysis to verify that greater values of h(0) do not fall near the fitted curve and behave in a different manner, closer to a linear behavior (Figure 15).

Fig. 15 Average altitude after 200 days as a function of the initial altitude h(0) for Strategy 2 (extended view). Equation 7 function parameters are also displayed. The color figure can be viewed online. 

4.2. Spacecraft Mass Sweep

A fundamental parameter for solar sailing is the sail loading σ (equation 3). As the spacecraft’s mass m increases, so does σ, and consequently the sail’s capacity to maneuver the spacecraft decreases (equation 1). Given that Strategy 2 was able to increase h av(t) over time, an effort was made to determine how much m could be increased while keeping a constant sail size (and surface area A) and still maintaining a positive slope for the average altitude regression line (Figure 16). Given the nature of equation 1, the data is fitted to the function presented in equation 8

fx=c1x-c2+c3. (8)

Fig. 16 Average altitude after 200 days as a function of the spacecraft’s mass for Strategy 2. The parameters of equation 8 are also displayed. The color figure can be viewed online. 

The figure shows that, even for a massive spacecraft of nearly half a ton, the regression line’s slope still maintains a positive value. Despite being a fitted parameter, the fact that c 3 is greater than zero indicates a limit, given this mission’s initial conditions and the employment of Strategy 2, which would guarantee no decay for the spacecraft’s average altitude. It means that, in theory, any spacecraft can use this technique.

In the search for negative hav(200)t, the initial altitude was changed to smaller values and the same curve fit procedures were performed (Table 3). From Figure 17 it is possible to see that hav(200)t is more sensitive to variations of h(0) for lower m, as expected. In fact, a negative hav(200)t starts to be observed for negative values of c 1, which happens with h(0) lower than 700km. This indicates a minimum altitude where Strategy 2 can be employed to keep a positive hav(200)t for a wide range of spacecraft masses. The present study showed that the initial altitude of the spacecraft is the determinant variable to allow this technique to be used to keep the spacecraft in orbit, and the mass does not influence the results.

TABLE 3 CURVE FIT PARAMETERS FOR THE AVERAGE ALTITUDE* 

h(0) c1 c2 c3
720 0.185 0.530 0.000087
715 0.147 0.538 0.000099
710 0.107 0.526 0.000112
705 0.063 0.330 0.000118
700 0.019 -3.681 0.000104
695 -0.035 1.671 0.000133
690 -0.090 1.398 0.000128

* Slope of the regression line after 200 days as a function of the initial altitude h(0) for Strategy 2.

Fig. 17 Curve fit for the average altitude regression line’s slope after 200 days as a function of the spacecraft’s mass, for Strategy 2 and different initial altitudes h(0). The color figure can be viewed online. 

4.3. Orbital Plane Inclination Sweep

Solar sails are an intriguing propelling method with a unique peculiarity. Considering an heliocentric inertial frame, their resulting thrust will always have a radially positive component. As long as a sail is open, it will reflect the incident sunlight, transfer linear momentum to or from the spacecraft to produce a resulting thrust with a positive radial component, whether this is a desired effect or not. The only alternative to avoid the aforementioned limitation is to reduce the sail’s exposed surface area to zero. From a mission control perspective, this would mean increasing either the azimuth angle α or the elevation angle δ to 90.

When considering a problem with the objective of increasing the average altitude over time, the desired direction in which the resulting SRP acceleration component needs to be maximized is different for various values of the orbital plane inclination. This results in slightly different desired azimuth angles over the span of an entire orbital revolution. In addition, it would be interesting to investigate different mission scenarios, considering whatever engineering limitations or obligations an organization might have in determining the orbital plane inclination of its mission. Inevitably, the variation of this orbital parameter is another interesting analysis to be made (Figure 18). Given the trigonometric relation of this parameter to the system’s dynamics, the data is fitted to the form presented in equation 9

f(x) = c1 sin(c2(x - c3)) + c4. (9)

Fig. 18 Average altitude after 200 days as a function of the orbital plane inclination i for Strategy 2. The parameters of equation 9 are also displayed. The color figure can be viewed online. 

As is well known, a smaller exposed surface area results in a lower SRP net thrust. In turn, this means lower specific orbital energy gains or losses and, theoretically, a reduced solar sailing maneuverability. But this fact could be used to the mission advantage. In the situation of higher inclinations, the use of Strategy 2 is specially capable of making the best use of the minimization of losses. In other words, whenever the spacecraft is traveling in the direction of the Sun, the energy losses can be reduced even further, compensating the inferior energy gains. This justifies the increasing average altitude regression line slope, as seen in Figure 18. It also shows that the strategy works for all ranges of inclinations studied [0,90]o but, since hav(200)t increases with the inclination, higher inclinations can be used if the goal is to increase the altitude of the spacecraft. This is particularly interesting if an escape from the Earth is desired.

Once again, the initial altitude was changed to lower values, followed by the same curve fit procedures, in the search for a minimum orbital plane inclination where Strategy 2 is capable of maintaining a positive hav(200)t Table (4). These curves are displayed in Figure 19, which also indicates that, for each one, there is an interception point with the hav(200)t=0 axis, as well as an inclination value with a maximum hav(200)t. In turn, these particular points are presented in Figure 20.

TABLE 4 CURVE FIT PARAMETERS FOR THE AVERAGE ALTITUDE* 

h(0) c1 c2 c3 c4
720 0.080 0.022 8.799 0.015
710 0.097 0. 022 8.672 -0.009
700 0.118 0. 022 8.982 -0.037
690 0.152 0. 022 7.482 -0.079
680 0.208 0.021 4.076 -0.144
670 -0.365 0.018 -9.752 -0.312
660 2.664 0.007 -132.831 -2.626

* Slope of the regression line after 200 days as a function of the initial altitude h(0) for Strategy 2.

Fig. 19 Curve fit for the average altitude regression line’s slope after 200 days as a function of the orbital plane inclination i for Strategy 2 and different initial altitudes h(0). The color figure can be viewed online. 

Fig. 20 Average altitude regression line’s slope after 200 days maximum and null values as a function of the orbital plane inclination i and the initial altitude h(0). The color figure can be viewed online. 

Figure 20 also displays a dashed red line that represents the mean value of i for max hav(200)t from the fitted curves, which is approximately 79.8. In spite of this value, the verified best case scenario of i for the employment of Strategy 2 is an inclination of 90. In the bottom graph, the blue area represents regions of h(0) and i where Strategy 2 is able to maintain a positive hav(200)t, while the red area represents the opposite region, of negative hav(200)t. Two remarks have to be made in this case scenario. First, for h(0) > 720km, Strategy 2 is able to maintain an average altitude gain for any i, including equatorial orbits. Second, for lower h(0), the drag forces are greatly superior, reducing the linear behavior of h av(t) and making it very difficult to fit parameters into the same function used for the other h(0) cases. This can already be verified with the parameter values of h(0) = 660km in Table 4, which are different from the rest by an order of magnitude. Nevertheless, it was possible to verify that for h(0) < 650km there are no positive hav(200)t for any i, which indicates that this is a limiting h(0) where Strategy 2 cannot maintain an average altitude gain, even for the best value of the inclination.

5. CONCLUSION

This study had the purpose of investigating different attitude strategies for the LightSail-2 mission and how they affect the solar sail’s capacity of maintaining the spacecraft’s altitude over time. Strategy 1 effectively replicated the original approach from the LightSail-2 mission and its implementation served as a base of comparison. Strategy 2, as expected, achieved a maximum performance result, at the cost of constantly maneuvering the sail with small attitude corrections. This strategy made it even possible to increase the spacecraft’s average altitude over time. Strategy 3 presents a simpler alternative to implement a smaller number of attitude maneuvers to keep the solar sail at a desired performance. Its implementation could mean a decrease of more than 70 times the number of attitude maneuvers implemented in LightSail-2’s mission without losing its solar sail’s maneuverability. For cases of up to one maneuver every 10 orbital revolutions, the sail also proved to be able to increase the spacecraft’s average altitude over time.

Making use of Strategy 2’s better performance, the values of the spacecraft’s mass, initial altitude and orbital plane inclination were changed in desired intervals to investigate regions where the sail could still maintain an average altitude gain over time. It was determined that, for the same inclination and for a wide range of spacecraft’s masses of up to half a ton, the initial altitude could be reduced by 20km to a value of 700km. In addition, the sail’s performance could be increased with a raise in the inclination, with a best case scenario for the highest value of 90. In turn, with this inclination, the initial altitude could be reduced a further 50km to a value of 650km.

Having examined the final results achieved and given the initial objectives, this study was successful in analyzing alternative attitude strategies able to keep a desired solar sail performance, while still keeping a simple enough implementation for the attitude control system employed in the LightSail-2 mission.

Acknowledgements

First and foremost, it is important to acknowledge the work and efforts of The Planetary Society, an organization concerned not only with the popularization of space technology and science, an already noble task, but also concerned with the practical advancements of the field, making missions such as LightSail-2 possible, which depend so much on the visionary collaboration of so many scientists, engineers and enthusiasts. The authors would also wish to express their appreciation for the support provided by Grants # 406841/2016-0 and 301338/2016-7 from the National Council for Scientific and Technological Development (CNPq), Grants # 2018/19959-0 and 2016/24561-0 from São Paulo Research Foundation (FAPESP) and to the financial support from the National Council for the Improvement of Higher Education (CAPES).

REFERENCES

Angrilli, F. & Bortolami, S. 1990, ESA Journal, 14, 431 [ Links ]

Bate, R. R., Mueller, D. D., & White, J. E. 1971, Fundamentals of astrodynamics (New York, NY: Dover) [ Links ]

Betts, B., Nye, B., Vaughn, J., et al. 2017a, International Symposium on Solar Sailing, Kyoto, Japan [ Links ]

Betts, B., Spencer, D., Bellardo, J., et al. 2019, IAF Space Propulsion Symposium (Washington, DC: IAF), 709 [ Links ]

Betts, B., Spencer, D., Nye, B., et al. 2017b, International Symposium on Solar Sailing, Kyoto, Japan [ Links ]

D’Ambrosio, A., Circi, C., & Zeng, X. 2019, AdSpR, 63, 3691 [ Links ]

Everhart, E. 1985, IAUC 83, Dynamics of Comets: Their Origin and Evolution, ed. A. Carusi & G. B. Valsecchi (Dordrecht: ASSL), 115, 185 [LINK] [ Links ]

Fernandez, J. M., Visagie, L., Schenk, M., et al. 2014, AcAau, 103, 204 [ Links ]

Forward, R. L. 1990, JSpRo, 27, 411 [ Links ]

Funase, R., Kawaguchi, J., Mori, O., Sawada, H., & Tsuda, Y. 2012, 53rd AIAA/ASME/ASCE/AHS/ASC, Structures, Structural Dynamics and Materials Conference (Honolulu, HI: AIAA), https://doi.org/10.2514/6.2012-1748 [ Links ]

Guerman, A., Smirnov, G., & Pereira, M. C. 2008, Advances in the Astronautical Sciences, 129, 607 [ Links ]

Heiligers, J., Fernandez, J., Stohlman, O., & Wilkie, W. 2019, Astrodynamics, 3, 231 [ Links ]

Johnson, L., Whorton, M., Heaton, A., et al. 2011, AcAau, 68, 571, https://doi.org/10.1016/j.actaastro.2010.02.008 [ Links ]

Keaton, P. W. 1986, STIN 86 (Washington: NASA) [ Links ]

Mansell, J., Spencer, D. A., Plante, B., et al. 2020, AIAA Scitech 2020 Forum (Orlando, FL: AIAA), https://doi.org/10.2514/6.2020-2177 [ Links ]

McInnes, C. R. 1999, AcAau, 45, 567, https://doi.org/10.1016/S0094-5765(99)00177-0 [ Links ]

______. 2003a, RSPTA, 361, 2989, https://doi.org/10.1098/rsta.2003.1280 [ Links ]

______. 2003b, AdSpR, 31, 1971, https://doi.org/10.1016/S0273-1177(03)00172-8 [ Links ]

______. 2004, Solar Sailing (Technology, Dynamics and Missions Applications) (London, UK: Praxis Publishing) [ Links ]

McNutt, L., Johnson, L., Kahn, P., Castillo-Rogez, J., & Frick, A. 2014, AIAA Space 2014 Conference and Exposition (San Diego, CA: AIAA), https://doi.org/10.2514/6.2014-4435 [ Links ]

Meireles, L. G. 2019, Análise de modelos nãoideais da dinâmica de atitude de velas solares e seus efeitos no movimento orbita, Master Thesis, INPE, São José dos Campos, http://urlib.net/rep/8JMKD3MGP3W34R/3SUPEUSLinks ]

Mori, O., Matsumoto, J., Chujo, T., et al. 2020, Astrodynamics, 4, 233 [ Links ]

Mori, O., Sawada, H., Funase, R., et al. 2010, JSAST, 8, 4, https://doi.org/10.2322/tastj.8.To_4_25 [ Links ]

Mori, O., Tsuda, Y., Sawada, H., et al. 2012, JSAST, 10, 413, https://doi.org/10.2322/tastj.10.Po_4_13 [ Links ]

National Oceanic and Atmospheric Administration (NOAA), National Aeronautics and Space Administration (NASA), United States Air Force (USAF). 1976, US Standard Atmosphere, 1976 (Washington: NOAA-S/T 76-1562) [ Links ]

Palla, C., Kingston, J., & Hobbs, S. 2017, 7th European Conference on Space Debris (Darmstadt: ESA Space Debris Office) [ Links ]

Plante, B., Spencer, D., Betts, B., et al. 2017, Fourth International Symposium on Solar Sailing, Kyoto, Japan [ Links ]

Spencer, D. A., Betts, B., Bellardo, J. M., et al. 2021, AdSpR, 67, 2878, https://doi.org/10.1016/j.asr.2020. 06.029 [ Links ]

The Planetary Society. 2020, LightSail 2 Mission Control, https://www.planetary.org/explore/projects/lightsail-solar-sailing/lightsail-mission-control.htmlLinks ]

Tsuda, Y., Mori, O., Funase, R., et al. 2011, AcAau, 69, 833, https://doi.org/10.1016/j.actaastro.2011.06.005 [ Links ]

Tsuda, Y., Mori, O., Funase, R., et al. 2013, AcAau, 82, 183, https://doi.org/10.1016/j.actaastro.2012.03.032 [ Links ]

Viquerat, A., Schenk, M., Lappas, V., & Sanders, B. 2015, 2nd AIAA Spacecraft Structures Conference (Kissimmee, FL: AIAA), https://doi.org/10.2514/6.2015-1627 [ Links ]

Vulpetti, G. & Scaglione, S. 1999, AcAau, 44, 123, https://doi.org/10.1016/S0094-5765(99)00038-7 [ Links ]

Vulpetti, G., Johnson, L., & Matloff, G. L. 2015, Solar sails (New York, NY: Copernicus Books) [ Links ]

Wie, B. 2007, JSpRo, 44, 545, https://doi.org/10.2514/1.23084 [ Links ]

Zeng, X., Gong, S., Li, J., & Alfriend, K. T. 2016, JGCD, 39, 1223, https://doi.org/10.2514/1.G001061 [ Links ]

Zeng, X., Jiang, F., & Junfeng, L. 2015, RAA, 15, 597, https://doi.org/10.1088/1674-4527/15/4/011 [ Links ]

Zhang, X. & Zhou, C. 2017, ChJA, 30, 1719, https://doi.org/10.1016/j.cja.2017.08.006 [ Links ]

Received: April 22, 2021; Accepted: September 27, 2021

Cristiano F. de Melo and Maria Cecília Pereira: Federal University of Minas Gerais, Av. Presidente Antônio Carlos, 6627 - Pampulha, Belo Horizonte, Minas Gerais, 31270-901, Brazil (cristiano.fiorilo@demec.ufmg.br, cecilia@demec.ufmg.br).

Lucas G. Meireles and Antônio F. B. A. Prado: National Institute for Space Research, Av. dos Astronautas, 1758 - Jardim da Granja, São José dos Campos, São Paulo, 12227-010, Brazil (meireleslg@gmail.com, antonio.prado@inpe.br).

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License