Received 31 March 2016; accepted 27 May 2016; published 31 May 2016

ABSTRACT

We investigate solar sail in the circular restricted three-body problem, where the larger primary is a source of radiation and the smaller primary is an oblate spheroid in the system. Firstly, the differential equations of motion for solar sail in the system combined effects of radiation and oblateness of celestial bodies are built. Then the positions of the solar sail collinear Lagrange points are calculated as mass ratio or oblateness changes in certain extent. Linearization near the collinear equilibria of the system is applied. A linear quadratic regulator is used to stabilize the nonlinear system. Three different cases of solar sail equilibrium orbits are studied each with different choices for the weight matrices. The simulations reveal that solar sail equilibrium orbits can be stable under active control by changing three angles, incident angle, cone angle and clock angle of the solar sail.

Keywords:

Solar Sail, Restricted Three-Body Problem, Equilibrium Orbits, Oblateness

Subject Areas: Aerospace Engineering

1. Introduction

The solar sail has become a hot topic of research during the last few decades [1] - [5] . It uses sunlight to generate propulsion in space by reflecting solar photon flux from a large, mirror-like sail made of lightweight, highly reflective polyimide film material. Especially after the successful flight of the Interplanetary Kite-craft Accelerated by Radiation of the Sun (IKAROS) in 2010 [6] , a number of new space missions using a solar sail have been proposed, such as NanoSail-D [7] , Cubesail [8] , DeorbitSail [9] and Sunjammer [10] . The reason why the solar sail has increasingly attracted the world’s attention is that, when compared with conventional spacecraft, it can accelerate continually without a propellant. It is this very characteristic that makes the solar sail to execute a few novel space missions that conventional spacecraft could not execute very well; for instance: Solar Polar Imager, Heliostorm, Interstellar Probe and Near Earth Asteroid Rendezvous [11] .

The restricted three-body problem (RTBP) has been considered a basic dynamic model ever since the scientists have studied the solar sail [12] [13] . The RTBP describes the motion of an infinitesimal mass, usually a satellite, solar sail or asteroid, moving under the gravitational effect of two finite masses, normally called primaries, which move in circular orbits about their center of mass under their mutual attraction [5] . Usually we do not take into account of the effect of the infinitesimal mass of the motion of the primaries. However, an absolute spherical celestial body is very rare in space, most of the planets are oblate. Sharma and Rao [14] investigated the locations of the five equilibrium points by considering the effect of oblateness of the more massive primary. Sharma [15] [16] discussed the existence of periodic orbits in the RTBP when the more massive primary was an oblate spheroid. Douskos [17] [18] focused on the equilibrium points and their stability in the Hill’s problem with oblateness. Singh [19] analyzed the combined effects of perturbations, oblateness, and radiation of the primaries on the nonlinear stability of the Lagrange points. However, the above literatures did not take account of the impacts of oblateness on solar sail orbits in RTBP. Solar sail as one of new important deep space probes must be considered the influences of oblateness of celestial bodies if it could observe and study well in the long-term space missions. Therefore, the effects of the oblateness of the primaries should be included when we establish the model of solar sail RTBP and investigate the motion of the solar sail in the RTBP. Thus we achieve better simulation results, which can provide an optimal design for the flight path of the solar sail.

This paper is organized as follows. Section 2 introduces the equations of motion of the system in the circular restricted three-body problem with the larger primary, a source of radiation and the smaller primary, an oblate spheroid. In Section 3, we investigate the equilibrium points of the system with the variations of lightness number of solar sail oroblateness of the smaller primary. Then, in Section 4, linearization near the collinear Lagrange points is taken into account, and the linear quadratic regulator (LQR) is developed to stabilize the nonlinear system. The simulation is given in Section 5. The conclusions are discussed in Section 6.

2. Equations of Motion

The restricted three-body problem with the larger primary a source of radiation and the smaller primary an oblate spheroid is investigated. We use a barycentric, rotating and dimensionless coordinate system Oxyz; the origin is at the barycenter of the primaries; the axis x is along the line joining with the primaries; the direction of the orbital angular velocity ω of the smaller primary defines the axis z; and the axis y completes the right-handed triad. We describe the circular restricted three-body problem in Figure 1. For convenience the dimensionless form is often used [20] . The two primaries have masses m₁ and m₂ respectively, the mass of the infinitesimal body, the solar sail, is m₃. The distance between the primaries is, and the gravitational constant is chosen to be unity. When it comes to the RTBP, the unit mass, length, and time of the system are defined as

(1)

Then, in this system the masses of two primaries are, , where μ is the mass ratio for the system,. The distances of two primaries to the barycenter are,.

Considering the oblateness of the smaller primary, the equations of motion of the solar sail in the rotating coordinate system can be written as

(2)

(3)

(4)

where Ω is the pseudo-potential function [18]

(5)

Ω_x, Ω_y, Ω_z are the components of the partial derivative of the pseudo-potential function Ω on each coordinate

axis; and are the distances of the solar sail from the primaries respectively, and is the perturbed mean motion of the primaries; β is the lightness num-

ber of the solar sail. is the oblateness coefficient of the smaller primary, where R_E and R_P are the equatorial and polar radii of the smaller primary, and R is the distance between the two primaries [21] . a_x, a_y, a_z are the projections of the acceleration produced by the solar radiation pressure force on the axis Ox, Oy, Oz. Figure 2 is a flat solar sail model, where is the angle between n (vectors are expressed by black italicized letters in this article) and incident light r₁ from larger primary; φ is the cone angle between the projection of n in xy plane; γ is the clock angle between the projection of n in xy plane and xz plane. The acceleration produced by solar radiation pressure force can be expressed [5]

