#### Abstract

This paper is dedicated to the development of a model of the attitude dynamics for a nonideal Simple Solar Photon Thruster (SSPT) and to the analysis of sailcraft motions with respect to their centre of mass. Derivation of the expressions for force and torque due to solar radiation that is valid for the case, when there is a misalignment of the SSPT axis with the sun direction, is followed by study of sailcraft dynamics and stability properties. Analysis of stability shows that an ideally reflecting sail is unstable, while for a sailcraft with nonideal collector, the symmetry axis is stable with respect to the Sun direction for large variety of system parameters. The motion around symmetry axis is always unstable and requires an active stabilizer.

#### 1. Introduction

Solar sails are attractive spacecrafts for long missions, for they create a propulsion effort using freely available solar radiation pressure. The initial concepts of such vehicle have been proposed in the 1920s [1, 2] as an alternative for minimizing on-board energy consumption. With the recent development of new films, materials, and structures its implementation has become possible making solar sails a subject of detailed research regarding eventual missions, the respective control issues, the sailcraft structure, lab testing, and so forth (see, e.g., [3–12], just to name a few), and leading to the implementation of the first sailcraft mission by JAXA in 2010 [10].

The large majority of these studies are focused on the analysis of the Flat Solar Sail (FSS); its principal element is a large flat surface formed by a reflecting membrane. Control of FSS orbital motion is performed by turning the entire sail with respect to the sunlight direction. As a result, the radiation pressure on the surface changes, enabling the control of the vehicle trajectory. Attitude dynamics of such sailcraft are studied, for example, in [4, 11, 12].

Another possible scheme for a sailcraft is the Solar Photon Thruster (SPT). Its main idea is to separate the functions of collecting and reflecting the incoming light. Even though the SSPT's general scheme has been proposed long ago [1, 2, 13], the detailed research has begun quite recently [14–22]. SPT consists of a parabolic reflecting surface (collector) that is always pointing towards the sun, taking advantage of the totality of the sail area. The incoming light is then reflected on a system of smaller mirrors capable of controling the direction of the outgoing light; the sailcraft is, therefore, controlled more easily by the rotation of one smaller mirror.

The SPT can be either Simple SPT (SSPT) or Dual Reflector SPT (DR SPT) depending on the system of smaller mirrors used. The SSPT is provided with two main elements, the Collector C and the Director mirror D [1, 13–15, 18, 19]. Ideally, the symmetry axis of the collector C remains aligned to the sun-sailcraft direction, and all incoming light is concentrated on the collector focus. Director D is placed at this focus; it can rotate changing the direction of the outgoing light in order to control the force acting on the sailcraft and therefore its motion.

Several versions of DR SPT structure (see, e.g., [1, 13, 19–21]) include one more parabolic mirror, called Reflector R. The incoming light from C hits R, becoming a parallel beam that is reflected on D. The director can rotate controlling the resulting propulsion.

Previous studies [17, 19] show that SSPT has many advantages compared to DR SPT. Although its heat tolerance requirements are higher than DR SPT (all the incoming light is concentrated at one point of the director), the advantages are more numerous. When compared to DR SPT, SSPT is much simpler, being therefore easier to deploy and maintain during the flight, and more compact and lighter, thus saving space for the payload. SSPT is more robust than DR SPT regarding the deformations of the collector shape and misalignments with respect to sun-sailcraft direction. For the above reasons, SSPT has been chosen to be the subject of more detailed studies.

In [22], the first model of force acting on SSPT which is assumed to be perfectly aligned with the sun-sailcraft direction is presented and some trajectory simulations are performed. Here, elaboration of a force and torque model for SSPT attending to the eventual axis misalignment permits one to study the sailcraft attitude dynamics and stability properties.

#### 2. Force and Torque Model for a Nonideal Sailcraft

Figure 1 shows the general scheme of SSPT. The Collector C is a big parabolic reflecting surface which concentrates the parallel beam of the incoming sunlight at its focus . Director D is a relatively small plane mirror located at the focus of the collector. The force acting on the sail is controlled by rotating D which changes the direction of outgoing light varying the magnitude and direction of the total propulsion force.

