Some New Locally Optimal Control Laws for Sailcraft Dynamics in Heliocentric Orbits

The concept of solar sailing and its developing spacecraft is presented. The gravitational and solar radiation forces are considered. The effect of source of radiation pressure and the force due to coronal mass ejections and solar wind on the sailcraft configurations is modeled. Some analytical control laws with some mentioned input constraints for optimizing sailcraft dynamics in heliocentric orbit using lagrange’s planetary equations are obtained. Optimum force vector in a required direction is maximized by deriving optimal sail cone angle. Ignoring the absorbed and diffusely reflected parts of the radiation, some special cases are obtained. New control laws that maximize thrust to obtain certain required maximization in some particular orbital element are obtained.


Introduction
The theoretical concept of solar sails is old-standing and dates back to Johannes Kepler when he observed that comet tails point away from the Sun and he suggested that the Sun caused this propulsion effect.After that it transferred into practice by the Soviet pioneers of astronautics.Solar sails are a form of spacecraft propulsion using a combination of light (radiation pressure) and high speed ejected gasses from a star (e.g., solar wind and coronal mass ejection) to push large ultra-thin mirrors to high speeds.This concept of using solar radiation pressure as a means of propulsion for space vehicles was first introduced in the 1920s by the father of Russian astronautics, Tsiolkovsky [1] and Tsander [2].It is also used in the spacecraft attitude dynamics.Rizvi [3] he developed a control method for the solar sail normal vector to trace a desired circular coning trajectory at orbit rate.He finally concludes that the control torques can be applied to the sailcraft to enable orbit rate cone tracing of the sail normal and yield the desired orbital effects.
In his book, Wright [4] presents a detailed analysis on some possible solar sail applications.During his time at the Jet Propulsion Laboratory (JPL), Wright was actively involved in the planning of a rendezvous mission to comet Halley using solar sail technology.In 1977, a solar electric propulsion concept was selected instead, primarily because of technology maturity.Not long thereafter, the Halley rendezvous mission was dropped by NASA.All these concepts were first experienced with the Mariner 10 mission to Mercury and Venus.Mariner 10 was also the first spacecraft to use a gravity assist trajectory, accelerating as it entered the gravitational influence of Venus, then using the planet's gravity field to move onto a slightly different course to reach Mercury.Since then, there have been several attempts to realize a solar sail mission.
The sail concept is performed by gaining momentum from an ambient source, solar electromagnetic radiation.Using momentum gained only byreflecting ambient sunlight, the sail slowly but continuously accelerated to accomplish a wide-range of potential missions.Light sails could also be driven by energy beams to extend their range of operations, which is strictly beam sailing rather than solar sailing.Solar sailcrafts offer the possibility of low-cost operations combined with long operating lifetimes.Since they have few moving parts and use no propellant, they can potentially be used numerous times for delivery of payloads.
The mechanism of this concept refering to the momentum carried by individual photons is extremely small.Thus, to provide a suitably large momentum transfer, we require the sail to have a large surface, while maintaining mass as low a as possible.At best, a solar sail will experience only 9 N of force per square km of sail located in Earth orbit.Adding the impulse due to incident and reflected photons, it is found that the thrust vector is directed normally to the surface of the sail; hence, by controlling the orientation of the sail relative to the Sun, we can gain or lose orbital angular momentum.
A true solar sail has yet to fly; however, significant steps have been taken since the beginning of the 1990s.In February 1993, under the guidance of Vladimir Syromiatnikov, the Russian Space Regatta Consortium deployed a 20 m spinning reflector, Znamya, from a progress supply vehicle.Observed from the Mir station, this deployment showed that spin deployment could be controlled by passive means.In May 1996, a large deployable reflector was demonstrated during the shuttle mission STS-77.The 14 m inflatable antenna experiment was primarily designed as a radio frequency reflector, but the promise of inflatable technology towards solar sailing was clearly demonstrated despite mission anomalies.In Köln, in December 1999, the German space agency, DLR, in association with ESA deployed a square, 20 m solar sail.This deployment now forms the basis for a future DLR/ESA in-orbit deployment demonstration perhaps in early 2006.Furthermore, in August 2004, the Institute of Space and Astronautical Science in Japan, ISAS, deployed two solar sails in space from a S-310 sounding rocket.NASA and several private enterprises all seeking to advance solar sailing are also conducting further work.It is thus clear that the technology is currently undergoing a revival in interest and a renewed drive towards flight status.
The literature on sailing dynamics, development, and attitude control is wealthy and it is of great interest to sketch some important ideas on these problems.The first solar sail trajectories were calculated by Tsu [5] and London [6].Tsu investigated various means of propulsion and showed that in many cases solar sails show superior performance when compared to chemical and ion propulsion systems.He used approximated heliocentric equations of motion to obtain spiraling trajectories.London presented similar spiral solutions for Earth-Mars transfers with constant sail orientation using the exact equations of motion.Optimal solar sail trajectories were first computed by Zhukov and Lebedev [7] for interplanetary missions between coplanar circular orbits.In 1980, Jayaraman [8] published similar minimumtime trajectories for transfers between the Earth and Mars.Two years later, Wood et al. [9] presented an analytical proof to show that the orbital transfer times obtained by Jayaraman [8] were incorrect due to the incorrect application of a transversality condition of variational calculus and an erroneous control law.Powers et al. [10] and Powers and Coverstone [11] obtained results similar to those reported in Wood's paper and obtained solutions for transfers to synchronous orbits.The more general time-optimal control problem of three-dimensional, inclined and elliptic departure and rendezvous planet orbits was discussed by Sauer [12].Hughes and McInnes [13] used genetic algorithms and sequential quadratic programming to obtain interplanetary trajectories via a direct method.Dachwald [14] presented a novel approach based on evolutionary neurocontrollers (ENC) to calculate optimal solar sail trajectories for interplanetary missions.
Wie [15] considered the orbital dynamics of a solar sail in the Earth-Sun circular restricted 3-body problem.The equations of motion of the sail are given by a set of nonlinear autonomous ordinary differential equations, which are nonconservative due to the noncentral nature of the force on the sail.They found there are equilibria admitting homoclinic paths, where the stable and unstable invariant manifolds are identical.As well as these periodic orbits about these equilibria also admit homoclinic paths.
Waters and McInnes [16] developed an attitude control systems for solar sail spacecraft are presented.He analysed a sailcraft in an Earth centered elliptic orbit, with particular emphasis on the significant effect of a solar-pressure disturbance torque (caused by an uncertain center-of-mass and center-of-pressure offset).
Gong et al. [17] investigated the time-optimal interplanetary transfer trajectories to a circular orbit of given inclination and radius.They derived optimal control law from the principle of maximization.An indirect method is used to solve the optimal control problem by selecting values for the initial adjoint variables.The conditions for the existence of the time-optimal transfer are dependent on the lightness number of the sail and the inclination and radius of the target orbit.A numerical method is used to obtain the boundary values for the time-optimal transfer trajectories.For the cases where no time-optimal transfer trajectories exist, first-order necessary conditions of the optimal control are proposed to obtain feasible solutions.The results show that the transfer time decreases as the minimum distance from the Sun decreases during the transfer duration.For a solar sail with a small lightness number, the transfer time may be evaluated analytically for a three-phase transfer trajectory.The analytical results are compared with previous results and the associated numerical results.The transfer time of the numerical result is smaller than the transfer time from previous results and is larger than the analytical result.
In contrast to ordinary ballistic navigation, the mission analysis for solar sail trajectories is not a simple task.The difference is that ballistic navigation is determined by a finite sequence of manoeuvres, with each of them being fully defined by Keplerian theory; 5 orbital parameters (the time they are executed, the intensity of the pulse, and the 3 angles determining its direction).But solar sailing, on the contrary, is a continuous process and the mathematics for reaching final conditions is far more complex, even more than in the case of electric propulsion, since the thrust intensity is related to the sail orientation: the orientation of the force vector applying to a perfectly reflective solar sail is normal to the sail, in the antisolar direction.Its intensity is proportional to the square cosine of the angle between the normal to the sail and the sunline.