(6)

3. Collinear Equilibrium Points

The equilibrium points of the system are the solutions of the equations

(7)

(8)

(9)

where

(10)

(11)

. (12)

We suppose that three collinear equilibrium points lie on the axis x, then we use to denote the collinear equilibria. These collinear Lagrange points satisfy Equations (7), (8) and (9)

(13)

(14)

. (15)

According to Equations (13)-(15), we find that collinear equilibrium points will be get when φ and γ equal zero, that is to say, when the surface of solar sail is perpendicular to the radiate light, Equation (13) has effective solution. From Equation (13) we see that the position of collinear equilibrium points varies with the magnitude of the mass ratio μ, oblateness A and lightness number β. In the following, we set μ = 0.001, then the positions of L₁, L₂, L₃ varying with the variations of A and β are shown in Figure 3. It is clear that these collinear equilibrium points are unstable. The effects of oblateness of smaller primary on L₁, L₂ are obvious, but it has little impact on L₃.

4. Linearization near the Collinear Equilibria

To further investigate the characteristics of the solar sail orbit in the circular restricted three-body problem with oblateness, we need to linearize the system because the differential equations are nonlinear. Given that the collinear Lagrange points of the nonlinear system are, we introduce small perturbations such that we define Equations (16), (17), and (18)

(16)

(17)

. (18)

Substitute Equations (16), (17), and (18) into Equations (2), (3) and (4), and assume that the sail acceleration is constant under the small perturbation from the collinear equilibrium point [22] [23] ; then we obtain the variational equations

(19)

(20)

(21)

where

(22)

(23)

(24)

. (25)

Therein is the evaluation of the second order partial derivative of the potential function at the equilibrium points. With this method we can get the linear dynamic model and establish its state-space equation expressed in matrix notation as Equation (26)

(26)

where the six-dimensional state vector is defined, and

. (27)

The LQR controller is developed to stabilize the nonlinear system in the neighborhood of the collinear libration point. We apply a linear feedback control to Equation (26) that minimizes the quadratic cost function

(28)

where the matrices Q and R represent the weights of the state and control, which are symmetric positive semidefinite and free to be chosen. We obtain the gain matrix by solving the algebraic Riccati equation [24]

. (29)

The closed-loop system is then obtained as

. (30)

A necessary and sufficient condition for the collinear equilibrium points to be linearly stable is that the real part of the eigenvalues of the matrix A − BK are all less than or equal to zero [25] .

5. Simulation

In this section, we choose μ = 0.001, A = 0.001, β = 0.05 as basic parameters, and as an example, we set initial conditions as, , , , , , the time of simulation is 15 units time. Three different cases of solar sail equilibrium point orbits are displayed in Figure 4(a), Figure 6(a) and Figure 8(a), and each case has a different choice of weight matrices Q and R, which are elaborated in Table 1.

Figure 4(a) is a spiral orbit of the sail in the vicinity of L₁. L₁ is a sink of this system, small perturbation near the collinear Lagrange point L₁ will be asymptotic stable. In Figure 4(b), the left parts are the time curve graphs of projections of solar radiation pressure acceleration on each axis. The projections on ξ and η change mildly and approach zeroquickly. But axis ζ’sprojection changes as a sine curve; the amplitude gets smaller and smaller. The right parts are the variations of three angles, where is about between 56.3˚ and 90˚, −90˚ ≤ γ, f ≤ 90˚.

Different from Case A, Figure 6(a1) is the Quasiperiodic orbits of the sail, which are also asymptotic stable. But these kinds of orbits will not come close to L₁, just like Figure 6(a2) shown, these quasiperiodic motion will last forever. In Figure 6(b), the left parts, the projections on ξ and η change like those in Case A; the projection on axis ζ changes as a sine curve and the amplitude remains about the same. The right parts are the variations of three angles, where changes in a small scale about from 89.8˚ to 90˚, −90˚ ≤ γ, f ≤ 90˚. Figure 8(a) displays a direct perturbation trajectory of the sail derived from L₁ point. The maximum acceleration of solar radiation pressure is about 0.009906 units acceleration. The minimum is about 56.5˚, the maximum is 90˚. Meanwhile, the position and velocity of the sail in the vicinity of L₁ varying with time are shown in Figure 5, Figure 7

and (k)(l). The minimum and maximum quantitative values of them are displayed in Table 2. In case A the position and velocity swing back and forth and the amplitudes become smaller. In case B cyclical variations are reflected in the position and velocity, which show that the quasiperiodic orbits of the sail are asymptotically stable. However, in case C position and velocity fluctuate dramatically during about the first four units time; after that they change only slightly near the collinear equilibrium point (Figures 4-9).

In Table 1, we have chosen the weighting matrices Q and R and evaluated the gain matrix K and the eigenvalue of matrix A − BK in three different cases. It is quite clear that the real parts of the six eigenvalues of matrix A − BK in all three cases are negative numbers.

6. Conclusion

In this paper, we investigate the solar sail equilibrium orbits in the circular restricted three-body problem with oblateness. We find that oblateness has little impact on the position of L₃. An LQR controller is used to obtain the numerical solution of the components of solar radiation pressure acceleration. We solve Equation (25) to get the changing laws of angles of the sail, which can make the system stable near the collinear Lagrange points. We can choose different weight matrices Q and R to obtain different solar sail orbits for different space mission requirements.

Cite this paper

Ming Song,Xingsuo He,Yehao Yan,Dongsheng He, (2016) Solar Sail Equilibrium Orbits in the Circular Restricted Three-Body Problem with Oblateness. Open Access Library Journal,03,1-10. doi: 10.4236/oalib.1102620

References

NOTES

^*Corresponding author.