Advanced
Ground Base Laser Torque Applied on LEO Satellites of Various Geometries
Ground Base Laser Torque Applied on LEO Satellites of Various Geometries
International Journal of Aeronautical and Space Sciences. 2012. Dec, 13(4): 484-490
Copyright ©2012, The Korean Society for Aeronautical Space Science
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
Author
Metrics
Cited by
TagCloud
About the Authors
N. S. Khalifa
asmaa_2000_2000@yahoo.com

Abstract
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.
Keywords
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
Worst-case expected disturbance torques for a 1-U CubeSat at 700 km [11]
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
where, u is the satellite altitude, P is the laser power, θ is the laser divergence, ϕ is the elevation angle,
Lager Image
is the molecule scattering coefficient at sea level,
Lager Image
is the aerosols scattering coefficient and is the sea level altitude.
- 2.2 The laser force.
The total radiant force exerted on a flat non-perfectly reflecting surface is given by [9] :
Lager Image
where,
Lager Image
where
Lager Image
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, Bf and Bb 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,
Lager Image
, is not directed normal to the surface. But, it inclines by an angle ϑ , known as the cone angle, to the incident direction,
Lager Image
as depicted in figure 1 . So, the force components in the incident direction will be:
Lager Image
- 2.3 The coordinate systems.
The geocentric equatorial system is used with the unit vectors;
Lager Image
directed parallel to the Earth equatorial plane,
Lager Image
directed in the plane that contains the meridian of the sub-satellite point and
Lager Image
directed normal to the equatorial plane. As shown in Fig. 1 , the incident radiation vector,
Lager Image
given by:
Lager Image
where
Lager Image
is the satellite position vector and
Lager Image
is the Earth radius vector (station coordinates) . For an oblate Earth, the Earth radius vector is given by [4] :
Lager Image
Lager Image
where ae 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,
Lager Image
in the geocentric coordinate system, is given by [4] :
Lager Image
with
Lager Image
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:
Lager Image
with
Lager Image
Lager Image
Lager Image
- 2.4 The Laser Torque
The radiation torque,
Lager Image
acting on a spacecraft is given by
Lager Image
The laser intensity delivered to the satellite surface
the general expression [1] , [14] and [15] :
Lager Image
where
Lager Image
is the vector from the spacecraft'scenter of mass to the element of the satellite projected area dA . The geocentric components of laser torque,
Lager Image
acting on a spacecraft is given by:
Lager Image
Lager Image
Lager Image
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=t1-t 0 , can be written as:
Lager Image
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
Lager Image
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.
- 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
where its origin
Lager Image
is at the geometric center of the satellite, as follows:
Lager Image
where p is the radius of the circular base and ϕ is the azimuthal angle. The position vector,
Lager Image
can be transformed into the following geocentric coordinate:
Lager Image
with
Lager Image
Lager Image
Lager Image
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.
Lager Image
The illuminated surfaces of the circular cylindrical satellite
- 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
where its origin
Lager Image
is at the geometric center of the satellite, as follows:
Lager Image
where ρ is the radius of the sphere, ϑ is the polar angle and ϕ is the azimuthal angle. The position vector,
Lager Image
can be transformed into the following geocentric coordinate:
Lager Image
Lager Image
The illuminated surfaces of the spherical satellite
with
Lager Image
Lager Image
Lager Image
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.
- 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
and
Lager Image
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:
Lager Image
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.
Lager Image
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
Lager Image
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.
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 -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.
Acknowledgements
The author is indebted for the helpful discussions andcomments of Prof. Kh. I. Khalil and Prof. El-Saftawy, M.I .
View Fulltext  
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