MODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION
|
|
- Katrina Barker
- 6 years ago
- Views:
Transcription
1 MODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION 1 SHEN YU JUN, 2 TAN YAN QUAN, 3 TAN GUOXIAN 1,2,3 Raffles Science Institute, Raffles Institution, 1 Raffles Institution Lane, Singapore yujunauthor@gmail.com Abstract- An orbit determination and propagation package is developed in Python to predict satellite trajectories accurately with position differences of within 1 km after one orbit when compared to simulated satellite true position. The effect of Earth s non-uniform gravitational fields, caused by equatorial bulge and non-sphericity, is modelled using zonal spherical harmonics from J2 to J6 and the orbit is propagated numerically via Runge-Kutta-4 (RK4) method. The full J2-6 gravitational coefficients model yields improvement in position accuracy over the J2-only model. Furthermore, the use of additional state vector observations via differential correction allow for refinement of initial orbital parameters, hence achieving good agreement with simulated satellite positions. The orbit determination and propagation software package is tested using two real-life applications: satellite collision avoidance and Earth observation. The first application is a simulation of the orbital hypervelocity collision between Iridium 33 and Cosmos 2251 using publicly available NORAD Two Line Elements (TLE) data. The two satellites simulated time of closest approach is 16:55:40 UTC, 10 February 2009, accurate to within one minute. The second application is a geo-registration of the International Space Station (ISS) using photographs taken by ISS Expedition 35 crew and corresponding ISS TLEs, calculatingits position over a rotating Earth coordinate frame and generating a ground-track. The location of ISS over an unknown longitude and latitude location on 15:01:42 UTC 27 March 2013 was interpolated to be 1.00 N, E, matching actual ISS position to within 0.1 degrees. Key words- Orbit Propagation, Orbit Determination, Differential Correction, Conjunction Analysis, Geo-registration. I. INTRODUCTION Accurate knowledge of a satellite s orbit is required for mission planning and geo-registration of ground positions. Hence, orbit determination must return reliable orbital elements data, which can provide more accurate positions of the satellite across time. Orbit determination improves initial orbital elements from a set of observation data, such as state vectors (position, velocity) or angles (azimuth, elevation) [1]. These initial elements uniquely define an orbit and thus a satellite s position in time [2]. Fig. 1 shows the flow from observation or sensor-collected data to refine position and velocity information. For an Earth orbit satellite, dominant perturbations are caused by the non-uniform gravitational field around an oblate Earth [3, 4]. In this paper, we modelled J2 to J6 coefficients using numerical propagation by Runge- Kutta-4 (RK4) method to develop an orbit propagation and determination model in Python. To test the orbit propagation and determination model s effectiveness in actual orbital collision avoidance (conjunction analysis) and photo georegistration applications, we calculated the 2009 satellite collision between Iridium 33 and Cosmos 2251 based on open-source Two Line Elements (TLE) data [5]. This was the first, catastrophic hypervelocity orbital collision [6]. The accuracy of the orbit propagator determines the confidence of two objects closest approach distance, ensuring safe space operations. Additionally, ground identification was performed using actual International Space Station (ISS) photos taken during Expedition 35 in 2013 [7] and corresponding TLE sets. II. PROPAGATION MODEL WITH J2-6 COEFFICIENTS Classical two body interaction results in perfect elliptical motion with only a central force from a point-like mass. In reality, for an Earth-orbiting satellite propagation, dominant perturbation effects are the J2 attraction (due to Earth s equatorial bulge) and zonal spherical harmonics from J3 to J6 (higher order non-sphericity terms) [2]. The effects of higher order J- gravitational terms can be solved numerically as no analytical solution exists [8]. The classic RK4 numerical method is used [2]. A detailed propagator flowchart is presented in Fig. 2. The equation of motion is Newton s Second Law (Eq. 1) with gravitational harmonics attached. Fig.1. Overview of orbit determination steps 198
2 Fig.2. RK4 orbit propagation function In Eq. 1, the J x and J y terms are symmetric due to equatorial bulge. All three J terms depend on inclination and are described by Eqs. 2& 3. using STK [10] to simulate satellites position and velocity. Gaussian noise was then added to one set to simulate GPS data ( GPS data), while the other data set was unchanged ( real data). The position error between the propagated result and real set is computed using at every time step by Eq. 4. Fig. 3 shows the position error for the propagator compared to true satellite positions using real simulated data without any noise, as a test of propagator accuracy only. TheJ2-6 model (the entirety of Eqs. 2 & 3) results in 30 km position error over 1 day. The J2onlymodel (up to J 2 in Eqs. 2 &3) gives up to 40 km in position difference, with large periodic fluctuations. In particular, between 500 min to 800 min, the J2 error is larger than J2-6 case by up to 50%. Hence, a thorough numerical solution for J2-6 model yields much better consistency and reliability of propagation. The J2-6 model will be used subsequently throughout this paper. Orbit propagation was developed in numerical and scientific Python libraries, using constants for Earth equatorial radius r E, gravitational coefficients J 2 to J 6 [2].The inputs are initial position and velocity and output is the satellite s full trajectory for a specified mission duration. III. PROPAGATION MODEL ACCURACY WITH SIMULATED GPS DATA To verify the accuracy of the RK4 propagator section, results were compared with simulated satellite position and the error plotted over time. For this paper, two HPOP ephemeris data were generated Fig.3. J2 and J2-6 models position error with real position over 1 day 199
3 Fig.4. Differential correction to improve initial orbital elements IV. ORBIT DETERMINATION WITH DIFFERENTIAL CORRECTION and velocity data, since observations currently come from the same source. (Eq. 7) Based on preliminary elements and subsequent propagation, we can improve our knowledge of the initial orbital elements by considering the rest of the observation data, such as GPS readings [1, 9], to obtain more accurate propagation. The improvement function used in this paper is differential correction as illustrated in Fig. 4. It is coded in numerical Python, incorporating the earlier propagation function. For any propagated quantity in the state vector at time n, the total error comes from initial error in each of the orbital elements used to propagate, in this case classical orbital elements a, e, i, ω,ω and θ. Eq. 5 shows the error in x displacement; the y and z errors are calculated similarly. Each partial derivative at each time step can be linearly approximated with the difference in propagated quantity (x, y, z) when only one initial element is perturbed (by 1% variation) [1]. (Eq. 6) Hence, we can form the residuals matrix B and derivatives matrix A as required in Fig. 4 to solve for the correction matrix C. An identity matrix W is used to weigh the reliability of the observed position To compute A, only the x, y, z components of state vector and their dependence on initial orbital 200
4 elements(classical or Cartesian) are considered. To simplify calculations, we neglect velocity terms. In B, additional observation data is introduced. Thus, A is 3n by 6 and B is 3n by 1, where n is the number of observations. x, y, z components must be calculated individually to obtain a unique trajectory in 3D space [9]. If we have more observation readings than elements, then a Root Mean Square fit is used. Finally, the differential correction process is iterated, using the improved orbital elements. In essence, initial orbital elements are refined. V. DIFFERENTIAL CORRECTION EFFECTIVENESS TEST propagation. Yet, with increasing number of additional readings, the gradient tends towards a minimum, which is bounded by the error between the J2-6 propagator and the commercial HPOP propagator. This minimum occurs when the initial positional error has been fully rectified by orbit determination. Practically, orbit determination is performed using only a limited number of data sets to propagate forward for a much longer epoch (eg. Use 1 hour of readings for 1 day of propagation) [2,4]. Additionally, the oscillatory fluctuations of error has a period of 1 orbit (~90 minutes for Low Earth Orbit), caused by non-symmetric spherical harmonics of Earth s gravitational field [4]. One of the HPOP ephemeris generated has Gaussian error added (Section 3). This is to simulate operational GPS data with it accuracy of about 100m (at high altitudes and velocity) [2, 4].Hence, using a single GPS state vector at t= 0 for forward propagation (using J2-6 propagator) will result in significant deviations from the actual orbit. To refine this state vector at t=0, orbit determination is conducted using subsequent state vectors at t=0, 1, 2 to obtain more accurate propagator data ( prop data). We investigate how the number of subsequent state vectors used affects the overall accuracy of propagated data, by considering the positional error at every time step between prop and GPS data (Eq. 8). To study the error of the overall package, the positional error over time is plotted between propagated data and true HPOP data for various corrected (with orbit determination) and uncorrected propagators (Fig. 5). The number of readings used to orbit determine is varied from 1 to 50, with 1 being a direct uncorrected propagation, at a frequency of 1 position vector per minute. The orbit is propagated for 1 day (1440 min). Fig.5. Total system error with and without orbit determination From the graph, orbit determination results in lower gradient of error, indicating more accurate orbit VI. CASE STUDY - I: CONJUNCTION ANALYSIS FOR COLLISION AVOIDANCE An important application of accurate orbit determination is in conjunction analysis: predicting possible collisions to take evasive action. To analyse satellites (especially foreign or defunct ones) and debris orbits, publicly available Two Line Element (TLE) data sets can be used [2, 6]. These data sets, published by NORAD, provide the 6 classical orbital elements for a satellite at regular intervals (~daily). In this study, TLE sets for both Iridium 33 (satellite 1 ) and Cosmos 2251 (satellite 2 ) were downloaded from CelesTrak [3] and propagated using the Python code developed to find the positions of each satellite in an Earth Centered Inertial (ECI) frame. The collision between Iridium 33 (US commercial) and Cosmos 2251 (Soviet military) on February 10, 2009, 1656 UTC time was catastrophic [6], generating thousands of space debris and highlighting the importance of space situational awareness [9]. The distance between the two satellites ( separation ) can be calculated by Eq. 9: Two satellites have the highest risk of collision at their closest approach (when separation is at a minimum) [6]. Using only one TLE data set results in less accurate propagations; the use of multiple data sets with orbit determination results in more accurate closest approach prediction. Using the J2-6 propagator only, when the last TLEs for each satellite (1 day prior) were used to propagate directly, the time of closest approach was 20:12:03 UTC (Fig. 6) and predicted separation was 12 km inbetween. The satellites are periodically closer and further apart due to their intersecting elliptical orbits around Earth. From the propagated data, there is a global minimum in separation between the two satellites as expected. However, note that the collision timing is inaccurate by 4 hours, suggesting significant error growth from initial TLE inaccuracy. 201
5 Fig.6. Closest approaches with last available TLE only Given this limitation, we conduct orbit determination with five prior TLE sets for both satellites to refine the accuracy of initial propagator input. This results in the precise time of closest approach at 16:55:40, accurate to the minute, albeit at a slightly larger separation of 9.8 km (Fig. 7), which is caused by the relative inaccuracy (about 1.5 km) and infrequency (about once daily) of TLE data sets [2]. The significant improvement in time suggests a much more reliable knowledge of position when orbit determination is conducted, as initial input is refined. This paper uses two photos, and the corresponding spacecraft nadir point (point on Earth directly beneath satellite), from Gateway to Astronaut Photography of Earth, a catalogue of NASA astronaut photography [5, 7]. One photo was the unknown target (Singapore, Table 1) and another as a known calibration (Miami, Table 1), taken 1 day apart during ISS Expedition 35. Using TLE data sets from CelesTrak from March 25-30, the motion of the ISS is propagated with orbit determination in the Python code developed. Using known times of image acquisition, ISS nadir point is calculated and corroborated with NASA data. As photos may be taken with the camera not pointing radially down, camera orientation is required to calculate the displacement of the centre point of the photo from the nadir point. Hence, with an accurate orbit determination and propagation system, operators can then predict the motion of satellites and subsequently geo-register ground positions of photos taken [4]. Results of the TLE-based study for ISS location are shown in Table 1, together with photograph details. Table 1 ISS position prediction over Singapore and actual value Fig.7. Improved closest approach with TLE orbit determination VII. CASE STUDY - II: GROUND IDENTIFICATION FOR SATELLITE TRACKING In addition to the ECI frame, the Earth Centered Earth Fixed (ECEF) frame is used to determine a satellite s position relative to a rotating Earth. Considering the rotation rate of Earth ( rad s [1, 2]), and the position of a known ground calibration point (Latitude: θ, Longitude: φ ), a satellite s position can then be converted from ECEF frame to a more familiar notation: longitude, latitude and altitude. Satellite position in ECI frame (x, y, z coordinates) is calculated at calibration ( cali ) and image acquisition ( targ ), and used to determine relative change in latitude and longitude (Eqs. 10 & 11). Latitude = sin sin + θ (10) Longitude = tan y x tan y x +φ ω (t t ) (11) The position of the ISS has been accurately predicted on 27 March 2013 over Singapore with orbit determination. The latitude and longitude values of ISS match completely. With further knowledge of camera orientation, the coordinates of the photo centre can be calculated towithin an error of 1km. This is more accurate because interpolation is conducted with TLE readings both before and after time of image acquisition. The motion of the ISS across Earth can be plotted in a groundtrack (Fig. 8) over a Mercator projection [2, 4]. 14 orbits are shown, from the time ISS passes over Singapore to Miami. The groundtrack shows the satellite s nadir over time as it orbits Earth. Each orbit (curve) is Eastwards of the previous orbit due to Earth rotation, and so the projections move rightwards [2]. These projections indicate the regions where image acquisition is possible. For the ISS, its orbital plane inclination bounds the latitude regions from 56.5 N to 56.5 S. 202
6 Fig.8. Groundtrack of ISS from 27 th to 28 th March 2013 CONCLUSION This study has fully developed and verified a reliable orbit determination Python code using the method of differential correction and a J2-6 gravitational perturbation propagator. The differential correction and propagation methods are collectively tested, producing a good agreement with simulated GPS data with 20 km positional difference after one day. Practical applications considered are orbital collision avoidance based on TLE data propagated and satellite image acquisition. The developed propagator algorithm is shown to predict accurately the real-life case studies of the 2009 Iridium-Cosmos satellite collision and ISS position relative to a rotating Earth. ACKNOWLEDGEMENTS The authors are grateful for the guidance by Dr Lou Kok Yong and Mr Kenny Chen from DSO National Laboratories during the study, and would like to acknowledge the support of Raffles Science Institute in providing the research attachment opportunity and support for attending the conference. REFERENCES [1] R. Bate, D. Mueller, J. White (1971) Fundamentals of Astrodynamics, Dover Publications Inc [2] D. A. Vallado (2007) Fundamentals of Astrodynamics and Applications, 3rd Edition, Springer [3] H. D. Curtis (2010) Orbital Mechanics for Engineering Students, Elsevier [4] O. Montenbruck (2000) Satellite Orbits: Models, Methods, and Applications, Springer [5] CelesTrak. Retrieved from [6] T. S. Kelso (2009) Analysis of the Iridium 33-Cosmos 2251 collision, Advances in the Astronautical Sciences, AAS [7] Gateway to Astronaut Photography of Earth, Earth Science and Remote Sensing Unit, NASA Johnson Space Center. Retrieved from [8] N. Ayat, M. Mehdipour (2006) Accurate Doppler Prediction Scheme for Satellite Orbits [9] D. A. Danielson, D. Canright, D. N. Perini, P. W. Schumacher, Jr. (1999) The Naval Space Command automatic Differential Correction process [10] Systems Toolkit (2015). Analytical Graphics Inc. Retrieved from 203
An Analysis of N-Body Trajectory Propagation. Senior Project. In Partial Fulfillment. of the Requirements for the Degree
An Analysis of N-Body Trajectory Propagation Senior Project In Partial Fulfillment of the Requirements for the Degree Bachelor of Science in Aerospace Engineering by Emerson Frees June, 2011 An Analysis
More informationSpace Surveillance with Star Trackers. Part II: Orbit Estimation
AAS -3 Space Surveillance with Star Trackers. Part II: Orbit Estimation Ossama Abdelkhalik, Daniele Mortari, and John L. Junkins Texas A&M University, College Station, Texas 7783-3 Abstract The problem
More informationACCURACY ASSESSMENT OF GEOSTATIONARY-EARTH-ORBIT WITH SIMPLIFIED PERTURBATIONS MODELS
ARTIFICIAL SATELLITES, Vol. 51, No. 2 2016 DOI: 10.1515/arsa-2016-0005 ACCURACY ASSESSMENT OF GEOSTATIONARY-EARTH-ORBIT WITH SIMPLIFIED PERTURBATIONS MODELS Lihua Ma, Xiaojun Xu, Feng Pang National Astronomical
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 3. The Orbit in Space Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation: Space We need means of describing orbits in three-dimensional space. Example: Earth s oblateness
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 3B. The Orbit in Space and Time Gaëtan Kerschen Space Structures & Systems Lab (S3L) Previous Lecture: The Orbit in Time 3.1 ORBITAL POSITION AS A FUNCTION OF TIME 3.1.1 Kepler
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5. Dominant Perturbations Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation Assumption of a two-body system in which the central body acts gravitationally as a point
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5. Dominant Perturbations Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation Assumption of a two-body system in which the central body acts gravitationally as a point
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5. Numerical Methods Gaëtan Kerschen Space Structures & Systems Lab (S3L) Why Different Propagators? Analytic propagation: Better understanding of the perturbing forces. Useful
More informationOptimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming
Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology
More informationOn Sun-Synchronous Orbits and Associated Constellations
On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843,
More informationAN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT
AAS 16-366 AN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT Jason A. Reiter * and David B. Spencer INTRODUCTION Collision avoidance maneuvers to prevent orbital
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationCOLLISION RISK ASSESSMENT AND MITIGATION STRATEGY FOR THE GSOC GEO SATELLITES
COLLISION RISK ASSESSMENT AND MITIGATION STRATEGY FOR THE GSOC GEO SATELLITES Saika Aida (1), Michael Kirschner (2), Florian Meissner (3) (1) DLR German Space Operations Center (GSOC), Münchner Str.20,
More informationEarth-Centered, Earth-Fixed Coordinate System
Fundamentals of Global Positioning System Receivers: A Software Approach James Bao-Yen Tsui Copyright 2000 John Wiley & Sons, Inc. Print ISBN 0-471-38154-3 Electronic ISBN 0-471-20054-9 CHAPTER FOUR Earth-Centered,
More informationExperimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations
Experimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations Aman Saluja #1, Manish Bansal #2, M Raja #3, Mohd Maaz #4 #Aerospace Department, University of Petroleum and Energy
More informationIncreasing the Accuracy of Orbital Position Information from NORAD SGP4 Using Intermittent GPS Readings
Increasing the Accuracy of Orbital Position Information from NORAD SGP4 Using Intermittent GPS Readings Michael Greene Robert E. Zee Space Flight Laboratory University of Toronto 12 August 2009 23 rd Annual
More informationOn Directional Measurement Representation in Orbit Determination
On Directional Measurement Representation in Orbit Determination Evan M. Ward and Greg Carbott U.S. Naval Research Lab, Washington, DC 20375 Precision Orbit Determination (OD) is often critical for successful
More informationMULTI PURPOSE MISSION ANALYSIS DEVELOPMENT FRAMEWORK MUPUMA
MULTI PURPOSE MISSION ANALYSIS DEVELOPMENT FRAMEWORK MUPUMA Felipe Jiménez (1), Francisco Javier Atapuerca (2), José María de Juana (3) (1) GMV AD., Isaac Newton 11, 28760 Tres Cantos, Spain, e-mail: fjimenez@gmv.com
More informationOrekit at the U.S. Naval Research Laboratory. Evan Ward
Orekit at the U.S. Naval Research Laboratory Evan Ward U.S. Naval Research Laboratory Astrodynamics and Navigation Section, Washington DC November 16, 2017 Outline Introduction Geolocation with Orekit
More informationAttitude Determination using Infrared Earth Horizon Sensors
SSC14-VIII-3 Attitude Determination using Infrared Earth Horizon Sensors Tam Nguyen Department of Aeronautics and Astronautics, Massachusetts Institute of Technology 77 Massachusetts Avenue, Cambridge,
More informationAccuracy Assessment of SGP4 Orbit Information Conversion into Osculating Elements
Accuracy Assessment of SGP4 Orbit Information Conversion into Osculating Elements Saika Aida (1), Michael Kirschner (2) (1) DLR German Space Operations Center (GSOC), Oberpfaffenhofen, 82234 Weßling, Germany,
More informationIMPROVED ORBITAL DEBRIS TRAJECTORY ESTIMATION BASED ON SEQUENTIAL TLE PROCESSING ABSTRACT
IAC-09.A6.2.9 IMPROVED ORBITAL DEBRIS TRAJECTORY ESTIMATION BASED ON SEQUENTIAL TLE PROCESSING Alana R. Muldoon, Gabriel H. Elkaim Department of Computer Engineering, University of California, Santa Cruz
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E 6 Part B, Lecture 6 17 March, 2017 C O N T E N T S In this lecture, we will look at various existing satellite tracking techniques. Recall that we
More informationCOMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING
COMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING Fran Martinez Fadrique, Alberto Águeda Maté, Joan Jorquera Grau, Jaime Fernández Sánchez, Laura Aivar García
More informationFundamentals of Astrodynamics and Applications
Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA
More informationInfrared Earth Horizon Sensors for CubeSat Attitude Determination
Infrared Earth Horizon Sensors for CubeSat Attitude Determination Tam Nguyen Department of Aeronautics and Astronautics Massachusetts Institute of Technology Outline Background and objectives Nadir vector
More informationAnalytical Method for Space Debris propagation under perturbations in the geostationary ring
Analytical Method for Space Debris propagation under perturbations in the geostationary ring July 21-23, 2016 Berlin, Germany 2nd International Conference and Exhibition on Satellite & Space Missions Daniel
More informationWhat Happened to BLITS? An Analysis of the 2013 Jan 22 Event
What Happened to BLITS? An Analysis of the 2013 Jan 22 Event T.S. Kelso Center for Space Standards & Innovation N.N. Parkhomenko, V.D. Shargorodsky, V.P. Vasiliev, V.S. Yurasov Open Joint-Stock Company
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 4B. Non-Keplerian Motion Gaëtan Kerschen Space Structures & Systems Lab (S3L) 2. Two-body problem 4.1 Dominant perturbations Orbital elements (a,e,i,ω,ω) are constant Real satellites
More informationPROBABILITY OF COLLISION WITH SPECIAL PERTURBATIONS DYNAMICS USING THE MONTE CARLO METHOD
AAS 11-435 PROBABILITY OF COLLISION WITH SPECIAL PERTURBATIONS DYNAMICS USING THE MONTE CARLO METHOD Chris Sabol, Christopher Binz and Alan Segerman, Kevin Roe, and Paul W. Schumacher, Jr. INTRODUCTION
More informationA Fast Method for Embattling Optimization of Ground-Based Radar Surveillance Network
A Fast Method for Embattling Optimization of Ground-Based Radar Surveillance Network JIANG Hai University of Chinese Academy of Sciences National Astronomical Observatories, Chinese Academy of Sciences
More informationcloseap: GMV S SOLUTION FOR COLLISION RISK ASSESSMENT
closeap: GMV S SOLUTION FOR COLLISION RISK ASSESSMENT D. Escobar Antón (1), A. Pérez Cambriles (1), F. M. Martínez Fadrique (1), A. Águeda Maté (1) (1) GMV S.A., Calle Isaac Newton 11, P.T.M. Tres Cantos,
More informationInfrared Earth Horizon Sensors for CubeSat Attitude Determination
Infrared Earth Horizon Sensors for CubeSat Attitude Determination Tam Nguyen Department of Aeronautics and Astronautics Massachusetts Institute of Technology Outline Background and objectives Nadir vector
More informationBenefits of a Geosynchronous Orbit (GEO) Observation Point for Orbit Determination
Benefits of a Geosynchronous Orbit (GEO) Observation Point for Orbit Determination Ray Byrne, Michael Griesmeyer, Ron Schmidt, Jeff Shaddix, and Dave Bodette Sandia National Laboratories ABSTRACT Determining
More informationCollision Risk Assessment for Spacecraft with CASS
Collision Risk Assessment for Spacecraft with CASS Ma Chaowei 2011, 17-18 Oct 2011 Beijing Institute of Tracking and Telecommunications Technology Ma Chaowei 2 Contents Introduction Basic Information of
More informationOrbit Propagatorr and Geomagnetic Field Estimator for NanoSatellite: The ICUBE Mission
Vol:7, No:7, 23 Orbit Propagatorr and Geomagnetic Field Estimator for NanoSatellite: The ICUBE Mission Lv Meibo, Naqvi Najam Abbas, Hina Arshad, and Li YanJun International Science Index, Physical and
More informationUSA Space Debris Environment, Operations, and Modeling Updates
USA Space Debris Environment, Operations, and Modeling Updates Presentation to the 51 st Session of the Scientific and Technical Subcommittee Committee on the Peaceful Uses of Outer Space United Nations
More informationorbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG
A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits
More informationAPPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING
APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types
More informationImproving LEO prediction precision with TLEs
Improving LEO prediction precision with TLEs LIU Wei,WANG Ronglan,YAN Ruidong,GONG Jiancun (1. Center for Space Science and Applied Research, CAS,Beijing 119, China) Abstract: TLE is the only publicly
More informationThird Body Perturbation
Third Body Perturbation p. 1/30 Third Body Perturbation Modeling the Space Environment Manuel Ruiz Delgado European Masters in Aeronautics and Space E.T.S.I. Aeronáuticos Universidad Politécnica de Madrid
More informationFigure 1. View of ALSAT-2A spacecraft
ALSAT-2A TRANSFER AND FIRST YEAR OPERATIONS M. Kameche (1), A.H. Gicquel (2), D. Joalland (3) (1) CTS/ASAL, 1 Avenue de la Palestine, BP 13, Arzew 31200 Oran, Algérie, email:mo_kameche@netcourrier.com
More informationSince the first orbital launch in 1957, the number of artificial objects in Earth orbit has been steadily increasing. This
SpaceOps Conferences 28 May - 1 June 2018, 2018, Marseille, France 2018 SpaceOps Conference 10.2514/6.2018-2720 Collision probability through time integration Implementation and operational results Vincent
More informationTHE INITIAL ORBIT DETERMINATION OF LEO SATELLITES OBSERVED NEAR THE LOCAL ZENITH
THE INITIAL ORBIT DETERMINATION OF LEO SATELLITES OBSERVED NEAR THE LOCAL ZENITH M.A. Earl (a) (a) Royal Military College of Canada, Kingston ON, earl-m@castor.ca Abstract Our large satellite population
More informationAIR FORCE INSTITUTE OF TECHNOLOGY
Space Based Satellite Tracking and Characterization Utilizing Non-Imaging Passive Sensors THESIS Bradley R. Townsend, Captain, USA AFIT/GA/ENY/08-M06 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE
More informationThe Space-based Telescopes for Actionable Refinement of Ephemeris (STARE) mission
LLNL-PRES-641541 Performance Measures x.x, x.x, and x.x SSC13-XI-11 The Space-based Telescopes for Actionable Refinement of Ephemeris (STARE) mission Vincent Riot, Willem de Vries, Lance Simms, Brian Bauman,
More informationFRAGMENTATION EVENT MODEL AND ASSESSMENT TOOL (FREMAT) SUPPORTING ON-ORBIT FRAGMENTATION ANALYSIS
FRAGMENTATION EVENT MODEL AND ASSESSMENT TOOL (FREMAT) SUPPORTING ON-ORBIT FRAGMENTATION ANALYSIS Roxana Larisa Andrișan (1), Alina Georgia Ioniță (1), Raúl Domínguez González (2), Noelia Sánchez Ortiz
More informationJimmy Dale Hicks Jr. Auburn, Alabama December 8, Keywords: orbit debris, state estimation, kalman filtering
Performance Comparison of an Extended Kalman Filter and an Iterated Extended Kalman Filter for Orbit Determination of Space Debris with Poor Apriori Information and Intermittent Observations by Jimmy Dale
More informationAIR FORCE INSTITUTE OF TECHNOLOGY
OPTIMAL ORBITAL COVERAGE OF THEATER OPERATIONS AND TARGETS THESIS Kimberly A. Sugrue, Captain, USAF AFIT/GA/ENY/07-M17 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE INSTITUTE OF TECHNOLOGY Wright-Patterson
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationPeriodicity Characterization of Orbital Prediction Error and Poisson Series Fitting
Periodicity Characterization of Orbital Prediction Error and Poisson Series Fitting Bai Xianzong *, Chen Lei, Tang Guoin College of Aerospace and Materials Engineering, National University of Defense Technology,
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5B. Orbital Maneuvers Gaëtan Kerschen Space Structures & Systems Lab (S3L) Previous Lecture: Coplanar Maneuvers 5.1 INTRODUCTION 5.1.1 Why? 5.1.2 How? 5.1.3 How much? 5.1.4 When?
More informationValidation of Astrodynamic Formation Flying Models Against SPACE-SI Experiments with Prisma Satellites
Validation of Astrodynamic Formation Flying Models Against SPACE-SI Experiments with Prisma Satellites Drago Matko, Tomaž Rodič, Sašo Blažič, Aleš Marsetič, Krištof Oštir, Gašper Mušič, Luka Teslić, Gregor
More informationCreating Satellite Orbits
Exercises using Satellite ToolKit (STK) vivarad@ait.ac.th Creating Satellite Orbits 1. What You Will Do Create a low-earth orbit (LEO) satellite Create a medium-earth orbit (MEO) satellite Create a highly
More informationAIR FORCE INSTITUTE OF TECHNOLOGY
INVESTIGATING ANALYTICAL AND NUMERICAL METHODS TO PREDICT SATELLITE ORBITS USING TWO-LINE ELEMENT SETS THESIS Adam T. Rich, Captain, USAF AFIT-ENY-MS-17-M-286 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY
More informationObservation of Light Curves of Space Objects. Hirohisa Kurosaki Japan Aerospace Exploration Agency Toshifumi Yanagisawa.
Observation of Light Curves of Space Objects Hirohisa Kurosaki Japan Aerospace Exploration Agency Toshifumi Yanagisawa Japan Aerospace Exploration Agency Atsushi Nakajima Japan Aerospace Exploration Agency
More informationOrbit Representation
7.1 Fundamentals 223 For this purpose, code-pseudorange and carrier observations are made of all visible satellites at all monitor stations. The data are corrected for ionospheric and tropospheric delays,
More informationResults and Analyses of Debris Tracking from Mt Stromlo
Results and Analyses of Debris Tracking from Mt Stromlo Jizhang Sang 1, Ian Ritchie, Matt Pearson, and Craig Smith EOS Space Systems Pty Ltd, Australia Abstract In the last 2 years, EOS Space Systems has
More informationORBITAL DECAY PREDICTION AND SPACE DEBRIS IMPACT ON NANO-SATELLITES
Journal of Science and Arts Year 16, No. 1(34), pp. 67-76, 2016 ORIGINAL PAPER ORBITAL DECAY PREDICTION AND SPACE DEBRIS IMPACT ON NANO-SATELLITES MOHAMMED CHESSAB MAHDI 1 Manuscript received: 22.02.2016;
More informationRELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS
(Preprint) AAS 12-202 RELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS Hemanshu Patel 1, T. Alan Lovell 2, Ryan Russell 3, Andrew Sinclair 4 "Relative navigation using
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More informationGPS Results for the Radio Aurora Explorer II CubeSat Mission
GPS Results for the Radio Aurora Explorer II CubeSat Mission Jessica Arlas, Sara Spangelo, This paper presents the performance of the Global Positioning System (GPS) subsystem for the Radio Aurora explorer
More informationOptimization of Eccentricity during a Two Burn Station Acquisition Sequence of a Geosynchronous Orbit
Optimization of Eccentricity during a Two Burn Station Acquisition Sequence of a Geosynchronous Orbit A project present to The Faculty of the Department of Aerospace Engineering San Jose State University
More informationA Comparison of Low Cost Transfer Orbits from GEO to LLO for a Lunar CubeSat Mission
A Comparison of Low Cost Transfer Orbits from GEO to LLO for a Lunar CubeSat Mission A presentation for the New Trends in Astrodynamics conference Michael Reardon 1, Jun Yu 2, and Carl Brandon 3 1 PhD
More informationDEFINITION OF A REFERENCE ORBIT FOR THE SKYBRIDGE CONSTELLATION SATELLITES
DEFINITION OF A REFERENCE ORBIT FOR THE SKYBRIDGE CONSTELLATION SATELLITES Pierre Rozanès (pierre.rozanes@cnes.fr), Pascal Brousse (pascal.brousse@cnes.fr), Sophie Geffroy (sophie.geffroy@cnes.fr) CNES,
More informationAlgorithms for inverting radio occultation signals in the ionosphere
Algorithms for inverting radio occultation signals in the ionosphere This document describes the algorithms for inverting ionospheric radio occultation data using the Fortran 77 code gmrion.f and related
More informationYahsat s Approach to GEO Collision Avoidance A Case Stduy with STTW-4
Yahsat s Approach to GEO Collision Avoidance A Case Stduy with STTW-4 John Baker Senior Manager - Flight Dynamics Yahsat Agenda Objective Collision Risk Know your enemy Situational Awareness Yahsat Stationkeeping
More informationCelestial Mechanics III. Time and reference frames Orbital elements Calculation of ephemerides Orbit determination
Celestial Mechanics III Time and reference frames Orbital elements Calculation of ephemerides Orbit determination Orbital position versus time: The choice of units Gravitational constant: SI units ([m],[kg],[s])
More informationA.I. Nazarenko 1
Model Study of the Possibilities of Space Debris Cataloging A.I. Nazarenko 1 anazarenko32@mail.ru The objective of this paper is the estimation of the characteristics of the measuring means and the software,
More informationHybrid (Ion and Chemical) GEO Stationkeeping Maneuver Planning Software
Hybrid (Ion and Chemical) GEO Stationkeeping Maneuver Planning Software J. K. Skipper, D. Racicot, S. Li, R. Provencher and J. Palimaka Telesat Canada, Ottawa, Ontario, Canada. Abstract In the geochronous
More informationIMPROVED CONJUNCTION ANALYSIS VIA COLLABORATIVE SPACE SITUATIONAL AWARENESS
IMPROVED CONJUNCTION ANALYSIS VIA COLLABORATIVE SPACE SITUATIONAL AWARENESS T.S. Kelso (1), D. Vallado (2), J. Chan (3), and B. Buckwalter (4) (1) Center for Space Standards & Innovation, 7150 Campus Drive,
More informationImprovement of Orbits of Geostationary Satellites from Observations Over a Time Interval of Days
Improvement of Orbits of Geostationary Satellites from Observations Over a Time Interval of 15-2 Days * Konstantin V.Grigoriev, ** Felix R. Hoots * Pulkovo Observatory, St-Petersburg, Russia ** GRC International,
More informationConjunction Risk Assessment and Avoidance Maneuver Planning Tools
Conjunction Risk Assessment and Avoidance Maneuver Planning Tools Saika Aida DLR German Space Operations Center (GSOC), Oberpfaffenhofen, 834 Weßling, Germany +49 8153 8-158, saika.aida@dlr.de ASTRACT
More informationSpacecraft De-Orbit Point Targeting using Aerodynamic Drag
AIAA SciTech Forum 9-13 January 2017, Grapevine, Texas AIAA Guidance, Navigation, and Control Conference AIAA 2017-1268 Spacecraft De-Orbit Point Targeting using Aerodynamic Drag Sanny R. Omar 1 and Riccardo
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) The Orbit in Space and Time Gaëtan Kerschen Space Structures & Systems Lab (S3L) Newton s laws F = ma F g mm = G r 1 2 uˆ 2 r Relative motion μ r = r 3 r Energy conserv. 2 v μ
More informationUnderstanding the Differences between LS Algorithms and Sequential Filters
Understanding the Differences between LS Algorithms and Sequential Filters In order to perform meaningful comparisons between outputs from a least squares (LS) orbit determination algorithm and orbit determination
More informationPrediction accuracy analysis from orbital elements generated for a new space object catalogue
Prediction accuracy analysis from orbital elements generated for a new space object catalogue James Bennett, Michael Lachut, Ben Greene SERC Limited, Mt Stromlo Observatory, Cotter Road, Weston Creek,
More informationCIRCULAR MOTION AND SHM : Higher Level Long Questions.
CIRCULAR MOTION AND SHM : Higher Level Long Questions. ***ALL QUESTIONS ARE HIGHER LEVEL**** Circular Motion 2012 Question 12 (a) (Higher Level ) An Olympic hammer thrower swings a mass of 7.26 kg at the
More informationCOVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE. European Space Surveillance Conference 7-9 June 2011
COVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE European Space Surveillance Conference 7-9 June 2011 Pablo García (1), Diego Escobar (1), Alberto Águeda (1), Francisco
More informationIAC-16.A Jason A. Reiter a *, David B. Spencer b
IAC-16.A6.7.5 Trading Spacecraft Propellant Use and Mission Performance to Determine the Optimal Collision Probability in Emergency Collision Avoidance Scenarios Jason A. Reiter a *, David B. Spencer b
More informationAstrodynamics 103 Part 1: Gauss/Laplace+ Algorithm
Astrodynamics 103 Part 1: Gauss/Laplace+ Algorithm unknown orbit Observer 1b Observer 1a r 3 object at t 3 x = r v object at t r 1 Angles-only Problem Given: 1a. Ground Observer coordinates: ( latitude,
More informationAnalysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints
Analysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints F.J. de Bruijn, E. Gill Abstract This paper investigates the design of the relative motion orbits of geostationary
More informationOssama Abdelkhalik and Daniele Mortari Department of Aerospace Engineering, Texas A&M University,College Station, TX 77843, USA
Two-Way Orbits Ossama Abdelkhalik and Daniele Mortari Department of Aerospace Engineering, Texas A&M University,College Station, TX 77843, USA November 17, 24 Abstract. This paper introduces a new set
More informationCALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION
ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University
More informationSolar Glint for Optical Ranging
Solar Glint for Optical Ranging Anja Brokaw, Senior Software Engineer Cosmic AES, anja.brokaw@cosmicaes.com Co-authors Adrian Johnson, Principal Engineer Cosmic AES, adrian.johnson@cosmicaes.com Phillip
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) L04: Non-Keplerian Motion Gaëtan Kerschen Space Structures & Systems Lab (S3L) Non-Keplerian Motion 4 Dominant Perturbations Analytic Treatment Numerical Methods Concluding Remarks
More informationAttitude Determination using Infrared Earth Horizon Sensors
Attitude Determination using Infrared Earth Horizon Sensors Tam N. T. Nguyen Department of Aeronautics and Astronautics Massachusetts Institute of Technology 28 th Annual AIAA/USU Conference on Small Satellites
More informationORBIT DETERMINATION ANALYSIS FOR A JOINT UK-AUSTRALIAN SPACE SURVEILLANCE EXPERIMENT
ORBIT DETERMINATION ANALYSIS FOR A JOINT UK-AUSTRALIAN SPACE SURVEILLANCE EXPERIMENT M. Rutten, * T. Bessell, * N. Gordon, *, * NSID, Defence Science and Technology Organisation, Australia School of ITEE,
More informationAccuracy of EOCFI Software functions in the Lib and Orbit libraries
estec European Space Research and Technology Centre Keplerlaan 1 2201 AZ Noordwijk The Netherlands T +31 (0)71 565 6565 F +31 (0)71 565 6040 www.esa.int Accuracy of EOCFI Software functions in the Lib
More informationRADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION
AAS 15-440 RADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION Austin B. Probe, * Brent Macomber,* Julie I. Read,* Robyn M. Woollands,* and John L. Junkins
More informationLaser de-spin maneuver for an active debris removal mission - a realistic scenario for Envisat
Laser de-spin maneuver for an active debris removal mission - a realistic scenario for Envisat Daniel Kucharski Space Environment Research Centre, Mt Stromlo Observatory, Weston Creek ACT, 2611, AUSTRALIA
More informationMAE 180A: Spacecraft Guidance I, Summer 2009 Homework 2 Due Tuesday, July 14, in class.
MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 2 Due Tuesday, July 14, in class. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details
More informationHomework #1 Solution: March 8, 2006
12.540 Homework #1 Solution: March 8, 2006 Question 1: (a) Convert geodetic position 290 deg Long 42 deg latitude ellipsoidal height 0 m into Cartesian and geocentric coordinates. (b) How far apart on
More informationMEETING ORBIT DETERMINATION REQUIREMENTS FOR A SMALL SATELLITE MISSION
MEETING ORBIT DETERMINATION REQUIREMENTS FOR A SMALL SATELLITE MISSION Adonis Pimienta-Peñalver, Richard Linares, and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY,
More informationOrbital Data Applications for Space Objects
Orbital Data Applications for Space Objects Lei Chen Xian-Zong Bai Yan-Gang Liang Ke-Bo Li Orbital Data Applications for Space Objects Conjunction Assessment and Situation Analysis 123 Lei Chen National
More informationOrbital Resonances in the Vinti Solution
Air Force Institute of Technology AFIT Scholar Theses and Dissertations 3-24-2016 Orbital Resonances in the Vinti Solution Laura M. Duffy Follow this and additional works at: https://scholar.afit.edu/etd
More informationThe Simplest Tether Control Law in a Small Satellite
The Simplest Tether Control Law in a Small Satellite Yosuke Nakamura (Kyushu University) Designing of a small satellite based on the Satellite Design Contest is progressing in Kyushu University. This is
More informationA TLE-BASED REPRESENTATION OF PRECISE ORBIT PREDICTION RESULTS. Jizhang Sang and Bin Li. School of Geodesy and Geomatics, Wuhan University
A TLE-BASED REPRESENTATION OF PRECISE ORBIT PREDICTION RESULTS Jizhang Sang and Bin Li School of Geodesy and Geomatics, Wuhan University ABSTRACT The computing efficiency of the analytic SGP4 algorithm
More informationSpacecraft Orbit Anomaly Representation Using Thrust-Fourier-Coefficients with Orbit Determination Toolbox
Spacecraft Orbit Anomaly Representation Using Thrust-Fourier-Coefficients with Orbit Determination Toolbox Hyun Chul Ko and Daniel J. Scheeres University of Colorado - Boulder, Boulder, CO, USA ABSTRACT
More informationOptElec: an Optimisation Software for Low-Thrust Orbit Transfer Including Satellite and Operation Constraints
OptElec: an Optimisation Software for Low-Thrust Orbit Transfer Including Satellite and Operation Constraints 7th International Conference on Astrodynamics Tools and Techniques, DLR, Oberpfaffenhofen Nov
More information