- Open Access
- Total Downloads : 227
- Authors : Pravin. V. P, T. Velmurugan, S. Balamurugan
- Paper ID : IJERTV3IS050603
- Volume & Issue : Volume 03, Issue 05 (May 2014)
- Published (First Online): 15-05-2014
- ISSN (Online) : 2278-0181
- Publisher Name : IJERT
- License: This work is licensed under a Creative Commons Attribution 4.0 International License
Orbital Maneuvering using Electrodynamic Tethers
Pravin.V. P1, T. Velmurugan2
Department of Mechanical Engineering, Regional Centre of Anna University, Tirunelveli Region,
Tamil Nadu, India
Scientist/Engineer, Flight Mechanics Division, Vikram Sarabhai Space Centre (VSSC), Trivandrum, Kerala, India
Abstract Orbital debris has influence in future space missions considerably. A number of catastrophic collisions are expected in the following years. But to remove these threats from low earth orbit, the conventional techniques using propellant are not economically efficient. This innovative topic shows the method to remove debris using electrodynamic tethers, which is a non-conventional technique. The main objective is to use electrodynamic tethers for effective transfer of a debris removal satellite to different orbits for the collection of debris or attaching some deorbiting mechanism on the debris in the most optimal way. In brief, the paper defines optimal path for the debris removing vehicle to visit a predetermined number of debris. Matlab software is used as an optimization tool in this work. For this purpose various kinds of maneuvers that a satellite can perform are analyzed and the coding in Matlab for each maneuvers are performed. This coding is done in order to derive the delta-v required for various types of maneuvers. The Lamberts problem was used to determine the velocity vectors of the vehicle in orbit on defining the position vectors. Also the optimal path, which the vehicle can travel using electrodynamic tethers is solved using Matlab code. The method is chosen to reduce the payload required and thereby reduce the cost of maneuvering.
Keywords Electrodynamic Tethers (EDT), Travelling Salesman Problem (TSP), Lamberts Problem, Lorentz force, Orbital debris, Low Earth Orbit (LEO), Oblate Earth effect, RAAN rate
Space explorations started several centuries back. In the beginning, these explorations helped in formulation of assumptions and theories regarding the evolution of earth. The Low Earth Orbit (LEO) extending from 160km to 2000km is the main operational area for active satellite, but due to the presence of space debris it is highly unstable. There are 300000 debris of the size around 1cm and 19000 debris of size more than 5cm including the inactive satellites and junks existing in the Low Earth Orbit. A number of collisions between artificial satellites have occurred in the recent past. One of the major events is the anti-satellite test by Chinese researchers in 2007 which completely destroyed the satellite
Feng Yun and created 2900 pieces of debris with size greater than 10cm. In 2009, the biggest hypervelocity collision ever
happened when the satellite Iridium 33 collided with another satellite named Kosmos 2251 . Iridum 33 was active during the collision event, but the inactive Russian satellite Kosmos 2251 (Inactivated in 1995) was left unnoticed which resulted in the collision. The collision created 1600 fragments which added threat to other active satellites. EDT can play a key role in changing the potential risks in the removal of debris. The method is cost effective because it doesnt need any fuel, fuel tank or additional propulsion accessories. All it requires is an electron emitter, balancing mass and long tether. Some experiments were already carried out by some agencies and organizations regarding the use of EDT in space.
Kessler Syndrome is the phenomena in which, the collision of one debris with another debris or an active satellite causes a series of collision process. Even if no launches are done from now on it cant be avoided. Only possibility is to remove the potentially risky debris.
The Lamberts problem is solved to find the velocity vectors. Gravity field variation for oblate earth effect also determined. The computation of velocity increment required for various orbital maneuvers are also coded in matlab. In order to determine the orbital elements and define the path of the vehicle, through different orbits Ode45 is used. The transfer from one debris orbit to another is found for a given current, tether length and initial conditions the orbit at any time can be defined.
As a future work travelling sales man problem is preferred to solve the logical concept involved in path optimization. However the problem solving using Branch and Bound algorithm is under progress.
Orbital maneuvers are performed in space for various purposes such as altitude raising, path changing, collision avoidance, etc. The v determinations for the following orbital maneuvers are completed .
A. Hohmann Transfer
In hohmann transfer, Va is provided at point A. This allows the satellite to enter an elliptical trajectory. After it reaches the point B in final orbit another firing is done in opposite direction as Vb.
The total V required is the sum of above velocity increments.
Vtotal = Va + Vb (1)
D. Phasing Maneuver
Phasing maneuver is used to change the position of a satellite in an orbit in order to rendezvous with another space object in the same orbit. From figure 5, it can be seen that the new orbit formed will be very much close to the initial orbit and dock with target, which allows proceeding the maintenance work, refueling, etc.
Time Period, T
2 a 32
Fig 1: Hohmann Transfer
Bi-elliptic Hohmann Transfer
In the case of bi-elliptic hohmann transfer, two semielliptical transfer paths used for maneuvering as shown in the diagram from A to B and B to C. The transfer orbits must be coaxial.
Where, Gravitational parameter, Âµ GM 398600 km/s4 G-Gravitational constant of earth
M- Mass of earth; a- Semimajor axis
Vtotal = Va + Vb + Vc (2)
Fig 2: Bi-elliptic Hohmann transfer
Apse Line Rotation
Fig 4: Phasing Maneuver
Non-Hohmann transfer with common Apse Line
This transfer strategy is used to transfer between orbits in same plane and sharing a common apse line. Since the non- hohmann transfer is preferred, the transfer trajectory is not necessarily a tangent to either orbit. v is same as apse line rotation.
ra , rd Radius at initial and final orbit.
a , d True anomaly at initial and final orbit.
Apse line can be said as an imaginary line which is defined by an orbits eccentricity vector. The satellite can be shifted from one orbit to another at the point at which both the
orbits intersect each other (A and B).
ra cos a
rd cos d
2 V2 2 V V cos
cos a cos d ra cos a rd cos d
2 1 2
V1,V2 Velocity of satellite at orbits 1 and 2.
1,2 True anomalies of satellite at orbits 1 and 2. Apse line rotation angle.
Fig 3: Apse line rotation
Fig 5: Non-Hohmann transfer with common apse line
Chase maneuver is used to chase down a target space object. The chase maneuver is performed to catch up with an object and dock with it. Docking mechanism is used for maintenance, refueling, deorbiting, etc.
Optimization techniques usually include an iterative process which determines the solution for a given problem.
Time Period, T
Lamberts problem is aimed at determining the velocity vectors from the given position vectors (r1, r2) of a vehicle in
1 e 2
e cos b
any orbit .
Through Newtons iterative process, the value of z can be accurately determined. This z value can be used for further calculations where A is a constant determined by the true anomaly and the magnitude of position vectors. The functions
a , b True anomaly at initial and final orbit.
h- Specific Angular momentum of the orbit
Âµ- Gravitational parameter
2 Sz A
A z Sz 1
C z 1 z
The Newtons iterative process is,
z i 1
Fig 6: Chase Maneuver
Plane Changing Maneuver
The plane changing maneuver is a complex one because the velocity increment should be given at an angle so that the
After iterative process, the z value is introduced in the Lagrange variables given by,
satellite moves into an orbit in another plane. In this case, both the orbits will share a common focus even though they dont have a common plane.
y z z S z 1
V V 2
V2 V2 2V V
r 1 r 2
t 1 t 2
t 1 t 2
Vr ,Vt Radial and Tangential velocity at an orbit
Now these Lagrange variables can be directly introduced to determine the velocity vector.
v1 1 r f r
g 2 1
v2 1 g r r
g 2 1
Fig 7: Plane changing maneuver
After determining the velocity vector the orbital elements can be easily determined. The orbital elements are Semimajor Axis (a), Eccentricity (e), RAAN (), Argument of Perigee (), Inclination (i) and True anomaly () as indicated in the figure 8.
Fig 8: Orbital Elements
B. TRAVELLING SALESMAN PROBLEM
TSP is a logical approach towards the complex problem. The minimal path can be found using this technique. For example, to visit 5 cities in the most optimal way as illustrated in figure 9, we need to choose the most efficient path from the combination of all the paths. Branch and bound algorithm is the method used for this purpose.
The branch and bound algorithm consists of a number of constraints. It eliminates the chance of sub tours between the cities which creates more complexity to the problem. The following diagram shows one of the possible paths to travel between the cities.
Fig 9: Travelling Salesman Problem
In the case of space object, the condition becomes more complex, because the distance between space objects vary with time increment.
TABLE I. Velocity increment for various maneuvers
Velocity Increment Required
First orbit apogee radius = 11000 km First orbit perigee radius = 9000 km Second orbit radius = 20000 km
Bi-elliptic hohmann transfer
First orbit radius = 10000 km Second orbit radius = 20000 km First elliptical transfer orbit apogee radius = 15000 km
Apse line rotation
First orbit apogee radius =15000 km First orbit perigee radius = 9000 km Second orbit apogee radius =12000 km Second orbit perigee radius =8000 km Apse line rotation angle=45 degree
First orbit apogee radius =11000 km First orbit perigee radius = 9000 km Second orbit true anomaly= 90 degree
Apogee radius =20000 km Perigee radius =10000 km
True anomaly of chaser = 45 degree True anomaly of target = 150degree Time required for the chase = 1 hr
Non-hohmann transfer with common apse line
Enter the apogee of initial orbit 20000 Enter the perigee of initial orbit 10000 True anomaly of initial orbit with respect to apse line 45
Electrodynamic tethers are long conducting wires which can be used in satellites for an economic maneuvering process as shown in figure 10.
The basic principle of tethered satellite is that of a motor. That is the conversion of electric potential into kinetic energy in the presence of a magnetic field. In real case the electric potential across the tether in the presence of earths magnetic field generates a Lorentz force.
The magnetic field around the earth is affected by so many factors. Among them the main factors are,
Magnetic and rotational tilt of earth
Solar wind effect
To determine the magnetic field action Euler-Hill reference frame is considered and the magnetic field components are noted .
In order to set the orbital elements, the basic parameters required are, the semi lactus rectum, specific angular momentum and the radius of the orbit. The parameters are given by,
p a1 e2
1 e cos
In order to find the trajectory of satellite after the tether based perturbative acceleration, the rate at which the orbital elements change over a period of time is found.
By using the Gauss form of variational equation,
The components of magnetic field vector can be taken as,
da 2 a 2 p
B 2 m sin sin i
e sin f r f t
dt h r
B m cos sin i
de 1 p sin f r p rcos r ef t
B m cos i
1 p cosf
k r 3
dt h e
r sin cos i
h sin i
d r sin
h sin i
di r cos f
d h 1
p cos f r p rsin f t
dt r 2 h e
Fig 10: Electrodynamic Tethers
Components of perturbative acceleration act along the radial, tangential and normal directions of the orbit.
Thus for a non-liberating EDT, the perturbative accelerations are,
f r IL Bk sin cos B j sin
The rate of change in orbital parameters can be found out using variational method. All the six orbital elements of initial orbit are essential in further calculation.
The rate of orbital element change with time can be used over a period of time to get the exact position of the spacecraft.
IL Bi sin Bk
Fig 11: Tether Deployment
f h IL B j cos cos Bi sin cos
Fig 12: Orbital element variation generated in Matlab
For a given, mass of satellite = 420kg ; Current applied
= 5 A ; Liberation angle, = = 0 deg ;Length of tether = 15km; Semiajor axis = 6878km ; Eccentricity = 0.02;RAAN= 30 deg ; Inclination= 98 deg ; Argument of perigee = 50 deg; True anomaly = 0 deg; Time = 10 days
The result shown in the figure 12 was generated using ode45 in Matlab. Ode45 is used to solve the ordinary equations with given initial values. Orbital element variation over the orbits with time increment is determined from the graph obtained.
The Semimajor axis shows huge variation considering other orbital elements. Over a period of 10 days, the semimajor axis varies by a distance of above 1000km.
The eccentricity shows slight variation, which shows that the satellite will still be in an almost circular orbit.
The true anomaly shows much higher variation because of the position change over time.
RAAN rate has based on natural perturbations are not considered here. Only effect of the EDT is considered. The inclination is the angle created between the equatorial plane and the satellite orbital plane which also shows slight change during the process.
Orbital debris removal in the most cost effective way.
Refueling of active satellites in the earth orbit.
For maintenance of active satellites in space.
Kessler Syndrome can be avoided.
To remove the potentially risk debris from the path of active satellites to avoid collision.
The space debris has become a huge threat to the future space missions. Removing the potentially risk debris will avoid the risk of collision with active satellites. NASA is now researching the scope of electrodynamic tethers in this field. In this work, the travel between one orbit to another using electrodynamic tethers is determined. The result shows that within a short span of time tethers can maneuver the satellite. Huge velocity increment required in conventional method can be avoided. As a future work, TSP inclusion into the work is under process. By including TSP any number of debris can be visited by the spacecraft in optimal way.
The authors gratefully acknowledge the support of FMD division in VSSC Trivandrum for their constant encouragement and support in course of the activity.
V.Adimurthy, A.S.Ganeshan, Space Debris Mitigation Measures In India, 54th International Astronautical Congress of the International Astronautical Federation, the International Academy of Astronautics and the International institute of Space Law, Bremen, Germany, 29 Sept-3 Oct, 2003, Paper Number: IAC-03-IAA.5.4.09
Brent W Barbee, Salvatore Alfano, Elfego Pinon, Kenn Gold, David Gaylor, Design of spacecraft missions to remove multiple orbital debris objects, Advances in the Astronautical Sciences, Vol. 144, Univelt, Inc., San Diego, CA, 2012, pp. 93-110
Alexander Karl, Active Removal of space debris-discussing technical and economical issues, International Astronautical Congress, 2006,
Paper Number: IAC-06-B6.4.04
Paul Williams, Optimal Orbital Transfer with Electrodynamic Tether,
Journal of Guidance, Control, and Dynamics, Vol.28, No.2, March-
Howard.D.Curtis, Orbital Mechanics for Engineering Students, Second Edition, Elsevier Publication, 2010
R.Zhong, Z.H.Zhu, Long term dynamics and optimal control of nano- satellite deorbit using a short electrodynamic tether, Advances in Space Research, Oct 2013, pp.1530-1544