Modeling the Source of Radiation Pressure
The absorbed energy heats the sail, which reradiates that energy from the front and rear surfaces.A solar sail is an oriented surface such that the acceleration experienced by the surface is a function of the surfaces attitude.The variation of the solar radiation pressure with solar distance can be approximated by an inverse square variation relationship; this assumption breaks down at low solar radius, when the finite angular size of the solar disc must be considered.The modeling of the source of radiation pressure is distinct and independent from the modeling of solar radiation pressure force, which is dependent on the optical properties of the sail surface.
Initially, we assume the solar disc has uniform brightness, such that an element of the solar disc will appear equally bright when viewed from any aspect angle.Thus, the specific intensity   is time independent and isotropic across the solar disc.Therefore, the solar radiation pressure exerted on a radially oriented, perfectly reflecting sail at a heliocentric distance , can be written as where is the velocity of light and  0 is the angular radius of the solar disc that is given by  0 = sin −1 ( Sun /) as shown in Figure 1.The geometric symmetry about the azimuth angle  leads to  independent of .Thus, (1) reduces to where  0 is the frequency integrated-specific intensity.

The Force Model on the Sailcraft Configurations
where  = 4.563 × 10 −6 N/m 2 is the nominal solar-radiation pressure constant 1 AU from the Sun,  is the surface area, n is a unit vector normal to the sail film, and ê is a unit vector pointing from the Sun to the sail, as shown in Figure 2. The solar radiation pressure varies inversely with the square of the distance from the Sun.The unit vector ê , as is clear from the geometry depicted in Figure 1, can be written as a linear combination of n and a unit vector transverse to it, t; that is, ê = cos n + cos  t.Therefore, (3) can be written as Now, we can rewrite (4) in the coordinate system ê ê ê , where ê is a unit vector perpendicular to ê , and ê completes the right-handed coordinate system (STW).Thus, the thrust unit vector (sail normal) in the three-dimensional space is defined by Now, the force due solar radiation pressure becomes Carrying out the dot products, we get

