Astrodynamics (AERO0024)

Transcription

1 Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L)

2 Motivation 2

3 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem 6.3 Gravity assist 3

4 6. Interplanetary Trajectories 6.1 Patched conic method Problem statement Sphere of influence Planetary departure Hohmann transfer Planetary arrival Sensitivity analysis and launch windows Examples 4

5 ? Hint #1: design the Earth-Mars transfer using known concepts Hint #2: division into simpler problems Hint #3: patched conic method 5

6 What Transfer Orbit? Constraints? Dep. V 1 Arr. Arr. V 2 Dep. Dep. Arr. Motion in the heliocentric reference frame 6

7 Planetary Departure? Constraints? v Earth / Sun v? V 1 Motion in the planetary reference frame 7

8 Planetary Arrival? Similar Reasoning Transfer ellipse SOI SOI Arrival hyperbola Departure hyperbola 8

9 Patched Conic Method Three conics to patch: 1. Outbound hyperbola (departure) 2. The Hohmann transfer ellipse (interplanetary travel) 3. Inbound hyperbola (arrival) Problem statement 9

10 Patched Conic Method Approximate method that analyzes a mission as a sequence of 2-body problems, with one body always being the spacecraft. If the spacecraft is close enough to one celestial body, the gravitational forces due to other planets can be neglected. The region inside of which the approximation is valid is called the sphere of influence (SOI) of the celestial body. If the spacecraft is not inside the SOI of a planet, it is considered to be in orbit around the sun Problem statement 10

11 Patched Conic Method Very useful for preliminary mission design (delta-v requirements and flight times). But actual mission design and execution employ the most accurate possible numerical integration techniques Problem statement 11

12 Sphere of Influence (SOI)? Let s assume that a spacecraft is within the Earth s SOI if the gravitational force due to Earth is larger than the gravitational force due to the sun. Gm m r Gm m E sat S sat 2 2 E, sat rs, sat r 5 E, sat km Sphere of influence 12

13 Sphere of Influence (SOI) Third body: sun or planet Spacecraft m 2 d m j r ρ r d ρ r Gm 3 j 3 3 r d O m 1 Central body: sun or planet Disturbing function (L04) Sphere of influence 13

14 If the Spacecraft Orbits the Planet r G m m r r p v sv sp pv r 3 pv Gms 3 3 r pv rsv rsp p:planet v: vehicle s:sun r A P pv p s Primary gravitational acceleration due to the planet Perturbation acceleration due to the sun Sphere of influence 14

15 If the Spacecraft Orbits the Sun r Gm m r Gm r r s v pv sp sv 3 sv p 3 3 r sv rpv rsp p:planet v: vehicle s:sun r A P sv s p Primary gravitational acceleration due to the sun Perturbation acceleration due to the planet Sphere of influence 15

16 SOI: Correct Definition due to Laplace It is the surface along which the magnitudes of the acceleration satisfy: P A p s P A s p Measure of the planet s influence on the orbit of the vehicle relative to the sun Measure of the deviation of the vehicle s orbit from the Keplerian orbit arising from the planet acting by itself r SOI m m p s 2 5 r sp Sphere of influence 16

17 SOI: Correct Definition due to Laplace If P A p s P A s p the spacecraft is inside the SOI of the planet. A A p The previous (incorrect) definition was 1 The moon lied outside the SOI and was in orbit about the sun like an asteroid! s Sphere of influence 17

18 SOI Radii Planet SOI Radius (km) SOI radius (body radii) Mercury 1.13x Venus 6.17x OK! Earth 9.24x Mars 5.74x Jupiter 4.83x Neptune 8.67x Sphere of influence 18

19 Validity of the Patched Conic Method The Earth s SOI is 145 Earth radii. This is extremely large compared to the size of the Earth: The velocity relative to the planet on an escape hyperbola is considered to be the hyperbolic excess velocity vector. vsoi v This is extremely small with respect to 1AU: During the elliptic transfer, the spacecraft is considered to be under the influence of the Sun s gravity only. In other words, it follows an unperturbed Keplerian orbit around the Sun Sphere of influence 19

