Ground Base Laser Torque Applied on LEO Satellites of Various Geometries

International Journal of Aeronautical and Space Sciences.
2012.
Dec,
13(4):
484-490

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/bync/ 3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

- Published : December 30, 2012

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

This paper is devoted to investigate the feasibility of using a medium power ground-based laser to produce a torque on LEO satellites of various shapes. The laser intensity delivered to a satellite is calculated using a simple model of laser propagation in which a standard atmospheric condition and linear atmospheric interaction mechanism is assumed. The laser force is formulated using a geocentric equatorial system in which the Earth is an oblate spheroid. The torque is formulated for a cylindrical satellite, spherical satellites and for satellites of complex shape. The torque algorithm is implemented for some sun synchronous low Earth orbit cubesats. Based on satellites perigee height, the results demonstrate that laser torque affecting on a cubesat has a maximum value in the order of 10
^{-9}
which is comparable with that of solar radiation. However, it has a minimum value in the order of 10
^{-10}
which is comparable with that of gravity gradient. Moreover, the results clarify the dependency of the laser torque on the orbital eccentricity. As the orbit becomes more circular it will experience less torque. So, we can conclude that the ground based laser torque has a significant contribution on the low Earth orbit cubesats. It can be adjusted to obtain the required control torque and it can be used as an active attitude control system for cubesats.
Worst-case expected disturbance torques for a 1-U CubeSat at 700 km [11]
where, u is the satellite altitude,
P
is the laser power,
θ
is the laser divergence, ϕ is the elevation angle,
is the molecule scattering coefficient at sea level,
is the aerosols scattering coefficient and is the sea level altitude.
where,
where
is a unit vector directed through the force direction, c is the speed of light and S is the radiation irradiance at the debris surface,
η
is the radiation incident angle to the debris surface normal,
β
is the debris surface specularity,
ρ'
is the debris surface reflectivity,
B_{f}
and
B_{b}
are the non-Lambartian coefficient of the front and back surfaces of the spacecraft, respectively, α' is the spacecraft absorption coefficient, and ε
_{f}
and ε
_{b}
are the front and back debris surface emissivity, respectively.
For non-perfectly reflecting surfaces, The force vector,
, is not directed normal to the surface. But, it inclines by an angle
ϑ
, known as the cone angle, to the incident direction,
as depicted in
figure 1
. So, the force components in the incident direction will be:
directed parallel to the Earth equatorial plane,
directed in the plane that contains the meridian of the sub-satellite point and
directed normal to the equatorial plane. As shown in
Fig. 1
, the incident radiation vector,
given by:
where
is the satellite position vector and
is the Earth radius vector (station coordinates) . For an oblate Earth, the Earth radius vector is given by
[4]
:
where
a_{e}
is the Earth’s Equatorial radius,
f'_{e}
is the Earth’s flattening,
φ_{g}
is the geodetic latitude,
h
is the height above sea level and
θ
is the sidereal time.
The satellite position vector,
in the geocentric coordinate system, is given by
[4]
:
with
where Ω is the longitude of the ascending node, ω is the argument of perigee,
v
is the true anomaly,
i
the inclination,
a
is the semi-major axis and
e
is the eccentricity. Using the geocentric coordinate system, the incident radiation vector is:
with
acting on a spacecraft is given by
The laser intensity delivered to the satellite surface
the general expression
[1]
,
[14]
and
[15]
:
where
is the vector from the spacecraft'scenter of mass to the element of the satellite projected area
dA
. The geocentric components of laser torque,
acting on a spacecraft is given by:
For a high laser repetition rate, many laser shots fire toward the satellite over an interval of time each one of intensity S. Consequently, the laser force can be considered as a continuous function. Thus, the total laser torque over an interval of time,
I=t_{1}-t
_{0}
, can be written as:
where
t
_{0}
is the time of the starting laser firing and
t
_{1}
is the time of its stopping.
In order to evaluate the previous integrals, the vector
must be determined. So, satellite geometry must be considered. In the next sections, three particular cases; circular cylindrical satellite, spherical satellite and satellite of complex shape will be studied.
where its origin
is at the geometric center of the satellite, as follows:
where
p
is the radius of the circular base and ϕ is the azimuthal angle. The position vector,
can be transformed into the following geocentric coordinate:
with
The illuminated surfaces of the cylinder satellite are a circular flat surface of radius ρ and an area
A
_{1}
= πρ
^{2}
in addition to a portion σ of the cylinder side of height H approximated and an area,
A
_{2}
= 2σπρ
H
, as in
fig. 2
.
Substituting into eqs. (14-16) and integrating over the areas
A
_{1}
and
A
_{2}
, laser torque applied on a cylindrical satellite is obtained.
The illuminated surfaces of the circular cylindrical satellite
where its origin
is at the geometric center of the satellite, as follows:
where ρ is the radius of the sphere, ϑ is the polar angle and ϕ is the azimuthal angle. The position vector,
can be transformed into the following geocentric coordinate:
The illuminated surfaces of the spherical satellite
with
The illuminated surfaces of the spherical satellite can be considered as a hemisphere with an area
A
= 2πρ
^{2}
, as shown in
fig. 3
.
Substituting into eqs. (14-16) and integrating over the area A, the laser torque applied on a spherical satellite is obtained.
and
In case of clear weather conditions, the propagation of a medium-powered 5 kW ground-based laser with a divergence of .1 mrad is illustrated in the following figure:
laser intensity delivered to LEO satellites of different altitudes
As illustrated in
fig 4
, the laser intensity has a significant dependency on the altitude. It has a maximum value ~1.4W/ m
^{2}
at the altitude of 314.4 km. However, it has a minimum value ~0.2W/m
^{2}
at an altitude of 823 km.
In the present work, the laser torque applied is calculated under the following postulates:
1-The laser source is located in the Helwan SLR station in Egypt with latitude of 29.86° N, longitude of 31.34° E and elevation of 145.46 m.
2-The radiation falls normal to the satellite surface (i.e.
η
= 0 ) and the satellite projected area is half of the total area.
3-The satellite surface is made of Aluminium, so it can be considered as a perfectly reflecting surface. Therefore, the laser force will be duplicated, where the satellite surface specularity
β
=1 and the satellite surface reflectivity
ρ'
= 1 While the spacecraft absorption coefficient α' = 0.
Based on the previous postulates and for a given time, the torque produced by a 5 kw laser beam is calculated for some
laser torque affecting on LEO cubesats at different altitudes using 5 kw laser beam at a given time.
Dependency of the laser torque on satellite eccentricity
low Earth orbit cubesats and the results are illustrated in the following table:
As illustrated in
table 2
, the satellites are arranged in ascending order according to their perigee height as it is clear in the table that the laser torque depends on the latitude, since the radiation force is inversely depended on the satellite altitude. The laser torque has a maximum value in the order of for a cubesat with a perigee height close to 600 km, however, it has a minimum value in the order for a cubesat perigee height close to 800 km.
In addition to the satellite height, the dependency of the ground based laser torque on the orbital eccentricity is studied and illustrated in
table (3)
.
As illustrated in table 3, the satellites are chosen to be nearly of the same altitude and inclination and they are arranged in descending order according to their eccentricity as it clear in the table that the laser torque directly depends on the orbital eccentricity.
^{-9}
which is comparable with that of solar radiation. However, it has a minimum value in the order of 10
^{-10}
which is comparable with that of the gravity gradient. Moreover, the laser torque depends directly on the orbital eccentricity whereas when the orbit becomes more circular it will experience less torque. Therefore, we can conclude that the laser torque can be used as an active attitude control system for cubesats.