The Forces due to Coronal Mass Ejections and Solar Wind.
Mass is carried away from the Sun into interplanetary space, primarily in the form of protons, through several mechanisms: solar wind, solar flares, and coronal mass ejections.
The majority of the solar mass carried away is in the solar wind, a persistent stream of plasma particles generated in the photosphere of the Sun that flows out in all directions.At the Earth, the speed of the solar wind is typically ∼400 km/s with a particle density of about 10 protons/cm 3 , Breen et al. [18].When Mariner 10 flew by Mercury in 1974 it measured a solar wind speed of 423 km/s and a quiescent solar wind particle density of about 60 protons/cm 3 .Assuming the following simplifying assumptions: (a) spherical envelope for the angular span of ejections, and (b) all particle collisions with the sailcraft are inelastic.The perturbing accelerations on the sailcraft, due to the coronal mass ejections, and due to solar wind respectively are given by Abd El-Salam [19]: where (  /) SW ,  SW and V SW are the th area-mass ratio exposed to the solar wind stream, particle density, and the velocity of the solar wind and   is the mass of the proton, while (  /) CME ,  CME ,  CME and V CME are the area-mass ratio exposed to the coronal mass ejections, the overall volume, the mass, and the velocity of the coronal mass ejections, respectively.Adding ( 8) and ( 7), we obtain the force model due to the solar radiation pressure, solar wind, and the coronal mass ejection:

Optimum Force Vector
Referring to Figure 2, we see the orientation of the solar sail, and so the thrust force vector is described relative to the sunline by the sail pitch angle, , and clock angle, .We note that the sail control angles  and  can be defined as either Now, the sail thrust vector is defined by the cone and clock angles in the radial-transverse-normal (STW) frame.In order to optimize the sail control angles, we define a required direction, ξ, along which the component of the sail thrust is to be maximized: The force in this required direction, namely, ξ-direction is Using (11), we can obtain Arranging the terms and using the trigonometric identities, we get Consider the following definitions: Equation ( 14) can be simplified to Differentiating ( 16) with respect to the cone angle and finding the turning points, (  /), then we can obtain the optimal sail cone angle, which maximizes the force in the required direction as Dividing on 3cos 2  sin  yields which can be simplified to Equating the numerator to zero and rearranging the terms yield tan α cos ( − δ) A 2 cos α Multiplying both sides by 1/cos 2  = 1 + tan 2  yields Rearranging the terms yields Dividing on cos α yields This is a transcendental equation which is so difficult to be solved analytically, but in the following, we will consider some special cases.

Special Case I.
Ignoring the absorbed and diffusely reflected parts of the radiation, of the coefficients (15) become Dividing on −A 1 yields Let us now define these new coefficients: Now, the cubic equation (25) becomes which has the solution in the maximizing direction where 2 ) Since tan  * 2,3 given by (29) are imaginary solutions, therefore, these two roots are rejected.And thus the only considered root is the first root tan  * 1 .

