CHARACTERISTICS and DESIGN STRATEGIES for NEAR RECTILINEAR HALO ORBITS WITHIN the EARTH-MOON SYSTEM a Thesis Submitted to the Fa
Total Page:16
File Type:pdf, Size:1020Kb
CHARACTERISTICS AND DESIGN STRATEGIES FOR NEAR RECTILINEAR HALO ORBITS WITHIN THE EARTH-MOON SYSTEM A Thesis Submitted to the Faculty of Purdue University by Emily M. Zimovan In Partial Fulfillment of the Requirements for the Degree of Master of Science in Aeronautics and Astronautics August 2017 Purdue University West Lafayette, Indiana ii iii THE PURDUE UNIVERSITY GRADUATE SCHOOL STATEMENT OF THESIS APPROVAL Professor Kathleen C. Howell, Chair School of Aeronautics and Astronautics Professor James M. Longuski School of Aeronautics and Astronautics Professor Carolin E. Frueh School of Aeronautics and Astronautics Approved by: Professor Weinong Wayne Chen Aeronautics and Astronautics Associate Head for Graduate Education iv To the women in science, technology, engineering, and mathematics that have come before me, and to those that will come after me. v ACKNOWLEDGMENTS First, a huge thank you for my parents. You instilled in me from a young age an invaluable love of math and science. Thank you for not just pushing me to question the world around me, but to find the answers. Without your encouragement and support, I would not be where I am today. From the bottom of my heart, I want to thank my fianc´e,Chris. When it felt like no one else in the world could understand where I was coming from or what I was going through, I knew that you actually did. I truly couldn't have done this without you. Your constant love and support mean the world to me. I'm grateful for you in so many ways. I want to extend my immense gratitude to my advisor, Professor Kathleen C. Howell. Thank you for giving me so many opportunities throughout my time in the research group. Your thoughts, encouragement, advice, and support have been invaluable to my research and time at Purdue. The members of my research group deserve a thank you, as well. You've listened to me present countless times over the past two years and have given me advice and thoughts that have led my research to where it is now. Thanks for helping me to learn all that I have. I would also like to thank my committee members, Professors James M. Longuski and Carolin Frueh. I appreciate the time you took to review my thesis and provide me with advice and feedback. Additionally, for helping fund my studies at both Purdue University and NASA Johnson Space Center, I would like to thank the School of Aeronautics and Astronautics, the College of Engineering, and the Purdue University vi Graduate School for awarding me with the Charles C. Chappelle Fellowship and the Ross Fellowship. Also, support provided by the NASA JSC NNX13AK60A grant is greatly appreciated. I am also grateful for the use of the outstanding facilities provided by the Rune and Barbara Eliasen Visualization Laboratory. vii TABLE OF CONTENTS Page LIST OF TABLES :::::::::::::::::::::::::::::::::: x LIST OF FIGURES ::::::::::::::::::::::::::::::::: xi ABSTRACT ::::::::::::::::::::::::::::::::::::: xvi 1 INTRODUCTION :::::::::::::::::::::::::::::::: 1 1.1 Previous Contributions ::::::::::::::::::::::::::: 2 1.1.1 Multi-Body Dynamics ::::::::::::::::::::::: 2 1.1.2 Differential Corrections ::::::::::::::::::::::: 3 1.1.3 Periodic Orbits ::::::::::::::::::::::::::: 4 1.2 Current Work :::::::::::::::::::::::::::::::: 6 2 DYNAMICAL MODELS ::::::::::::::::::::::::::::: 9 2.1 The Circular Restricted Three-Body Model :::::::::::::::: 10 2.1.1 Assumptions for the Circular Restricted Three-Body Model :: 10 2.1.2 Equations of Motion in the Circular Restricted Three-Body Model11 2.1.3 Integrals of the Motion ::::::::::::::::::::::: 17 2.2 The N-Body Ephemeris Model :::::::::::::::::::::: 25 2.2.1 Assumptions for the N-Body Ephemeris Model ::::::::: 25 2.2.2 Equations of Motion in the N-Body Ephemeris Model :::::: 27 viii Page 3 DIFFERENTIAL CORRECTIONS ::::::::::::::::::::::: 31 3.1 State Transition Matrix :::::::::::::::::::::::::: 31 3.1.1 STM for the Circular Restricted Three-Body Model ::::::: 33 3.1.2 STM for the N-Body Ephemeris Model :::::::::::::: 35 3.1.3 Monodromy Matrix ::::::::::::::::::::::::: 37 3.2 Differential Corrections Strategies - Shooting Methods ::::::::: 38 3.2.1 Simple Targeting Scheme Formulation :::::::::::::: 40 3.2.2 Multi-Segment Corrections Algorithm Formulation ::::::: 43 3.2.3 Additional Constraints ::::::::::::::::::::::: 46 3.3 Coordinate Frames ::::::::::::::::::::::::::::: 52 3.3.1 Coordinate Frame Transformation: Rotating Frame to Arbi- trary Inertial Frame :::::::::::::::::::::::: 52 3.3.2 Coordinate Frame Transformation: Rotating Frame to J2000 Inertial Frame :::::::::::::::::::::::::::: 55 4 TRAJECTORY CONSTRUCTION: CR3BP AND THE EPHEMERIS MODEL61 4.1 Periodic Orbit Computation in the CR3BP :::::::::::::::: 61 4.1.1 Continuation of an Orbit Family ::::::::::::::::: 68 4.1.2 Examples of Periodic Orbit Families ::::::::::::::: 71 4.2 Poincar´eMapping and Surfaces of Section :::::::::::::::: 84 4.3 Computation of Periodic Orbit Invariant Manifolds ::::::::::: 85 4.4 Transitioning CR3BP Solutions into the Ephemeris Model ::::::: 87 4.4.1 Apse Angle ::::::::::::::::::::::::::::: 88 5 MISSION APPLICATIONS: ORBITAL CHARACTERISTICS AND TRANS- FER DESIGN ::::::::::::::::::::::::::::::::::: 93 ix Page 5.1 Stability ::::::::::::::::::::::::::::::::::: 94 5.1.1 CR3BP: Stability Index :::::::::::::::::::::: 96 5.2 Eclipsing :::::::::::::::::::::::::::::::::: 104 5.2.1 Resonance ::::::::::::::::::::::::::::: 105 5.2.2 Eclipse Duration ::::::::::::::::::::::::: 112 5.3 Transfers :::::::::::::::::::::::::::::::::: 121 5.3.1 Accessing NRHOs from Low Earth Orbit :::::::::::: 121 5.3.2 Accessing Distant Retrograde Orbits from NRHOs :::::: 126 6 CONCLUDING REMARKS :::::::::::::::::::::::::: 145 6.1 Summary of the Present Investigation :::::::::::::::::: 145 6.2 Concluding Remarks ::::::::::::::::::::::::::: 146 6.3 Recommendations for Future Work ::::::::::::::::::: 148 REFERENCES ::::::::::::::::::::::::::::::::::: 149 x LIST OF TABLES Table Page 2.1 Characteristic Quantities for Nondimensionalization ::::::::::::: 15 2.2 Rotating Coordinates of the Libration Points in the Earth-Moon and Sun- Earth Systems ::::::::::::::::::::::::::::::::: 22 2.3 Value of the Jacobi Constant for each Libration Point in the Earth-Moon and Sun-Earth Systems :::::::::::::::::::::::::::: 22 4.1 Initial Conditions for the Distant Retrograde Orbits ::::::::::::: 73 4.2 Initial Conditions for the Period-3 Distant Retrograde Orbits :::::::: 76 4.3 Initial Conditions for the L1 Halo Orbits ::::::::::::::::::: 79 4.4 Initial Conditions for the L2 Halo Orbits ::::::::::::::::::: 82 5.1 Shadow events for a January 25, 2025 apoapsis insertion into a 4:1 synodic resonant NRHO. ::::::::::::::::::::::::::::::: 117 5.2 Shadow events for a January 25, 2025 apoapsis insertion into a 9:2 synodic resonant NRHO. ::::::::::::::::::::::::::::::: 118 5.3 Shadow events for a January 1, 2025 apoapsis insertion into a 9:2 synodic resonant NRHO. ::::::::::::::::::::::::::::::: 118 xi LIST OF FIGURES Figure Page 2.1 Geometry of a Three-Body System Defined in an Inertial Frame :::::: 12 2.2 Inertial and Rotating Coordinate Frames Definitions :::::::::::: 13 2.3 Location of the Libration Points as Viewed in the Rotating Frame ::::: 21 2.4 The ZVS and corresponding ZVC for JCL3 < JC < JCL2. The forbidden regions are shaded in grey. ::::::::::::::::::::::::::: 23 2.5 The ZVCs for decreasing values of Jacobi constant in the Earth-Moon three-body system. Earth and Moon are plotted twice their actual size. The forbidden regions are shaded in grey. ::::::::::::::::::: 26 2.6 Inertial Frame Geometry of an N-Body System ::::::::::::::: 28 3.1 Single Shooter Corrections Algorithm Schematic ::::::::::::::: 41 3.2 Multiple Shooter Corrections Algorithm Schematic ::::::::::::: 43 4.1 A distant retrograde orbit about the Moon in the Earth-Moon rotating frame, computed using a single shooting algorithm. ::::::::::::: 66 4.2 An L2 southern halo orbit computed using a single shooting method. ::: 67 4.3 Family of distant retrograde orbits about the Moon computed in the Earth- Moon system, colored by Jacobi constant value. ::::::::::::::: 74 4.4 Jacobi constant value and orbital period as a function of perilune radius over a range of periodic orbits in the distant retrograde orbit family in the Earth-Moon system. :::::::::::::::::::::::::::::: 74 xii Figure Page 4.5 Family of period-3 distant retrograde orbits computed in the Earth-Moon system, colored by Jacobi constant value. ::::::::::::::::::: 77 4.6 Jacobi constant value and orbital period as a function of maximum ex- cursion in thex ^ direction for representative orbits in the period-3 distant retrograde periodic orbit family in the Earth-Moon system. ::::::::: 77 4.7 Family of L1 halo orbits computed in the Earth-Moon system, colored by Jacobi constant value. ::::::::::::::::::::::::::::: 80 4.8 Jacobi constant value and orbital period as a function of perilune radius for the periodic L1 halo orbit family in the Earth-Moon system. :::::: 80 4.9 Family of L2 halo orbits computed in the Earth-Moon system, colored by Jacobi constant value. ::::::::::::::::::::::::::::: 83 4.10 Jacobi constant value and orbital period as a function of perilune radius for the periodic L2 halo orbit family in the Earth-Moon system. :::::: 83 4.11 Poincar´emap schematic, adapted from Perko. [46] :::::::::::::: 85 4.12 Transition of four revolutions of a CR3BP L2 halo orbit (red dashed) to the EMSJ ephemeris model (blue solid). ::::::::::::::::::: 89 4.13 An ephemeris trajectory with a wider spread relative to a corresponding CR3BP orbit. ::::::::::::::::::::::::::::::::::