DOCSLIB.ORG

Analysis of Electric Propulsion Systems for Drag Compensation of Small Satellites in Low Earth Orbits the Universtiy of Manchester

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Analysis of Electric Propulsion Systems for Drag Compensation of Small Satellites in Low Earth Orbits the Universtiy of Manchester Analysis of Electric Propulsion Systems for Drag Compensation of Small Satellites in Low Earth Orbits The Universtiy of Manchester Teodor Bozhanov, ID: 9023890 Supervisor: Dr. Peter C.E. Roberts Final Report 02 May 2017 Abstract Small satellites, in particular CubeSats, have been the study focus of many major space agencies, universities and private organisations. Recent studies have shown that small satellites can reach up to 95% of the operational capabilities of large satellites, for only 5% of the cost. One drawback however, is that CubeSats have no onboard propulsion system and therefore no orbit maintenance and manoeuvre capabilities. This, coupled with the fact that they are usually inserted at very low Earth orbits (VLEO (under 400 km)), means that orbital lifetime is extremely limited. Consequently, small satellites and CubeSats have reduced lifetime and operational capabilities, limiting the range of their missions. Electric propulsion systems can generate the required thrust for drag compensation, while being extremely efficient. This results in a relatively low propellant fraction, reducing the negative impact on the available payload. This study focuses on using various Electric Propulsion (EP) and Atmosphere Breath- ing Electric Propulsion (ABEP) systems to increase the lifetime and usefulness of the satellites. The scope of the study is limited to VLEO ranging from 100 to 300 km. This zone is ideal for Earth observation and reconnaissance missions. In addition, it falls within the range of the ABEP system, where higher atmospheric density is more favourable. The computer model generated for this study accounts for 6 different types of per- turbations and is able to model the change of orbital parameters with high degree of accuracy. Preliminary results from several different types of EP show an increase of orbital lifetime between 200 and 600%. In some cases the associated velocity change is sufficient for performing small orbital transfers, rendezvous and docking manoeuvres. Orbit raising of about 900 - 1200 km was observed from a couple of thrusters. It was seen that a trade-off zone, between EP and ABEP, starts to form from about 250 km below, where ABEP thrusters become more effective. 1 Contents 1 Introduction 9 1.1 History of Electric Propulsion . .9 1.2 Motivation . 12 1.3 Aims and Objectives . 15 1.4 Methodology . 16 2 Physical Background 17 2.1 Space Environment . 17 2.1.1 Atmosphere and Vacuum Environment . 17 2.1.2 Neutral Environment . 19 2.2 Electric Propulsion Systems . 20 2.2.1 Electrothermal Acceleration . 21 2.2.2 Electrostatic Acceleration . 23 2.2.3 Electromagnetic Acceleration . 24 2.3 Atmosphere Breathing Electric Propulsion . 25 3 Literature Review 27 3.1 Introduction . 27 3.2 Relevant Research . 27 3.2.1 Electric Propulsion . 27 3.2.2 Atmosphere Breathing Electric Propulsion (ABEP) . 30 3.2.3 Drag and Atmosphere Models . 33 3.3 Summary . 36 4 Methodology 37 4.1 Mission Setup . 37 4.2 Methodology . 37 4.3 Assumptions and Limitations . 38 5 Preliminary Analysis 41 5.1 Atmosphere Models . 41 5.1.1 Drag Variation With Altitude . 43 2 xCONTENTS Teodor Bozhanov 5.1.2 Drag Variation With Latitude and Longitude . 43 5.1.3 Drag Variation With Solar Activity . 46 5.2 Orbital Propagators . 48 5.2.1 Energy Methods . 48 5.2.2 Gauss's Planetary Equations (GPE) . 50 5.3 EP Systems . 53 5.3.1 Pulsed Plasma Thrusters . 53 5.3.2 Ion Thrusters . 54 5.3.3 Field Emission Electric Propulsion . 56 5.4 ABEP systems . 57 5.4.1 Intake . 57 5.4.2 Thrusters . 60 6 Results and Discussion 62 6.1 Drag Compensation With EP . 62 6.1.1 ARC PPT Thruster . 64 6.1.2 RIT Ion Thruster . 65 6.1.3 IFM Nano FEEP Thruster . 67 6.1.4 Summary . 69 6.2 Drag Compensation With ABEP . 74 6.2.1 Theoretical Predictions . 74 6.2.2 Numerical Simulations . 77 6.3 Future Work . 79 7 Conclusion 81 8 Project Management 83 8.1 Semester I . 83 8.2 Semester II . 85 3 List of Figures 1.1 Number of Commercial Satellites with on-board EP system (Hoskins, W.A. et al. 2013). 12 1.2 Small satellites trends (2016-2022) (Doncaster and Shulman 2016). 14 1.3 Space Debris Population at different altitudes. Active Debris Re- moval (ADR05) of 5 objects per year needs to be conducted to main- tain current debris density (Liou, Johnson, and Hill 2010). 15 2.1 Variation of atmospheric layers with altitude (km) (Tewari 2007) . 18 2.2 Single particle momentum transfer. (Tribble 1995) . 19 2.3 Oxidised silver which is flaked off, exposing the underlying fresh ma- terial which is oxidised again (Rooij 2010). 20 2.4 1D Schematic of an Electrothermal Thruster (Jahn 2006) . 22 2.5 1D Schematic of Electrostatic Acceleration (Sforza 2016) . 23 2.6 Schematic of MPD accelerator (Sforza 2016) . 24 2.7 Schematic of a typical PPT (Sforza 2016) . 25 2.8 Atmosphere Breathing Ion Engine (Nishiyama 2003) . 26 3.1 Two Air Intake concepts: Funnel Concept a); Bypass Concept b). 31 3.2 Variation of the different atmospheric constituents with altitude (Schon- herr et al. 2015) . 31 3.3 Variation of Relative Density with Solar Activity (Tribble 1995) . 34 4.1 Analysis of Electric Propulsion Systems: Flowchart . 39 5.1 Variation of Density With Altitude for 5 Different Atmosphere Models 42 5.2 Variation of Drag Force With Altitude for NRLMSISE and MSIS. 43 5.3 Drag Variation With Latitude, at Constant Longitude and Solar Flux, for Every Hour at 300 km Altitude. 44 5.4 Drag Variation With Longitude, at Constant Latitude and Solar Flux, for Every Hour at 300 km Altitude. 45 5.5 Maximum Drag Variation With Solar Flux, at Constant Longitude, Latitude and Time. 46 4 xLIST OF FIGURES Teodor Bozhanov 5.6 Solar Flux Drag Variation at 116 km altitude. 47 5.7 Orbital Decay for 1U, 2 kg CubeSat . 49 5.8 Orbital Decay for 1U, 2 kg CubeSat using GPE . 51 5.9 Orbital Decay for 1U, 2 kg CubeSat Using Extended GPE . 52 5.10 Extended GPE Propagator: Orbital Decay of a 1U Cubesat (different mission types) . 53 5.11 NASA S-iEPS thruster: tank configuration (left) and thruster circuit arrangement (Krejci, Mier-Hicks, et al. 2015). 55 5.12 IFM Nano Thruster Family (Reissner et al. 2015). 56 5.13 BUSEK long annular intake (Hohman 2012). 58 5.14 JAXA long annular bypass intake with a conical diffusion region (Fu- jita 2004). 58 5.15 DSMC analysis with integrated turbomolecular pump, where = 0.95 and P are the molecular flow transmission probabilities of incident particles. High speed case: 3 stages turbo (Fujita 2004). 59 5.16 Schematic of Inductive Plasma Thruster (Romano, Binder, Herdrich, Fasoulas, et al. 2016). 60 2 5.17 IPT thrust variation with altitude for Air at Af = 1m , ηc = 0.35 (Romano, Binder, Herdrich, Fasoulas, et al. 2016). 61 6.1 Orbital Lifetime with no propulsion system in place. 63 6.2 PPT Drag Compensation for 1U CubeSat. 64 6.3 PPT Drag Compensation for 2U CubeSat. 65 6.4 RIT-µX Drag Compensation for 1U CubeSat. 66 6.5 RIT-µX Drag Compensation for 1U CubeSat 200 km. 66 6.6 RIT-µX Drag Compensation for 2 and 3U CubeSats. 67 6.7 IFM Nano Drag Compensation for 2U CubeSat. 68 6.8 IFM Nano Drag Compensation for 2U CubeSat 200 km. 68 6.9 IFM Nano Drag Compensation for 3U CubeSat. 69 6.10 IFM Nano Drag Compensation for 3U CubeSat 200 km. 69 6.11 Various Thruster Arrangements (Reissner et al. 2015) . 71 6.12 Lifetime percentage increase with 30 g of propellant. 72 6.13 Change of ∆V with Propellant Mass. 73 6.14 Mass flow variation with altitude (Romano, Massuti, and Herdrich 2014) . 74 6.15 Variation of Drag force for 3U CubeSat . 76 6.16 ABEP orbital lifetime at 180 km, operational for 200 h. 77 6.17 ABEP orbital lifetime at 250 km, operational for 200 h. 78 8.1 Semester 1 Gantt Chart . 84 5 xLIST OF FIGURES Teodor Bozhanov 8.2 Gantt Chart for Semester 2. 86 6 List of Tables 5.1 Percentage Difference Between 4 Atmosphere Models and the Control Model(MSIS)............................... 42 5.2 Percentage Difference Between Energy Method and Commercially Available Software . 49 5.3 Percentage Difference Between Different Orbit Propagators . 51 5.4 Comparison of different PPT systems (Colleti, Ciaralli, and Gabriel 2015), (Krejci, Seifert, and Scharlemann 2013). 54 5.5 Comparison of various Ion Thrusters (Wright and Ferrer 2015), (Kre- jci, Mier-Hicks, et al. 2015). 55 5.6 Comparison of different FEEP systems (Wright and Ferrer 2015), (Reissner et al. 2015) . 56 5.7 BUSEK results of DSMC simulation (Romano, Binder, Herdrich, and S. 2015) . 58 5.8 JAXA results of DSMC simulation (Romano, Binder, Herdrich, and S. 2015) . 59 6.1 Electric thruster selection. 62 6.2 Extended GPE Propagator: Simulation Starting Parameters . 63 6.3 Simulation data for 1U CubeSats. 70 6.4 Simulation data for 2U CubeSats. 70 6.5 Simulation data for 3U CubeSats. 71 6.6 Atmosphere Breathing Electric Propulsion, Theoretical Predictions. 75 6.7 Theoretical Predictions for ABEP Systems. 76 7 xLIST OF TABLES Teodor Bozhanov Acronyms ABEP Atmosphere Breathing Electric Propulsion ABIE Atmosphere Breathing Ion Engine AFRL Air Force Research Laboratory AO Atomic Oxygen CAT Cathode-Arc Thruster CIRA COSPAR International Reference Model COSPAR Committee on Space Research DSMC Direct Simulation Monte Carlo DTM Drag Temperature Model EP Electric Propulsion EPRB Electric Propulsion Research Building ESA European Space Agency EUV Extreme Ultraviolet FEEP Field Emission Electric Propulsion GOCE Gravity field and steady-state Ocean Circulation Explorer IPG Inductive Plasma Generators Isp Specific Impulse ISRU In Situ Resource Utilisation JAXA Japan Aerospace eXploration Agency LEO Low Earth Orbit LISA Laser Interferometer Space Antenna MPD Magnetoplasmadynamic Thrusters MIT Massachusetts Institute of Technology MSIS Mass Spectrometer Incoherent Scatter NASA National Aeronautics and Space Administration NEP Nuclear Electric Propulsion PPT Pulsed Plasma Thruster SERT Space Electric Rocket Test SMART Small Missions for Advanced Research and Technology TLE Two-Line Element VASIMIR Variable Specific Impulse Magnetoplasma Rocket VLEO Very Low Earth.
Load more

