# Some new locally optimal control laws for sailcraft dynamics in heliocentric orbits.

1. IntroductionThe 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 by reflecting 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 Koln, 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 minimum-time 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.

2. 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 [I.sub.r] 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 r, can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where is the velocity of light and [[theta].sub.0] is the angular radius of the solar disc that is given by [[theta].sub.0] = [sin.sup.-1]([R.sub.Sun]/r) as shown in Figure 1. The geometric symmetry about the azimuth angle [phi] leads to P independent of r. Thus, (1) reduces to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)

where [I.sub.0] is the frequency integrated-specific intensity.

3. The Force Model on the Sailcraft Configurations

3.1. The Force due to Solar Radiation Pressure. Solar radiation exerts a pressure on the ith surface of area [A.sub.i] of the solar sail due to photons impinging on a surface in space. If a fraction [[rho].sub.a] of the impinging photons is absorbed, a fraction [[rho].sub.s] is specularly reflected, and a fraction [[rho].sub.d] is diffusely reflected by a surface, then we have [[rho].sub.a] + [[rho].sub.s] + [[rho].sub.d] = l. The force due to solar radiation pressure (in brief SRP) acting on such a flat Lambertian surface located one astronomical unit (1AU) from the Sun is modeled as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)

where P = 4.563 x [10.sup.-6] N/[m.sup.2] is the nominal solar-radiation pressure constant 1 AU from the Sun, A is the surface area, [??] is a unit vector normal to the sail film, and [[??].sub.s] 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 [[??].sub.s], as is clear from the geometry depicted in Figure 1, can be written as a linear combination of [??] and a unit vector transverse to it, [??]; that is, [[??].sub.s] = cos [alpha][??] + cos [alpha][??]. Therefore, (3) can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

Now, we can rewrite (4) in the coordinate system [[??].sub.s][[??].sub.t][[??].sub.w], where [[??].sub.t] is a unit vector perpendicular to [[??].sub.s], and [[??].sub.w] completes the right-handed coordinate system (STW). Thus, the thrust unit vector (sail normal) in the three-dimensional space is defined by

[??] = cos [alpha][[??].sub.s] + sin [alpha] cos [delta][[??].sub.t] + sin [alpha] sin [delta][[??].sub.w]. (5)

Now, the force due solar radiation pressure becomes

[F.sub.SRP] = [PA.sub.i] cos [alpha]{([[rho].sub.a] + [[rho].sub.s])[[??].sub.s] + (2[[rho].sub.s]cos [alpha] + 2/3 [[rho].sub.d]) x [cos [alpha][[??].sub.s] + sin [alpha] cos [delta][[??].sub.t] + sin [alpha] sin [delta][[??].sub.w]]}. (6)

Carrying out the dot products, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)

3.2. 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.sup.3], Breen et al. [18]. When Mariner10flew by Mercuryin1974itmeasuredasolar wind speed of 423 km/s and a quiescent solar wind particle density of about 60 protons/[cm.sup.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]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (8)

where [([A.sub.i]/m).sub.SW], [n.sub.SW] and [v.sub.SW] are the ith area-mass ratio exposed to the solar wind stream, particle density, and the velocity of the solar wind and [M.sub.P] is the mass of the proton, while [([A.sub.i]/m).sub.CME], [[tau].sub.CME], [M.sub.CME] and [v.sub.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:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)

4. 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 sun-line by the sail pitch angle, [alpha], and clock angle, [delta]. We note that the sail control angles [alpha] and [delta] can be defined as either

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Now, thesailthrustvectorisdefinedbythe cone andclock angles in the radial-transverse-normal (ST W) 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:

[??] = cos [??][[??].sub.s] + sin [??] cos [??][[??].sub.t] + sin [??] sin [??][[??].sub.w]. (11)

The force in this required direction, namely, [??]-direction is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

Using (11), we can obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)

Arranging the terms and using the trigonometric identities, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