In [22], the general characteristics of a compound solar sail have been analysed. The results show that for minimizing collector film overheating in order to lower the thermal requirements, the light beam should not reflect on the collector more than once. To avoid secondary reflections, the angle between the collector’s symmetry axis and the normal to the director’s surface must satisfy the following condition: Here, and are the focal distance and the effective radius of the collector, respectively. If the control angle is limited by condition (1), the components of the force and torque due to solar radiation pressure can be obtained analytically.

A force and torque model for SSPT has been developed in [22] assuming that all sailcraft elements possess ideally reflecting surfaces. But as the collector is a large surface, it should be made of a thin film and it cannot be an ideal reflector. The director is much smaller, so it can have better optical properties. So, here a model for force and torque considering a nonideal collector is developed. The director is supposed to be ideal and small enough so that the effect of its shadow can be disregarded. Consider a SSPT whose symmetry axis is misaligned in relation to sunlight direction.

To calculate the vectors of force and torque due to solar radiation on the SSPT, consider a parabolic surface (collector) shown in Figure 2. Its equation in the reference frame , in which the center of mass coincides with the focus of the paraboloid and is its symmetry axis, is given by or

Suppose that the incoming sunlight has the direction (Figure 2). It is absorbed by the collector at the point of the element of the parabolic surface with the position vector , producing the force where is the intensity of the light flow at the current point of the orbit. The respective torque about the center of the mass of the SSPT is

A fraction of the incident light is reflected on the element of the surface , and a fraction of the reflected light is reflected specularly, while a fraction is scattered from the collector surface in a nonspecular way [1]. After hitting the collector, specularly reflected light has direction as follows: where is the internal normal unit vector to the surface at point given by

In (7), , , and are the coordinate orts. This fraction of light that is reflected specularly generates a force while the nonspecularly reflected light generates a force given by Moreover, a portion of the absorbed light is reemitted backwards creating the force In (10), and are the surface front and back non-Lambertian coefficients and and are the front and back emissivities. The total force applied is described as It produces the torque with respect to the center of mass given by

Then, the light is reflected by the director whose orientation can be described by the normal to its surface . One can assume that the light hitting D is only the fraction specularly reflected on the collector: as the director is small compared to the collector, the respective spatial angle observed from the collector is small and the influence of the scattered light may be disregarded. The incoming light from point hits the director at point with the position vector that can be found from condition (13) as follows: so,

Force and torque transmitted by the light hitting the director are, respectively, as follows:

The reflection of the light at the director produces force and torque that can be written, respectively, as where is the direction of the light reflected on the control mirror D.

The interaction of the light beam caught by the element of the collector C with the sail structure generates the following elementary force:

The total force actuating on the sailcraft due to sunlight radiation pressure can be obtained by integration as follows:

After integration, one gets

Here, with

For the case when the ideally reflecting SSPT is completely aligned with the sun direction ( and ) and the director rotates in the plane of the orbit (, , and ), expressions (20) are reduced to

These expressions for the solar radiation force acting on SSPT have been used previously to study orbital transfer problems for SSPT and FSS [16, 22]. These studies show several advantages of SSPT scheme compared to FSS that are more evident for large-scale missions (i.e., for Solar System escape).

Calculation of the total torque due to the light pressure results in the following expressions: where is defined in (22) and

These models for force and torque are used later on to study the attitude dynamics of SSPT.

#### 3. SSPT Attitude Dynamics

Consider the problem of a SSPT attitude dynamics on a heliocentric orbit. This analysis is done using two right Cartesian reference frames located at the SSPT center of mass . is the reference frame connected with the sailcraft; its axes are SSPT central principle axes of inertia. is the frame related to the orbital motion of the SSPT; its unit vectors along the axes are where is the radius vector of point with respect to the center of the sun and is SSPT’s center of mass velocity. Obviously, the sunlight follows in the direction of axis.

