2010:048 MASTER'S THESIS
Behavior and Relative Velocity of Debris near Geostationary Orbit
Lin Gao
Luleå University of Technology Master Thesis, Continuation Courses Space Science and Technology Department of Space Science, Kiruna
2010:048 - ISSN: 1653-0187 - ISRN: LTU-PB-EX--10/048--SE
CRANFIELD UNIVERSITY
SCHOOL OF ENGINEERING
MSc THESIS
Academic Year 2009-10
Lin Gao
Behavior and Relative Velocity of Debris near Geostationary Orbit
Supervisor: Dr. S.E.Hobbs
May 2010
This thesis is submitted in partial (45%) fulfillment of the requirements for the degree of Master of Science
©Cranfield University 2010. All rights reserved. No part of this publication may be reproduced without the written permission of the copyright owner.
i
Abstract
A general model is developed describing third-body gravity perturbation to debris’ orbit. Applying this model to debris released from geostationary orbit tells their motion in both short and long term. Without considering the Moon’s precession around the solar pole, the relative velocity between GEO debris can be calculated. This is an important coefficient for simulating the GEO debris environment and can serve as an input to break up models. ii iii
Acknowledgements
To my parents and dear friend M.Z. for her support. Sincerely thank you to Dr. S.E.Hobbs for his supervision. iv CONTENTS v
Contents
Contents v
List of figures viii
Abbreviations x
1 Introduction 1
1.1Background...... 1
1.2AimoftheThesis...... 1
1.3DocumentStructure...... 2
2 Literature Review 3
2.1GeosynchronousRegion...... 3
2.2SpaceDebrisModels...... 4
2.3PerturbationSource...... 5
2.3.1 NonhomogeneityandOblatenessoftheEarth...... 6
2.3.2 AtmosphereDrag...... 6
2.3.3 GravitationalPerturbation...... 7
2.3.4 SolarRadiationandSolarWind...... 7
2.4MotionoftheMoon...... 8
2.4.1 Cassini’sLaws...... 8
2.4.2 Relative Motion between the Moon, Sun, and Earth ...... 8
2.4.3 LunarStandstill...... 9 vi CONTENTS
3 Gravity Perturbation Model Development 11
3.1GeneralAnalysis...... 11
3.2GravityPerturbationModel...... 12
3.2.1 AnalyticalApproach...... 12
3.2.2 QuantitativeApproach...... 13
3.2.3 Validation...... 17
3.2.4 AmendmentofTheModel...... 17
4 Relative Velocity Between GEO Debris 23
4.1VerificationoftheModel...... 23
4.1.1 ComparisonBetweenOtherModels...... 23
4.1.2 Application of the Model to the Precession of Lunar Orbit . . 26
4.2RelativeVelocityoftheGEODebris...... 28
4.2.1 Motion of GEO Debris with Sun’s Perturbation ...... 28
4.2.2 CombinationofLunarandSolarGravityPerturbation.... 33
4.3TheRelativeVelocity...... 35
5 Discussion of the Results 39
5.1AboutTheModel...... 39
5.1.1 TheValidityoftheModel...... 39
5.1.2 TheMeaningofDevelopingsuchaModel...... 40
5.2AboutTheResult...... 40
6 Conclusions 42
References 43
A Calculation of iR 46
B Matlab code 48
B.1Numericallysolvethegroupofdifferentialfunctions...... 48 Contents vii
B.2PlotthePrecessionofLunarOrbitalPole...... 50
B.3 Plot Inclination and the Precession of GEO Debris’ Orbital Pole . . . 51
B.4PlotoftheRevolvingAscendingNode...... 53
B.5SimulationoftheDistributionofDebris...... 54
B.6SimulationoftheDistributionofRelativeVelocities...... 54
C Application of Matlab cftool toolbox 57 viii LIST OF FIGURES
List of Figures
2.1DistributionofGEODebris[1]...... 3
2.2PayloadsandupperstageslaunchedintoGEO[2]...... 4
2.3GEOregion...... 5
2.4ApparentpathsoftheSunandMoononthecelestialsphere...... 9
2.5ApparentmotionoftheSunandtheMoon.[3]...... 9
2.6Parametersofthelunarorbit...... 10
3.1MaximumrelativevelocitybetweenGEOdebris...... 12
3.2TypicalpositionbetweenSunandEarth...... 13
3.3MomentappliedbySuntoGEOdebrisatSolstice...... 13
3.4GeneralsituationforSolstice...... 14
3.5 Position of M and m in (x,y,z)and(x ,y ,z ) ...... 14
3.6 iR is in the direction of OA when the inclination of orbital plane i =0 18
4.1 The New Defined i ...... 24
4.2ComparisonofCoordinateSystems...... 25
4.3ThePrecessionofLunarOrbitalPole...... 27
4.4Moon’sOrbitalInclinationtoEclipticPlane...... 27
4.5TrackoftheOrbitalPoleofGEODebris...... 29
4.6TheInclinationoftheOrbit...... 29
4.7Debris’RelativeMotiontotheEarthSurface...... 30 LIST OF FIGURES ix
4.8 The Ascending Node is perpendicular to the Projection of the Orbital Pole...... 31
4.9 The Ascending Node Evolves with the Same Period as ω ...... 31
4.10 Curves have different amplitude and phase. They are all moving westwardsbecauseoftheEarthrotation...... 33
4.11DistributionofGEOdebris...... 33
4.12 The simulation of debris distribution is quite similar to the observation. 34
4.13TheDistributionofRelativeVelocities...... 37
A.1 Direction of iR changswithL ...... 46
C.1 Curve Fitting Toolbox Interface. Thetaz vs Time is plotted as blue andcurvefittingresultistheredline...... 57
C.2CurveFittingResult...... 58 x Abbreviations
Acronyms and definitions
GEO Geostationary Earth Orbit IADC Inter-Agency Space Debris Coordination Committee SHM Simple-Harmonic-Motion Lunar Orbital Pole the pole perpendicular to the lunar orbit and it points to the north Ecliptic Pole the pole perpendicular to the ecliptic plane and it points to the north Introduction 1
Chapter 1
Introduction
1.1 Background
Geostationary Earth Orbit (GEO) is widely used because of its unique characteris- tics. Its orbital period is exactly one sidereal day with the altitude of 35,786 km. Satellites using this orbit can thus stay above a fixed point relative to the Earth, which is very important for some satellites like communication satellites.
But after their lifetimes these satellites are no longer actively controlled and become orbital debris. Break-ups and explosions would further contribute to the debris population. Since there is no nature removal mechanism in GEO region, these debris would have lifetimes exceeding a million years. Their orbit would drift from 15 degrees north to 15 degrees south of the equatorial plane and 52 km above and below the geosynchronous arc because of the gravity perturbation from the Sun and the Moon [4]. The nonhomogeneity and oblateness of the Earth cause migration west and east around the Earth. Net effect of these motions leads to a torus around the Earth, where only 32% of the 1,124 known objects (2004) are under active control [1].
Because of the threat placed by these debris, risk and damage assessment are in- dispensable in spacecrafts design. In breakup models the initial condition, e.g., the relative velocity, is of extreme importance in the case of collision (rather than explo- sion or rupture) [5]. And the distribution of relative velocity between GEO debris is exactly what this thesis is pursuing.
1.2 Aim of the Thesis
The speed of any GEO objects is around 3.07 km/s with no more than 0.48% vari- ation. So the relative velocity is decided mainly by the direction, i.e., the orbit inclination. This places the necessity of a thorough understanding of third-body 2 Introduction gravity perturbation on debris’ orbit. After developing a model of gravity pertur- bation, the motion of GEO debris can be studied and relative velocity would be a direct result.
1.3 Document Structure
The key to know the motion of GEO debris is to understand the mechanism how the Sun and the Moon apply gravity perturbation to them. Most of the theories describing third-body gravity perturbation were done under the assumption that the orbital inclination i ≈ 0. Chapter 3 is the development of a perturbation model without this assumption.
In Section 4.1, solution of the model was compared to the work by Alby.F. It was further applied to the lunar orbit and results were quite optimistic. These prove that my model is reasonable and correct.
The motion of GEO debris in short and long term considering only the Sun is discussed in Section 4.2.1. SHM for GEO debris in short term is assumed to be a good approximation because of the similarity between the actual and computed result from the model considering only the Sun. In Section 4.2.2 detailed discussion about the approximation is provided. Calculation of the relative velocity between debris is then listed in Section 4.3. Literature Review 3
Chapter 2
Literature Review
Geosynchronous region is fairly crowded with retired satellites and fragments that were generated by break-ups, explosions and collisions. And the number of launches each year is generally increasing (Figure 2.2). In all the cataloged GEO objects only 31% are under control (Figure 2.1 (a)). Break-ups of spacecrafts contribute about 43% of catalogued objects, and 85% of all space debris larger than 5 cm in diameter [6]. Explosions and collisions happen less frequently but contribute about 50 percent of all tracked objects [1]. Figure 2.1 (b) is a plot of the GEO debris torus.
(a) cataloged GEO objects (b) debris torus around GEO
Figure 2.1: Distribution of GEO Debris [1].
2.1 Geosynchronous Region
Because of the uniqueness of this orbit, [7] defines the protected geosynchronous region as a segment of the spherical shell with:
• lower altitude = geostationary altitude minus 200 km • upper altitude = geostationary altitude plus 200 km 4 Literature Review
Figure 2.2: Payloads and upper stages launched into GEO [2]
• -15 degrees ≤ latitude ≤ +15 degrees
• geostationary altitude (ZGEO) = 35,786 km (the altitude of the geostationary Earth orbit)
This region is where this thesis focuses at. [7] requires that spacecrafts be pushed away from geostationary orbit at the end of the disposal to avoid interference with active spacecrafts. According to their guide, the satellites must meet the following two conditions:
1. A minimum increase in perigee altitude of: 235km+(1000·CR · A/m)
2. An eccentricity less than or equal to 0.003.
Here CR is the solar radiation pressure coefficient. A/m is the aspect area to dry mass ratio (m2kg−1).
2.2 Space Debris Models
NASA’s breakup model of EVOLVE 4.0, also called the NASA Standard Breakup Model, is quite helpful in designing satellites by providing results of fragmentations Literature Review 5
Figure 2.3: GEO region in the debris environment. Initial conditions of the breakup such as total mass of the parent object and the collision velocity are basic inputs of the model. [5]
There are some well developed models about debris environment, such as EVOLVE (NASA) and MASTER-2005 (ESA). They both emphasize at estimating the future evolution of the debris environment. Take EVOLVE as an example, typical time frame is 100 years and its premise is to predict future debris environment according to recent launch and debris degradation rates [8]. MASTER (Meteoroid and Space Debris Terrestrial Environment Reference) is a software used to analyze space debris flux and spatial densities [9].
2.3 Perturbation Source
Under the ideal case that GEO debris experience no other force than Earth gravity, they will continuously move in the geostationary orbit with the same speed and relatively still to neighboring objects. Then there is no need to worry about them. Unfortunately that does not exist in practice. In practice perturbing accelerations acting on space objects result from nonhomogeneity and oblateness of the Earth, atmosphere drag, third-body gravitational perturbation, solar radiation and solar wind. 6 Literature Review
2.3.1 Nonhomogeneity and Oblateness of the Earth
The Earth is oblate and its mass is not distributed homogeneously. This results in the difference between reality and the ideal model of the Earth gravity potential function as:
μ U(r)=− (2.1) r
[10] gives a simplified form of U: ∞ μ Re n U ≈− 1 − ( ) JnPn(sφ) r n=2 r μ = [U0 + UJ2 + UJ3 + ···] (2.2) r
− Re 2 1 2 − here U0 = 1, UJ2 =( r ) J2 2 (3 sin φ 1). J2 is at least 400 times larger than any other Jn coefficients. As a result the orbit of GEO debris can be viewed as getting perturbation from UJ2. The Lagrange’s planetary equations can be used: da 2 ∂U = , (2.3) dt na ∂M √ − 2 − 2 de 1 e ∂U − 1 e ∂U = 2 2 , (2.4) dt na e ∂M na e ∂ω di −1 ∂U ∂U = √ +cos(i) , (2.5) dt na2 1 − e2 sin(i) ∂Ω ∂ω dΩ 1 ∂U = √ , (2.6) dt 2 − 2 ∂i √na 1 e sin(i) dω 1 − e2 ∂U cos(i) ∂U = − √ , (2.7) dt na2e ∂e na2 1 − e2 sin(i) ∂i dM 2 ∂U 1 − e2 ∂U = n − − (2.8) dt na ∂a na2e ∂e
Substitute UJ2 into Eq. 2.3, Eq. 2.4 and Eq. 2.5: da de di = = = 0 (2.9) dt dt dt
This means the change of da, de and di is zero for each period.
2.3.2 Atmosphere Drag
GEO orbit locates in the exosphere of the Earth atmosphere. Here particles are so far away thus hardly collide with each other. The atmosphere drag can be ignored. Literature Review 7
2.3.3 Gravitational Perturbation
Theoretically all planets and the Sun apply gravitational perturbation on GEO debris. According to [10] its magnitude is determined by Kp :
Gmp Kp = 3 (2.10) nrp where G is Gravitational Constant, mp is the mass of the perturbing body, rp is its distance from Earth center, n is the mean motion of the GEO object. And
−3 Kmoon =5.844 × 10 deg/day (2.11) −3 Ksun =2.69 × 10 deg/day (2.12)
Even the closet planet to the Earth, Venus, its distance from Earth is at least (depends on its relative position with the Earth) 108 times that of the Moon while itsmassisonlyabout4.49 times the Moon’s mass. The planet with highest mass in the solar system is Jupiter, which is 2.58 × 104 times heavier but also 1.64 × 103 times further than Moon from the center of Earth. As a result, none of these planets will apply influence to the orbit of GEO debris by more than 10−5 times of Moon’s influence. And only gravitational perturbation from the Sun and the Moon need to be considered.
2.3.4 Solar Radiation and Solar Wind
The Sun radiates electromagnetic waves from X-rays to radio waves, which cause pressure on any objects in the GEO orbit and change the eccentricity of the their orbit. Solar wind contains emitted ionized nuclei and electrons from the Sun. The momentum flux caused by solar wind is much less than that caused by the solar radiation by a factor of 100 to 1,000 [10]. As a result in this thesis only solar radiation pressure is considered.
The solar radiation pressure |FR| is the product of the momentum flux P , the cross- sectional area A of the debris perpendicular to the sun line and the absorption characteristic coefficient Cp: |FR| = PACp (2.13)
Cp range from 0 to 2; Cp = 1 for a black body and Cp = 2 for a body reflecting all light.
In the orbits around the Earth, P is:
452.7 2 P = × 10−8(N/m ) (2.14) 1.0004 + 0.334 cos(D) 8 Literature Review
D is the “phase” of the year starting on July 4, the day of earth aphelion [11].
From Eq. 2.13 and Eq. 2.14:
452.7ACp −8 |FR| = × 10 N (2.15) 1.0004 + 0.334 cos(D)
Eq. 2.15 shows that the solar radiation pressure relates directly to the “phase” of the year.
2.4 Motion of the Moon
2.4.1 Cassini’s Laws
The Cassini’s laws describe the motion of the rotation of the Moon based on obser- vation. [12]
1. The Moon rotates counterclockwise with an uniform motion around an axis fixed with respect to its surface; the period of rotation is the same as the Moon’s sidereal period of revolution around the Earth.
2. The angle between the rotation axis and the ecliptic pole of the Moon is constant ≈ 1◦35 .
3. The ecliptic pole of the Moon, the axis normal to the Moon’s orbit and its rotation axis are always in the same plane.
2.4.2 Relative Motion between the Moon, Sun, and Earth
The Moon revolves around the Earth in a month (27.3 days), the Earth revolves around the Sun in a year (365.25 days). The Sun’s track in the celestial sphere lies in a fixed plane tilted 23.44◦ to the equatorial plane, while the lunar orbital inclination varies from 23.736◦ +5.133◦ =28.869◦ to 23.736◦ − 5.133◦ =18.603◦ (Figure 2.4). This is caused by the precession of the lunar pole around the ecliptic pole with an cone angle of 5.13◦ (Figure 2.5).
Orbital parameters of the Moon is defined in Figure 2.6: ism, im and is are the inclination between the ecliptic orbit and the lunar orbit, the lunar orbit and the equatorial orbit, the ecliptic orbit the equatorial orbit respec- tively. Ω is the right ascension angle of the lunar orbit in the ecliptic orbit from the vernal equinox axis. [10] gives an approximate expression of im as:
im[deg] = 23.736 + 5.133 sin [ωm(t − 2, 001.7433)] (2.16) Literature Review 9
Figure 2.4: Apparent paths of the Sun and Moon on the celestial sphere.
Figure 2.5: Apparent motion of the Sun and the Moon.[3]
here ωm is the frequency of the precession of lunar pole with a period of 18.6 Years, t is time expressed in years.
2.4.3 Lunar Standstill
Due to the gravity of the Sun, the Moon’s orbit is not fixed relative to the Earth but precessing around the solar pole with period of 18.6 Years. With constant obliquity of the ecliptic of 5.133◦, the inclination of lunar orbit to the equatorial plane varies from 18.603◦ to 28.869◦. Moving in the celestial sphere, the orbit inclination is the extreme declination that the Moon reaches in a month. When the inclination is 10 Literature Review
Figure 2.6: Parameters of the lunar orbit.
18.603◦ (called minor standstill) or 28.869◦ (called major standstill), the change of inclination is very slow. The most recent lunar major standstill happened in 2006 with consequences be able to be observed for more than 3 years. Gravity Perturbation Model Development 11
Chapter 3
Gravity Perturbation Model Development
3.1 General Analysis
According to the definition of GEO region [7], the radius of concerned orbits range from (rGEO-200) km to (rGEO+200) km, rGEO = 42164 (km) is the GEO orbital radius. So the minimum velocity of GEO objects is: μ vmin = rGEO + 200
μ is the geocentric gravitational constant, μ = 398600.4415 ± 0.0008 km3s−2 [13].