Out-of-plane orbital maneuvers using swing-bys with the Moon

Similar documents
SIMULATION OF THREE PHASE THREE LEG TRANSFORMER BEHAVIOR UNDER DIFFERENT VOLTAGE SAG TYPES

Analytical Modeling of Natural Convection in Horizontal Annuli

Chapter 3, Solution 1C.

Design of Analog Integrated Circuits

Feedback Principle :-

Physic 231 Lecture 33

Transient Conduction: Spatial Effects and the Role of Analytical Solutions

Shell Stiffness for Diffe ent Modes

A New Method for Solving Integer Linear. Programming Problems with Fuzzy Variables

Chapter 7. Systems 7.1 INTRODUCTION 7.2 MATHEMATICAL MODELING OF LIQUID LEVEL SYSTEMS. Steady State Flow. A. Bazoune

Advances in Engineering Research (AER), volume 102 Second International Conference on Mechanics, Materials and Structural Engineering (ICMMSE 2017)

Grade 12 Physics Exam Review

Exploiting vector space properties for the global optimization of process networks

Chem 204A, Fall 2004, Mid-term (II)

Wp/Lmin. Wn/Lmin 2.5V

Spring 2002 Lecture #17

Conduction Heat Transfer

CHAPTER 3 ANALYSIS OF KY BOOST CONVERTER

Lecture 12. Heat Exchangers. Heat Exchangers Chee 318 1

Chapter 6 : Gibbs Free Energy

Introduction to Electronic circuits.

Bi-level Optimization Method of Air-conditioning System Based on Office Building Energy Storage Characteristics

The APS Transfer Line from Linac to Injector Synchrotron.

Approach: (Equilibrium) TD analysis, i.e., conservation eqns., state equations Issues: how to deal with

Physics 107 HOMEWORK ASSIGNMENT #20

A Proposal of Heating Load Calculation considering Stack Effect in High-rise Buildings

55:041 Electronic Circuits

IGEE 401 Power Electronic Systems. Solution to Midterm Examination Fall 2004

Center of Mass and Momentum. See animation An Object Tossed Along a Parabolic Path.

CIRCLE YOUR DIVISION: Div. 1 (9:30 am) Div. 2 (11:30 am) Div. 3 (2:30 pm) Prof. Ruan Prof. Naik Mr. Singh

4DVAR, according to the name, is a four-dimensional variational method.

_J _J J J J J J J J _. 7 particles in the blue state; 3 particles in the red state: 720 configurations _J J J _J J J J J J J J _

Section 3: Detailed Solutions of Word Problems Unit 1: Solving Word Problems by Modeling with Formulas

element k Using FEM to Solve Truss Problems

ME2142/ME2142E Feedback Control Systems. Modelling of Physical Systems The Transfer Function

6. ELUTRIATION OF PARTICLES FROM FLUIDIZED BEDS

Conservation of Energy

Chapter 5: Force and Motion I-a

WYSE Academic Challenge 2004 Sectional Physics Solution Set

Water vapour balance in a building moisture exposure for timber structures

CONVEX COMBINATIONS OF ANALYTIC FUNCTIONS

Differentiation Applications 1: Related Rates

Downscaling Geopotential Height Using Lapse Rate

Phys101 Final Code: 1 Term: 132 Wednesday, May 21, 2014 Page: 1

Energy & Work

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

AP Physics Kinematic Wrap Up

Comparison of Building Codes and Insulation in China and Iceland

PHYSICS 536 Experiment 12: Applications of the Golden Rules for Negative Feedback

Integrating Certified Lengths to Strengthen Metrology Network Uncertainty

EE 204 Lecture 25 More Examples on Power Factor and the Reactive Power

ECE 2100 Circuit Analysis

lecture 5: Nucleophilic Substitution Reactions

A method of constructing rock-analysis diagrams a statistical basks.

Cambridge Assessment International Education Cambridge Ordinary Level. Published

EE 221 Practice Problems for the Final Exam

Analysis The characteristic length of the junction and the Biot number are

CTN 2/23/16. EE 247B/ME 218: Introduction to MEMS Design Lecture 11m2: Mechanics of Materials. Copyright 2016 Regents of the University of California

PT326 PROCESS TRAINER

and the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:

Phys 344 Ch 5 Lect 4 Feb 28 th,

Problem Set 5 Solutions - McQuarrie Problems 3.20 MIT Dr. Anton Van Der Ven

Module B3. VLoad = = V S V LN