Special Case II.
In addition to the simplification addressed in the special Case I, If we consider also the solar wind and the coronal mass ejection does not depend on the cone angle .Then, (27) can be reduced to the following quadratic equation: The solution of this equation yields the required cone angle  =  * to maximize the force vector in the required direction: 4 tan α cos ( − δ) . (32)

Heliocentric Orbital Dynamics
In the first instance, the Sun-centered orbital dynamics of solar sail spacecraft will be described, since this is traditionally where solar sailing is applied and where the bulk of the mission concepts exists, due to the large ΔV requirements.The vector equation of motion of a solar sail spacecraft moving in a heliocentric orbit is defined by (33), where  is the position vector of the spacecraft with respect to the Sun at time .The gravitational parameter of the Sun is defined by .For an ideal sail, the thrust vector is aligned along the sail normal direction, : The parameter, , is known as the sail lightness number and is defined as the ratio of the local solar radiation pressure acceleration produced by the sail to the local solar gravitational acceleration.This number is a useful parameter since it is independent of the solar distance and defines the performance of the sail.Another useful performance parameter is the sail characteristic acceleration, ac, which is the solar radiation pressure-induced acceleration of an ideal sail, pitched face on to the Sun ( = 0) at 1 AU from the Sun.It is easy to calculate that a characteristic acceleration of 1.0 mm s-2 corresponds to a dimensionless sail lightness number of 0.1686.As has been discussed, the solar sail characteristic acceleration is dependent on the surface area and reflectivity of the sail assembly and the mass of the spacecraft.

Lagrange's Planetary Equations
The state of a spacecraft can be described by a vector of 6 orbital elements, namely, semimajor axis, , eccentricity, , inclination, , argument of perihelion, , right ascension of the ascending node, Ω, true anomaly, , or any other time elements.Theses 6 elements are equivalent to 6 Cartesian position and velocity components.To measure the rate of change in these elements, we use the very famous system of 6 first-order differential equations known as Lagrange planetary equations; see (35).They are used to solve the equations motion of the sailcraft.These equations are particularly useful when we want to maximize the rate of change of a particular orbital element.This is useful when we want to modify one orbital element, while leaving the other timeaveraged elements unchanged.One form of this system is given by where , , , Ω, , and  are the usual Keplerian orbital elements, semimajor axis, eccentricity, inclination, longitude of the ascending node, and true anomaly, respectively,  is the semilatus rectum,  the gravitational parameter, and  is the orbital mean motion.The components of the solar sail thrust are denoted by , ,  radial, transverse, and normal, respectively.All Lagrange planetary equations can be written in compact form as only one equation as follows: where X(≡, , , Ω, , ) denotes any element of the Keplerian orbital elements.F, Ξ X are the force vector and primer vector of optimization, respectively.These vectors can be conveniently written as where The primer vector in the direction of maximizing the semimajor axis is The primer vector in the direction of maximizing the eccentricity is The primer vector in the direction of maximizing the inclination is The primer vector in the direction of maximizing the longitude of the ascending node is The primer vector in the direction of maximizing the argument of periapsis is Using this system, we can maximize directly any orbital elements or any other dynamical orbital parameters such as, the radius of periapsis and apoapsis.For example, if the radius of apoapsis,   = (1 + ) is differentiated, we obtain

Some New Locally Optimal Control Laws
To maximize the thrust to obtain certain required maximization in some particular orbital element, we set (/)(X/) = 0 and thus obtain a new set of locally optimal control laws.These control laws cannot guarantee global optimality, and they are often termed closed-loop methods.Global optimality requires the use of numerical methods, and even then, the true optimum solution is hard to attain.

Optimal Control Law for the Rate of Change of Semimajor Axis. The primer vector components maximize the first Lagrange planetary equation for the rate of change of semimajor axis.
In what follows, we will find an optimum cone angle that maximizes the rate of change of semimajor axis: The force in this required direction is given by Consider the following definitions: Equation ( 46) can be simplified to If we differentiate (49) with respect to the cone angle and find the turning points , (  /), then we can obtain the optimal sail cone angle, which maximizes the force in the required direction as follows: Dividing on 3cos 2  sin  yields which can be written as Equating the numerator to zero and rearranging the terms yield Multiplying both sides by 1/cos Finally, we can obtain This is a transcendental equation which is so difficult to be solved analytically, but in the following, we will consider some special cases.

