University of Cincinnati

Total Page:16

File Type:pdf, Size:1020Kb

University of Cincinnati UNIVERSITY OF CINCINNATI _____________ , 20 _____ I,______________________________________________, hereby submit this as part of the requirements for the degree of: ________________________________________________ in: ________________________________________________ It is entitled: ________________________________________________ ________________________________________________ ________________________________________________ ________________________________________________ Approved by: ________________________ ________________________ ________________________ ________________________ ________________________ Use of Near-Frozen Orbits for Satellite Formation Flying A thesis submitted to the Division of Research and Advanced Studies of the University of Cincinnati in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE (M.S.) in the Department of Aerospace Engineering and Engineering Mechanics of the College of Engineering 2001 by Heidi L. Davidz B.S., The Ohio State University, 1997 Committee Chair: Dr. Trevor Williams Abstract There is growing interest in flying coordinated clusters of small spacecraft to perform missions once accomplished by single, larger spacecraft. Using these satellite clusters reduces cost, improves survivability, and increases the flexibility of the mission. One challenge in implementing these satellite clusters is maintaining the formation as it experiences orbital perturbations, most notably due to the non-spherical Earth. Certain aspects of the orbital geometry can remain virtually fixed over extended periods of time due to a natural phenomenon called a frozen orbit. Specifically, the elements of the orbital geometry that can remain fixed are the argument of perigee (a measure of where the orbit is closest to the Earth) and the eccentricity (a measure of how circular or elliptical the orbit is). For satellite formations, using this frozen orbit phenomenon results in considerable propellant savings. In this study, a discussion of current literature on this topic is given. Some examples of formations at near-frozen conditions are shown. There is also a discussion of the propellant impact of near-frozen conditions. If two orbits meet a certain set of initial conditions, the orbits will naturally stay in the vicinity of each other. These orbits are sometimes called “Hill’s orbits”. An algorithm is developed here that determines if the Hill’s orbit conditions will be met given the initial differences in eccentricity and argument of perigee for two satellites. Acknowledgments I would like to thank Professor Trevor Williams for his guidance in conducting this research. Also, I would like to thank my husband Michael Cluff for his patience and support as I completed this work. Table of Contents 1 Introduction ............................................................................................................. 9 1.1 Introduction to Satellite Formation Flight.........................................................9 1.2 Propulsion for Formation Flight......................................................................11 1.3 Proposed Work................................................................................................ 13 2 Satellite Formation Flight Dynamics .................................................................... 13 2.1 Hill’s Orbits..................................................................................................... 13 2.2 Perturbations.................................................................................................... 17 3 Frozen Orbits......................................................................................................... 23 3.1 Definition of Frozen Orbit............................................................................... 23 3.2 Equation Development.................................................................................... 27 3.3 Experience with Frozen Orbits........................................................................ 28 4 Near-Frozen Orbit Test Cases............................................................................... 32 4.1 Introduction ..................................................................................................... 32 4.2 Development of Parameter Values.................................................................. 33 4.3 Development of Difference Equations............................................................ 34 4.4 Simulation Results...........................................................................................36 4.5 Case 4 Motion ................................................................................................. 49 5 Application of Cases to the Clohessy-Wiltshire/Hill’s Equations ........................ 50 5.1 Introduction ..................................................................................................... 50 5.2 Review of the CW/Hill’s Equations................................................................ 50 5.3 Coordinate Systems......................................................................................... 51 5.4 “Algorithm 6” from Vallado and McClain ..................................................... 51 Page 1 5.5 Application of “Algorithm 6” to Find PQW Radius and Velocity Vectors.... 52 5.6 Find IJK Radius and Velocity Vectors............................................................ 54 5.7 Find RSW Radius and Velocity Vectors......................................................... 54 5.8 Identify the CW/Hill’s Initial Conditions .......................................................56 5.9 Numerical Example......................................................................................... 57 6 Propellant Calculations......................................................................................... 64 6.1 Derivation of the Tangential Thrust ∆V Equation .......................................... 65 6.2 Propellant Calculation Results ........................................................................70 7 Summary ............................................................................................................... 78 8 Suggestions for Future Work ................................................................................ 79 9 Bibliography.......................................................................................................... 80 10 Appendix ............................................................................................................... 83 10.1 Matlab Code to Generate Frozen Orbit Circulations .................................... 83 10.2 Matlab Code for Running Cases, Pseudo State Space Program ................... 84 10.3 Matlab Code for Running Cases, ODE Program .......................................... 85 10.4 Matlab Code to Generate Propellant Calculation Graphs ............................. 86 Page 2 Table of Figures Figure 1: Zonal Harmonics ............................................................................................... 17 Figure 2: Sectorial Harmonics........................................................................................... 17 Figure 3: Tesseral Harmonics ........................................................................................... 17 Figure 4: Circulations Around the Frozen Orbit Conditions............................................ 26 Figure 5: Case 1 Eccentricity ............................................................................................ 37 Figure 6: Case 2 Eccentricity ............................................................................................ 37 Figure 7: Case 3 Eccentricity ............................................................................................ 38 Figure 8: Case 4 Eccentricity ............................................................................................ 38 Figure 9: Case 1 Argument of Perigee.............................................................................. 39 Figure 10: Case 2 Argument of Perigee............................................................................ 40 Figure 11: Case 3 Argument of Perigee............................................................................ 40 Figure 12: Case 4 Argument of Perigee............................................................................ 41 Figure 13: Case 1 Eccentricity vs. Argument of Perigee .................................................. 42 Figure 14: Case 2 Eccentricity vs. Argument of Perigee .................................................. 42 Figure 15: Case 3 Eccentricity vs. Argument of Perigee .................................................. 43 Figure 16: Case 4 Eccentricity vs. Argument of Perigee .................................................. 43 Figure 17: Case 1 Eccentricity Difference ........................................................................ 44 Figure 18: Case 2 Eccentricity Difference ........................................................................ 45 Figure 19: Case 3 Eccentricity Difference ........................................................................ 45 Figure 20: Case 4 Eccentricity Difference ........................................................................ 46 Figure 21: Case 1 Argument of Perigee Difference.......................................................... 47 Figure 22: Case 2 Argument of Perigee Difference.......................................................... 47 Page 3 Figure 23: Case 3 Argument of
Recommended publications
  • Analysis of Perturbations and Station-Keeping Requirements in Highly-Inclined Geosynchronous Orbits
    ANALYSIS OF PERTURBATIONS AND STATION-KEEPING REQUIREMENTS IN HIGHLY-INCLINED GEOSYNCHRONOUS ORBITS Elena Fantino(1), Roberto Flores(2), Alessio Di Salvo(3), and Marilena Di Carlo(4) (1)Space Studies Institute of Catalonia (IEEC), Polytechnic University of Catalonia (UPC), E.T.S.E.I.A.T., Colom 11, 08222 Terrassa (Spain), [email protected] (2)International Center for Numerical Methods in Engineering (CIMNE), Polytechnic University of Catalonia (UPC), Building C1, Campus Norte, UPC, Gran Capitan,´ s/n, 08034 Barcelona (Spain) (3)NEXT Ingegneria dei Sistemi S.p.A., Space Innovation System Unit, Via A. Noale 345/b, 00155 Roma (Italy), [email protected] (4)Department of Mechanical and Aerospace Engineering, University of Strathclyde, 75 Montrose Street, Glasgow G1 1XJ (United Kingdom), [email protected] Abstract: There is a demand for communications services at high latitudes that is not well served by conventional geostationary satellites. Alternatives using low-altitude orbits require too large constellations. Other options are the Molniya and Tundra families (critically-inclined, eccentric orbits with the apogee at high latitudes). In this work we have considered derivatives of the Tundra type with different inclinations and eccentricities. By means of a high-precision model of the terrestrial gravity field and the most relevant environmental perturbations, we have studied the evolution of these orbits during a period of two years. The effects of the different perturbations on the constellation ground track (which is more important for coverage than the orbital elements themselves) have been identified. We show that, in order to maintain the ground track unchanged, the most important parameters are the orbital period and the argument of the perigee.
    [Show full text]
  • Astrodynamics
    Politecnico di Torino SEEDS SpacE Exploration and Development Systems Astrodynamics II Edition 2006 - 07 - Ver. 2.0.1 Author: Guido Colasurdo Dipartimento di Energetica Teacher: Giulio Avanzini Dipartimento di Ingegneria Aeronautica e Spaziale e-mail: [email protected] Contents 1 Two–Body Orbital Mechanics 1 1.1 BirthofAstrodynamics: Kepler’sLaws. ......... 1 1.2 Newton’sLawsofMotion ............................ ... 2 1.3 Newton’s Law of Universal Gravitation . ......... 3 1.4 The n–BodyProblem ................................. 4 1.5 Equation of Motion in the Two-Body Problem . ....... 5 1.6 PotentialEnergy ................................. ... 6 1.7 ConstantsoftheMotion . .. .. .. .. .. .. .. .. .... 7 1.8 TrajectoryEquation .............................. .... 8 1.9 ConicSections ................................... 8 1.10 Relating Energy and Semi-major Axis . ........ 9 2 Two-Dimensional Analysis of Motion 11 2.1 ReferenceFrames................................. 11 2.2 Velocity and acceleration components . ......... 12 2.3 First-Order Scalar Equations of Motion . ......... 12 2.4 PerifocalReferenceFrame . ...... 13 2.5 FlightPathAngle ................................. 14 2.6 EllipticalOrbits................................ ..... 15 2.6.1 Geometry of an Elliptical Orbit . ..... 15 2.6.2 Period of an Elliptical Orbit . ..... 16 2.7 Time–of–Flight on the Elliptical Orbit . .......... 16 2.8 Extensiontohyperbolaandparabola. ........ 18 2.9 Circular and Escape Velocity, Hyperbolic Excess Speed . .............. 18 2.10 CosmicVelocities
    [Show full text]
  • Orbit Options for an Orion-Class Spacecraft Mission to a Near-Earth Object
    Orbit Options for an Orion-Class Spacecraft Mission to a Near-Earth Object by Nathan C. Shupe B.A., Swarthmore College, 2005 A thesis submitted to the Faculty of the Graduate School of the University of Colorado in partial fulfillment of the requirements for the degree of Master of Science Department of Aerospace Engineering Sciences 2010 This thesis entitled: Orbit Options for an Orion-Class Spacecraft Mission to a Near-Earth Object written by Nathan C. Shupe has been approved for the Department of Aerospace Engineering Sciences Daniel Scheeres Prof. George Born Assoc. Prof. Hanspeter Schaub Date The final copy of this thesis has been examined by the signatories, and we find that both the content and the form meet acceptable presentation standards of scholarly work in the above mentioned discipline. iii Shupe, Nathan C. (M.S., Aerospace Engineering Sciences) Orbit Options for an Orion-Class Spacecraft Mission to a Near-Earth Object Thesis directed by Prof. Daniel Scheeres Based on the recommendations of the Augustine Commission, President Obama has pro- posed a vision for U.S. human spaceflight in the post-Shuttle era which includes a manned mission to a Near-Earth Object (NEO). A 2006-2007 study commissioned by the Constellation Program Advanced Projects Office investigated the feasibility of sending a crewed Orion spacecraft to a NEO using different combinations of elements from the latest launch system architecture at that time. The study found a number of suitable mission targets in the database of known NEOs, and pre- dicted that the number of candidate NEOs will continue to increase as more advanced observatories come online and execute more detailed surveys of the NEO population.
    [Show full text]
  • Arxiv:1505.07033V2 [Astro-Ph.IM] 4 Sep 2015
    September 7, 2015 0:46 manuscript Journal of Astronomical Instrumentation c World Scientific Publishing Company A CUBESAT FOR CALIBRATING GROUND-BASED AND SUB-ORBITAL MILLIMETER-WAVE POLARIMETERS (CALSAT) Bradley R. Johnson1, Clement J. Vourch2, Timothy D. Drysdale2, Andrew Kalman3, Steve Fujikawa4, Brian Keating5 and Jon Kaufman5 1Department of Physics, Columbia University, New York, NY 10027, USA 2School of Engineering, University of Glasgow, Glasgow, Scotland G12 8QQ, UK 3Pumpkin, Inc., San Francisco, CA 94112, USA 4Maryland Aerospace Inc., Crofton, MD 21114, USA 5Department of Physics, University of California, San Diego, CA 92093-0424, USA Received (to be inserted by publisher); Revised (to be inserted by publisher); Accepted (to be inserted by publisher); We describe a low-cost, open-access, CubeSat-based calibration instrument that is designed to support ground- based and sub-orbital experiments searching for various polarization signals in the cosmic microwave background (CMB). All modern CMB polarization experiments require a robust calibration program that will allow the effects of instrument-induced signals to be mitigated during data analysis. A bright, compact, and linearly polarized astrophysical source with polarization properties known to adequate precision does not exist. Therefore, we designed a space-based millimeter-wave calibration instrument, called CalSat, to serve as an open-access calibrator, and this paper describes the results of our design study. The calibration source on board CalSat is composed of five \tones" with one each at 47.1, 80.0, 140, 249 and 309 GHz. The five tones we chose are well matched to (i) the observation windows in the atmospheric transmittance spectra, (ii) the spectral bands commonly used in polarimeters by the CMB community, and (iii) The Amateur Satellite Service bands in the Table of Frequency Allocations used by the Federal Communications Commission.
    [Show full text]
  • Up, Up, and Away by James J
    www.astrosociety.org/uitc No. 34 - Spring 1996 © 1996, Astronomical Society of the Pacific, 390 Ashton Avenue, San Francisco, CA 94112. Up, Up, and Away by James J. Secosky, Bloomfield Central School and George Musser, Astronomical Society of the Pacific Want to take a tour of space? Then just flip around the channels on cable TV. Weather Channel forecasts, CNN newscasts, ESPN sportscasts: They all depend on satellites in Earth orbit. Or call your friends on Mauritius, Madagascar, or Maui: A satellite will relay your voice. Worried about the ozone hole over Antarctica or mass graves in Bosnia? Orbital outposts are keeping watch. The challenge these days is finding something that doesn't involve satellites in one way or other. And satellites are just one perk of the Space Age. Farther afield, robotic space probes have examined all the planets except Pluto, leading to a revolution in the Earth sciences -- from studies of plate tectonics to models of global warming -- now that scientists can compare our world to its planetary siblings. Over 300 people from 26 countries have gone into space, including the 24 astronauts who went on or near the Moon. Who knows how many will go in the next hundred years? In short, space travel has become a part of our lives. But what goes on behind the scenes? It turns out that satellites and spaceships depend on some of the most basic concepts of physics. So space travel isn't just fun to think about; it is a firm grounding in many of the principles that govern our world and our universe.
    [Show full text]
  • Mission Design for the Lunar Reconnaissance Orbiter
    AAS 07-057 Mission Design for the Lunar Reconnaissance Orbiter Mark Beckman Goddard Space Flight Center, Code 595 29th ANNUAL AAS GUIDANCE AND CONTROL CONFERENCE February 4-8, 2006 Sponsored by Breckenridge, Colorado Rocky Mountain Section AAS Publications Office, P.O. Box 28130 - San Diego, California 92198 AAS-07-057 MISSION DESIGN FOR THE LUNAR RECONNAISSANCE ORBITER † Mark Beckman The Lunar Reconnaissance Orbiter (LRO) will be the first mission under NASA’s Vision for Space Exploration. LRO will fly in a low 50 km mean altitude lunar polar orbit. LRO will utilize a direct minimum energy lunar transfer and have a launch window of three days every two weeks. The launch window is defined by lunar orbit beta angle at times of extreme lighting conditions. This paper will define the LRO launch window and the science and engineering constraints that drive it. After lunar orbit insertion, LRO will be placed into a commissioning orbit for up to 60 days. This commissioning orbit will be a low altitude quasi-frozen orbit that minimizes stationkeeping costs during commissioning phase. LRO will use a repeating stationkeeping cycle with a pair of maneuvers every lunar sidereal period. The stationkeeping algorithm will bound LRO altitude, maintain ground station contact during maneuvers, and equally distribute periselene between northern and southern hemispheres. Orbit determination for LRO will be at the 50 m level with updated lunar gravity models. This paper will address the quasi-frozen orbit design, stationkeeping algorithms and low lunar orbit determination. INTRODUCTION The Lunar Reconnaissance Orbiter (LRO) is the first of the Lunar Precursor Robotic Program’s (LPRP) missions to the moon.
    [Show full text]
  • Small Satellites in Inclined Orbits to Increase Observation Capability Feasibility Analysis
    International Journal of Pure and Applied Mathematics Volume 118 No. 17 2018, 273-287 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu Special Issue ijpam.eu Small Satellites in Inclined Orbits to Increase Observation Capability Feasibility Analysis DVA Raghava Murthy1, V Kesava Raju2, T.Ramanjappa3, Ritu Karidhal4, Vijayasree Mallikarjuna Kande5, G Ravi chandra Babu6 and A Arunachalam7 1 7 Earth Observations System, ISRO Headquarters, Bengaluru, Karnataka India [email protected] 2 4 5 6ISRO Satellite Centre, Bengaluru, Karnataka, India 3SK University, Ananthapur, Andhra Pradesh, India January 6, 2018 Abstract Over the period of past four decades, Remote Sensing data products and services have been effectively utilized to showcase varieties of applications in many areas of re- sources inventory and monitoring. There are several satel- lite systems in operation today and they operate from Po- lar Orbit and Geosynchronous Orbit, collect the imagery and non-imagery data and provide them to user community for various applications in the areas of natural resources management, urban planning and infrastructure develop- ment, weather forecasting and disaster management sup- port. Quality of information derived from Remote Sensing imagery are strongly influenced by spatial, spectral, radio- metric & temporal resolutions as well as by angular & po- larimetric signatures. As per the conventional approach of 1 273 International Journal of Pure and Applied Mathematics Special Issue having Remote Sensing satellites in near Polar Sun syn- chronous orbit, the temporal resolution, i.e. the frequency with which an area can be frequently observed, is basically defined by the swath of the sensor and distance between the paths.
    [Show full text]
  • Five Special Types of Orbits Around Mars
    JOURNAL OF GUIDANCE,CONTROL, AND DYNAMICS Vol. 33, No. 4, July–August 2010 Engineering Notes Five Special Types of Orbits sensing satellites. Lagrangian planetary equations, which give the mean precession rate of the line of node, provide the theoretical Around Mars foundation of sun-synchronous orbits. Boain described basic mechanical characteristics of Earth’s sun-synchronous orbits and ∗ † ‡ provided several design algorithms to satisfy practical mission Xiaodong Liu, Hexi Baoyin, and Xingrui Ma requirements [4]. Several Martian satellites have been placed into Tsinghua University, sun-synchronous orbits, such as MGS [1] and MRO [2]. 100084 Beijing, People’s Republic of China Orbits at the critical inclination can keep eccentricity and argument of perigee invariable on average. For example, the Molniya DOI: 10.2514/1.48706 and Tundra orbits applied such conditions to stop the rotation of argument of perigee. Orlov introduced the concept of the Earth Nomenclature critical inclination [5], and Brouwer eliminated short-period terms based on canonical transformations [6]. Later, Coffey et al. analyzed a = semimajor axis, km the averaged Hamilton system reduced by the Delaunay normal- e = orbital eccentricity ization and provided geometrical interpretation of the critical i = orbital inclination, deg. inclination for Earth’s artificial satellites [7]. Some research has J2 = second-order zonal harmonic regarded orbits at the critical inclination as types of frozen orbits J3 = third-order zonal harmonic [8,9]. J4 = fourth-order zonal harmonic Frozen orbits have always been an important focus of orbit design, J22 = second-degree and -order tesseral harmonic because their average eccentricity and argument of perigee remain M = mean anomaly, deg.
    [Show full text]
  • Open Rosen Thesis.Pdf
    THE PENNSYLVANIA STATE UNIVERSITY SCHREYER HONORS COLLEGE DEPARTMENT OF AEROSPACE ENGINEERING END OF LIFE DISPOSAL OF SATELLITES IN HIGHLY ELLIPTICAL ORBITS MITCHELL ROSEN SPRING 2019 A thesis submitted in partial fulfillment of the requirements for a baccalaureate degree in Aerospace Engineering with honors in Aerospace Engineering Reviewed and approved* by the following: Dr. David Spencer Professor of Aerospace Engineering Thesis Supervisor Dr. Mark Maughmer Professor of Aerospace Engineering Honors Adviser * Signatures are on file in the Schreyer Honors College. i ABSTRACT Highly elliptical orbits allow for coverage of large parts of the Earth through a single satellite, simplifying communications in the globe’s northern reaches. These orbits are able to avoid drastic changes to the argument of periapse by using a critical inclination (63.4°) that cancels out the first level of the geopotential forces. However, this allows the next level of geopotential forces to take over, quickly de-orbiting satellites. Thus, a balance between the rate of change of the argument of periapse and the lifetime of the orbit is necessitated. This thesis sets out to find that balance. It is determined that an orbit with an inclination of 62.5° strikes that balance best. While this orbit is optimal off of the critical inclination, it is still near enough that to allow for potential use of inclination changes as a deorbiting method. Satellites are deorbited when the propellant remaining is enough to perform such a maneuver, and nothing more; therefore, the less change in velocity necessary for to deorbit, the better. Following the determination of an ideal highly elliptical orbit, the different methods of inclination change is tested against the usual method for deorbiting a satellite, an apoapse burn to lower the periapse, to find the most propellant- efficient method.
    [Show full text]
  • Single Axis Tracking with the RC1500B the Apparent Motion Of
    Single Axis Tracking with the RC1500B The apparent motion of an inclined orbit satellite appears as a narrow figure 8 pattern aligned perpendicular to the geo-stationary satellite arc. As the inclination of the satellite increases both the height and the width of the figure 8 pattern increase. The single axis tracker can follow the long dimension of the figure 8 but cannot compensate for the width of the figure 8 pattern. A paper (available on our web site - http://www.researchconcepts.com/Files/track_wp.pdf) describes the height and width of the figure 8 pattern as a function of the inclination of the satellite’s orbit. Note that the inclination of the satellite increases with time. The maximum rate of increase is approximately 0.9 degrees per year. As the inclination increases, the width of the figure 8 pattern will also increase. This has two implications for system performance. One, the maximum signal loss due to antenna misalignment will increase with time, and two, the antenna must have range a of motion sufficient to accommodate the greatest satellite inclination that will be encountered. Many people feel that single axis tracking is viable for antenna’s up to 3.8 meters at C band and 2.4 meters at Ku band. Implicit in this is the fact that the inclination of most commercial satellites is not allowed to exceed 5 degrees. This assumption should be verified before a system is fielded. A single axis tracking system must be in precise mechanical alignment to minimize loss due to the mount’s inability to compensate for the width of the figure eight pattern.
    [Show full text]
  • Measurement Techniques and New Technologies for Satellite Monitoring
    Report ITU-R SM.2424-0 (06/2018) Measurement techniques and new technologies for satellite monitoring SM Series Spectrum management ii Rep. ITU-R SM.2424-0 Foreword The role of the Radiocommunication Sector is to ensure the rational, equitable, efficient and economical use of the radio- frequency spectrum by all radiocommunication services, including satellite services, and carry out studies without limit of frequency range on the basis of which Recommendations are adopted. The regulatory and policy functions of the Radiocommunication Sector are performed by World and Regional Radiocommunication Conferences and Radiocommunication Assemblies supported by Study Groups. Policy on Intellectual Property Right (IPR) ITU-R policy on IPR is described in the Common Patent Policy for ITU-T/ITU-R/ISO/IEC referenced in Annex 1 of Resolution ITU-R 1. Forms to be used for the submission of patent statements and licensing declarations by patent holders are available from http://www.itu.int/ITU-R/go/patents/en where the Guidelines for Implementation of the Common Patent Policy for ITU-T/ITU-R/ISO/IEC and the ITU-R patent information database can also be found. Series of ITU-R Reports (Also available online at http://www.itu.int/publ/R-REP/en) Series Title BO Satellite delivery BR Recording for production, archival and play-out; film for television BS Broadcasting service (sound) BT Broadcasting service (television) F Fixed service M Mobile, radiodetermination, amateur and related satellite services P Radiowave propagation RA Radio astronomy RS Remote sensing systems S Fixed-satellite service SA Space applications and meteorology SF Frequency sharing and coordination between fixed-satellite and fixed service systems SM Spectrum management Note: This ITU-R Report was approved in English by the Study Group under the procedure detailed in Resolution ITU-R 1.
    [Show full text]
  • Aas 12-269 Frozen Orbits for Scientific Missions Using
    AAS 12-269 FROZEN ORBITS FOR SCIENTIFIC MISSIONS USING ROTATING TETHERS Hodei Urrutxua,∗ Jes ´usPelaez´ y and Martin Laraz We derive a semi-analytic formulation that permits to study the long-term dynam- ics of fast-rotating inert tethers around planetary satellites. Since space tethers are extensive bodies they generate non-keplerian gravitational forces which depend solely on their mass geometry and attitude, that can be exploited for controlling science orbits. We conclude that rotating tethers modify the geometry of frozen orbits, allowing for lower eccentricity frozen orbits for a wide range of orbital in- clination, where the length of the tether becomes a new parameter that the mission analyst may use to shape frozen orbits to tighter operational constraints. INTRODUCTION Science orbits for missions to planetary satellites have, in general, low altitudes and near-polar inclinations so that the entire surface can be mapped, and the science requirements of the mission can be accomplished.1,2 However, designing such an orbit can be difficult because the dynamical environment of many planetary satellites is highly perturbed due to their proximity to the central planet.3 One such example is the Moon, on which we focus on the current study. High-inclination orbits around the Moon are known to be unstable.4,5,6 However, the onset of instability can be delayed by a proper orbit design, hence maximizing the orbital lifetime.7 In this regard, the so called frozen orbits offer an interesting starting point for the design of science orbits,1,8,9 since all their orbital elements, with the only exception of the mean anomaly and the longitude of the ascending node, remain constant.
    [Show full text]