Recommended publications
  • PDF (03Ragozzine Exo-Interiors.Pdf)
    20 Chapter 2 Probing the Interiors of Very Hot Jupiters Using Transit Light Curves This chapter will be published in its entirety under the same title by authors D. Ragozzine and A. S. Wolf in the Astrophysical Journal, 2009. Reproduced by permission of the American Astro- nomical Society. 21 Abstract Accurately understanding the interior structure of extra-solar planets is critical for inferring their formation and evolution. The internal density distribution of a planet has a direct effect on the star-planet orbit through the gravitational quadrupole field created by the rotational and tidal bulges. These quadrupoles induce apsidal precession that is proportional to the planetary Love number (k2p, twice the apsidal motion constant), a bulk physical characteristic of the planet that depends on the internal density distribution, including the presence or absence of a massive solid core. We find that the quadrupole of the planetary tidal bulge is the dominant source of apsidal precession for very hot Jupiters (a . 0:025 AU), exceeding the effects of general relativity and the stellar quadrupole by more than an order of magnitude. For the shortest-period planets, the planetary interior induces precession of a few degrees per year. By investigating the full photometric signal of apsidal precession, we find that changes in transit shapes are much more important than transit timing variations. With its long baseline of ultra-precise photometry, the space-based Kepler mission can realistically detect apsidal precession with the accuracy necessary to infer the presence or absence of a massive core in very hot Jupiters with orbital eccentricities as low as e ' 0:003.
    [Show full text]
  • Orbital Lifetime Predictions
    Orbital LIFETIME PREDICTIONS An ASSESSMENT OF model-based BALLISTIC COEFfiCIENT ESTIMATIONS AND ADJUSTMENT FOR TEMPORAL DRAG co- EFfiCIENT VARIATIONS M.R. HaneVEER MSc Thesis Aerospace Engineering Orbital lifetime predictions An assessment of model-based ballistic coecient estimations and adjustment for temporal drag coecient variations by M.R. Haneveer to obtain the degree of Master of Science at the Delft University of Technology, to be defended publicly on Thursday June 1, 2017 at 14:00 PM. Student number: 4077334 Project duration: September 1, 2016 – June 1, 2017 Thesis committee: Dr. ir. E. N. Doornbos, TU Delft, supervisor Dr. ir. E. J. O. Schrama, TU Delft ir. K. J. Cowan MBA TU Delft An electronic version of this thesis is available at http://repository.tudelft.nl/. Summary Objects in Low Earth Orbit (LEO) experience low levels of drag due to the interaction with the outer layers of Earth’s atmosphere. The atmospheric drag reduces the velocity of the object, resulting in a gradual decrease in altitude. With each decayed kilometer the object enters denser portions of the atmosphere accelerating the orbit decay until eventually the object cannot sustain a stable orbit anymore and either crashes onto Earth’s surface or burns up in its atmosphere. The capability of predicting the time an object stays in orbit, whether that object is space junk or a satellite, allows for an estimate of its orbital lifetime - an estimate satellite op- erators work with to schedule science missions and commercial services, as well as use to prove compliance with international agreements stating no passively controlled object is to orbit in LEO longer than 25 years.
    [Show full text]
  • Aerodynamic Characteristics of Naca 0012 Airfoil Section at Different Angles of Attack
    AERODYNAMIC CHARACTERISTICS OF NACA 0012 AIRFOIL SECTION AT DIFFERENT ANGLES OF ATTACK SUPREETH NARASIMHAMURTHY GRADUATE STUDENT 1327291 Table of Contents 1) Introduction………………………………………………………………………………………………………………………………………...1 2) Methodology……………………………………………………………………………………………………………………………………….3 3) Results……………………………………………………………………………………………………………………………………………......5 4) Conclusion …………………………………………………………………………………………………………………………………………..9 5) References…………………………………………………………………………………………………………………………………………10 List of Figures Figure 1: Basic nomenclature of an airfoil………………………………………………………………………………………………...1 Figure 2: Computational domain………………………………………………………………………………………………………………4 Figure 3: Static Pressure Contours for different angles of attack……………………………………………………………..5 Figure 4: Velocity Magnitude Contours for different angles of attack………………………………………………………………………7 Fig 5: Variation of Cl and Cd with alpha……………………………………………………………………………………………………8 Figure 6: Lift Coefficient and Drag Coefficient Ratio for Re = 50000…………………………………………………………8 List of Tables Table 1: Lift and Drag coefficients as calculated from lift and drag forces from formulae given above……7 Introduction It is a fact of common experience that a body in motion through a fluid experience a resultant force which, in most cases is mainly a resistance to the motion. A class of body exists, However for which the component of the resultant force normal to the direction to the motion is many time greater than the component resisting the motion, and the possibility of the flight of an airplane depends on the use of the body of this class for wing structure. Airfoil is such an aerodynamic shape that when it moves through air, the air is split and passes above and below the wing. The wing’s upper surface is shaped so the air rushing over the top speeds up and stretches out. This decreases the air pressure above the wing. The air flowing below the wing moves in a comparatively straighter line, so its speed and air pressure remain the same.
    [Show full text]
  • NASA Process for Limiting Orbital Debris
    NASA-HANDBOOK NASA HANDBOOK 8719.14 National Aeronautics and Space Administration Approved: 2008-07-30 Washington, DC 20546 Expiration Date: 2013-07-30 HANDBOOK FOR LIMITING ORBITAL DEBRIS Measurement System Identification: Metric APPROVED FOR PUBLIC RELEASE – DISTRIBUTION IS UNLIMITED NASA-Handbook 8719.14 This page intentionally left blank. Page 2 of 174 NASA-Handbook 8719.14 DOCUMENT HISTORY LOG Status Document Approval Date Description Revision Baseline 2008-07-30 Initial Release Page 3 of 174 NASA-Handbook 8719.14 This page intentionally left blank. Page 4 of 174 NASA-Handbook 8719.14 This page intentionally left blank. Page 6 of 174 NASA-Handbook 8719.14 TABLE OF CONTENTS 1 SCOPE...........................................................................................................................13 1.1 Purpose................................................................................................................................ 13 1.2 Applicability ....................................................................................................................... 13 2 APPLICABLE AND REFERENCE DOCUMENTS................................................14 3 ACRONYMS AND DEFINITIONS ...........................................................................15 3.1 Acronyms............................................................................................................................ 15 3.2 Definitions .........................................................................................................................
    [Show full text]
  • NEAR of the 21 Lunar Landings, 19—All of the U.S
    Copyrights Prof Marko Popovic 2021 NEAR Of the 21 lunar landings, 19—all of the U.S. and Russian landings—occurred between 1966 and 1976. Then humanity took a 37-year break from landing on the moon before China achieved its first lunar touchdown in 2013. Take a look at the first 21 successful lunar landings on this interactive map https://www.smithsonianmag.com/science-nature/interactive-map-shows-all-21-successful-moon-landings-180972687/ Moon 1 The near side of the Moon with the major maria (singular mare, vocalized mar-ray) and lunar craters identified. Maria means "seas" in Latin. The maria are basaltic lava plains: i.e., the frozen seas of lava from lava flows. The maria cover ∼ 16% (30%) of the lunar surface (near side). Light areas are Lunar Highlands exhibiting more impact craters than Maria. The far side is pocked by ancient craters, mountains and rugged terrain, largely devoid of the smooth maria we see on the near side. The Lunar Reconnaissance Orbiter Moon 2 is a NASA robotic spacecraft currently orbiting the Moon in an eccentric polar mapping orbit. LRO data is essential for planning NASA's future human and robotic missions to the Moon. Launch date: June 18, 2009 Orbital period: 2 hours Orbit height: 31 mi Speed on orbit: 0.9942 miles/s Cost: 504 million USD (2009) The Moon is covered with a gently rolling layer of powdery soil with scattered rocks called the regolith; it is made from debris blasted out of the Lunar craters by the meteor impacts that created them.
    [Show full text]
  • Evaluation of Drag and Lift in the Internal Flow Field of a Dual Rotor Spinning Unit Via CFD
    MATEC Web of Conferences 104, 02005 (2017) DOI: 10.1051/ matecconf/201710402005 IC4M & ICDES 2017 Evaluation of drag and lift in the internal flow field of a dual rotor spinning unit via CFD Nicholus Tayari Akankwasa1, Huiting Lin and Jun Wang1,2,a 1College of Textile, Donghua University, Shanghai, 201620, China 2The Key Lab of Textile Science and Technology, Ministry of Education, Shanghai 201620, China Abstract. In the present study, we evaluate the drag and lift magnitude in the new dual-feed rotor spinning unit using computational fluid dynamics technique. We adopt theoretical and numerical approach based on FLUENT to investigate the influence of new design on the drag and lift in the rotor interior. Results reveal that the drag and lift inside the rotor of the proposed model are reduced by 60-80% and 50-66% respectively as compared to the conventional unit. The velocity and pressure profiles become evenly distributed in the dual- feed rotor interior as opposed to the conventional rotor spinning unit and this modification is anticipated to improve fiber configurations. This phenomenon can be utilized to further optimize the rotor spinning unit and other wall-bounded engineering problems. 1 Introduction In the fluid dynamics concept, previous studies have focused more on the turbulent viscosity, flocculation, eddies and vortices among other flow properties. Growing interest in drag reduction in flows has increased in the recent years. Principally, the drag and lift phenomenon has been widely applied in solving aerodynamics and aeronautics problems. For wall-bounded flows, it is important to note that drag and lift has a significant impact on the processing of fibers and polymers.
    [Show full text]
  • Study of Extraterrestrial Disposal of Radioactive Wastes
    NASA TECHNICAL NASA TM X-71557 MEMORANDUM N--- XI NASA-TM-X-71557) STUDY OF N74-22776 EXTRATERRESTRIAL DISPOSAL OF READIOACTIVE WASTES. PART 1: SPACE TRANSPORTATION. AN.. DESTINATION CONSIDERATIONS. FOR, (NASA) Unclas 64 p HC $6.25 CSCL 18G G3/05 38494 STUDY OF EXTRATERRESTRIAL DISPOSAL OF RADIOACTIVE WASTES by R. L. Thompson, J. R. Ramler, and S. M. Stevenson Lewis Research Center Cleveland, Ohio 44135 May 1974 This information is being published in prelimi- nary form in order to expedite its early release. I STUDY OF EXTRATERRESTRIAL DISPOSAL OF RADIOACTIVE WASTES Part I Space Transportation and Destination Considerations for Extraterrestrial Disposal of Radioactive W4astes by R. L. Thompson, J. R. Ramler, and S. M. Stevenson Lewis Research Center CT% CSUMMARY I NASA has been requested by the AEC to conduct a feasibility study of extraterrestrial (space) disposal of radioactive waste. This report covers the initial work done on only one part of the NASA study, the evaluation and comparison of possible space destinations and space transportation systems. Only current or planned space transportation systems have been considered thus far. The currently planned Space Shuttle was found to be more cost-effective than current expendable launch vehicles by about a factor of 2. The Space Shuttle requires a third stage to perform the waste disposal missions. Depending on the particular mission, this third stage could be either a reusable space tug or an expendable stage such as a Centaur. Of the destinations considered, high Earth orbits (between geostationary and lunar orbit altitudes), solar orbits (such as a 0.90 AU circular solar orbit) or a direct injection to solar system escape appear to be the best candidates.
    [Show full text]
  • New Satellite Drag Modeling Capabilities
    44th AIAA Aerospace Sciences Meeting and Exhibit AIAA 2006-470 9 - 12 January 2006, Reno, Nevada New Satellite Drag M odeling Capabilities Frank A. Marcos * Air Force Research Laboratory , Hanscom AFB, MA 01731 -3010 This paper reviews the operational impacts of satellite drag, the historical and current capabilities, and requirements to deal with evo lving higher accuracy requirements. Modeling of satellite drag variations showed little improvement from the 1960’s to the late 1990’s. After three decades of essentially no quantitative progress, the problem is being vigorously and fruitfully attacked on several fronts. This century has already shown significant advances in measurements, models, solar and geomagnetic proxies and the application of data assimilation techniques to operational applications. While thermospheric measurements have been historica lly extremely sparse, new data sets are now available from intense ground -based radar tracking of satellite orbital decay and from satellite -borne accelerometers and remote sensors. These data provide global coverage over a wide range of thermospheric alti tudes. Operational assimilative empirical models, utilizing the orbital drag data, have reduced model errors by almost a factor of two. Together with evolving new solar and geomagnetic inputs, the satellite -borne sensors support development of advanced ope rational assimilative first principles forecast models. We look forward to the time when satellite drag is no longer the largest error source in determining or bits of low altitude satellites. I. Introduction Aerodynamic drag continues to be the larg est uncertainty in precision orbit determin ation for satellites operating below about 600 km. Drag errors impact many aerospace missions including satellite orbit location and prediction, collision avoidance warnings, reentry prediction, lifetime estimates and attitude dynamics.
    [Show full text]
  • Analysis of Orbital Decay Time for the Classical Hydrogen Atom Interacting with Circularly Polarized Electromagnetic Radiation
    Analysis of Orbital Decay Time for the Classical Hydrogen Atom Interacting with Circularly Polarized Electromagnetic Radiation Daniel C. Cole and Yi Zou Dept. of Manufacturing Engineering, 15 St. Mary’s Street, Boston University, Brookline, MA 02446 Abstract Here we show that a wide range of states of phases and amplitudes exist for a circularly polarized (CP) plane wave to act on a classical hydrogen model to achieve infinite times of stability (i.e., no orbital decay due to radiation reaction effects). An analytic solution is first deduced to show this effect for circular orbits in the nonrelativistic approximation. We then use this analytic result to help provide insight into detailed simulation investigations of deviations from these idealistic conditions. By changing the phase of the CP wave, the time td when orbital decay sets in can be made to vary enormously. The patterns of this behavior are examined here and analyzed in physical terms for the underlying but rather unintuitive reasons for these nonlinear effects. We speculate that most of these effects can be generalized to analogous elliptical orbital conditions with a specific infinite set of CP waves present. The article ends by briefly considering multiple CP plane waves acting on the classical hydrogen atom in an initial circular orbital state, resulting in “jump- and diffusion-like” orbital motions for this highly nonlinear system. These simple examples reveal the possibility of very rich and complex patterns that occur when a wide spectrum of radiation acts on this classical hydrogen system. 1 I. INTRODUCTION The hydrogen atom has received renewed attention in the past decade or so, due to studies involved with Rydberg analysis, chaos, and scarring [1], [2], [3], [4].
    [Show full text]
  • Collisions & Encounters I
    Collisions & Encounters I A rAS ¡ r r AB b 0 rBS B Let A encounter B with an initial velocity v and an impact parameter b. 1 A star S (red dot) in A gains energy wrt the center of A due to the fact that the center of A and S feel a different gravitational force due to B. Let ~v be the velocity of S wrt A then dES = ~v ~g[~r (t)] ~v ~ Φ [~r (t)] ~ Φ [~r (t)] dt · BS ≡ · −r B AB − r B BS We define ~r0 as the position vector ~rAB of closest approach, which occurs at time t0. Collisions & Encounters II If we increase v then ~r0 b and the energy increase 1 j j ! t0 ∆ES(t0) ~v ~g[~rBS(t)] dt ≡ 0 · R dimishes, simply because t0 becomes smaller. Thus, for a larger impact velocity v the star S withdraws less energy from the relative orbit between 1 A and B. This implies that we can define a critical velocity vcrit, such that for v > vcrit galaxy A reaches ~r0 with sufficient energy to escape to infinity. 1 If, on the other hand, v < vcrit then systems A and B will merge. 1 If v vcrit then we can use the impulse approximation to analytically 1 calculate the effect of the encounter. < In most cases of astrophysical interest, however, v vcrit and we have 1 to resort to numerical simulations to compute the outcome∼ of the encounter. However, in the special case where MA MB or MA MB we can describe the evolution with dynamical friction , for which analytical estimates are available.
    [Show full text]
  • MATH 240: HOMEWORK #4 1. the Drag Equation Suppose That An
    MATH 240: HOMEWORK #4 DUE IN FLORA'S MAILBOX BY NOON ON NOV. 25. 1. The drag equation Suppose that an object of mass m is falling toward the earth, and that it has height h(t) above the surface at time t. Newton's law says that (1.1) Force = mass × acceleration: Two kinds of forces act on the object: A. A downward force of gravity having magnitude mg, where g is the gravitation constant. B. The drag force FD due to wind resistance. This acts in the opposite direction to the velocity h0(t), and its magnitude is 0 2 (1.2) jFDj = cA(t)h (t) where c > 0 is constant and A(t) is the cross sectional area of the object at time t. Notice that if A(t) is constant, the drag force goes up as the square of the velocity. This is due to the fact that the energy imparted by each molecule of air hit during the fall is proportion to velocity, and the number of molecules hit per second is also proportional to velocity. 1. Show that if we assume h(t) is decreasing with time (corresponding to falling toward the earth), (1.1) becomes (1.3) −mg + cA(t)h0(t)2 = mh00(t) In particular, explain the signs of the terms on the left. 2. Suppose now that the object is a parachutist with a circular parachute. They open the parachute so that it's radius is b=pjh0(t)j at time t for some constant b > 0. Are they making the parachute larger or smaller as the velocity decreases in magnitude? What is A(t) in this case, and what differental equation does (1.3) become? Note: Be careful to make sure that the signs of terms in the differential equation agree with the fact that the drag force acts in the opposite direction to the velocity of the parachutist.
    [Show full text]
  • Low-Speed Aerodynamic Characteristics of a Delta Wing With
    Low-Speed Aerodynamic Characteristics of a Delta Wing with Deflected Wing Tips Thesis Presented in Partial Fulfillment of the Requirements for the Degree of Master of Science in the Graduate School of The Ohio State University By Colin Weidner Trussa Graduate Program in Aeronautical and Astronautical Engineering The Ohio State University 2020 Master’s Examination Committee Dr. Clifford Whitfield, Advisor Dr. Rick Freuler Dr. Matthew McCrink Copyrighted by Colin Weidner Trussa 2020 2 ABSTRACT The purpose of this work was to investigate the low-speed aerodynamic characteristics of a novel delta wing layout with deflected wing tips. This project is motivated by the ongoing unmanned aerial vehicle research and development at The Ohio State University Aerospace Research Center. The model under test for this study had four main design requirements: (1) high-speed, (2) highly maneuverable, (3) aerodynamically interesting, and (4) multi-configurable. The last three requirements are addressed directly in this report with specific emphasis on requirements two and three. A modular fuselage design satisfied requirement four, and the novel delta wing addressed requirements two and three. The novel delta wing has a leading-edge sweep of 60 degrees, a high-speed airfoil with a rounded leading-edge, and wing tips that can rotate a full 180 degrees about a hinge, located at 2/3rds of the half-span parallel to fuselage centerline. Three different wing tip deflection configurations were analyzed: positive, negative, and asymmetric. Positive wing tip deflection corresponds to the wing tips being deflected up towards the vertical tail. Negative wing tip deflection is when the wing tips are deflected down, away from the vertical tail.
    [Show full text]