High-Accuracy Orbital Dynamics Simulation Through Keplerian and Equinoctial Parameters
Total Page:16
File Type:pdf, Size:1020Kb
Load more
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. -
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 -
A Study of Ancient Khmer Ephemerides
A study of ancient Khmer ephemerides François Vernotte∗ and Satyanad Kichenassamy** November 5, 2018 Abstract – We study ancient Khmer ephemerides described in 1910 by the French engineer Faraut, in order to determine whether they rely on observations carried out in Cambodia. These ephemerides were found to be of Indian origin and have been adapted for another longitude, most likely in Burma. A method for estimating the date and place where the ephemerides were developed or adapted is described and applied. 1 Introduction Our colleague Prof. Olivier de Bernon, from the École Française d’Extrême Orient in Paris, pointed out to us the need to understand astronomical systems in Cambo- dia, as he surmised that astronomical and mathematical ideas from India may have developed there in unexpected ways.1 A proper discussion of this problem requires an interdisciplinary approach where history, philology and archeology must be sup- plemented, as we shall see, by an understanding of the evolution of Astronomy and Mathematics up to modern times. This line of thought meets other recent lines of research, on the conceptual evolution of Mathematics, and on the definition and measurement of time, the latter being the main motivation of Indian Astronomy. In 1910 [1], the French engineer Félix Gaspard Faraut (1846–1911) described with great care the method of computing ephemerides in Cambodia used by the horas, i.e., the Khmer astronomers/astrologers.2 The names for the astronomical luminaries as well as the astronomical quantities [1] clearly show the Indian origin ∗F. Vernotte is with UTINAM, Observatory THETA of Franche Comté-Bourgogne, University of Franche Comté/UBFC/CNRS, 41 bis avenue de l’observatoire - B.P. -
A Hot Subdwarf-White Dwarf Super-Chandrasekhar Candidate
A hot subdwarf–white dwarf super-Chandrasekhar candidate supernova Ia progenitor Ingrid Pelisoli1,2*, P. Neunteufel3, S. Geier1, T. Kupfer4,5, U. Heber6, A. Irrgang6, D. Schneider6, A. Bastian1, J. van Roestel7, V. Schaffenroth1, and B. N. Barlow8 1Institut fur¨ Physik und Astronomie, Universitat¨ Potsdam, Haus 28, Karl-Liebknecht-Str. 24/25, D-14476 Potsdam-Golm, Germany 2Department of Physics, University of Warwick, Coventry, CV4 7AL, UK 3Max Planck Institut fur¨ Astrophysik, Karl-Schwarzschild-Straße 1, 85748 Garching bei Munchen¨ 4Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 5Texas Tech University, Department of Physics & Astronomy, Box 41051, 79409, Lubbock, TX, USA 6Dr. Karl Remeis-Observatory & ECAP, Astronomical Institute, Friedrich-Alexander University Erlangen-Nuremberg (FAU), Sternwartstr. 7, 96049 Bamberg, Germany 7Division of Physics, Mathematics and Astronomy, California Institute of Technology, Pasadena, CA 91125, USA 8Department of Physics and Astronomy, High Point University, High Point, NC 27268, USA *[email protected] ABSTRACT Supernova Ia are bright explosive events that can be used to estimate cosmological distances, allowing us to study the expansion of the Universe. They are understood to result from a thermonuclear detonation in a white dwarf that formed from the exhausted core of a star more massive than the Sun. However, the possible progenitor channels leading to an explosion are a long-standing debate, limiting the precision and accuracy of supernova Ia as distance indicators. Here we present HD 265435, a binary system with an orbital period of less than a hundred minutes, consisting of a white dwarf and a hot subdwarf — a stripped core-helium burning star. -
Electric Propulsion System Scaling for Asteroid Capture-And-Return Missions
Electric propulsion system scaling for asteroid capture-and-return missions Justin M. Little⇤ and Edgar Y. Choueiri† Electric Propulsion and Plasma Dynamics Laboratory, Princeton University, Princeton, NJ, 08544 The requirements for an electric propulsion system needed to maximize the return mass of asteroid capture-and-return (ACR) missions are investigated in detail. An analytical model is presented for the mission time and mass balance of an ACR mission based on the propellant requirements of each mission phase. Edelbaum’s approximation is used for the Earth-escape phase. The asteroid rendezvous and return phases of the mission are modeled as a low-thrust optimal control problem with a lunar assist. The numerical solution to this problem is used to derive scaling laws for the propellant requirements based on the maneuver time, asteroid orbit, and propulsion system parameters. Constraining the rendezvous and return phases by the synodic period of the target asteroid, a semi- empirical equation is obtained for the optimum specific impulse and power supply. It was found analytically that the optimum power supply is one such that the mass of the propulsion system and power supply are approximately equal to the total mass of propellant used during the entire mission. Finally, it is shown that ACR missions, in general, are optimized using propulsion systems capable of processing 100 kW – 1 MW of power with specific impulses in the range 5,000 – 10,000 s, and have the potential to return asteroids on the order of 103 104 tons. − Nomenclature -
AFSPC-CO TERMINOLOGY Revised: 12 Jan 2019
AFSPC-CO TERMINOLOGY Revised: 12 Jan 2019 Term Description AEHF Advanced Extremely High Frequency AFB / AFS Air Force Base / Air Force Station AOC Air Operations Center AOI Area of Interest The point in the orbit of a heavenly body, specifically the moon, or of a man-made satellite Apogee at which it is farthest from the earth. Even CAP rockets experience apogee. Either of two points in an eccentric orbit, one (higher apsis) farthest from the center of Apsis attraction, the other (lower apsis) nearest to the center of attraction Argument of Perigee the angle in a satellites' orbit plane that is measured from the Ascending Node to the (ω) perigee along the satellite direction of travel CGO Company Grade Officer CLV Calculated Load Value, Crew Launch Vehicle COP Common Operating Picture DCO Defensive Cyber Operations DHS Department of Homeland Security DoD Department of Defense DOP Dilution of Precision Defense Satellite Communications Systems - wideband communications spacecraft for DSCS the USAF DSP Defense Satellite Program or Defense Support Program - "Eyes in the Sky" EHF Extremely High Frequency (30-300 GHz; 1mm-1cm) ELF Extremely Low Frequency (3-30 Hz; 100,000km-10,000km) EMS Electromagnetic Spectrum Equitorial Plane the plane passing through the equator EWR Early Warning Radar and Electromagnetic Wave Resistivity GBR Ground-Based Radar and Global Broadband Roaming GBS Global Broadcast Service GEO Geosynchronous Earth Orbit or Geostationary Orbit ( ~22,300 miles above Earth) GEODSS Ground-Based Electro-Optical Deep Space Surveillance -
A Highly Elliptical Orbit Space System for Hydrometeorological Monitoring of the Arctic Region by V
A highly elliptical orbit space system for hydrometeorological monitoring of the Arctic region by V. V. Asmus1, V. N. Dyadyuchenko2, Y. I. Nosenko3, G. M. Polishchuk4 and V. A. Selin3 The lack of reliable, frequently high latitudes. It has therefore been • Monitoring of climate change updated information on the Earth’s suggested that demonstration of polar ice caps is a signifi cant problem a hydrometeorological system of • Data collection and relay for weather forecasting, affecting satellites on highly elliptical orbit from land-, sea- and air-based forecast skill for the entire planet. The (HEO), called the “Arctica” system, observing platforms poor numerical weather prediction should be created to provide the (NWP) skill for the Arctic region necessary complex information for the • Exchange and dissemination of and the Earth’s northern territories diffi cult tasks involved in developing processed hydrometeorological is caused primarily by errors in the whole Arctic region. and heliogeophysical data. determining initial conditions, which depend on the quality of initial Signifi cantly, the hydrometeorological Further progress in global and data. Until now, initial data have observations carried out in the regional numerical weather prediction been received from meteorological Arctic within the framework of the depends to a large extent on: geostationary satellites, which are International Polar Year 2007-2008 not very effective in scanning high are not provided with remote-sensing • Quasi-continuous reception latitudes and polar-orbiting -
GPS Applications in Space
Space Situational Awareness 2015: GPS Applications in Space James J. Miller, Deputy Director Policy & Strategic Communications Division May 13, 2015 GPS Extends the Reach of NASA Networks to Enable New Space Ops, Science, and Exploration Apps GPS Relative Navigation is used for Rendezvous to ISS GPS PNT Services Enable: • Attitude Determination: Use of GPS enables some missions to meet their attitude determination requirements, such as ISS • Real-time On-Board Navigation: Enables new methods of spaceflight ops such as rendezvous & docking, station- keeping, precision formation flying, and GEO satellite servicing • Earth Sciences: GPS used as a remote sensing tool supports atmospheric and ionospheric sciences, geodesy, and geodynamics -- from monitoring sea levels and ice melt to measuring the gravity field ESA ATV 1st mission JAXA’s HTV 1st mission Commercial Cargo Resupply to ISS in 2008 to ISS in 2009 (Space-X & Cygnus), 2012+ 2 Growing GPS Uses in Space: Space Operations & Science • NASA strategic navigation requirements for science and 20-Year Worldwide Space Mission space ops continue to grow, especially as higher Projections by Orbit Type* precisions are needed for more complex operations in all space domains 1% 5% Low Earth Orbit • Nearly 60%* of projected worldwide space missions 27% Medium Earth Orbit over the next 20 years will operate in LEO 59% GeoSynchronous Orbit – That is, inside the Terrestrial Service Volume (TSV) 8% Highly Elliptical Orbit Cislunar / Interplanetary • An additional 35%* of these space missions that will operate at higher altitudes will remain at or below GEO – That is, inside the GPS/GNSS Space Service Volume (SSV) Highly Elliptical Orbits**: • In summary, approximately 95% of projected Example: NASA MMS 4- worldwide space missions over the next 20 years will satellite constellation. -
Elliptical Orbits
1 Ellipse-geometry 1.1 Parameterization • Functional characterization:(a: semi major axis, b ≤ a: semi minor axis) x2 y 2 b p + = 1 ⇐⇒ y(x) = · ± a2 − x2 (1) a b a • Parameterization in cartesian coordinates, which follows directly from Eq. (1): x a · cos t = with 0 ≤ t < 2π (2) y b · sin t – The origin (0, 0) is the center of the ellipse and the auxilliary circle with radius a. √ – The focal points are located at (±a · e, 0) with the eccentricity e = a2 − b2/a. • Parameterization in polar coordinates:(p: parameter, 0 ≤ < 1: eccentricity) p r(ϕ) = (3) 1 + e cos ϕ – The origin (0, 0) is the right focal point of the ellipse. – The major axis is given by 2a = r(0) − r(π), thus a = p/(1 − e2), the center is therefore at − pe/(1 − e2), 0. – ϕ = 0 corresponds to the periapsis (the point closest to the focal point; which is also called perigee/perihelion/periastron in case of an orbit around the Earth/sun/star). The relation between t and ϕ of the parameterizations in Eqs. (2) and (3) is the following: t r1 − e ϕ tan = · tan (4) 2 1 + e 2 1.2 Area of an elliptic sector As an ellipse is a circle with radius a scaled by a factor b/a in y-direction (Eq. 1), the area of an elliptic sector PFS (Fig. ??) is just this fraction of the area PFQ in the auxiliary circle. b t 2 1 APFS = · · πa − · ae · a sin t a 2π 2 (5) 1 = (t − e sin t) · a b 2 The area of the full ellipse (t = 2π) is then, of course, Aellipse = π a b (6) Figure 1: Ellipse and auxilliary circle. -
Interplanetary Trajectories in STK in a Few Hundred Easy Steps*
Interplanetary Trajectories in STK in a Few Hundred Easy Steps* (*and to think that last year’s students thought some guidance would be helpful!) Satellite ToolKit Interplanetary Tutorial STK Version 9 INITIAL SETUP 1) Open STK. Choose the “Create a New Scenario” button. 2) Name your scenario and, if you would like, enter a description for it. The scenario time is not too critical – it will be updated automatically as we add segments to our mission. 3) STK will give you the opportunity to insert a satellite. (If it does not, or you would like to add another satellite later, you can click on the Insert menu at the top and choose New…) The Orbit Wizard is an easy way to add satellites, but we will choose Define Properties instead. We choose Define Properties directly because we want to use a maneuver-based tool called the Astrogator, which will undo any initial orbit set using the Orbit Wizard. Make sure Satellite is selected in the left pane of the Insert window, then choose Define Properties in the right-hand pane and click the Insert…button. 4) The Properties window for the Satellite appears. You can access this window later by right-clicking or double-clicking on the satellite’s name in the Object Browser (the left side of the STK window). When you open the Properties window, it will default to the Basic Orbit screen, which happens to be where we want to be. The Basic Orbit screen allows you to choose what kind of numerical propagator STK should use to move the satellite. -
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. -
Perturbation Theory in Celestial Mechanics
Perturbation Theory in Celestial Mechanics Alessandra Celletti Dipartimento di Matematica Universit`adi Roma Tor Vergata Via della Ricerca Scientifica 1, I-00133 Roma (Italy) ([email protected]) December 8, 2007 Contents 1 Glossary 2 2 Definition 2 3 Introduction 2 4 Classical perturbation theory 4 4.1 The classical theory . 4 4.2 The precession of the perihelion of Mercury . 6 4.2.1 Delaunay action–angle variables . 6 4.2.2 The restricted, planar, circular, three–body problem . 7 4.2.3 Expansion of the perturbing function . 7 4.2.4 Computation of the precession of the perihelion . 8 5 Resonant perturbation theory 9 5.1 The resonant theory . 9 5.2 Three–body resonance . 10 5.3 Degenerate perturbation theory . 11 5.4 The precession of the equinoxes . 12 6 Invariant tori 14 6.1 Invariant KAM surfaces . 14 6.2 Rotational tori for the spin–orbit problem . 15 6.3 Librational tori for the spin–orbit problem . 16 6.4 Rotational tori for the restricted three–body problem . 17 6.5 Planetary problem . 18 7 Periodic orbits 18 7.1 Construction of periodic orbits . 18 7.2 The libration in longitude of the Moon . 20 1 8 Future directions 20 9 Bibliography 21 9.1 Books and Reviews . 21 9.2 Primary Literature . 22 1 Glossary KAM theory: it provides the persistence of quasi–periodic motions under a small perturbation of an integrable system. KAM theory can be applied under quite general assumptions, i.e. a non– degeneracy of the integrable system and a diophantine condition of the frequency of motion.