Mutual orientation of these frames can be described by the orthogonal matrix , where the elements are the direction cosines of the axes of the system with respect to . The transition from system to system can be done by three Euler’s rotations about axes 1, 2, and 3 through the angles (roll), (pitch), and (yaw). The direction cosines are expressed through these angles as

At the heliocentric trajectories, the characteristic values of the angular acceleration due to orbital motion and gravity-gradient torque are much smaller than those created by cm-cp offset. Taking into account only the latter effect, the equations of attitude motion can be written as where , , and are the principal moments of inertia of the sailcraft and , , and are the components of the solar radiation torque given by (24). The projections of the sailcraft absolute angular velocity onto the axes , , , and can be expressed as

The direction of sunlight is . Equations of motion (28) and (29) admit the following trivial solution:

Solution (30) corresponds to the nominal equilibrium orientation. According to the SSPT force and torque model, the solar radiation torque acting on an SPT in the vicinity of this equilibrium has the components where

Assuming that the attitude motion of SSPT occurs in the vicinity of the equilibrium orientation, (28) and (29) can be linearized, resulting in

Analyzing linearized system (33), one can notice that the motion of the axis does not depend on the yaw angle . The rotation of the sailcraft about the collector symmetry axis is unstable. Considering that the sailcraft is an axisymmetric body, that is, , the stability conditions for motion of the axis can be written using the characteristic equation

For stability of zero equilibrium position of system (33), characteristic equation (34) should have purely imaginary roots, which happens if and only if

As one can see, , and so conditions (35) are equivalent to

axis should be stable for all possible orientations of the director. Considering that director rotation is limited by condition (1), the components of the normal to the director surface and must be limited by

This condition is satisfied for all possible values of and if and only if

Writing condition (38) in terms of the physical characteristics of the collector film, one gets where .

##### 3.1. Stability of a Solar Sail with Ideal Collector

For a sailcraft with an ideal collector, the surface optical parameters are [1] as follows: In this case, condition (39) becomes

So, for all geometries of the sail, the axis is unstable.

##### 3.2. Stability of a Solar Sail with Black Collector

For a sailcraft with an ideally black film in collector, the surface optical parameters are [1] In this case, condition (39) becomes

Condition (43) means that for an ideal black collector the motion of the axis is stable for . So, for black sails with , SSPT axis is stable.

##### 3.3. Stability of a General Solar Sail

For sailcraft with general properties of the film, condition (39) holds for some . Considering the case described in [1], the film properties are

Condition (39) for this case means that the sailcraft axis is stable if .

The stability of axis does not mean that the sailcraft is stable. Its motion about axis is always unstable, being a combination of harmonic oscillation and rotation with constant angular velocity. Meanwhile, the orientation of a SSPT about axis should be stabilized since it determines the director position, so a stabilizer for is always necessary.

#### 4. Conclusion

In this paper, a new model for force and torque actuating on a compound solar sail is developed and used to study sailcraft motions. The attitude dynamics of a SSPT is studied for the case when the collector film possesses nonideal optical properties and the sailcraft symmetry axis is not aligned to the sun-sailcraft direction. Analysing the sailcraft stability, one can see that the solar radiation pressure destabilizes an ideal sail, but non-ideal collector film contributes to stabilization, resulting in stability, of its symmetry axis for all orientations of director when collector focal distance is large enough. However, an active control system is always necessary to guarantee stabilization of the sailcraft orientation with respect to the collector’s symmetry axis.

#### Conflict of Interests

The authors declare that they have no conflict of interests.

#### Acknowledgments

This research is supported by the Brazilian program “Science without Borders” and Project ODySSea (PTDC/CTESPA/098030/2008) funded by the Portuguese Foundation for Science and Technologies (FCT), the Portuguese Operational Program for Competitiveness Factors (COMPETE), the Portuguese Strategic Reference Framework (QREN), and the European Regional Development Fund (FEDER).