Activity Guide Loops and Random Numbers

ENGI 4430 Parametric Vector Functions Page 2-01

A Note on the Linear Programming Sensitivity. Analysis of Specification Constraints. in Blending Problems

Verification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Circuits Op-Amp. Interaction of Circuit Elements. Quick Check How does closing the switch affect V o and I o?

GENERATION OF ROBUST ERROR RECOVERY LOGIC IN ASSEMBLY SYSTEMS USING MULTI-LEVEL OPTIMIZATION AND GENETIC PROGRAMMING

Kinetics of Particles. Chapter 3

Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Professor and Chair Mechanical Engineering Department Christian Brothers University 650 East Parkway South Memphis, TN

SUMMER REV: Half-Life DUE DATE: JULY 2 nd

Theoretical study of third virial coefficient with Kihara potential

Phys. 344 Ch 7 Lecture 8 Fri., April. 10 th,

Chapter 3: Cluster Analysis

Kepler's Laws of Planetary Motion

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.

Van der Waals-coupled electronic states in incommensurate double-walled carbon nanotubes

ASTRODYNAMICS. o o o. Early Space Exploration. Kepler's Laws. Nicolaus Copernicus ( ) Placed Sun at center of solar system

SELECTION OF MODEL PARAMETERS OF BIOGAS IC ENGINE. Karol Cupiał, Grzegorz Katolik

XXIX CILAMCE November 4 th to 7 th, 2008 Maceió - Brazil

A HYDRAULIC OPEN LOOP SYSTEM FOR CONTROLLED EXCAVATION ALONG PRESCRIBED PATH. E. Bundy, W. Gutkowski

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Journal of Intelligent and Robotic Systems Vol. 35, pp

Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff

Intersection of an Ellipsoid and a Plane

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Where a licence is displayed above, please note the terms and conditions of the licence govern your use of this document.

SVD Ladder Mock-up Assembly for the Belle II Experiment

A Comparative Study of Resonances in the Dynamics of the Mimas-Tethys System

1 PreCalculus AP Unit G Rotational Trig (MCR) Name:

1.2.1 Vectors. 1 P age. Examples What is the reference vector angle for a vector that points 50 degrees east of south?

Chapter 3 Kinematics in Two Dimensions; Vectors

The influence of a semi-infinite atmosphere on solar oscillations

Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff

Transcription:

Jurnal f Physcs: Cnference Seres PAPER OPEN ACCESS Out-f-plane rbtal maneuvers usng swng-bys wth the Mn Related cntent - Pwered Swng-By Maneuvers arund the Mn A F Slva, A F B A Prad and O C Wnter cte ths artcle: J B S Net et al 5 J. Phys.: Cnf. Ser. 64 4 ew the artcle nlne fr updates and enhancements. hs cntent was dwnladed frm IP address 48.5.3.83 n 9/8/8 at :3

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 Out-f-plane rbtal maneuvers usng swng-bys wth the Mn J B S Net, A F B A Prad and J K S Frmga Natnal Insttute fr Space Research INPE, Av ds Astrnautas 758, SJC-SP, Brazl. E-mal: net.jbs9@gmal.cm.br Abstract. hs paper has the gal f shwng sme cases f plane change maneuvers usng a swng-by wth the Mn t decrease the magntude f the mpulses used, when cmpared t a classcal Hhmann maneuver. he analytcal mdel s based n the "patched-cncs" apprach, where a seres f "tw-bdy" prblems s cnsdered t buld the whle maneuver. A study f the effects f the sem-majr axes f the transfer rbts were made, t cmplete sme prevus studes made n the lterature. he results shw that, fr sme fnal nclnatns, the use f the swng-by n the Mn s really advantageus.. Intrductn Durng the desgn phase f artfcal satellte mssns, ne the mst mprtant factrs s the cst. Optmzatn methds are usually appled t search fr lwer cst slutns frm the cnceptn t the end f the useful lfe f the equpment. Mre recently, studes are als made cnsderng the later dspsal f the spacecraft. Lkng fr slutns t reduce the magntude f the mpulses requred n rbtal maneuvers, several methds have been develped and used, such as the classcal transfers f Hhmann and B-ellptcal [], bth f them cnsderng cplanar rbts. Nw, f the bjectve s t make a plane change, t s necessary t apply sme knd f prpulsn wth a cmpnent perpendcular t the rbtal plane. hs type f maneuver are usually very expensve n terms f fuel cnsumptn and dfferent ptns have been cnsdered t reduce the fuel cnsumptn. An alternatve slutn, already cnsdered n the lterature, s t use a gravty asssted maneuver (swng-by) wth the Mn t make the ut-f-plane part f the maneuver, as shwn n []. he swng-by s a type f maneuver that des nt use prpulsn based n engnes, but takes advantage f the gravtatnal nfluence f a celestal bdy when the spacecraft passes thrugh ts zne f nfluence [3-6]. he smaller bdy that s passng by the larger ne s accelerated r braked by the gravty feld f the celestal bdy and s ts trajectry s mdfed. Successful mssns, such as yager and [7], shwed the effcency f the swng-by n savng fuel. An example f the use f ths maneuver t change the rbtal nclnatn f a spacecraft s the Ulysses mssn, presented n [8], where the spacecraft was placed n a near plar rbt arund the Sun, n rder t bserve the ples f the Sun. hs paper ams t shw that, fr sme desred fnal nclnatns, a swng-by wth the Mn (n the Earth-Mn system) can be mre ecnmcal then usng plane change classcal maneuvers. he Cntent frm ths wrk may be used under the terms f the Creatve Cmmns Attrbutn 3. lcence. Any further dstrbutn f ths wrk must mantan attrbutn t the authr(s) and the ttle f the wrk, jurnal ctatn and DOI. Publshed under lcence by IOP Publshng Ltd

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 cntrbutn f the present paper s t analyze the nfluence f the sem-majr axs f the transfer rbt f the spacecraft, whch was a pnt nt cnsdered n prevus studes [].. Methdlgy he cmplete maneuver s dvded n three parts. In the frst phase an mpulse s used t put the satellte n a clse apprach trajectry wth the mn. In the secnd stage the three-dmensnal swngby s made and, n the last phase, tw mpulses are appled t take the satellte back t an rbt wth the same sem-majr axs and eccentrcty f the ntal rbt, but nw wth the fnal desred nclnatn. he analytcal mdel used s presented n [9], and t cnsders a system cnsstng f three bdes: a prmary (M, the Earth), a secndary (M, the Mn) and a thrd bdy f neglgble mass wth ts mtn gverned by the ther tw prevus cted bdes (M 3, the artfcal satellte). he three phases f the prblem s dmnated by the celestal mechancs f the tw bdes, s t s pssble t use the patched-cncs apprach. Intally, t s cnsdered that the satellte s n an rbt arund the Earth that s cplanar wth the rbt f the Mn. Fr the calculatns, t was assumed a canncal system f unts, where: () the gravtatnal cnstant s equal t ne; () the unt f dstances s the Earth-Mn dstance (3844.5 km); () the angular velcty f the prmares s equal t ne; (v) the mass f the mn s gven by µ L =m /(m +m ) (where m and m are the masses f the Earth and Mn, respectvely), and the mass f the Earth (µ ) s gven by µ = - µ L, s the ttal mass f the system s equal t ne; (v) the tme unt s set such that the perd f the prmares s π... Intal mpulse he mpulse used n the frst step t generate the transfer rbt (Δ ) and the eccentrcty f ths rbt (e ) are gven by: a e a a e e a e a () a where a s the sem-majr axs f the transfer rbt and ts mnmum value s gven by a = ( + a )/, where a s the sem-majr axs f the ntal rbt f the satellte arund the Earth and e s the eccentrcty f ths rbt... he Swng-by maneuver In the secnd step the calculatn f the swng-by starts when the satellte appraches the Mn and ts mtn s gverned by ths bdy. Fgure shws the gemetry f the swng-by, where: r ap s the dstance f the clsest apprach f the satellte wth the mn (perlune) and p s the velcty at ths pnt; the angles α and β specfy the perapss pstn n the three-dmensnal space and λ s the angle frmed between the vectr and the hrzntal plane passng by the perlune. Fgure shws the sum f the velcty vectrs durng the swng-by, where: and are the velcty vectrs wth respect t the mn befre and after the swng-by, respectvely, and the angle between them s δ, where δ s the s called deflectn angle; and are the satellte velcty vectrs befre and after the swng-by, respectvely, wth respect t the Earth; s the velcty vectr f the Mn wth respect t the Earth; γ s called flght-path-angle, whch s the angle between and ; ϕ s the angle between and. p ()

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 Fgure. Gemetry f the three-dmensnal Swngby. Fgure. he sum f vectrs durng the swng-by. Intally, t s calculated and the true anmaly (θ): (3) a cs a e e It s nw pssble t calculate γ and the magntude f the velcty f the satellte wth respect t the mn ( ): (4) tan e sn e cs cs (5) (6) Here ϕ and δ are calculated by: cs (7) sn r ap L (8) he apprach angle α and the angle λ are gven belw, respectvely, by: 3

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 he vectrs and are gven by [9]: (9) sn tan tan () sn cs cs,cs sn,sn cs ( sn sn cs cs sn, sn sn sn cs cs,cs sn ),, sn cs cs,cs sn,sn cs cs,cs sn cs ( sn sn cs cs sn, sn sn sn,, and allws the calculatn f the mpulse generated by the swng-by ( ), gven by [9]: sn cs cs,cs sn, sn sn he calculatn f the energy change due t the swng-by s gven by: () () (3) (4) E cs sn sn (5) he angular mmentum vectrs befre ( C ) and after the swng-by ( C ) are gven by [9]: C R d L cs sn sn cs sn sn sn ) C R dl,sn sn cs cs sn, (, sn sn cs cs sn, ( cs cs cs (6) cs cs cs (7) cs sn sn cs sn sn sn ) where R d L,, s the pstn vectr f the Mn wth respect t the Earth, and d L s the Earth- Mn dstance. hen, the magntude f the angular mmentum change (ΔC) s gven by [9]: C C d C sn d L cs sn sn sn,sn, cs sn (8) 4

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 [9]: After the swng-by maneuver s cmpleted t s pssble t calculate the fnal nclnatn ( ) by Cs C C Z sn sn cs cs sn cs cs cs cs sn sn cs sn sn sn (9).3. Fnal mpulses In the thrd step, after the swng-by s cmpleted, tw mpulses are used t return the satellte t ts ntal rbt (wth a and e ), but keepng the mdfed nclnatn. Frst, t s calculated the sem-majr axs f the satellte after the swng-by (a ): a () Next, t s calculated the sem-majr axs f the transfer rbt (a 3 ) and the frst mpulse (Δ ): a a 3 e a e a () () a3 e a a e When the satellte reaches the pergee f the transfer rbt ne last mpulse s used t crcularze the rbt (Δ 3 ): 3 a (3) 3 a e a a e After the cmplete maneuver, all mpulses are added: (4) t 3.4. Classcal Maneuver test the swng-by as an alternatve t the mre cnventnal types f plane change maneuver, ths technque s cmpared t a mn-mpulsve maneuver, gven by []: m a (5) sn e a 5

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 3. Results Accrdng t the ntal bjectves f the present paper and usng the prpsed technque, tw smulatns are made fr each value f the sem-majr axs f the transfer rbt: a (mnmum semmajr axs) and 5a. Fr bth smulatns, the ntal eccentrcty was e =, but n the frst smulatn (smulatn ) the ntal sem-majr axs was a =.7 and the perapss dstance r ap =.46 []. In the secnd smulatn (smulatn ) the ntal value f the sem-majr axs was a =.69 and fr the perapss dstance t was used agan r ap =.46. he ntal sem-majr axs and the perapss dstance are shwn n canncal unts (the Earth-Mn dstance). he smulatns are made fr β rangng frm t π and the results are shwn n fgures 3 t, dsplayng the fnal nclnatn ( ) n the vertcal axs as a functn f β n the hrzntal axs, bth n radans. Sme graphcs shw Δ t - Δ m, the dfference between the ttal mpulse used by the maneuver that ncluded the swng-by and the magntude f the mn-mpulsve maneuver, gven n canncal unts f velcty ( canncal unt f velcty s apprxmately.3 km/s), als as a functn f β n the hrzntal axs. Fgure 3. vs. β, fr smulatn cnsderng the sem-majr axs a. Fgure 4. Δ t - Δ m vs. β, fr smulatn cnsderng the sem-majr axs a. Fgure 5. vs. β, fr smulatn cnsderng the sem-majr axs 5a. Fgure 6. Δ t - Δ m vs. β, fr smulatn cnsderng the sem-majr axs 5a he results fr the slutn, cnsderng the sem-majr axs a, shwn n fgures 3 and 4, shw that the use f the swng-by s justfed nly fr the fnal nclnatns between.9 rad and.95 rad, whch crrespnds t β between.63 rad and.94 rad. Insde ths range the savngs n the magntude f the mpulses reaches.5 canncal unts. Hwever, the results gven by Slutn cnsderng the sem-majr axs 5a, as shwn n fgures 5 and 6, ndcates that the use f the swng-s justfed by the frst famly f slutns fr the fnal nclnatns between.9 rad and.9 rad. Fr values f β between.8 rad and.99 rad, the savngs n the magntude f the mpulse reaches values up t.86. 6

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 Nw, fr the secnd famly f slutns, the use f the swng-by s nly justfed fr values f β between.44 rad and.99 rad, wth fnal nclnatns between.9 rad and.3 rad, wth the savngs n the magntude f the mpulses up t. canncal unts. Fgure 7. vs. β, fr smulatn cnsderng the sem-majr axs a. Fgure 8. Δ t - Δ m vs. β, fr smulatn cnsderng the sem-majr axs a. Fgure 9. by β, smulatn fr 5a. Fgure. Δ t - Δ m by β, smulatn fr 5a Fgures 7 and 8 shw the results fr the slutn cnsderng the sem-majr axs as a. In ths case the use f the swng-by s nt justfed. he results f the slutn fr a sem-majr axs f 5a are shwn n fgures 9 and. It ndcates that the use f the swng-by s nly justfed fr the frst famly f slutns and fr the fnal nclnatns between.93 rad and.5 rad, whch crrespnds t β between.86 rad and.94 rad. he ecnmy n the magntude f the mpulse can reach up t.98 canncal unts. In general, the results are accrdng t the expected fr ths type f maneuver. here s a cst t send the spacecraft t the Mn and back, even n planar maneuvers. take advantage f ths maneuver, t s necessary t have savngs n the ut-f-plane maneuver that cmpensates ths extra cnsumptn. hs s the physcal reasns f the lmtatns n the regn f advantages f the prpsed technque. It can be appled when n the stuatns where the nclnatn change s larger, s the savngs n ths part f the maneuver cmpensates the extra csts nvlved. Althugh expected, the results shwn n the present paper quantfes the exact lcatns f the regns where the maneuver has gans ver the standard nes and hw much are the gans. In partcular, the present paper extends studes made befre by cnsderng the effects f the sem-majr f the transfer rbt n the savngs btaned. hs s a free mprtant parameter when desgnng a mssn. 7

XII Clóqu Brasler de Dnâmca Orbtal CBDO Jurnal f Physcs: Cnference Seres 64 (5) 4 IOP Publshng d:.88/74-6596/64//4 4. Cnclusn After analyzng the results, t s cncluded that, fr sme fnal nclnatns desred fr the rbts, the use f the swng-by maneuver as an alternatve t the classcal maneuvers can be a very gd ptn, generatng savngs even greater than canncal unts f velcty. Regardng the analyss f the nfluence f the sem-majr axs f the transfer rbt, the result shw that larger sem-majr axs allws the exstence f new famles f slutns, wth hgher fnal nclnatns. hs type f analyss can help mssn desgners t desgn mssns that can beneft frm the technque shwn here. 5. References [] Prad A F B A and Rs-Net A 993 Um estud bblgráfc sbre prblema de transferêncas de órbtas Revsta braslera de cêncas mecâncas 5 pp 65-78 [] rre K S and Prad A F B A 6 Changng nclnatn f earth satelltes usng the gravty f the mn Mat. Prb. Eng. 6 pp -3 [3] Gmes M and Prad A F B A A study f the mpact f the ntal energy n a clse apprach f a clud f partcles WSEAS rans. Mat. 9 pp 8-8 [4] Prad A F B A and Brucke R A 995 Effects f atmspherc drag n swng-by trajectry Acta Astr. 36 pp 85-9 [5] Gmes M and Prad A F B A 8 Swng-by maneuvers fr a clud f partcles wth planets f the slar system WSEAS rans. Ap. her. Mec. 3 pp 869-878 [6] Prad A F B A and Brucke R A 995 A classfcatn f swng-by trajectres usng the mn App. Mec. Rev. 48 pp 38-4 [7] Khlhase C E and Penz P A 977 yager mssn descrptn Sp. Sc. Rev. pp77- [8] Wenzel K P, Marsden R G, Page D E and Smth E J 99 he ULYSSES mssn Astr. and Astrp. Suppl. 9 pp 7-9 [9] Prad A F B A An analytcal descrptn f the clse apprach maneuver n three dmensns Prc. Int. Astr. Cngress (R de Janer) [] Ry A E 988 Orbtal mtn (Brstl) 8