20 Outbound Hyperbola The spacecraft necessarily escapes using a hyperbolic trajectory relative to the planet. Hyperbolic excess speed When this velocity vector is added to the planet s heliocentric velocity, the result is the spacecraft s heliocentric velocity on the interplanetary elliptic transfer orbit at the SOI in the solar system. Lecture 02: v v v v v esc SOI esc Is v SOI the velocity on the transfer orbit? Planetary departure 20

21 Magnitude of V SOI The velocity v D of the spacecraft relative to the sun is imposed by the Hohmann transfer (i.e., velocity on the transfer orbit). H. Curtis, Orbital Mechanics for Engineering Students, Elsevier Planetary departure 21

22 Magnitude of V SOI By subtracting the known value of the velocity v 1 of the planet relative to the sun, one obtains the hyperbolic excess speed on the Earth escape hyperbola. 2R sun 2 vsoi vd v1 1 v R 1 R1 R 2 Imposed Known Lecture Planetary departure 22

23 Direction of V SOI What should be the direction of v SOI? For a Hohmann transfer, it should be parallel to v Planetary departure 23

24 Parking Orbit A spacecraft is ordinary launched into an interplanetary trajectory from a circular parking orbit. Its radius equals the periapse radius r p of the departure hyperbola. H. Curtis, Orbital Mechanics for Engineering Students, Elsevier Planetary departure 24

25 ΔV Magnitude and Location a Lecture 02: known v 2 h e a r p known h 2 h (1 e) e v 2 1 e 1 h r v p rv p rp h v v p v c r 1 cos 1 p rp e Planetary departure 25

26 Planetary Departure: Graphically Departure to outer or inner planet? H. Curtis, Orbital Mechanics for Engineering Students, Elsevier Planetary departure 26

27 Circular, Coplanar Orbits for Most Planets Planet Inclination of the orbit to the ecliptic plane Eccentricity Mercury 7.00º Venus 3.39º Earth 0.00º Mars 1.85º Jupiter 1.30º Saturn 2.48º Uranus 0.77º Neptune 1.77º Pluto 17.16º Hohmann transfer 27

28 Governing Equations v D 2R sun 2 v1 R 1 R1 R2 1 v D v 1 v 2 sun 2R1 va 1 R R R v 2 v A Signs? v 2 - v A, v D - v 1 >0 for transfer to an outer planet v 2 - v A, v D - v 1 <0 for transfer to an inner planet Hohmann transfer 28

29 Schematically Transfer to outer planet Transfer to inner planet Hohmann transfer 29

30 Arrival at an Outer Planet For an outer planet, the spacecraft s heliocentric approach velocity v A is smaller in magnitude than that of the planet v 2. v v v 2 A v v 2 A v 2 v and have opposite signs Planetary arrival 30

31 Arrival at an Outer Planet The spacecraft crosses the forward portion of the SOI Planetary arrival 31

32 Enter into an Elliptic Orbit If the intent is to go into orbit around the planet, then must be chosen so that the v burn at periapse will occur at the correct altitude above the planet. r p 1 2 rv p 2 h (1 e) 2 (1 e) v vp, hyp vp, capture v r r r r 2 p p p p Planetary arrival 32

33 Planetary Flyby Otherwise, the specacraft will simply continue past periapse on a flyby trajectory exiting the SOI with the same relative speed v it entered but with the velocity vector rotated through the turn angle. e 1 rv p 2 2sin 1 1 e Planetary arrival 33

34 Sensitivity Analysis: Departure The maneuver occurs well within the SOI, which is just a point on the scale of the solar system. One may therefore ask what effects small errors in position and velocity (r p and v p ) at the maneuver point have on the trajectory (target radius R 2 of the heliocentric Hohmann transfer ellipse). 21 v R2 2 r 1 p rp v p 2 R2 Rv 1 D vdvrp rp vd vp 1 2 sun Sensitivity analysis and launch windows 34