Consider the following definitions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

Equation (14) can be simplified to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

Differentiating (16) with respect to the cone angle and finding the turning points, ([partial derivative][F.sub.[xi]]/[partial derivative][alpha]), then we can obtain the optimal sail cone angle, which maximizes the force in the required direction as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

Dividing on 3[cos.sup.2][alpha] sin [alpha] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (18)

which can be simplified to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

Equating the numerator to zero and rearranging the terms yield

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)

Multiplying both sides by l/[cos.sup.2][alpha] = 1 + [tan.sup.2][alpha] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)

Rearranging the terms yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

Dividing on cos [??] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

This is a transcendental equation which is so difficult to be solved analytically, but in the following, we will consider some special cases.

4.1. Special Case I. Ignoring the absorbed and diffusely reflected parts of the radiation, of the coefficients (15) become

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

Dividing on -[A.sub.1] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

Let us now define these new coefficients:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)

Now, the cubic equation (25) becomes

[tan.sup.3][alpha] + [B.sub.1][tan.sup.2][alpha] + [B.sub.2] tan [alpha] + [B.sub.3] = 0, (27)

which has the solution in the maximizing direction [xi]

[(tan [[alpha].sup.*.sub.1]).sub.[xi]] = - 1/3 [B.sub.1] + B + [B.sup.2.sub.1] - 3[B.sub.2]/9B, (28)

[(tan [[alpha].sup.*.sub.2, 3]).sub.[xi]] = - 1/3 [B.sub.1] - 1/2 B - [B.sup.2.sub.1] - 3[B.sub.2]/18B [+ or -] i[square root of 3]/2 [B.sub.1] [B - [B.sup.2.sub.1] - 3[B.sub.2]/9B], (29)

where

B = [(1/6 [B.sub.1][B.sub.2] - 1/2 [B.sub.3] - 1/27 [B.sup.3.sub.1] + [(1/27 [B.sup.3.sub.1][B.sub.3] - 1/108 [B.sup.2.sub.1][B.sup.2.sub.2] - 1/6 [B.sub.1][B.sub.2][B.sub.3] + 1/27 [B.sup.3.sub.2] + 1/4 [B.sup.2.sub.2]).sup.1/2]).sup.1/3]. (30)