linear atmospheric mechanism of laser interactions
;
radiation torque
;
The Attitude Determination and Control System (ADCS)
;
geocentric equatorial coordinate systems
;
spherical and cylindrical coordinate systems
;
LEO satellites and cubesats.

1. Introduction

The Attitude Determination and Control System (ADCS) is to stabilize the spacecraft against attitude disturbing influences resulting from the environment in the Earth’s orbit in order to orient it in the desired fixed pointing. Active control of the attitude of a satellite can be achieved by a number of different actuators: reaction wheels, thrusters, control moment gyroscopes or magnetic torquers. However, the magnetic torquers are widely used actuators for geostationary satellites, small satellites, and microsatellites. These high-tech devices interact with the Earth’s magnetic field and create a control torque, which can be adjusted to the required value. Combined with one or more reaction wheels, they provide all the control you need to maintain your spacecraft’s attitude
[8]
and
[14]
.
Unlike thrusters, magnetorquers are lightweight, reliable, and energy-efficient. A further advantage over momentum wheels and control moment gyroscopes is the absence of moving parts and therefore significantly higher reliability. However, they require a thoughtful design and careful assembly.
The main disadvantage of magnetorquers is that the strength of the magnets should be chosen to be strong enough to overcome the greatest expected disturbances (as given in
table 1
). A broader disadvantage is the dependence on Earth's magnetic field strength, making this approach unsuitable for deep space missions, and also more suitable for low Earth orbits as opposed to higher ones like the geosynchronous. The dependence on the highly variable intensity of Earth's magnetic field is also problematic since the attitude control problem becomes highly nonlinear. It is also impossible to control attitude in all three axes even if the full three coils are used, since torque can only be generated perpendicular to the Earth's magnetic field vector
[11]
and
[13]
.
The main issue of this work is to investigate the feasibility of using a medium power (5kw) to produce a torque on low Earth orbit cubesats in order to overcome some disadvantages of the traditional attitude actuators. The use of a ground based laser to produce a perturbing force goes back to the study of laser orbital perturbations by some scientists. The use of this medium power has the advantage of being inexpensive, and suffering a less amount of atmospheric attenuation. Moreover, it is safe in use where it cannot damage the satellite surface
[3]
,
[10]
and
[6]
.
Worst-case expected disturbance torques for a 1-U CubeSat at 700 km[11]