35 Sensitivity Analysis: Earth-Mars, 300km Orbit sun R km, R km, r 6678 km v km / s, v km / s D km / s, km / s 2 R r 2 p v R r v p A 0.01% variation in the burnout speed v p changes the target radius by 0.067% or km. A 0.01% variation in burnout radius r p (670 m!) produces an error over km. p p p Sensitivity analysis and launch windows 35

36 Sensitivity Analysis: Launch Errors Ariane V Trajectory correction maneuvers are clearly mandatory Sensitivity analysis and launch windows 36

37 Sensitivity Analysis: Arrival The heliocentric velocity of Mars in its orbit is roughly 24km/s. If an orbit injection were planned to occur at a 500 km periapsis height, a spacecraft arriving even 10s late at Mars would likely enter the atmosphere Sensitivity analysis and launch windows 37

38 Cassini-Huygens

39 Existence of Launch Windows Phasing maneuvers are not practical due to the large periods of the heliocentric orbits. The planet should arrive at the apse line of the transfer ellipse at the same time the spacecraft does Sensitivity analysis and launch windows 39

40 Rendez-vous Opportunities nt nt 2 n n Tsyn n n t T syn 2 n n 1 2 T syn TT 1 2 T T 1 2 Synodic period Sensitivity analysis and launch windows 40

41 Transfer Time Lecture 02 t sun R R 2 3/ 2 nt Sensitivity analysis and launch windows 41

42 Earth-Mars Example T syn days It takes 2.13 years for a given configuration of Mars relative to the Earth to occur again. t s days 44 The total time for a manned Mars mission is = days = 2.66 years Sensitivity analysis and launch windows 42

43 Earth-Mars Example 1. In 258 days, Mars travels 258/688*360=135 degrees. Mars should be ahead of 45 degrees. 2. In 258 days, the Earth travels 258/365*360=255 degrees. At Mars arrival, the Earth is 75 degrees ahead of Mars. 3. At Mars departure, the Earth should be behind Mars of 75 degrees. 4. A return is possible if the Earth wins =210 degrees w.r.t. Mars. The Earth wins 360/ /688=0.463 degrees per day. So one has to wait 210/0.46=453 days Sensitivity analysis and launch windows 43

44 Earth-Jupiter Example: Hohmann Galileo s original mission was designed to use a direct Hohmann transfer, but following the loss of Challenger Galileo's intended Centaur booster rocket was no longer allowed to fly on Shuttles. Using a lesspowerful solid booster rocket instead, Galileo used gravity assists instead Examples 44

45 Earth-Jupiter Example: Hohmann Velocity when leaving Earth's SOI: v v v D 2R E sun 2 1 R 1 R1 R km/s Velocity relative to Jupiter at Jupiter's SOI: J sun 2R v2 va v R 2 R1 R km/s Transfer time: years Examples 45

46 Earth-Jupiter Example: Departure Velocity on a circular parking orbit (300km): v c E R E 7.726km/s h 2 2 v v km/s km/s r p e 2 rv p Examples 46

47 Earth-Jupiter Example: Arrival Final orbit is circular with radius=6r 2 2 (1 e) v v km/s r r p p J e= Examples 47

48 Hohmann Transfer: Other Planets Planet Mercury Venus Mars Jupiter Saturn Pluto v departure (km/s) Transfer time (days) Assumption of circular, co-planar orbits and tangential burns Examples 48

49 Venus Express: A Hohmann-Like Transfer C 3 = 7.8 km 2 /s 2 Time: 154 days Real data Why? C 3 = 6.25 km 2 /s 2 Time: 146 days Hohmann Examples 49

50 6.1.7 Examples

51 6. Interplanetary Trajectories 6.2 Lambert s problem 51

52 Motivation Section 6.1 discussed Hohmann interplanetary transfers, which are optimal with respect to fuel consumption. Why should we consider nontangential burns (i.e., non- Hohmann transfer)? L Lambert s problem 52