Since tan [[alpha].sup.*.sub.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 [[alpha].sup.*.sub.1].

4.2. 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 [alpha]. Then, (27) can be reduced to the following quadratic equation:

-2[A.sub.2] tan [??] cos ([delta] - [??]) [tan.sup.2][alpha] -3[A.sub.2] tan [alpha] + [A.sub.2] tan [??] cos ([delta] - [??]) = 0. (31)

The solution of this equation yields the required cone angle [alpha] = [[alpha].sup.*] to maximize the force vector in the required direction:

[(tan [[alpha].sup.*]).sub.[xi]] = -3 [+ or -] [square root of (9 + 8[tan.sup.2][??][cos.sup.2]([delta] - [??]))]/4 tan [??] cos ([delta] - [??]). (32)

5. 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 [DELTA]v requirements. The vector equation of motion of a solar sail spacecraft moving in a heliocentric orbit is defined by (33), where r is the position vector of the spacecraft with respect to the Sun at time t. The gravitational parameter of the Sun is defined by [mu].For an ideal sail, the thrust vector is aligned along the sail normal direction, n:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (33)

For the simple analysis to follow, it is appropriate to show how this equation is represented in two-dimensional polar coordinates, where [theta] is the azimuth angle from the x-axis, is the distance of the spacecraft from the Sun. The pitch or cone angle, [alpha], is defined as the angle between the sail normal and the radial vector.

Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (34)

The parameter, [beta], 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 ([alpha] = 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.

6. Lagrange's Planetary Equations

The state of a spacecraft can be described by a vector of 6 orbital elements, namely, semimajor axis, a, eccentricity, e, inclination, i, argument of perihelion, [omega], right ascension of the ascending node, [OMEGA], true anomaly, f, 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 time-averaged elements unchanged. One form of this system is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (35)

where a, e, i, [OMEGA], [omega], and f are the usual Keplerian orbital elements, semimajor axis, eccentricity, inclination, longitude of the ascending node, and true anomaly, respectively, p is the semilatus rectum, [mu] the gravitational parameter, and n is the orbital mean motion. The components of the solar sail thrust are denoted by S, T, W radial, transverse, and normal, respectively.

All Lagrange planetary equations can be written in compact form as only one equation as follows:

dX/df = F x [[XI].sub.X], (36)

where X([equivalent to] a, e, i, [OMEGA], [omega],f) denotes any element of the Keplerian orbital elements. F, [[XI].sup.X] are the force vector and primer vector of optimization, respectively. These vectors can be conveniently written as

F = S[[??].sub.s] + T[[??].sub.t] + W[[??].sub.w], [[XI].sup.X] = [[XI].sup.X.sub.s][[??].sub.s] + [[XI].sup.X.sub.t][[??].sub.t] + [[XI].sup.X.sub.w][[??].sub.w], (37)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (38)

The primer vector in the direction of maximizing the semimajor axis is

[[XI].sup.a] = [[XI].sup.a.sub.s][[??].sub.s] + [[XI].sup.a.sub.t][[??].sub.t] = 2p[r.sup.2]/[mu][(1 - [e.sup.2]).sup.2] [e sin f[[??].sub.s] + p/r [[??].sub.t]]. (39)

The primer vector in the direction of maximizing the eccentricity is

[[XI].sup.e] = [[XI].sup.e.sub.s][[??].sub.s] + [[XI].sup.e.sub.t][[??].sub.t] = [r.sup.2]/[mu] [sin f[[??].sub.s] + (1 + r/p (e + cos f))[[??].sub.t]]. (40)

The primer vector in the direction of maximizing the inclination is

[[XI].sup.i] = [[XI].sup.i.sub.w][[??].sub.w] = [r.sup.3]/[mu]p cos (f + [omega])[[??].sub.w]. (41)

The primer vector in the direction of maximizing the longitude of the ascending node is

[[XI].sup.[OMEGA]] = [[XI].sup.[OMEGA].sub.w][[??].sub.w] = [r.sup.3]/[mu]p sin i sin (f + [omega])[[??].sub.w]. (42)

The primer vector in the direction of maximizing the argument of periapsis is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (43)

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, [r.sub.A] = a(1 + e) is differentiated, we obtain

d[r.sub.A]/df = da/df (1 + e) a de/df. (44)

7. Some New Locally Optimal Control Laws

To maximize the thrust to obtain certain required maximization in some particular orbital element, we set ([partial derivative]/[partial derivative][alpha])(dX/df) = 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.

7.1. 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:

[partial derivative]/[partial derivative][alpha] (da/df) = 2p[r.sup.2]/[mu][(1 - [e.sup.2]).sup.2] [e sin f([partial derivative]S/[partial derivative][alpha]) + p/r ([partial derivative]T/[partial derivative][alpha])] = 0 [right arrow] [e sin f([partial derivative]S/[partial derivative][alpha]) + p/r ([partial derivative]T/[partial derivative][alpha])] = 0. (45)

The force in this required direction is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (46)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (47)

Consider the following definitions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (48)

Equation (46) can be simplified to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (49)

If we differentiate (49) with respect to the cone angle and find the turning points, ([partial derivative][F.sub.[xi]]/[partial derivative][alpha]), then we can obtain the optimal sail cone angle, which maximizes the force in the required direction as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (50)

Dividing on 3[cos.sup.2][alpha] sin [alpha] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (51)

which can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (52)

Equating the numerator to zero and rearranging the terms yield

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (53)

Multiplying both sides by l/[cos.sup.2][alpha] = l + [tan.sup.2][alpha] yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (54)

Finally, we can obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (55)

This is a transcendental equation which is so difficult to be solved analytically, but in the following, we will consider some special cases.

7.2. Special Case I. Ignoring the absorbed and diffusely reflected parts of the radiation, of the coefficients (47) become

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (56)

Thus, (55) becomes

[partial derivative]/[partial derivative][alpha] (da/dt) = [tan.sup.3][alpha] + 2[A.sub.2]/[A.sub.1] p/r cos [delta][tan.sup.2][alpha] + (1 + 3[A.sub.2]/[A.sub.1]) tan [alpha] - [A.sub.2]/[A.sub.1] p/r cos[delta] = 0. (57)

Let us now define the new coefficients:

[K.sub.1] = 2[A.sub.2]/[A.sub.1] p/r cos [delta], [K.sub.2] = (1 + 3[A.sub.2]/[A.sub.1]), [K.sub.3] = - [A.sub.2]/[A.sub.1] p/r cos[delta]. (58)

Equation (57) becomes clear as a cubic equation:

[tan.sup.3][alpha] + [K.sub.1][tan.sup.2][alpha] + [K.sub.2] tan [alpha] + [K.sub.3] = 0, (59)

which has the same solution as given by (27), but upon replacing the coefficients B's with the newly defined coefficients K's we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (60)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (61)

Since tan [[alpha].sup.*.sub.2, 3] are imaginary solution, therefore, these two roots are rejected. And thus the only considered root is the first root tan [[alpha].sup.*.sub.1].

7.3. 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:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (62)

which can be simplified to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (63)

Consider the following definitions:

[A.sub.0] = [([A.sub.i]/m).sub.SW] [n.sub.SW][M.sub.P][v.sup.2.sub.SW] + [([A.sub.i]/m).sub.CME]([M.sub.CME]/[[tau].sub.CME])[v.sup.2.sub.CME], [A.sub.2] = 2[[rho].sub.s][A.sub.i]P. (64)

Equation (63) can be simplified to

da/dt = 2p[r.sup.2]/[mu][(1 - [e.sup.2]).sup.2] x [[A.sub.0] e sin f + [A.sub.2] (e sin f cos [alpha] + p/r sin [alpha] cos [delta])[cos.sup.2][alpha]]. (65)

Differentiating (65) with respect to a and equating the result to zero, we obtain

[partial derivative]/[partial derivative][alpha] (da/dt) = -3e sin f [cos.sup.2] [alpha] sin [alpha] + p/r cos [alpha] cos [delta] -3 p/r [sin.sup.2] [alpha] cos [alpha] cos [delta]= 0. (66)

Dividing on 3[cos.sup.2][alpha] sin [alpha] yields

[partial derivative]/[partial derivative][alpha] (da/dt) = -e sin f + 1/3 p/r sin [alpha] cos [alpha] cos [delta] - 3 p/r tan [alpha] cos [delta] = 0, (67)

which can be written as

(1/3 cos [alpha] sin [alpha]){-3e sin f cos [alpha] sin [alpha] + p/r cos [delta] -3 p/r tan [alpha] cos [alpha] sin [alpha] cos [delta]} = 0. (68)

Equating the numerator to zero and rearranging the terms yield

2(1 + e cos f/e sin f) cos [delta] [tan.sup.2][alpha] + 3 tan [alpha] - 1 + e cos f/e sin f cos[delta] = 0, (69)

which has the solution of quadratic equations

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (70)

This equation represents the control law that maximizes the rate of change of semimajor axis. It computes the actual sail pitch/cone angle profile necessary to increase the semimajor axis at a maximum rate (Figure 3).

7.4. Optimal Control Law for the Rate of Change of Eccentricity. when comparing (70) and (69), we observe that the required direction to maximize certain variable is given by tan [??] = (l + e cos f)/e sin f = (p/r)/e sin f, which represents the transverse component divided by the radial component of the primer vector. Thus, in case of eccentricity, set

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (71)

The maximum rate of change of the orbital eccentricity follows directly, taking into account the simplification addressed in special Case II (Figure 4), as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (72)

7.5. Optimal Control Law for the Rate of Change of Inclination. Changes in the out-of-plane orbital elements, such as inclination and right ascension of the ascending node, can be effected by the use of simple switching functions. For maximum rate of change of inclination, the solar sail thrust can be directed alternately above and below the orbit plane every half orbit by the sign function which has +1 or -1.

In what follows, we will find an optimum cone angle that maximizes the rate of change of inclination:

[partial derivative]/[partial derivative][alpha] (di/df) = [r.sup.3]/[mu]p sign [cos(f + [omega])] [partial derivative]/[partial derivative][alpha] W = 0, (73)

from which one can obtain

3 [[rho].sub.s]/[[rho].sub.d] [1 - 2[tan.sup.2][alpha]] = 1/cos [alpha] [[tan.sup.2][alpha] - 1]. (74)

Setting x = [tan.sup.2][alpha], the above equation can be written as

[x.sup.3] + (-1 - [[rho].sub.d]/3[[rho].sub.s])[x.sup.2] + (2[[rho].sub.d] - 3[[rho].sub.s]/3[[rho].sub.s]) x + (1 - [[rho].sub.d]/3[[rho].sub.s]) = 0, (75)

which has the solution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (76)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (77)

As a special case when ignoring the diffusive part of radiation, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (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:

[[DELTA].sub.i] = 88.2[beta](degrees per orbit). (79)

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[degrees], which enables us to maximize angular momentum.

7.6. Optimal Control Law for the Rate of Change of Ascending Node. The procedure is the same as in optimal control law for inclination

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (80)

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:

[DELTA][OMEGA] = (88.2[beta]/sin i)(degrees per orbit). (81)

However, closer orbits to the Sun have shorter orbit periods, and so the time to achieve an overall inclination change is shorter (Figure 5).

7.7. Optimal Control Law for the Rate of Change of Argument of Perihelion. when comparing (70) and (69), we observe that the required direction to maximize a certain variable is given by tan [??] =(1 + e cos f)/e sin f = (p/r)/e sin f, which represents the transverse component divided by the radial component of the primer vector. Thus, in case of argument of perihelion (Figure 6), we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (82)

8. 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.

http://dx.doi.org/ 10.1155/2013/353056

Acknowledgments

The author is deeply indebted and thankful to the deanship of the scientific research and to his helpful and distinct team of employees at Taibah University, Al-Madinah Al-Munawwarah, Saudia Arabia. This research work was supported by a Grant no. 3018/1434. The author would also wish to record here his indebtedness and thankfulness to the reviewers for their valuable and fruitful comments as well as for their potential reading and suggestions.

References

[1] K. E. Tsiolkovsky, Extension of Man Into Outer Space, United Scientific and Technical, Gaithersburg, Md, USA, 1921, Proceedings of the Symposium on Jet Propulsion, vol. 2, 1936.

[2] K. Tsander, "From a scientific heritage," NASA Technical Translation TTF-541, 1967.

[3] F. Rizvi, Solar sail attitude dynamics and coning control: on developing control methods for solar sail coning at orbit rate to attain desired orbital effects [M.S. thesis], Faculty of the Graduate School of the University of Colorado, Department of Aerospace Engineering Sciences, Boulder, Colo, USA, 2010.

[4] J. L. Wright, Space Sailing, Gordon and Breach Science, New York, NY, USA, 2nd edition, 1993.

[5] T. C. Tsu, "Interplanetary travel by solarsail," American Rocket Society Journal, vol. 29, pp. 422-427, 1959.

[6] H. S. London, "Some exact solutions of the equations of motion of a solar sail with constant sail setting," American Rocket Society Journal, vol. 30, pp. 198-200, 1960.

[7] A. N. Zhukov and V. N. Lebedev, "Variational problem of transfer between heliocentric orbits by means of a solar sail," Cosmic Research, vol. 2, pp. 45-50, 1964.

[8] T. S. Jayaraman, "Time optimal orbit transfer trajectory for solar sail spacecraft," Journal of Guidance and Control, vol. 3, no. 6, pp. 536-542, 1980.

[9] L. J. Wood, T. P. Bauer, and K. P. Zondervan, "Comment on timeoptimal orbit transfer trajectory for solar sail spacecraft," Journal of Guidance, Control, and Dynamics, vol. 5, no. 2, pp. 221-224, 1982.

[10] R. B. Powers, V. L. Coverstone, and J. E. Prussing, "Solar sail optimal orbit transfers to synchronous orbits," in Proceedings of the AIAA/AAS Astrodynamics Specialist Conference, Girdwood, Alaska, USA, August 1999.

[11] R. B. Powers and V. L. Coverstone, "Optimal solar sail orbit transfers to synchronous orbits," The Journal of the Astronautical Sciences, vol.49, no.2, pp. 269-281, 2001, Proceedings of the AIAA/AAS Astrodynamics Specialist Conference, Girdwood, Alaska, vol. 103 of Advances in the Astronautical Sciences, pages 523-538, 1999.

[12] C. G. Sauer, "Optimum solar sail interplanetary trajectories," in Proceedigs of the AIAA/AAS Astrodynamics Conference, vol. 792, p. 76, San Diego, Calif, USA, 1976.

[13] G. W. Hughes and C. R. McInnes, "Solar sail hybrid trajectory optimization," in Proceedings of the AAS/AIAA Astrodynamics Specialist Conference of Advances in the Astronautical Sciences, vol. 109, pp. 2369-2380, Quebec City, Canada, 2001.

[14] B. Dachwald, "Solar sail performance requirements for missions to the outer solar system and beyond," in Proceedings of the 55th International Astronautical Congress, IAC-04-S.P.11, pp. 7345-7353, October 2004.

[15] B. Wie, "Solar sail attitude control and dynamics, part 1," Journal of Guidance, Control, and Dynamics, vol.27, no.4, pp.526-535, 2004.

[16] T. J. Waters and C. R. McInnes, "Solar sail dynamics in the three-body problem: homoclinic paths of points and orbits," International Journal of Non-Linear Mechanics, vol. 43, no. 6, pp. 490-496, 2008.

[17] S.-P. Gong, Y.-F. Gao, and J.-F. Li, "Solar sail time-optimal interplanetary transfer trajectory design," Research in Astronomy and Astrophysics, vol. 11, no. 8, pp. 981-996, 2011.

[18] A. R. Breen, B. J. Thompson, M. Kojima et al., "Measurements of the solar wind over a wide range of heliocentric distances--a comparison of results from the first three whole sun months," Journal of Atmospheric and Solar-Terrestrial Physics, vol. 62, no. 16, pp. 1527-1543, 2000.

[19] F. A. Abd El-Salam, "Perturbative effects on a Mercurian orbiter due to the solar radiation pressure, solar wind and the coronal mass ejections," New Astronomy, vol. 12, no. 6, pp. 490-496, 2007.

F. A. Abd El-Salam (1, 2)

(1) Department of Mathematics, Faculty of Science, Taibah University, Al-Madinah 433, Saudi Arabia

(2) Department of Astronomy, Faculty of Science, Cairo University, Cairo 12613, Egypt

Correspondence should be addressed to F. A. Abd El-Salam; f.a.abdelsalam@gmail.com

Received 25 April 2013; Accepted 22 August 2013

Academic Editor: Song Cen

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Salam, F.A. Abd El- |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2013 |

Words: | 6434 |

Previous Article: | Asian option pricing with monotonous transaction costs under fractional Brownian motion. |

Next Article: | Impulsive state feedback control of cheese whey fermentation for single-cell protein production. |

Topics: |