Lager Image

2. Laser Torque Model.

- 2.1 Laser Intensity.

The laser intensity delivered to the satellite surface is determined using an analytical model of long range laser beam propagation through the atmosphere. For standard atmospheric conditions, this model considers only the linear mechanism of laser atmospheric interactions. Moreover, the atmospheric turbulences are countered by using the adaptive optics and technical capabilities of the laser system. Based on that model, the laser intensity is proportional to the laser power and inversely proportional to satellite altitude and laser divergence. The laser intensity is given by
[2]
and
[10]
:
Lager Image

Lager Image

Lager Image

- 2.2 The laser force.

The total radiant force exerted on a flat non-perfectly reflecting surface is given by
[9]
:
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 2.3 The coordinate systems.

The geocentric equatorial system is used with the unit vectors;
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 2.4 The Laser Torque

The radiation torque,
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 2.5 Laser Torque Applied on Circular Cylindrical Satellite.

The position vector of the surface elements can be expressed in terms of a coordinate system with unit vectors
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 2.6 Laser Torque Applied on Spherical Satellite.

The position vector of the surface elements can be expressed in terms of a coordinate system with unit vectors
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

- 2.7 Laser Torque Applied on Satellite of complex shape.

In order to compute the torque applied on spacecraft of complex shape, we can follow the following scheme
[5]
and
[7]
:
1-Approximate each surface by means of simple geometric shapes (planes, cylinders, cones, spheres,…, etc.).
2-Determined the torque applied on each surface independently.
3-Then, the total torque applied over the whole spacecraft is obtained by the vector sum of all torques applied on each elementary surface.
3. Numerical Application.

Based on the standard atmospheric conditions and taking into consideration that
Lager Image

Lager Image

Lager Image

laser torque affecting on LEO cubesats at different altitudes using 5 kw laser beam at a given time.

Lager Image

Dependency of the laser torque on satellite eccentricity

Lager Image

4. Conclusion

