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 F. 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.
Scheme of a Simple Solar Photon Thruster.
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:
(1)tanα~>4fcRc(4fc2-Rc2).
Here, fc and Rc 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 Oxyz, in which the center of mass O coincides with the focus F of the paraboloid and Oz is its symmetry axis, is given by
(2)x2+y2+4fc(z-fc)=0,
or
(3)z=14fc(4fc-x2-y2)=0.
Scheme of reflections on SSPT with axis misalignment.
Suppose that the incoming sunlight has the direction σ→ (Figure 2). It is absorbed by the collector at the point Ac of the element of the parabolic surface dS with the position vector ξ→=xi→+yj→+zk→, producing the force
(4)dP→1=-ρσ→(σ→·n→)dS,
where ρ=ΦrE2/rs2 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 O is
(5)dM→1=ξ→×dP→1.
A fraction r~ of the incident light is reflected on the element of the surface dS, and a fraction s of the reflected light is reflected specularly, while a fraction (1-s) is scattered from the collector surface in a nonspecular way [1]. After hitting the collector, specularly reflected light has direction σ→1 as follows:
(6)σ→1=σ→-2(σ→·n→)n→,
where n→ is the internal normal unit vector to the surface at point ξ→ given by
(7)n→=-xi→+yj→+2fck→x2+y2+4fc2.
In (7), i→, j→, and k→ are the Oxyz coordinate orts. This fraction of light that is reflected specularly generates a force
(8)dP→2r=ρr~sσ→1(σ→,n→)dS,
while the nonspecularly reflected light generates a force given by
(9)dP→2s=ρBf(1-s)r~n→(σ→,n→)dS.
Moreover, a portion of the absorbed light is reemitted backwards creating the force
(10)dP→2e=ρ(1-r~)εfBf-εbBbεf+εbn→(σ→,n→)dS.
In (10), Bf and Bb are the surface front and back non-Lambertian coefficients and εf and εb are the front and back emissivities. The total force applied is described as
(11)dP→2=dP→2r+dP→2s+dP→2e.
It produces the torque with respect to the center of mass O given by
(12)dM→2=ξ→×dP→2.
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 Ac hits the director at point Ad with the position vector that can be found from condition (13) as follows:
(13)0=(ξ→Ad·ν→)=(ξ→+τσ→1·ν→),
so,
(14)τ=-(ξ→·ν→)(σ→1·ν→),ξ→Ad=ξ→-(ξ→·ν→)(σ→1·ν→)σ→1.
Force and torque transmitted by the light hitting the director are, respectively, as follows:
(15)dP→3=-ρr~sσ→1(σ→·n→)dS,dM→3=ξ→Ad×dP→3=ρr~s(σ→1×ξ→)(σ→·n→)dS.
The reflection of the light at the director produces force and torque that can be written, respectively, as
(16)dP→4=ρr~sσ→2(σ→·n→)dS,dM→4=ξ→Ad×dP→4,
where
(17)σ→2=σ→1-2(σ→1·ν→)ν→
is the direction of the light reflected on the control mirror D.
The interaction of the light beam caught by the element dS of the collector C with the sail structure generates the following elementary force:
(18)dP→=dP→1+dP→2+dP→3+dP→4.
The total force actuating on the sailcraft due to sunlight radiation pressure can be obtained by integration as follows:
(19)P→=∫CdP→.
After integration, one gets
(20)Px=πρ[χx+ζxln(1+Rc24fc2)],Py=πρ[χy+ζyln(1+Rc24fc2)],Pz=πρ[χz+ζzln(1+Rc24fc2)].
Here,
(21)χx=16fc{+2νx(νz+νyσyσz-2νzσz2))]σxλΨ-6fcRc2×[+2νx(νz+νyσyσz-2νzσz2))-σxσz+r~s(-σxσz(1-2νx2)+2νx(νz+νyσyσz-2νzσz2))]},χy=16fc{+2νy(νz+νxσxσz-2νzσz2))]σyλΨ-6fcRc2×[+2νy(νz+νxσxσz-2νzσz2))-σyσz+r~s(-σyσz(1-2νy2)+2νy(νz+νxσxσz-2νzσz2))]},χz=4fcσz(-2fc+Rc2+4fc2)Ψ+Rc2[-2νz(νz+νxσxσz+νyσyσz))σz2+r~s(1-σz2(1-2νz2)-2νz(νz+νxσxσz+νyσyσz))],
with
(22)λ=16fc3+(Rc2-8fc2)Rc2+4fc2,Ψ=Bbεb(r~-1)+Bf[εf+r~εb-r~s(εf+εb)]εf+εb,ζx=8fc2r~s[-σxσz(1-2νx2)+νx(νz+2νyσyσz-3νzσz2)],ζy=8fc2r~s[-σyσz(1-2νy2)+νy(νz+2νxσxσz-3νzσz2)],ζz=4fc2r~s[1-3σz2(1-2νz2)-2νz(νz+2νxσxσz+2νyσyσz)].
For the case when the ideally reflecting SSPT is completely aligned with the sun direction (σx=σy=0 and σz=1) and the director rotates in the plane of the orbit (νx=sinα~, νy=0, and νz=cosα~), expressions (20) are reduced to
(23)PX=ΦrE2rs2πRc2[1-8fc2Rc2ln(1+Rc24fc2)]sin2α~,PY=0,PZ=2ΦrE2rs2πRc2[cos2α~-4fc2Rc2ln(1+Rc24fc2)cos2α~].
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:
(24)Msx=πρ{σyυΨ120fc2+Rc2σz4fc[(Rc2-4fc2)σy-4fc2r~s(2νxνyσx+(1-2νx2)σy)]σyυΨ120fc2},Msy=πρ{-σxυΨ120fc2+Rc2σz4fc[-(Rc2-4fc2)σx+4fc2r~s(2νxνyσy+(1-2νy2)σx)]σxυΨ120fc2},Msz=2πρfcRc2r~sνzσz(νxσy-νyσx),
where Ψ is defined in (22) and
(25)υ=64fc5+(3Rc4+4fc2Rc2-32fc4)Rc2+4fc2.
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 O. Oxyz is the reference frame connected with the sailcraft; its axes are SSPT central principle axes of inertia. OXYZ is the frame related to the orbital motion of the SSPT; its unit vectors along the axes are
(26)e→Z=r→OrO,e→Y=e→Z×v→OvO,e→X=e→Y×e→Z,
where r→O is the radius vector of point O with respect to the center of the sun and v→O is SSPT’s center of mass velocity. Obviously, the sunlight follows in the direction of OZ axis.
Mutual orientation of these frames can be described by the orthogonal matrix A=[aij], where the elements aij are the direction cosines of the axes of the Oxyz system with respect to OXYZ. The transition from system OX1X2X3 to system Ox1x2x3 can be done by three Euler’s rotations about axes 1, 2, and 3 through the angles γ (roll), α (pitch), and β (yaw). The direction cosines aij are expressed through these angles as
(27)a11=cos(x1,X1)=cosαcosβ,a12=cos(x2,X1)=-cosαsinβ,a13=cos(x3,X1)=sinα,a21=cos(x1,X2)=sinβcosγ+sinαcosβsinγ,a22=cos(x2,X2)=cosβcosγ-sinαsinβsinγ,a23=cos(x3,X2)=-cosαsinγ,a31=cos(x1,X3)=sinβsinγ-sinαcosβcosγ,a32=cos(x2,X3)=cosβsinγ+sinαsinβcosγ,a33=cos(x3,X3)=cosαcosγ.
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
(28)Ap˙+(C-B)qr=Msx,Bq˙+(A-C)rp=Msy,Cr˙+(B-A)pq=Msz,
where A, B, and C are the principal moments of inertia of the sailcraft and Msx, Msy, and Msz are the components of the solar radiation torque given by (24). The projections of the sailcraft absolute angular velocity onto the axes Oxyz, p, q, and r can be expressed as
(29)p=γ˙a11+α˙sinβ,q=γ˙a12+α˙cosβ,r=γ˙a13+β˙.
The direction of sunlight is σ→=(a31,a32,a33). Equations of motion (28) and (29) admit the following trivial solution:
(30)α=β=γ=0.
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
(31)Msx=Gνxνyα+(F+νx2G)γ,Msy=(F+νy2G)α+Gνxνyγ,Msz=Gνz(νyα+νxγ),
where
(32)F=πρ[υΨ120fc2+Rc24fc(Rc2-4fc2)-Rc2fcr~s],G=2πρRc2fcr~s.
Assuming that the attitude motion of SSPT occurs in the vicinity of the equilibrium orientation, (28) and (29) can be linearized, resulting in
(33)Bα¨-(F+νy2G)α-Gνxνyγ=0,Aγ¨-Gνxνyα-(F+νx2G)γ=0,Cβ¨=Gνz(νyα+νxγ).
Analyzing linearized system (33), one can notice that the motion of the Oz 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, A=B, the stability conditions for motion of the Oz axis can be written using the characteristic equation
(34)A2λ4-A[2F+G(νx2+νy2)]λ2+F[F+G(νx2+νy2)]=0.
For stability of zero equilibrium position of system (33), characteristic equation (34) should have purely imaginary roots, which happens if and only if
(35)2F+G(νx2+νy2)<0,F[F+G(νx2+νy2)]>0,[2F+G(νx2+νy2)]2-4F[F+G(νx2+νy2)]>0.
As one can see, G>0, and so conditions (35) are equivalent to
(36)F<-G(νx2+νy2).
Oz 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 νx and νy must be limited by
(37)16fc2Rc2(4fc2+Rc2)2<νx2+νy2≤1.
This condition is satisfied for all possible values of νx and νy if and only if
(38)F<-G.
Writing condition (38) in terms of the physical characteristics of the collector film, one gets
(39)15ψ(r~s+ψ-1)+1εb+εf[Bbεb(r~-1)+Bf(εf+r~εb-sr~fεb-sr~fεf)]×[2+(3ψ-2)(1+ψ)3/2]<0,
where ψ=(Rc/2fc)2.
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:
(40)r~=1,s=1,εb=εf,Bb=Bf=23.
In this case, condition (39) becomes
(41)4fc3πρψ2<0.
So, for all geometries of the sail, the Oz 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]
(42)r~=0,s=0,Bb=Bf=0.
In this case, condition (39) becomes
(43)4fc3πρψ(ψ-1)<0.
Condition (43) means that for an ideal black collector the motion of the Oz axis is stable for ψ<1. So, for black sails with Rc<2fc, 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 0<ψ<ψ0. Considering the case described in [1], the film properties are
(44)r~=0.88,s=0.94,εb=0.55,εf=0.05,Bb=0.55,Bf=0.79.
Condition (39) for this case means that the sailcraft Oz axis is stable if Rc<0.83fc.
The stability of Oz axis does not mean that the sailcraft is stable. Its motion about Oz axis is always unstable, being a combination of harmonic oscillation and rotation with constant angular velocity. Meanwhile, the orientation of a SSPT about Oz 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).
McInnesC. R.1999Berlin, GermanySpringer10.1007/978-1-4471-3992-8PolyahovaE. N.1986Moscow, RussiaMacDonaldM.McInnesC.Solar sail science mission applications and advancement20114811170217162-s2.0-8005407383710.1016/j.asr.2011.03.018WieB.Solar sail attitude control and dynamics20042745365442-s2.0-384315266710.2514/1.11133SleightD. W.MuheimD. M.Parametric studies of square solar sails using finite element analysisProceedings of the 45th AIAA/ASME/ASCE/AHS Structures, Structural Dynamics and Materials ConferenceApril 2004Palm Springs, California85972-s2.0-16244410125JinjinZ. G.ZhaoyuS.HongjunZ.WenhaoH.LianguanS.Structural analysis about a new solar sail20092249256LiQ.MaX.WangT.Reduced model for flexible solar sail dynamics20114834464532-s2.0-7995804584610.2514/1.48789CouceiroF.GamboaP.SilvaJ. M.GuermanA. D.Membrane-based space energy collector: a conceptual study20138781310.1016/j.actaastro.2013.01.012TaleghaniB. K.LivelyP. S.GasparJ. L.MurphyD. M.TrauttT. A.Dynamic and static shape test/analysis correlation of a 10 meter quadrant solar sailProceedings of the 46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials ConferenceApril 2005Austin, Tex, USA400340162-s2.0-28844496351TsudaY.MoriO.FunaseR.SawadaH.YamamotoT.SaikiT.EndoT.YonekuraK.HoshinoH.KawaguchiJ.Achievement of IKAROS—Japanese deep space solar sail demonstration mission20138221831882-s2.0-8486059149710.1016/j.actaastro.2012.03.032JinZ.TianshuW.Coupled attitude-orbit control of flexible solar sail for displaced solar orbit201350367568510.2514/1.A32369TsudaY.SaikiT.FunaseR.Generalized attitude model for spinning solar sail spacecraft20133696797410.2514/1.59516ForwardR. L.Solar photon thrustor19902744114162-s2.0-002546342010.2514/3.26158McInnesC. R.Payload mass fractions for minimum-time trajectories of flat and compound solar sails2000236107610782-s2.0-003431604510.2514/2.4651MengaliG.QuartaA. A.Time-optimal three-dimensional trajectories for solar photon thruster spacecraft20054223793812-s2.0-1814440039810.2514/1.9417GuermanA. D.SmirnovG.Attitude control of solar photon thruster vs plane solar sail2006123263522648GuermanA.SmirnovG.Attitude dynamics of a compound solar sailProceedings of the AIAA 57th International Astronautical Congress (IAC '06)October 2006Valencia, Spain409041002-s2.0-40749094050GuermanA. D.SmirnovG.Compound solar sail with optical properties: models and performance20074437327342-s2.0-3454729506110.2514/1.26121GuermanA. D.SmirnovG. V.PereiraM. C.Comparison of two compound solar sail schemes1281Proceedings of the International Conference on Numerical Analysis and Applied Mathematics (ICNAAM '10)September 20104854882-s2.0-7995446445310.1063/1.3498515DachwaldB.WurmP.Design concept and modeling of an advanced solar photon thruster2009134723740DachwaldB.WurmP.Mission analysis and performance comparison for an advanced solar photon thruster20114811185818682-s2.0-8005406691910.1016/j.asr.2011.01.030GuermanA. D.SmirnovG. V.PereiraM. C.Orbital dynamics of a simple solar photon thruster20092009112-s2.0-7764921982410.1155/2009/537256537256ZBL1188.70049