53 Non-Hohmann Trajectories Solution using Lambert s theorem (Lecture 05): If two position vectors and the time of flight are known, then the orbit can be fully determined. 6.2 Lambert s problem 53

54 Venus Express Example 6.2 Lambert s problem 54

55 Porkchop Plot: Visual Design Tool C3 contours C3 contours Arrival date Departure date In porkchop plots, orbits are considered to be non-coplanar and elliptic. 6.2 Lambert s problem 55

56 6.2 Lambert s problem

57 Type II transfer for cargo: the spacecraft travels more than a 180 true anomaly Type I transfer for piloted: the spacecraft travels less than a 180 true anomaly

58

59 6.2 Lambert s problem

60 6.2 Lambert s problem

61 6. Interplanetary Trajectories 6.3 Gravity assist Basic principle Real-life examples 61

62 ΔV Budget: Earth Departure SOYUZ Planet Mercury Venus Mars Jupiter Saturn Pluto C 3 (km 2 /s 2 ) [56.25] [139.24] Assumption of circular, co-planar orbits and tangential burns 62

63 ΔV Budget: Arrival at the Planet A spacecraft traveling to an inner planet is accelerated by the Sun's gravity to a speed notably greater than the orbital speed of that destination planet. If the spacecraft is to be inserted into orbit about that inner planet, then there must be a mechanism to slow the spacecraft. Likewise, a spacecraft traveling to an outer planet is decelerated by the Sun's gravity to a speed far less than the orbital speed of that outer planet. Thus there must be a mechanism to accelerate the spacecraft Basic principle 63

64 Prohibitive ΔV Budget? Use Gravity Assist Also known as planetary flyby trajectory, slingshot maneuver and swingby trajectory. Useful in interplanetary missions to obtain a velocity change without expending propellant. This free velocity change is provided by the gravitational field of the flyby planet and can be used to lower the Δv cost of a mission Basic principle 64

65 What Do We Gain? Spacecraft outbound velocity SOI Spacecraft inbound velocity V out = V in Basic principle 65

66 Gravity Assist in the Heliocentric Frame Resultant V out SOI Resultant V in Planet s sun relative velocity Δv v v, out, in Basic principle 66

67 A Gravity Assist Looks Like an Elastic Collision Inertial frame Frame attached to the train Frame attached to the train Inertial frame Basic principle 67

68 Leading-Side Planetary Flyby A leading-side flyby results in a decrease in the spacecraft s heliocentric speed (e.g., Mariner 10 and Messenger) Basic principle 68

69 Trailing-Side Planetary Flyby A trailing-side flyby results in an increase in the spacecraft s heliocentric speed (e.g., Voyager and Cassini-Huygens) Basic principle 69

70 What Are the Limitations? Launch windows may be rare (e.g., Voyager). Presence of an atmosphere (the closer the spacecraft can get, the more boost it gets). Encounter different planets with different (possibly harsh) environments. What about flight time? Basic principle 70

71 V V E J G A See Lecture 1 71

72 Rosetta Real-life examples 72

73 Messenger Real-life examples 73

74 6.3.2 Real-life examples

75 Hohmann Transfer vs. Gravity Assist Gravity assist Planet C3 (km 2 /s 2 ) Transfer time (days) Real mission C3 (km 2 /s 2 ) Transfer time (days) Mercury [56.25] 105 Messenger Cassini Saturn Huygens Remark: the comparison between the transfer times is difficult, because it depends on the target orbit. The transfer time for gravity assist mission is the time elapsed between departure at the Earth and first arrival at the planet Real-life examples 75

76 6. Interplanetary Trajectories 6.1 PATCHED CONIC METHOD Problem statement Sphere of influence Planetary departure Hohmann transfer Planetary arrival Sensitivity analysis and launch windows Examples 6.2 LAMBERT S PROBLEM 6.3 GRAVITY ASSIST Basic principle Real-life examples 76

77 Even More Complex Trajectories _Francais/Reveil_du_satellite_Rosetta_dans_moins_de_45_jours 77

78 Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L)