Special Case I.
Ignoring the absorbed and diffusely reflected parts of the radiation, of the coefficients (47) become Thus, (55) becomes Let us now define the new coefficients: Equation ( 57) becomes clear as a cubic equation: which has the same solution as given by ( 27), but upon replacing the coefficients B   with the newly defined coefficients K   we obtain where Since tan  * 2,3 are imaginary solution, therefore, these two roots are rejected.And thus the only considered root is the first root tan  * 1 .

Special Case II.
Ignoring the absorbed and diffusely reflected parts of the radiation and assuming that the solar wind and the coronal mass ejection do not depend on the cone angle, then (57) can be reduced to the following quadratic equation: Consider the following definitions: Equation ( 63) can be simplified to (71) In what follows, we will find an optimum cone angle that maximizes the rate of change of inclination: from which one can obtain Setting  = tan 2 , the above equation can be written as which has the solution where (77) As a special case when ignoring the diffusive part of radiation, we obtain However, closer orbits to the Sun have shorter orbit periods, and so the time to achieve an overall inclination change is shorter.This optimal angle is 35.26 ∘ , which enables us to maximize angular momentum.
which has the same solutions as those maximizing the inclination and is in good agreement with the previous results.If we substitute this control law into the ascending node in Lagrange planetary equations and integrate over one orbit, we can obtain the change in ascending node per orbit.We find that the change in ascending node per orbit is independent of orbit radius and only depends on the sail lightness number as follows: ΔΩ = ( 88.2 sin  ) (degrees per orbit) .
However, closer orbits to the Sun have shorter orbit periods, and so the time to achieve an overall inclination change is shorter (Figure 5).(82)

Conclusion
We can conclude the new contributions of the present work as follows.The force model on the sailcraft configurations is amended with the solar wind and coronal mass ejection forces.The solar radiation pressure model is modified so as to include the absorbed, the specularly reflected, and the diffusely reflected into the model.Some analytical control laws with some mentioned input constraints for optimizing sailcraft dynamics in Heliocentric orbit using Lagrange's planetary equations are mathematically explored.Optimum force vector in a required direction is maximized by deriving optimal sail cone angle.Ignoring the absorbed and diffusely reflected parts of the radiation, some special cases are obtained.New control laws that maximize thrust to obtain certain required maximization in some particular orbital element are obtained.

Figure 1 :
Figure 1: Solar radiation pressure due to a finite solar disc.

P l a n e p e r p e n d i c u l a r t o t h e s u n l iFigure 2 :
Figure 2: Definition of sail control angles.

Figure 3 : 2 cos 2 1 . 7 . 5 .
Figure 3: (a) The optimum cone angle that maximizes the rate of change of semimajor axis when taking the positive sign in the numerator.(b) The optimum cone angle that maximizes the rate of change of semimajor axis when taking the negative sign in the numerator.

Figure 4 :Figure 5 :
Figure 4: (a) The optimum cone angle that maximizes the rate of change of eccentricity when taking the positive sign in the numerator.(b)The optimum cone angle that maximizes the rate of change of eccentricity when taking the negative sign in the numerator.

3𝜌 𝑠 [ 1 −
2tan 2 ] = 0, ⇒ (tan  * )  = 1 √ 2 ⇒  *  = 35.26∘ .(78) This result is in agreement with the previous results.If we substitute the control law into the inclination equation in Lagrange planetary equations and integrate over one orbit, we can obtain the change in inclination per orbit.We find that the change in inclination per orbit is independent of orbit radius and only depends on the sail lightness number (straight number relation) as follows: Δ = 88.2(degrees per orbit) .

Figure 6 :
Figure 6: (a) The optimum cone angle that maximizes the rate of change of argument of perihelion when taking the positive sign in the numerator.(b) The optimum cone angle that maximizes the rate of change of argument of perihelion when taking the negative sign in the numerator.

7. 6 .
Optimal Control Law for the Rate of Change of AscendingNode.The procedure is the same as in optimal control law for inclination
Solar radiation exerts a pressure on the th surface of area   of the solar sail due to photons impinging on a surface in space.If a fraction   of the impinging photons is absorbed, a fraction   is specularly reflected, and a fraction   is diffusely reflected by a surface, then we have   +   +   = 1.The force due to solar radiation pressure (in brief SRP) acting on such a flat Lambertian surface located one astronomical unit (1 AU) from the Sun is modeled as 3.1.The Force due to Solar Radiation Pressure.