A simple model of radiation torque affecting on satellites of various shapes is formulated using the geocentric equatorial coordinate system. Based on the model, the intensity delivered to the target is proportional to the laser power and inversely proportional to the beam divergence as given by eq. 1. So, we can adjust the laser power and divergence (by choosing the beam wavelength) in order to obtain the required control torque.
The current numerical test confirms that the laser torque has significant contribution on the low Earth orbit cubesats. The results show that the laser torque has a maximum value in the order of 10
Acknowledgements

The author is indebted for the helpful discussions andcomments of Prof. Kh. I. Khalil and Prof. El-Saftawy, M.I .

Beletskii V. V.
1966
"motion of an artificial satellite about itscenter of mass"
NASA translation
NASA TF ? 429

El-Saftawy M. I.
,
Afaf M. A. E.
,
Khalifa N.S.
2007
"Analytical Studies of Laser Beam Propagation through theAtmosphere”
Proceeding of 6th International Conference on Laser Science and Applications(ICLAS-07)
Cairo, Egypt

El-saftawy M.I.
,
Makram I.
2004
"The laser shots as aperturbing force on spacecraft's orbit"
NRIAG Journal ofAstronomy And Astrophysics
5
(1)
00 -
00

Escobal P.R
1965
Methods of Orbit Determination
John Wiley and Sons, Inc.
New York, London, Sydney

Harris M.
,
Lyle R.
1969
“Spacecraft Radiation Torque”, NASA Space Vehicle Design Criteria (Guidance and Control)
NASA SP-8027

James Mason
,
Jan Stupl
,
William Marshall
,
Creon Levit
2011
“Orbital Debris-Debris Collision Avoidance”
Journal of Advances in Space Research
48
(10)
1643 -
1655

Ping J.
,
Sengoku A.
,
Nagaoka N.
,
Iwata T.
,
Matsumoto K.
,
Kawano N.
2001
"How solar radiation pressure acts onRSAT and VSAT with a small evolving tip-off in SELENE"
Earth Planets Space
53
919 -
925

Karla Patricia Vega
“Attitude Control System forCubeSat for Ions, Neutrals, Electrons and MAGnetic Field(CINEMA)”
Graduate Division of Engineering and Mechanical Engineering, University of California, Berkeley
M.S. Thesis

Mcinnes Colin Robert
1999
Solar Sailing: Technology,Dynamics and Mission Applications
Springer-Praxis Series in Space Science and Technology
Chichester

khalifa N.S.
2009
“Effect of an Artificial Radiant Force onthe Spacecraft's Orbit”
Department Astronomy and Meteorology, Cairo university
Ph. D. Thesis

Rawashdeh S.
2009
“PASSIVE ATTITUDE STABILIZATION FOR SMALL SATELLITES”
College of Engineering, University of Kentucky
Kentucky
M.S. Thesis

Prilutsky Oleg F.
,
Fomenkova M.N.
1990
“Laser beamscattering in the atmosphere”
Science & Global Security
2
(1)
79 -
86

Francois-Lavet V.
2010
“Study of passive and activeattitude control systems for the OUFTI nanosatellites”
Faculty of Applied Sciences, University of Liège
M.S. Thesis

Wertz James R
1978
”SpacecraftAttitude Determination and Control
D. Reidel
Dordrecht, Holland

Zanardi M.C.
,
Vilhena de Moraes
1999
“ABC Analyticaland semi-analytical analysis of an artificial satellite'srotational motion"
Celestial Mechanics and Dynamical Astronomy
75
(4)
227 -
250

Citing 'Ground Base Laser Torque Applied on LEO Satellites of Various Geometries
'

@article{ HGJHC0_2012_v13n4_484}
,title={Ground Base Laser Torque Applied on LEO Satellites of Various Geometries}
,volume={4}
, url={http://dx.doi.org/10.5139/IJASS.2012.13.4.484}, DOI={10.5139/IJASS.2012.13.4.484}
, number= {4}
, journal={International Journal of Aeronautical and Space Sciences}
, publisher={The Korean Society for Aeronautical & Space Sciences}
, author={Khalifa, N. S.}
, year={2012}
, month={Dec}