Interplanetary Mission Analysis

Size: px
Start display at page:

Download "Interplanetary Mission Analysis"

Transcription

1 Interplanetary Mission Analysis Stephen Kemble Senior Expert EADS Astrium Page 1

2 Contents 1. Conventional mission design. Advanced mission design options Page

3 1. Conventional mission design 1.1 Interplanetary transfers and Lambert s problem 1. Spheres of influence and patch conics 1.3 Planetary escape and capture Page 3

4 .1 Interplanetary transfers and Lambert s problem Page 4

5 Planet to planet transfers Large energy changes are required to reach another planet Energy is sum of kinetic energy (positive) and gravity potential (negative) Rendez-vous with Mercury requires the greatest energy change Venus and Mars are far less demanding Plot shows energy relative to Earth Specific energy is used (ie per unit spacecraft mass) Energy (J) 6.00E E+08.00E E E E E E+08 Mercury Venus Earth Mars Jupiter Saturn Uranus Neptune Pluto V µ E = r Page 5

6 Properties of planetary orbits sma (AU) sma (km) Energy Period (days) Eccentricity Mu Energy rel Earth Mercury E E E E+08 Venus E E E E+08 Earth E E E E+00 Mars E E E E+08 Jupiter E E E E+08 Saturn E E E E+08 Uranus E E E E+08 Neptune E E E E+08 Pluto E E E E+08 Page 6 Rmin V at Rmin Rmax V at Rmax (km) (km/s) (km) Mercury Venus Earth Mars Jupiter Saturn Uranus Neptune Pluto

7 Transfer between circular orbits Planets may be approximated to lie in circular, co-planar orbits,orbit 1, Earth, and orbit, the target The optimum transfer is then a Hohmann transfer Requires two impulses One leaving Earth Second impulse can be used to rendez-vous with the target. V ap r r 1 V pe E 1 = µ r 1 E µ = r E µ = r 1 + r V µ = + E r Page 7 µ µ V1 = + E + E 1 rp 1 rp 1 µ µ V + = + E E rp rp

8 General transfers In practice: Planet orbits are eccentric Each has a different inclination and so are not co-planar r c θ r 1 The general transfer problem is Lambert s problem Lambert proposed that the time of flight depends upon three quantities 1. The semi-major axis of the connecting ellipse. The chord length, c 3. The sum of the position radii from the focus or the connecting ellipse. Find a transfer between two positions in a given time. s a r1 + r + c = c = r1 + r r cosθ 1 r Page 8

9 Lambert s problem Difference in Mean anomalies of transfer orbit M µ a ( t t ) = ( t ) M 1 = n 1 3 t1 ( ) ( ( ) ) t t1 = E E1 e sin E sin E1 = e sin cos µ µ Substitute: Then can derive Substitute Can derive Gives the result Solution depends only on semimajor axis a Page 9 a 3 A = α = a 3 E E 1 E E A + B α sin = E E E cos cos + E B = e 1 3 a µ ( t t ) = ( A sin Acos B) 1 s a 3 a µ β = B A β sin = s c a E + E ( t t ) = (( α sin α ) ( β sin β )) 1

10 Lambert s problem () Can be applied to hyperbolic orbits Define further variables: α sinh = s a sinh β = s c a Gives the result for hyperbolic orbits Solution depends only on a 3 a µ ( t t ) = (( sinh α α ) ( sinh β β )) 1 In addition to semi-major axis, eccentricity, e can be found, for elliptical or hyperbolic cases p = ar r 4 1 β sin c θ sin α + p = a(1 e ) p = ar r 4 1 β sin c θ sinh α + Page 10

11 Solving Lamberts problem A transfer from Earth to Mars can be considered. Consider leaving Earth on 3rd October, 011 Transfer durations of 150, 180 and 00 days investigated. Departure and arrival epochs specify two radius vectors, from which c and s may be calculated. An estimate at the value of semimajor axis solving the equation can be made and substituted into the right hand side the equation: A plot of error between left and right hand sides versus semi-major axis can be obtained. The solution lies where the error equals zero, ie crosses the axis error (secs) E E+11.00E E E E sma (m) 3 a µ ( t t ) = (( α sinα ) ( β sin β )) Page 11

12 Solving Lamberts problem () Hyperbolic transfer The spacecraft again leaves Earth on 3rd October, 011 now the transfer duration is only 100 days. The departure and arrival epochs are once again specify the two radius vectors, from which the quantities c and s may be calculated. Now the transfer duration is so short that the connecting orbit must be hyperbolic In this case the semi-major axis at the solution is negative, indicating a solar system escape orbit. error (secs) E E E E E E sma (m) Page 1

13 Transfers to Mars The non-coplanar nature of the orbits is an important aspect The plane change requirement implied results in a ridge in the contour plots Results in the classical pork chop plot. Each synodic period now shows two local minima: short and long transfer types, typically 00 and 350 days durations Page Vinfinity Earth-Mars (m/s) Transfer time (days) Launch epoch (MJD)

14 Nature of the transfers A particular launch date can be selected The effect of varying transfer duration at that epoch can be examined Use 3 rd October 011. Two minima are expected as the transfer duration evolves. The short and long minima are seen separated by a local maximum at approximately 30 days Vinf total (m/s) Duration (days) Page 14

15 Nature of the transfers Examining the trajectories reveals the nature of the local maximum As the transfer approaches a conjunction type transfer, it becomes increasingly difficult to achieve the out of plane component of the rendez-vous position. The solution is to increase the heliocentric inclination of the transfer orbit. The solutions result in transfers switching from South to North around the local maximum Page 15

16 The synodic period The synodic period is that between the repetition of a particular relative orbit geometry between the planets in question, such as a particular difference in solar longitude. Such repetitions occur at fixed intervals for two circular, co-planar planet orbits. If the orbits are assumed circular, then the synodic period is calculated as: 360 τ = Where τ p1 is the orbital period of planet 1 and τ p is the orbital period of planet and ω p1 is the orbital period of planet 1 and ω p is the orbital period of planet Page 16 S = τ p1 τ p ( ω ω ) p1 360 p

17 The synodic period The synodic period for the planets, assuming circular orbits are: Mercury Venus Earth Mars Jupiter Saturn Uranus Neptune Mercury Venus Earth Mars Jupiter Saturn Uranus Neptune Pluto Page 17

18 Global repeat periods The absolute locations of the planets do not repeat at these intervals, only the relative locations repeat. The time between an exact repetition of an absolute transfer geometry is calculated by: τ = nτ = global S mτ p1 The time is the number of synodic periods required to generate a whole number of orbits of the departure planet, whose period is τ p1. This may in practice be a very long period for a precise repeat It can be approximated, if the repeat geometry departure longitude lies within a few degrees of the original departure longitude. Page 18

19 Vinfinity Earth-Mars (m/s) Examples of global repeats Transfer opportunities from Earth show a global repeat to the nearest planets as follows: Venus: Synodic period 1.6 years, Global repeat 8 years Mars: Synodic period.15 years, Global repeats near 13 and 15 years Jupiter: Synodic period 1.09 years, Global repeats near 1 years Mercury: Synodic period 0.3 years, Global repeats near 1 year Launch epoch (MJD) Page Transfer time (days)

20 Multi-revolution transfers Many locally optimal transfers are possible when the upper limit on transfer durations is extended. This allows an increased number of heliocentric revolutions. Instead of 0.5 revolutions, it is possible to use 1.5 or in principle, n.5 revolutions. For the example of co-planer, circular planet orbits, the minimum speed change required of the spacecraft is independent of the number of revolutions. The optimal transfer duration is simply incremented by n times the transfer orbit period. Page 0

21 Multi-revolution transfers () Lambert s problem now has two solutions for the semi-major axis, when n>0 Page 1

22 Multirevolution transfers to Mars Example of 1.5 revolution case from Earth to Mars Launch dates over two synodic periods are shown Minimum Vinfinities are comparable to 0.5 rev case Launch and arrival dates differ significantly to 0.5 revs Page Vinfinity Earth-Mars (m/s) Transfer time (days) Launch epoch (MJD)

23 . Spheres of influence and patch conics Page 3

24 Spheres of influence (1) Consider two major bodies (Sun and planet) Consider the ratio of the perturbing acceleration (from the second body) to the main acceleration (from the first body). Ratio can be evaluated for both the Sun s perturbation on motion about the planet and the planet s perturbation on motion about the Sun. r cs r ps a r pc A surface may then be found defined by r ps where these ratios are equal, in terms of the angular offset from the Sun s direction, a r pc is the distance from planet to Sun Page 4 r ps = r pc µ µ p c 1 ( 1+ 3cos a)

25 Spheres of influence () Hill s sphere of influence can also be considered: includes rotational effect similar to L1 and L distances The radii of the different spheres differ Values for the planets of the solar system differ significantly Classical sphere of Hill's sphere radius Influence Mercury Venus Earth Mars Jupiter Saturn Uranus Neptune Pluto Page 5

26 Using spheres of influence Consider a mission from Earth to Mars Keplerian orbits only consider a single gravity field Initially the spacecraft has a Keplerian orbit calculated relative to Earth: Earth s gravity dominates Generally hyperbolic After leaving the Earth s sphere of influence the spacecraft motion can be defined in terms of a Keplerian orbit about the Sun Generally elliptical Finally on reaching Mars the spacecraft has a Keplerian orbit calculated relative to Mars: Martian gravity dominates Generally hyperbolic 3 distinct orbits are used to describe the trajectory. The limitations of applicability are indicated by the sphere of influence of Earth and Mars Method is called patched conics Hyperbola wrt planet Velocity of planet Velocity of planet 1 Hyperbola wrt planet 1 Page 6

27 .3 Planetary escape and Capture Page 7

28 Escaping from a planet A hyperbolic orbit about a planet has a non-zero speed at infinite distance from the planet (if no other gravity field than the planet s were considered). When at several millions of kilometres, its speed relative to the planet tends towards the Excess hyperbolic velocity. The excess velocity targeted is that needed to achieve the interplanetary transfer, ie the heliocentric speed increment relative the planet This is approximately equal to the effective departure speed from the planet within the heliocentric domain. This quantity can be calculated from: V = µ a Page 8

29 Escaping from a planet The direction of the excess hyperbolic departure vector can be chosen by Choosing the appropriate orbit plane Choosing the appropriate pericentre of the initial hyperbolic orbit The asymptotic departure direction within the orbit plane is given by: where θ is the maximum true anomaly in the hyperbolic orbit, ie the velocity vector is asymptotically aligned with the radius vector. The DeltaV needed to reach this orbit is for a general elliptical orbit with apogee, θ cos 1 = 1 e Page 9 a p1 and perigee radius r pl1 µ 1 V = + V µ rpl rpl ( r 1 + r 1 1 pl1 apl1 )

30 Escaping from a planet The DeltaV needed to reach a given Vinfinity increases more quickly with target Vinfinity as Vinfinity increases In these examples the perigee altitude is 00km DV (m/s) Declinations that can be reached from a given inclination orbit are given by: sin Therefore maximum abs(declination) = inclination Apogee (km) DEC = sin isin( ω + θ ) Page 30

31 Approach to a planet The approach to a planet is characterised by three parameters: Excess hyperbolic speed (magnitude of the approaching, planet relative velocity vector,) Right Ascension and Declination of the of the approaching, planet relative velocity vector. The capture manoeuvre is a large retro-manoeuvre performed at planet pericentre It is the reverse of an escape manoeuvre. Evaluation of the location of the possible pericentre is required to assess the range of possible capture orbits that can be achieved without plane change. Page 31

32 Approach definitions Three planes are of primary importance Reference plane (eg ecliptic parallel) B plane: perpendicular to approach velocity vector Orbit plane: determined by selection of the Beta angle: locates the position of the pericentre of the orbit Approaching relative velocity vector at Beta angle > 90 deg B plane perpendicular to approach velocity vector Declination (negative as shown) Planet Relative orbit plane β angle measured from reference plane Plane parallel to Ecliptic plane: Reference plane Relative Approach plane Intersection of B plane and ecliptic parallel plane. Defines zero Beta angle Page 3

33 Capture pericentre The asymptotic approach direction within the orbit plane with respect to pericentre is given by the same relationship as the planetary departure case Two pericentre solutions exist for each approach orbit plane considered (as in the departure case) The approach plane is defined by the Beta angle A 360 degree range of Beta angles can be considered A limit exists on the inclination of the approach orbit that exists sin i sin DEC Page 33

34 Capture pericentre () Pericentre locations with different approach speeds at Mars 80 Pericentre locations with different approach speeds of 3 km/sec at Mars and a range of approach declinations Latitude Latitude Longitude -60 Page Longitude

35 Capture orbits example The family of possible capture orbits at Mars are shown, resulting from a typical interplanetary trajectory and approach to Mars Capture orbit has a high apocentre (circa km) Page 35

36 . Advanced mission design options.1 The three body problem. Designing missions to the Lagrange points.3 Gravity assist.4 Multiple Gravity Assists.5 Plane changing gravity assist.6 Gravitational escape and capture Page 36

37 .1 The three body problem Page 37

38 The three body problem The three body problem considers the motion of a spacecraft through a combined gravity field. The first body is the central body, the second is the orbiting, major body and the third is the spacecraft itself, whose mass is considered to be negligible in comparison to the other two. No analytical solutions exist that are comparable to Keplerian orbits for the two body case. Generalisations that describe the overall properties of the motion in certain circumstances exist. Also approximation methods enable certain classes of motion to be described by analytical expressions. The precise motion can only be obtained by numerical integration of the spacecraft state vector derivatives. Page 38

39 The Jacobi constant A constant of motion can be derived for a special case of the three body problem This is known as the circular, restricted three body problem The two major bodies are assumed to move in circular orbits about a common barycentre A good approximation for many planetary or moon orbits The constant is the Jacobi constant + U C Where V is the speed of the spacecraft relative to a rotating reference frame U is a generalised potential given by: U µ 1 µ r x r Y r1 r ω = + Where rx,ry, rz are the components of position with respect to the barycentre of the two major bodies, r1 and r are the distances from body 1 and body, is the angular velocity of the two bodies in circular orbit µ1 = Gm1, ie the gravitational parameter for mass m1 µ = Gm, ie the gravitational parameter for mass m V = ( ) Page 39

40 Zero velocity surfaces The limiting regions of the motion can be obtained by considering the case where speed is zero: µ 1 µ ω rx + ry + C = r1 r Zero velocity surfaces show limitations of the possible motion. Can evaluate contours of zero velocity, for different values of the constant, C. Only motion in the ecliptic is considered (ie rz=0) in the following examples Therefore this is now a case of a planar, circular, restricted three body problem. Motion is excluded within a particular shaded area that corresponds to a given value of C. X lies along the Sun to planet vector and Y (vertical axis, in km) lies perpendicular to X in the ecliptic. ( ) 0 Page 40

41 Earth zero velocity surface Earth lies on the X axis at 1.496e8 km. Series of constant velocity surfaces evaluated over a range of values of the Jacobi constant, C, starting with: CMaximum = (m/s)^ (Cmaximum) is a case where escape from Earth is not possible. Further surfaces may be generated with reducing C (ie C increasing): no longer a barrier to an initially Earth neighbouring trajectory transferring into the heliocentric domain. The contours in the figure represents the absolute value of the reduction in the Jacobi constant below the maximum (Cmaximum). Therefore a large value for the contour means that C now lies significantly below Cmaximum As C is reduced inaccessible regions shrink, until motion becomes possible through the locations of the first two colinear Lagrange points, L1 and L, at approximately +/- 1.5 million km along the X axis from Earth Page 41

42 The Lagrange libration points Points exist where a particle, if placed there, with no velocity in the rotating frame, will experience no resultant acceleration with respect to this rotating frame. Condition is: U = 0 5 points can be found in the orbit plane of the secondary body about the primary There are three, collinear, points, denoted L1, L and L3 Two further points complete an equilateral triangle with the baseline from primary to secondary, L4 and L5 L 3 L 5 d 1 d L 4 L 1 L Page 4

43 Lagrange points at the planets The L1 and L points are shown for each planet, using the approximation of the circular restricted three body problem Assumes that each planet lies in a circular orbit about the Sun at a value corresponding to its semi-major axis. The semi-major axis used is that from the JPL mean ephemeris at J000 Page 43 m/(m1+m) L1 (km) L (km) Mercury 1.660E-07.04E+05.10E+05 Venus.448E E E+06 Earth 3.003E E E+06 Earth-Moon 3.040E E E+06 Mars 3.7E E E+06 Jupiter 9.537E E E+07 Saturn.857E E E+07 Uranus 4.366E E E+07 Neptune 5.150E E E+08 Pluto 6.607E E E+06

44 Orbits at the Lagrange points The collinear points are unstable (stable in transverse direction, unstable in radial direction) Assuming a linearisation of the motion allows the derivation of differential equations describing possible motions (similar to Hill s equations for perturbed circular orbits about a single body) Local displacements relative to the Lagrange point are: δr = r Scaled by dividing by the separation between two bodies C11 is a constant d δr' dt' d δr' dt' d δr' dt' x y z dδr' dt' dδr' + dt' = δr' z y x c 11 = δr' = δr' x y ( 1+ c ) ( 1 c ) r Ln Page 44

45 Orbits at the Lagrange points () These differential equations may be solved: δr' δr' δr' x y z = A e 1 1 st' + A z e st' st' = aa1e aae = C cosω t' + C st' + A 3 ba sinω t' cosω t' + A z 4 xy 4 cosω t' + ba xy sinω t' 3 xy sinω t' xy A1 to 4 are constants determined by the initial conditions (also C1 and C) The two angular frequencies differ slightly, depend on c11 Ratios of motion amplitudes in x and y (a and b) also depend on c11 Page 45

46 Orbit properties at the Lagrange points Initial solutions may be found that remove the exponential components The relative amplitudes of the oscillating motion in the rotating x and y directions is determined by the constant, b. Orbit characteristics can be determined for each planet assuming circular orbits Units of angular velocity are relative to orbit angular frequency Page 46 µ wxy wz Txy Tz b (days) (days) Mercury 1.660E Venus.448E Earth 3.003E Mars 3.7E Jupiter 9.537E Saturn.857E Uranus 4.366E Neptune 5.150E Pluto 6.607E

47 Lissajous orbits The solution of these equations is a Lissajous figure. It is not a closed orbit Motion projection in xy is closed Initial conditions determine the initial phasing of the motions Figures show examples The origin is a collinear Lagrange point Z Y Y X Page 47

48 Halo orbits The full solution regarding motion about the Lagrange points shows a dependence of orbital frequency on amplitude The dependence differs between in and out of plane components The result is that critical ration of in/out of plane amplitudes yield solutions with identical frequencies This allows closed orbits to be achieved Choosing the appropriate phasing (via initial conditions) allows a Halo orbit to be achieved An example of such a Halo orbit was adopted for ISEE3, with in plane amplitude of km and out of plane amplitude km Az Page Ay

49 Stable and Unstable motions The solutions indicate the possible presence of exponential, time dependent terms. Can be suppressed by suitable selection of the initial states such that the constants A1 and A are zero. If the constant, A1 takes a non-zero value, then an exponentially increasing, time dependent term exists. The result is that a perturbation to an oscillatory solution yielding non-zero A1 will eventually lead to a complete departure from that orbit. Such motions are unstable. This type of perturbation process is one of stepping onto the unstable manifold of the orbit at the Lagrange point. If the constant, A takes a non-zero value, then an exponentially decreasing, time dependent term exists. The result is that a perturbation to an oscillatory solution yielding non-zero A will eventually lead to the perturbation reducing to zero and returning to the oscillatory solution. Such motions are stable. This type of perturbation process is one of stepping onto the stable manifold of the orbit at the Lagrange point Page 49

50 Stable and Unstable motions () The stable manifold may be used to execute transfers to reach the oscillatory solution. Start at a point on the stable manifold far removed from the target orbit. Evolution of the trajectory with time then results in the spacecraft approaching the orbit described by the oscillatory solution, as the exponential term tends to zero. The conditions for the constant A1 remaining at zero in the presence of a velocity perturbation to an oscillatory solution can be found in terms of a direction in which this velocity perturbation may be applied. This direction lies in the ecliptic plane. A perturbation that has a component that is perpendicular to this direction in the ecliptic leads to an unstable solution where A1 is nonzero. This fact may be used as a feature of orbit generation and maintenance strategies: Station keeping aims to remove the unstable component Page 50

51 . Designing missions to the Lagrange points Page 51

52 Transfers Transfers from an initial Earth elliptical orbit to a Lissajous or Halo orbit about the Earth-Sun L1 or L Lagrange points can be achieved without manoeuvre An example can be considered as follows: Initial perigee altitude of 500km and semi-major axis (at Earth perigee) of km. Osculating apogee (at perigee) of nearly 1.4 million km. Solar gravity perturbs the motion as move towards apogee. T Perigee is raised and the energy and angular momentum with respect to Earth are modified. The nature of this perturbation depends on the location of the apogee with respect to the Earth-Sun direction. Particular orientations can be found that enable freely reaching a Lissajous orbit about the L1 Lagrange point. Page 5

53 Searching for a transfer (1) The orientation of the initial orbit can be obtained by search in the longitude of the line of apses First choose apogee in 1.3 to 1.5 million km range Choose initial longitude to be within, for example, 0 degrees of Sun direction (for transfers to an L1 orbit) or anti-sun direction (for transfers to an L orbit) The figure, seen in a rotating reference frame, shows the effects of a 60 deg variation in perigee longitude. The behaviour includes orbits that return to Earth perigee; enter Lissajous orbits and also a series of orbits that escape Earth s influence Page 53

54 Searching for a transfer () Having established approximately the required longitude, a much more refined search is required Rule based methods can be used to automatically find such transfers Example shows a variation in longitude mof 0.05 degrees about the required transfer solution Sensitivity to initial conditions is therefore very high Page 54

55 Target orbit implications The minimum Lissajous orbit amplitude (in ecliptic) depends on the initial perigee altitude For a typical transfer orbit with perigee at 500km, minimum amplitude is km Greater amplitudes can be achieved by appropriate choice of apogee radius and apse longitude Higher enable smaller amplitudes, eg perigee at 36000km altitude can free inject to km amplitude orbit If a lower amplitude than the minimum free injection case is sought then a transfer manoeuvre can be performed Example assumes an initial 500km perigee altitude DeltaV (m/s) Halo SMA (km) Page 55

56 Transfer characteristics The transfer duration is typically 10 days until entering the Lissajous orbit However the vicinity of the Lagrange point is reached after approximately 30 days In the case of a free injection transfer, final dispersion correction manoeuvres can be performed after this time For injection to smaller amplitude orbits the manoeuvre is made after typically days Page 56

57 Transfer types (1) First consider the motion in the ecliptic Two longitudes give free injection transfers for a given choice of apogee radius (and the same perigee radius) Post injection orbits are similar in amplitude Transfers take slightly different durations Page 57

58 Transfer types () Out of plane motion of the target Lissajous orbit is determined by selecting an appropriate declination of the perigee departure vector (ie determined by the declination of the line of apses) 4 solutions exist for a given out of plane amplitude and left or right in ecliptic transfer type routes transfer to the North routes transfer to the South 1 North and 1 South transfer examples are shown in rotating frame, XZ plane Z rotating (km) Y rotating (km) Z rotating (km) E E E+09-1.E+09-1E+09-8E+08-6E+08-4E+08 -E+08 0 E X rotating (km) E E E+09-1.E+09-1E+09-8E+08-6E+08-4E+08 -E+08 0 E X rotating (km) E+09-8E+08-6E+08-4E+08 -E+08 0 E+08 4E+08 6E+08 8E+08 1E Type Type 1 Type Type 1 Type Type Page Y rotating (km)

59 .3 Gravity Assist Page 59

60 Gravity assist (1) Gravity assist manoeuvres enable an exchange of energy between a planet and the spacecraft A spacecraft approaches a planet from deep space. Hyperbolic orbit approach They rotate about their common barycentre The spacecraft departs with a modified orbit relative to the Sun Analysis is made using patched conics Hyperbolic orbit with respect to major body Sphere of influence of major body Departure orbit Approach orbit Page 60

61 Patch Conics Three conic sections are generated to describe the orbital phases. Phase 1: The approach orbit, Expressed in terms relative to the central body. Phase : The fly-by orbit Expressed in terms relative to the major body. This must be a hyperbolic orbit Phase 3: The departure orbit Expressed in terms relative to the central body Page 61

62 Patch and Link Conics Selection of the radius of the patch sphere is a key aspect The approximation generates a small velocity error compared to the full 3 body solution Typically the classical radius of the Laplace sphere of influence is selected A second approximation is that of Link conics With the link conic method, the fly-by, or gravity assist, is specified in a similar way, generally by ephemeris with respect to the fly-by body. From this ephemeris, it is possible to evaluate the asymptotic approach and departure velocity vectors, ie the relative velocity vector at infinite distance from the planet. Page 6 Limiting true anomaly,, depends only on eccentricity, e θ ± cos 1 = 1 e

63 Deflection The important feature of a fly-by is that the spacecraft s velocity w.r.t the planet is deflected The deflection only depends on eccentricity It is related to the limiting true anomaly The eccentricity is dependent in the choice of pericentre and the Excess hyperbolic speed α = *sin e = 1 r 1 perirel a 1 e = 1+ r perirel µ Deflection angle, α * V planet True anomaly of asymptotic departure direction, θ Page 63

64 Link approximations The link conic approximation therefore assumes that the velocity vector change is instantaneous The change takes place at the planet Therefore the time to traverse the sphere of influence is neglected Typically -3 days for the inner planets This is generally small compared to the periods of the interplanetary orbits It is equivalent to the approximation of using excess hyperbolic speed to model the departure or arrival to a planet Page 64

65 Gravity assist velocity change The gravity assist effect can be described by velocity vector triangles Consider the approach situation: Spacecraft velocity Planet velocity Relative velocity Relative velocity is deflected as the spacecraft swings by the planet Then add the velocity vectors V departure Approach V Γ V p Deflected V rel =V Γ departure V p α Body Approach V rel =V Γ relp Departure Γ relp Page 65

66 Alternative coplanar solutions Deflection can take place in one of two possible directions, depending on which side of the planet the spacecraft approaches. Determines the sign of alpha (shown as positive in the previous figure). Choice of solution can result in a considerable difference in the post gravity assist heliocentric state Page 66

67 Gravity assist effects Gravity assist can modify the interplanetary orbit in a variety of ways If co-planar heliocentric orbits are assumed (spacecraft and planet) with fly-by in the plane of the orbits then the spacecraft may fly around either side of the planet by making small targeting manoeuvres when distant from the planet. First the approach vector flight path angle is found V is the initial heliocentric velocity V *sin( Γ Γ planet ) V planet is the velocity of the planet Γ = tan 1 relp V * cos( Γ Γ planet ) V planet The departure flight path is then The departure heliocentric velocity is Page 67 V departure V = V + V V V cos( Γ Γ = V rel planet planet planet ) Γ departure = tan 1 V V planet *sin( α + Γ + V relp ) * cos( α + Γ ( V sin( α + Γ )) + ( V + V *cos( α + Γ )) ) = sqrt * relp planet relp relp )

68 Total system effect of a gravity assist The deflection of the planet (or major body) relative velocity vector means that the energy and angular momentum of the departure orbit relative to the central body are modified This energy and momentum exchange is matched by a change in these components for the major body The large mass of the major body means that its orbit is negligibly effected Many large spacecraft performing energy raising gravity assists at Earth could slow down the planet! Page 68

69 Gravity assist example: co-planar assumptions Consider a gravity assist at Venus with a spacecraft in an elitpical orbit Initial Ap Initial Pe Planet Planet Anom RV radius Hyper V Peri h (km) The possible orbit changes post fly-by are Post GA Ap Post GA Pe Post GA Ap Post GA Pe Page 69

70 General coplanar gravity assist cases Effect of fly-by altitude The previous example of a Venus crossing orbit (perihelion at 0.5AU) can be examined to determine the effect of fly-by altitude (ie pericentre radius) The aphelion raising solution is taken As the altitude increases the effectiveness reduces Ap-Pe (AU Ap Pe sma sma altitude (km) Page 70

71 .4 Multiple gravity assists Page 71

72 Multiple gravity assists Repeated gravity assists can be used to perform large modifications to the spacecraft orbit. Repeated gravity assists at the same planetary body or moon can be performed. This is generally a resonant gravity assist sequence. The spacecraft must reach, after the gravity assist, an orbit that is resonant with the body in question. This means that after an integer number of revolutions, the spacecraft will re-encounter the body. Resonance ratios can be: n:1 (n revolutions by the spacecraft, 1 revolution by the body) n:m (m revolutions by the body) 1:m (1 revolution by the spacecraft, m revolutions by the body). In these resonant cases, the excess hyperbolic speed is the same at successive fly-bys Page 7

73 Multiple gravity assist limits Such repeated gravity assists can be used to progressively raise apohelion Conversely, may be used to lower perihelion Eventually a limiting case is reached where the velocity with respect to the central body cannot be further increased/decreased Ie the relative velocity becomes aligned with fly-by body velocity Page 73 V V planet V + V max min V ( V ) = V flyby ( V ) = V flyby α V rel =V

74 Multiple gravity assist example Consider a series of resonant gravity assists in the Jovian system. The fly-by body is Ganymede. The initial orbit is chosen to be: Apocentre = 0 million km Pericentre = km Ganymede orbital radius = 1.07 million km At the end of the sequence the orbit is approaching the limiting case with apocentre just above Ganymede Apocentre (km) Pericentre (km) Page 74

75 Tisserand s Criterion Tisserand derived a criterion from which it is possible to compare orbits that appear significantly perturbed and deduce whether the relationship between them may be due to a gravity assist at a planet. An example is comets passing close to Jupiter The basis of the method is the circular, restricted three-body problem. The Jacobi constant is used V + U = V is the speed w.r.t the rotating frame When motion is considered at a sufficiently large distance from the fly-by planet, its gravity field may be neglected. Speed in the rotating frame may be converted to speed in an inertial frame, V IC. r c is the distance from the central body. µ VIC ω rc VI zˆ = r Page 75 C ( ) C C

76 Tisserand s Criterion () The constant may be expressed at distances far from the planet as: µ ωh cos i a In the co-planar case, then the orbit relationships before and after are given by: µ µ ω µ a1(1 e1 ) = ω µ a (1 e a a 1 For a given initial spacecraft orbit the evolution of the possible orbits under repeated gravity assists at the same planet may be evaluated. The locus of the evolution of the orbital elements is given by the above relationship, eg a relationship between a and e. = C ) Page 76

77 Tisserand s Criterion (3) The change in the orbit at a single gravity assist is determined by: Excess hyperbolic speed w.r.t. the planet (or moon) Gravity constant Fly-by pericentre radius This change therefore determines the size of the step along the locus. The locus can alternatively be expressed in terms of orbit period and pericentre, or apocentre and pericentre Period Step at gravity assist Locus for a given Vinfinity at target body Pericentre Page 77

78 Tisserand s Criterion Example Consider the example of an Earth crossing orbit and the locus of heliocentric orbit evolution for multiple gravity assists at Earth Both periodpericentre and apocentrepericentre plots are considered Apocentre (km) Period (days) 8.00E E E E E E+0.00E E E E E E E E E E E+08 Pericentre (km) 3.50E E+08.50E+08.00E E E E+07 Page E E E E E E E E E+08 Pericentre (km)

79 Designing sequences using Tisserand s Criterion Loci of possible orbit sequences at different fly-by bodies can be considered The intersection of loci indicates the possibility to switch from one sequence (at fly-by body 1) to another sequence, at fly-by body This assumes that the required phasing between bodies 1 and can be achieved to execute the transfer. Can always be achieved if wait for long enough (or perform a small phasing manoeuvre) This technique is very relevant to GA sequence design at a planetary moon system, where many possibilities exist for multiple fly-bys (orbit periods are shorter than heliocentric case) Page 79

80 Designing a gravity assist sequence in the Jovian system The overlapping loci for gravity assists at Callisto, Ganymede and Europa indicate the multitude of options for gravity assist combinations Apocentre (km) C-7000 C-6000 C-5000 G-7000 G-6000 G-5000 E-4000 E-3000 E-000 E-6000 G E E+05.00E E E E E E E E E+06 Page 80 Pericentre (km)

81 Designing a gravity assist sequence in the Jovian system () A sequence of gravity assists at Ganymede and Europa can be used to reach a low approach speed orbit relative to Europa from an initial capture orbit at Jupiter Manoeuvres may also be allowed in Jovian orbits Objective is to minimise the total transfer DeltaV: Manoeuvres in the Jovian system Orbit insertion manoeuvre at Europa V total = i V api + µ r Europa EuropaOrbit + V final r µ EuropaOrbit V api is the V applied at apocentre i to raise pericentre. V final is the excess hyperbolic speed at Europa from which insertion to Europa orbit is made, r Europaorbit is the radius of this circular orbit about Europa. Page 81

82 Example of gravity assist sequence in the Jovian system This tour uses an initial sequence of gravity assists at Ganymede (4) Then 1 at Europa and the next back to Ganymede to raise pericentre 3 more GA s at Europa then another at Ganymede to raise pericentre Then a sequence of Europa GA s with intermediate DV s to raise pericentre and so progressively reduce approach speed at Europa DeltaV to Europa circular orbit approx 1000 m/s Period (days) E E E E E E E E+06 Pericentre (km) G-7000 G-6000 G-5000 E-4000 E-3000 E-000 E-5000 G-4000 Route Page 8

83 Example of gravity assist sequence in the Jovian system () The sequence is Page 83 Event Apocentre (km) Pericentre (km) Period (days) V (m/s) Resonance DV (m/s) Capture :1 Gany :1 Gany :1 Gany : Gany Europa Gany Europa :4 Europa :4 Europa Gany Europa :3 DV Europa : DV Europa :4 DV Europa :10 DV

84 Example of use of gravity assists: Mission to Jupiter A direct transfer to Jupiter needs an Earth departing Vinfinity of approx 9 km/sec A multi-fly-by gravity assist route can substantially reduce this. The Galileo strategy : 1. Earth to Venus. Venus gravity assist to Earth (Vinfinity approx 6 km/sec) 3. Earth gravity assist to Earth via year resonant orbit (Vinfinity approx 9 km/sec) 4. Earth gravity assist to Jupiter Total transfer duration 6. years Page 84

85 .5 Plane changing gravity assists Page 85

86 3D Gravity Assist geometry and the B plane The general case considers non co-planar orbits of spacecraft and planet. The B plane is the plane perpendicular to the asymptotic relative velocity vector, or hyperbolic approach vector. For a given, targeted fly-by pericentre radius, the intersection of the forward projection of the approach hyperbola with the B plane can take place at any point on a circle around the planet. A very small deep space manoeuvre can modify the approach to lie anywhere on such a circle A plane may be defined, that contains the approaching asymptotic velocity vector and the velocity vector of the planet the Approach plane. This plane is therefore perpendicular to the B plane. The Beta angle, β, is the angle between the X axis (defined by plane intersection) and the location of the intersection of the approach relative velocity vector with the B plane. Cases 1 and in the figure are therefore at Beta angles of zero and 90 degrees.. Approaching relative velocity vector is offset to achieve an approach that is parallel to the approach plane B plane β = 90 deg 1. Approaching relative velocity vector contained in the approach plane Approach plane Z B Y B X B Page 86

87 3D deflection It is possible to define a further plane, the fly-by plane, as being the plane containing the approaching velocity vector relative to the major body (the asymptotic approach vector) and also the departing relative velocity vector. This plane is defined by the Beta angle, described previously. The deflection is in the fly-by plane Two deflection angle components can be defined: One in the approach plane and one perpendicular to it α φ α θ α θ is negative in this diagram β α B plane Approaching relative velocity vector Plane containing body velocity vector and spacecraft velocity vector, ie The approach plane Page 87

88 3D gravity assist example Consider the example of a gravity assist at the Jovian moon, Ganymede. Ganymede lies in a near circular orbit about Jupiter, with semi-major axis at 1.07 million km. Consider a range of initial orbits each crossing Ganymedes s orbit. The initial apocentre is 5 million km, and a range of pericentres from km to km are considered. The result is to generate a range of excess hyperbolic speeds with respect to Ganymede. Page 88 Pericentre (km) Speed at Ganymede (m/s) Vinfinity (m/s) Flight path angle (deg)

89 3D gravity assist example () The Beta angle of the fly-by is varied through a 360 degree range Post gravity assist orbits are then found. In each case it is assumed that the initial spacecraft orbit is coplanar with the orbit of Ganymede. Defines the approach plane The pericentre altitude at Ganymede is assumed to be 300km Page 89

90 3D gravity assist example (3) A β angle of 90 degrees maximises the inclination change It also yields a small change in the post gravity assist apocentre and pericentre. Inclination change is measured with respect to the initial approach plane Greatest inclination change occurs for the highest initial pericentre (or lowest Vinfinity) case Apocentre (m) Inclination (deg) Beta(deg) Beta(deg) Page 90

91 Multiple plane changing fly-bys It is possible to find a Beta angle that maintains the relative speed at the fly-by target, after the fly-by (ie the semi-major axis is unchanged). This means that the spacecraft may stay in a resonant orbit (ie same orbital period about the target). This angle is generally close to 90 degrees. Repeated gravity assists at the same planetary body or moon allows a large inclination change to be achieved. Each fly-by is designed such that the effect of the gravity assist is to increase inclination and also to achieve a velocity relative to the central body that yields a resonant orbit with respect to the major body. The spacecraft will then return to the major body after some integer number of revolutions about the central body. Subsequent fly-bys can further increase inclination. The post fly-by velocity relative to the central body an hence resonance can be maintained by choosing the appropriate Beta angle Page 91

92 Multiple plane changing fly-bys () Repeated gravity assists result in the plane containing planet and spacecraft velocities reaching 90 degrees relative to the initial approach plane (generally close to the fly-by body orbit plane). If the same resonant fly-bys are maintained then this plane will continue to rotate past 90 degrees. Maximum inclination is achieved close to the 90 degree rotation case. Inclination may be further increased by changing the resonant orbit Now the approach plane is unchanged and the approach vector is deflected in this 90 degree rotated plane A maximum inclination can be found Page 9

93 Maximum achievable inclination Repeated gravity assists result in the plane containing planet and spacecraft velocities reaching 90 degrees relative to the initial approach plane (generally close to the fly-by body orbit plane). Maximum inclination is achieved at the 90 degree rotation case. Plane containing spacecraft and fly-by body velocity vectors V limit V rel =V V 1 V rel =V The maximum inclination (relative to the initial approach plane) in the case of the same resonant orbit being maintained is given by: V i = tan 1 max V planet The maximum achievable inclination (by changing resonance) is given by: Page 93 i max = sinθ + V cosθ sin 1 V V planet i max V planet θ Fly-by body orbit plane ψ=90

94 Example of multiple plane changing gravity assists A heliocentric example of plane changing with multiple gravity assists cab be considered. Useful for missions that are required to perform out of ecliptic observations. A single gravity assist at Jupiter can achieve a 90 degree inclination with respect to the ecliptic. However considerable time and/or fuel is needed to reach Jupiter. An alternative strategy to reach high inclinations could be considered at planets within the inner solar system. Venus is a good choice Relatively high mass and easily reached from Earth. Page 94

95 Example of multiple plane changing gravity assists () The maximum, total inclination change that is achievable is dependent on the excess hyperbolic speed with respect to Venus Compare with previous equations Different excess hyperbolic speeds are achieved by orbits with different aphelion/perihelions High inclination requires high Vinfinity Repeated gravity assists will be required to achieve high inclinations. Considered as an option for the Solar Orbiter mission The objective of this mission is to observe the Sun from high latitudes when relatively close to the Sun. Therefore such a sequence is well suited to this type of Solar observing mission. Page 95

96 Example of multiple plane changing gravity assists (3) The initial orbit is: Apocentre = 13 million km, Pericentre = 33 million km. Excess hyperbolic speed with respect to Venus is 0. km/sec A 3: resonance may be reached with Venus. Fly-bys occur approximately every 450 days, orbital period is 150 days. 6 gravity assists are used in reaching a maximum inclination of nearly 35 degrees. This procedure will therefore take approximately 700 days, or approaching 7.5 years. Inclination (deg) Apocentre (km) Pericentre (km) Page

97 Example of multiple plane changing gravity assists (5) Inclination increases as the perihelion is progressively raised Orbits of Earth and Venus also shown Page 97

98 .6 Gravitational escape and capture Page 98

99 Natural Capture Temporary capture or transit from one heliocentric orbit to another by gravitational perturbation possible for comets and asteroids Jovian neighbourhood shows numerous irregular bodies possible captured via approaching through the Lagrange points. Example is Comet Shoemaker-Levy 9, which eventually impacted Jupiter in Jupiter studied regarding possible transitions between heliocentric resonant orbits, via Lagrange point passage. Example: Oterma and Gehrels 3 in :3 resonant orbit to 3: resonant orbit References: Carusi, Valsecchi, Koon, Lo, Marsden, Ross Page 99

100 Application to spacecraft transfers Technique applicable to transfers that include approach or escape speeds that are appropriate for the planet under consideration. Approach with significantly greater speed will not experience any significant advantage from the effect. For the inner planets, these speeds are not high enough to allow a direct transfer to another planet. Therefore methods may be sought that augment this effect to achieve the required transfers. Two techniques that may be used are: Use of additional deep space propulsive manoeuvres. These can be implemented by high specific impulse, low thrust systems. Use of gravity assist manoeuvres at the planet under consideration. This technique also requires intermediate propulsive manoeuvres. Page 100

101 Escape and capture mechanism Motion under multiple gravity fields can be exploited to achieve a beneficial orbit change In the three body problem initially bound orbits can achieve escape Conversely hyperbolic approach orbits can achieve capture The mechanism can be demonstrated via consideration of the 3 body problem and the Jacobi constant Figure shows energy surface. Assumes planet relative velocity in inertial oriented frame is always radial Earth is at origin, XY ecliptic rotating frame axes Page 101

102 Escape and capture mechanism () The energy is negative (bound orbit) close to Earth However achieves a positive value at distances >M km (in this example) The key is to find a trajectory that can achieve this Such trajectories require passages in the vicinity of the colinear Lagrange points Initial orbits (in escape case) are similar to free injection transfers to the Lagrange points Typical Earth escape vinfinity achievable 1000 m/s Page 10

103 Example of BepiColombo: Capture at Mercury A gravitational capture method is used at Mercury This requires a further reduced approach Vinfinity compared with the previous case (typically < 300 m/s). Then after flying close to one of the Mercury-Sun collinear Lagrange libration points, the spacecraft s osculating ephemeris on reaching Mercury pericentre is that of a high apocentre, captured orbit. This is a weakly bound orbit which may then remain captured in a high apocentre orbit for several revolutions. The chemical propulsion DeltaV to reduce apocentre from the weakly captured value ( to 00000km) to the target orbit at Mercury is now significantly reduced compared to hyperbolic approach. A net mass fuel mass reduction can be achieved by this strategy. Page 103

104 Example of BepiColombo: Capture at Mercury () Freely reaching weakly bound orbit significantly improves the robustness of the mission. The spacecraft can miss the nominal Mercury pericentre manoeuvre and remain in the high apocentre orbit and then re-attempt the apocentre reduction manoeuvre at subsequent Mercury pericentre passages. The use of a gravitational capture method means that the spacecraft must spend approximately half a Mercury orbit period longer before reaching pericentre Trajectory is seen in Mercury-Sun rotating frame Page 104

Astrodynamics (AERO0024)

Astrodynamics (AERO0024) Astrodynamics (AERO0024) L06: Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 Problem Statement? Hint #1: design the Earth-Mars transfer using known concepts

More information

Astrodynamics (AERO0024)

Astrodynamics (AERO0024) Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem

More information

Interplanetary Trajectory design for Rosetta and Solar Orbiter

Interplanetary Trajectory design for Rosetta and Solar Orbiter Jose Rodriguez-Canabal Memorial Interplanetary Trajectory design for Rosetta and Solar Orbiter Johannes Schoenmaekers H / Mission Analysis Section (HSO-GFA) Senior Adviser ESA / ESOC Email: johannes.schoenmaekers@esa.int

More information

L eaving Earth and arriving at another planet or asteroid requires

L eaving Earth and arriving at another planet or asteroid requires Designing Interplanetary Transfers L eaving Earth and arriving at another planet or asteroid requires a spacecraft to implement a sequence of manoeuvres. These include changes of velocity needed to escape

More information

Flight and Orbital Mechanics

Flight and Orbital Mechanics Flight and Orbital Mechanics Lecture slides Challenge the future 1 Flight and Orbital Mechanics AE-104, lecture hours 1-4: Interplanetary flight Ron Noomen October 5, 01 AE104 Flight and Orbital Mechanics

More information

Lecture D30 - Orbit Transfers

Lecture D30 - Orbit Transfers J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we

More information

AST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000).

AST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000). AST111 PROBLEM SET 2 SOLUTIONS Home work problems 1. Angles on the sky and asteroid motion An asteroid is observed at two different times. The asteroid is located at RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch

More information

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major

More information

Interplanetary Mission Opportunities

Interplanetary Mission Opportunities Interplanetary Mission Opportunities Introduction The quest for unravelling the mysteries of the universe is as old as human history. With the advent of new space technologies, exploration of space became

More information

Satellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN

Satellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,

More information

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS

ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.

More information

Interplanetary Travel

Interplanetary Travel Interplanetary Travel Interplanetary Travel Concept Patched Conic Hypothesis Departure & Arrival Manoeuvres Interplanetary Travel Concept Interplanetary travel is concerned with motion of manmade objects

More information

ASTRIUM. Interplanetary Path Early Design Tools at ASTRIUM Space Transportation. Nathalie DELATTRE ASTRIUM Space Transportation.

ASTRIUM. Interplanetary Path Early Design Tools at ASTRIUM Space Transportation. Nathalie DELATTRE ASTRIUM Space Transportation. Interplanetary Path Early Design Tools at Space Transportation Nathalie DELATTRE Space Transportation Page 1 Interplanetary missions Prime approach: -ST has developed tools for all phases Launch from Earth

More information

Solar Orbiter Ballistic Transfer Mission Analysis Synthesis

Solar Orbiter Ballistic Transfer Mission Analysis Synthesis European Space Agency Agence Spatiale Européenne directorate of operations and infrastructure ground systems engineering department mission analysis office MAO Working Paper No. 483 Issue 1, Rev. 0 Solar

More information

ASEN 5050 SPACEFLIGHT DYNAMICS Interplanetary

ASEN 5050 SPACEFLIGHT DYNAMICS Interplanetary ASEN 5050 SPACEFLIGHT DYNAMICS Interplanetary Prof. Jeffrey S. Parker University of Colorado Boulder Lecture 29: Interplanetary 1 HW 8 is out Due Wednesday, Nov 12. J2 effect Using VOPs Announcements Reading:

More information

Feasible Mission Designs for Solar Probe Plus to Launch in 2015, 2016, 2017, or November 19, 2008

Feasible Mission Designs for Solar Probe Plus to Launch in 2015, 2016, 2017, or November 19, 2008 Feasible Mission Designs for Solar Probe Plus to Launch in 2015, 2016, 2017, or 2018 2007 Solar Probe Study & Mission Requirements Trajectory study and mission design trades were conducted in the fall

More information

Patch Conics. Basic Approach

Patch Conics. Basic Approach Patch Conics Basic Approach Inside the sphere of influence: Planet is the perturbing body Outside the sphere of influence: Sun is the perturbing body (no extra-solar system trajectories in this class...)

More information

Chapter 14 Satellite Motion

Chapter 14 Satellite Motion 1 Academic Physics Mechanics Chapter 14 Satellite Motion The Mechanical Universe Kepler's Three Laws (Episode 21) The Kepler Problem (Episode 22) Energy and Eccentricity (Episode 23) Navigating in Space

More information

Fundamentals of Satellite technology

Fundamentals of Satellite technology Fundamentals of Satellite technology Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Orbital Plane All of the planets,

More information

Astromechanics. 6. Changing Orbits

Astromechanics. 6. Changing Orbits Astromechanics 6. Changing Orbits Once an orbit is established in the two body problem, it will remain the same size (semi major axis) and shape (eccentricity) in the original orbit plane. In order to

More information

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial

More information

OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES

OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES Gislaine de Felipe and Antonio Fernando Bertachini de Almeida Prado Instituto Nacional de Pesquisas Espaciais - São José dos Campos - SP - 12227-010

More information

Previous Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions

Previous Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions 2 / 48 Previous Lecture Orbital maneuvers: general framework Single-impulse maneuver: compatibility conditions closed form expression for the impulsive velocity vector magnitude interpretation coplanar

More information

Celestial Mechanics Lecture 10

Celestial Mechanics Lecture 10 Celestial Mechanics Lecture 10 ˆ This is the first of two topics which I have added to the curriculum for this term. ˆ We have a surprizing amount of firepower at our disposal to analyze some basic problems

More information

Massimiliano Vasile, Stefano Campagnola, Paolo Depascale, Stefano Pessina, Francesco Topputo

Massimiliano Vasile, Stefano Campagnola, Paolo Depascale, Stefano Pessina, Francesco Topputo A Toolbox for Preliminary Massimiliano Vasile, Stefano Campagnola, Paolo Depascale, Stefano Pessina, Francesco Topputo Mission Analysis and Design PAMSIT IMAGO ATOM-C EPIC Massimiliano Vasile, Stefano

More information

The Heliocentric Model of Copernicus

The Heliocentric Model of Copernicus Celestial Mechanics The Heliocentric Model of Copernicus Sun at the center and planets (including Earth) orbiting along circles. inferior planets - planets closer to Sun than Earth - Mercury, Venus superior

More information

INTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS

INTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS INTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS Francesco Topputo (), Massimiliano Vasile () and Franco Bernelli-Zazzera () () Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano,

More information

Design of a Multi-Moon Orbiter

Design of a Multi-Moon Orbiter C C Dynamical A L T E C S H Design of a Multi-Moon Orbiter Shane D. Ross Control and Dynamical Systems and JPL, Caltech W.S. Koon, M.W. Lo, J.E. Marsden AAS/AIAA Space Flight Mechanics Meeting Ponce, Puerto

More information

Extending the Patched-Conic Approximation to the Restricted Four-Body Problem

Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean

More information

APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING

APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types

More information

Chapter 8. Precise Lunar Gravity Assist Trajectories. to Geo-stationary Orbits

Chapter 8. Precise Lunar Gravity Assist Trajectories. to Geo-stationary Orbits Chapter 8 Precise Lunar Gravity Assist Trajectories to Geo-stationary Orbits Abstract A numerical search technique for designing a trajectory that transfers a spacecraft from a high inclination Earth orbit

More information

State Vector Reference Body Transformations Throughout The Solar System

State Vector Reference Body Transformations Throughout The Solar System Introduction One of the more mind-boggling programming challenges encountered during interplanetary trajectory design software development is how to transform a specified inertial position and velocity

More information

Study of the Transfer Between Libration Point Orbits and Lunar Orbits in Earth-Moon System

Study of the Transfer Between Libration Point Orbits and Lunar Orbits in Earth-Moon System Noname manuscript No. (will be inserted by the editor) Study of the Transfer Between Libration Point Orbits and Lunar Orbits in Earth-Moon System Yu Cheng, Gerard Gómez Josep J. Masdemont 3 Jianping Yuan

More information

A SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE

A SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE DOING PHYSICS WITH MATLAB A SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE Download Directory: Matlab mscripts mec_satellite_gui.m The [2D] motion of a satellite around the Earth is computed from

More information

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements

AST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends

More information

A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter

A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3: RESONANCE Érica C. Nogueira, Othon C. Winter Grupo de Dinâmica Orbital e Planetologia UNESP -- Guaratinguetá -- Brazil Antonio F.B. de

More information

Tides and Lagrange Points

Tides and Lagrange Points Ast111, Lecture 3a Tides and Lagrange Points Arial view of Tidal surge in Alaska (S. Sharnoff) Tides Tidal disruption Lagrange Points Tadpole Orbits and Trojans Tidal Bulges Tides Tidal Force The force

More information

Trajectories to the Moons (incl. a trajectory for an Enceladus orbiter)

Trajectories to the Moons (incl. a trajectory for an Enceladus orbiter) Trajectories to the Moons (incl. a trajectory for an Enceladus orbiter) Stefano Campagnola 1,Ryan Russell 2, and Nathan Strange 3 JPL-Caltech/CDS Meetings, August 3rd 2009 1 Aerospace and Mechanical Engineering,

More information

ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM

ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM JORGE K. S. FORMIGA 1,2 and ANTONIO F B A PRADO 2 National Institute for Space Research -INPE 1 Technology Faculty-FATEC-SJC

More information

Escape Trajectories from the L 2 Point of the Earth-Moon System

Escape Trajectories from the L 2 Point of the Earth-Moon System Trans. Japan Soc. Aero. Space Sci. Vol. 57, No. 4, pp. 238 244, 24 Escape Trajectories from the L 2 Point of the Earth-Moon System By Keita TANAKA Þ and Jun ichiro KAWAGUCHI 2Þ Þ Department of Aeronautics

More information

Interplanetary Trajectory Design using Dynamical Systems Theory

Interplanetary Trajectory Design using Dynamical Systems Theory Interplanetary Trajectory Design using Dynamical Systems Theory THESIS REPORT by Linda van der Ham 8 February 2012 The image on the front is an artist impression of the Interplanetary Superhighway [NASA,

More information

Patched Conic Interplanetary Trajectory Design Tool

Patched Conic Interplanetary Trajectory Design Tool Copyright by Martin James Brennan 2011 The Thesis committee for Martin James Brennan Certifies that this is the approved version of the following thesis: Patched Conic Interplanetary Trajectory Design

More information

Ulrich Walter. Astronautics. The Physics of Space Flight. 2nd, Enlarged and Improved Edition

Ulrich Walter. Astronautics. The Physics of Space Flight. 2nd, Enlarged and Improved Edition Ulrich Walter Astronautics The Physics of Space Flight 2nd, Enlarged and Improved Edition Preface to Second Edition Preface XVII Acknowledgments XIX List of Symbols XXI XV 1 Rocket Fundamentals 1 1.1 Rocket

More information

Libration Orbit Mission Design: Applications Of Numerical And Dynamical Methods

Libration Orbit Mission Design: Applications Of Numerical And Dynamical Methods Libration Orbit Mission Design: Applications Of Numerical And Dynamical Methods David Folta and Mark Beckman NASA - Goddard Space Flight Center Libration Point Orbits and Applications June 1-14, 22, Girona,

More information

The B-Plane Interplanetary Mission Design

The B-Plane Interplanetary Mission Design The B-Plane Interplanetary Mission Design Collin Bezrouk 2/11/2015 2/11/2015 1 Contents 1. Motivation for B-Plane Targeting 2. Deriving the B-Plane 3. Deriving Targetable B-Plane Elements 4. How to Target

More information

INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT ABOUT LIBRATION POINTS

INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT ABOUT LIBRATION POINTS 6 th International Conference on Astrodynamics Tools and Technique (ICATT) INDIRECT PLANETARY CAPTURE VIA PERIODIC ORBIT LI Xiangyu, Qiao Dong, Cui Pingyuan Beijing Institute of Technology Institute of

More information

MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design

MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time and flight path angle as a function of orbital position Relative orbital

More information

Chapter 8. Orbits. 8.1 Conics

Chapter 8. Orbits. 8.1 Conics Chapter 8 Orbits 8.1 Conics Conic sections first studied in the abstract by the Greeks are the curves formed by the intersection of a plane with a cone. Ignoring degenerate cases (such as a point, or pairs

More information

ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)

ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler

More information

Optimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission

Optimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission Optimal Gravity Assisted Orbit Insertion for Europa Orbiter Mission Deepak Gaur 1, M. S. Prasad 2 1 M. Tech. (Avionics), Amity Institute of Space Science and Technology, Amity University, Noida, U.P.,

More information

A Study of the Close Approach Between a Planet and a Cloud of Particles

A Study of the Close Approach Between a Planet and a Cloud of Particles A Study of the Close Approach Between a Planet a Cloud of Particles IIAN MARTINS GOMES, ANTONIO F. B. A. PRADO National Institute for Space Research INPE - DMC Av. Dos Astronautas 1758 São José dos Campos

More information

MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design

MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. Principles of Space Systems Design Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time and flight path angle as a function of orbital position Relative orbital

More information

Study of LEO to GEO transfers via the L1 Sun-Earth or Earth-Moon libration points

Study of LEO to GEO transfers via the L1 Sun-Earth or Earth-Moon libration points Study of LEO to GEO transfers via the L1 Sun-Earth or Earth-Moon libration points Author: Elisabet Herrera Sucarrat Advisor: Josep Masdemont Soler Projecte Final de Màster d'enginyeria Matemàtica Facultat

More information

Analysis for the Earth Escape Strategy Using Unstable Manifolds of Sun-Earth Halo Orbit and Lunar Gravity Assists

Analysis for the Earth Escape Strategy Using Unstable Manifolds of Sun-Earth Halo Orbit and Lunar Gravity Assists Analysis for the Earth Escape Strategy Using Unstable Manifolds of Sun-Earth Halo Orbit and Lunar Gravity Assists Hongru Chen ), Yasuhiro Kawakatsu ) and Toshiya Hanada ) ) Department of Aeronautics and

More information

Escape Trajectories from Sun Earth Distant Retrograde Orbits

Escape Trajectories from Sun Earth Distant Retrograde Orbits Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists30, pp. Pd_67-Pd_75, 06 Escape Trajectories from Sun Earth Distant Retrograde Orbits By Yusue OKI ) and Junichiro KAWAGUCHI ) ) Department of Aeronautics

More information

Today in Astronomy 111: multiple-body systems

Today in Astronomy 111: multiple-body systems Today in Astronomy 111: multiple-body systems Two simple topics in systems of three bodies: Transfer orbits, and how asteroids get knocked from orbit to orbit I The two-body potential without Coriolis

More information

Orbital Mechanics MARYLAND

Orbital Mechanics MARYLAND Orbital Mechanics Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit Interplanetary trajectories Planetary launch and

More information

BepiColombo Launch Window Design Based on a Phasing Loop Strategy

BepiColombo Launch Window Design Based on a Phasing Loop Strategy BepiColombo Launch Window Design Based on a Phasing Loop Strategy Carlos Corral van Damme * GMV, Tres Cantos, Madrid, 28760, Spain Rüdiger Jehn and Stefano Campagnola European Space Operations Centre,

More information

A space probe to Jupiter

A space probe to Jupiter Problem 3 Page 1 Problem 3 A space probe to Jupiter We consider in this problem a method frequently used to accelerate space probes in the desired direction. The space probe flies by a planet, and can

More information

arxiv: v1 [physics.pop-ph] 29 May 2009

arxiv: v1 [physics.pop-ph] 29 May 2009 Gravity assist in 3D like in Ulysses mission Herbert Morales Escuela de Física, Universidad de Costa Rica, San José, Costa Rica arxiv:0905.4788v1 [physics.pop-ph] 29 May 2009 (Dated: October 22, 2018)

More information

Gravitation. Luis Anchordoqui

Gravitation. Luis Anchordoqui Gravitation Kepler's law and Newton's Synthesis The nighttime sky with its myriad stars and shinning planets has always fascinated people on Earth. Towards the end of the XVI century the astronomer Tycho

More information

Dynamical Systems and Space Mission Design

Dynamical Systems and Space Mission Design Dynamical Systems and Space Mission Design Jerrold Marsden, Wang Sang Koon and Martin Lo Shane Ross Control and Dynamical Systems, Caltech shane@cds.caltech.edu Constructing Orbits of Prescribed Itineraries:

More information

arxiv:gr-qc/ v1 15 Nov 2004

arxiv:gr-qc/ v1 15 Nov 2004 Mission design for LISA Pathfinder arxiv:gr-qc/0411071v1 15 Nov 2004 M Landgraf, M Hechler, and S Kemble ESA/ESOC, Robert-Bosch-Straße 5, D-64293 Darmstadt, Germany E-mail: Markus.Landgraf@esa.int EADS

More information

Problem Description for the 6 th Global Trajectory Optimisation Competition. Anastassios E. Petropoulos 1

Problem Description for the 6 th Global Trajectory Optimisation Competition. Anastassios E. Petropoulos 1 Problem Description for the 6 th Global Trajectory Optimisation Competition Anastassios E. Petropoulos 1 Outer Planet Mission Analysis Group Jet Propulsion Laboratory California Institute of Technology

More information

Chapter 13. Gravitation

Chapter 13. Gravitation Chapter 13 Gravitation e = c/a A note about eccentricity For a circle c = 0 à e = 0 a Orbit Examples Mercury has the highest eccentricity of any planet (a) e Mercury = 0.21 Halley s comet has an orbit

More information

Lecture 13. Gravity in the Solar System

Lecture 13. Gravity in the Solar System Lecture 13 Gravity in the Solar System Guiding Questions 1. How was the heliocentric model established? What are monumental steps in the history of the heliocentric model? 2. How do Kepler s three laws

More information

Chapter 4. Integrated Algorithm for Impact Lunar Transfer Trajectories. using Pseudo state Technique

Chapter 4. Integrated Algorithm for Impact Lunar Transfer Trajectories. using Pseudo state Technique Chapter 4 Integrated Algorithm for Impact Lunar Transfer Trajectories using Pseudo state Technique Abstract A new integrated algorithm to generate the design of one-way direct transfer trajectory for moon

More information

arxiv: v1 [gr-qc] 7 Oct 2009

arxiv: v1 [gr-qc] 7 Oct 2009 Earth Flyby Anomalies Michael Martin Nieto 1 and John D. Anderson 2 Michael Martin Nieto is a fellow at Los Alamos National Laboratory and John D. Anderson is a senior research scientist emeritus at the

More information

ASEN 6008: Interplanetary Mission Design Lab Spring, 2015

ASEN 6008: Interplanetary Mission Design Lab Spring, 2015 ASEN 6008: Interplanetary Mission Design Lab Spring, 2015 Lab 4: Targeting Mars using the B-Plane Name: I d like to give credit to Scott Mitchell who developed this lab exercise. He is the lead Astrodynamicist

More information

Orbital Dynamics and Impact Probability Analysis

Orbital Dynamics and Impact Probability Analysis Orbital Dynamics and Impact Probability Analysis (ISAS/JAXA) 1 Overview This presentation mainly focuses on a following point regarding planetary protection. - How to prove that a mission satisfies the

More information

Orbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design

Orbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design Lecture #05 September 15, 2015 Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit

More information

Chapter 4. Motion and gravity

Chapter 4. Motion and gravity Chapter 4. Motion and gravity Announcements Labs open this week to finish. You may go to any lab section this week (most people done). Lab exercise 2 starts Oct 2. It's the long one!! Midterm exam likely

More information

Design of Low Energy Space Missions using Dynamical Systems Theory

Design of Low Energy Space Missions using Dynamical Systems Theory Design of Low Energy Space Missions using Dynamical Systems Theory Koon, Lo, Marsden, and Ross W.S. Koon (Caltech) and S.D. Ross (USC) CIMMS Workshop, October 7, 24 Acknowledgements H. Poincaré, J. Moser

More information

Identifying Safe Zones for Planetary Satellite Orbiters

Identifying Safe Zones for Planetary Satellite Orbiters AIAA/AAS Astrodynamics Specialist Conference and Exhibit 16-19 August 2004, Providence, Rhode Island AIAA 2004-4862 Identifying Safe Zones for Planetary Satellite Orbiters M.E. Paskowitz and D.J. Scheeres

More information

The Law of Ellipses (Kepler s First Law): all planets orbit the sun in a

The Law of Ellipses (Kepler s First Law): all planets orbit the sun in a Team Number Team Members Present Learning Objectives 1. Practice the Engineering Process a series of steps to follow to design a solution to a problem. 2. Practice the Five Dimensions of Being a Good Team

More information

Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the

Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the Geocentric-Heliocentric debate (Brahe s data and Kepler s analysis)

More information

Ellipses. Gravitation and Mechanics MOPS Ellipse Foci. Any Point on the Curve Distance. Distance B Foci

Ellipses. Gravitation and Mechanics MOPS Ellipse Foci. Any Point on the Curve Distance. Distance B Foci MOPS0513-02-00 Chapter 3. Objectives: Gravitation and Mechanics Upon completion of this chapter you will be able to describe the force of gravity, characteristics of ellipses, and the concepts of Newton's

More information

Celestial Mechanics and Satellite Orbits

Celestial Mechanics and Satellite Orbits Celestial Mechanics and Satellite Orbits Introduction to Space 2017 Slides: Jaan Praks, Hannu Koskinen, Zainab Saleem Lecture: Jaan Praks Assignment Draw Earth, and a satellite orbiting the Earth. Draw

More information

ESA s Juice: Mission Summary and Fact Sheet

ESA s Juice: Mission Summary and Fact Sheet ESA s Juice: Mission Summary and Fact Sheet JUICE - JUpiter ICy moons Explorer - is the first large-class mission in ESA's Cosmic Vision 2015-2025 programme. Planned for launch in 2022 and arrival at Jupiter

More information

COMPARISON OF LOW-ENERGY LUNAR TRANSFER TRAJECTORIES TO INVARIANT MANIFOLDS

COMPARISON OF LOW-ENERGY LUNAR TRANSFER TRAJECTORIES TO INVARIANT MANIFOLDS AAS 11-423 COMPARISON OF LOW-ENERGY LUNAR TRANSFER TRAJECTORIES TO INVARIANT MANIFOLDS Rodney L. Anderson and Jeffrey S. Parker INTRODUCTION In this study, transfer trajectories from the Earth to the Moon

More information

A Low-Cost Mission for LISA Markus Landgraf, Florian Renk, Pierre Joachim, Rüdiger Jehn HSO-GFA

A Low-Cost Mission for LISA Markus Landgraf, Florian Renk, Pierre Joachim, Rüdiger Jehn HSO-GFA A Low-Cost Mission for LISA Markus Landgraf, Florian Renk, Pierre Joachim, Rüdiger Jehn HSO-GFA LISA Internal Final Presentation July 8 th, 2011 CDF HSO-GFA Page 1 Overview Basic working assumptions Operational

More information

Chapter 13. Universal Gravitation

Chapter 13. Universal Gravitation Chapter 13 Universal Gravitation Planetary Motion A large amount of data had been collected by 1687. There was no clear understanding of the forces related to these motions. Isaac Newton provided the answer.

More information

ABOUT COMBINING TISSERAND GRAPH GRAVITY-ASSIST SEQUENCING WITH LOW-THRUST TRAJECTORY OPTIMIZATION

ABOUT COMBINING TISSERAND GRAPH GRAVITY-ASSIST SEQUENCING WITH LOW-THRUST TRAJECTORY OPTIMIZATION ABOUT COMBINING TISSERAND GRAPH GRAVITY-ASSIST SEQUENCING WITH LOW-THRUST TRAJECTORY OPTIMIZATION Volker Maiwald German Aerospace Center (DLR) Institute of Space Systems Department of System Analysis Space

More information

Orbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design

Orbital Mechanics MARYLAND U N I V E R S I T Y O F. Orbital Mechanics. ENAE 483/788D - Principles of Space Systems Design Planetary launch and entry overview Energy and velocity in orbit Elliptical orbit parameters Orbital elements Coplanar orbital transfers Noncoplanar transfers Time in orbit Interplanetary trajectories

More information

A = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc.

A = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc. Q13.1 The mass of the Moon is 1/81 of the mass of the Earth. Compared to the gravitational force that the Earth exerts on the Moon, the gravitational force that the Moon exerts on the Earth is A. 81 2

More information

Rigorous Global Optimization of Impulsive Space Trajectories

Rigorous Global Optimization of Impulsive Space Trajectories Rigorous Global Optimization of Impulsive Space Trajectories P. Di Lizia, R. Armellin, M. Lavagna K. Makino, M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16 19, 2006 Motivation

More information

Ay 1 Lecture 2. Starting the Exploration

Ay 1 Lecture 2. Starting the Exploration Ay 1 Lecture 2 Starting the Exploration 2.1 Distances and Scales Some Commonly Used Units Distance: Astronomical unit: the distance from the Earth to the Sun, 1 au = 1.496 10 13 cm ~ 1.5 10 13 cm Light

More information

SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS. Sandro da Silva Fernandes. Cleverson Maranhão Porto Marinho

SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS. Sandro da Silva Fernandes. Cleverson Maranhão Porto Marinho SUN INFLUENCE ON TWO-IMPULSIVE EARTH-TO-MOON TRANSFERS Sandro da Silva Fernandes Instituto Tecnológico de Aeronáutica, São José dos Campos - 12228-900 - SP-Brazil, (+55) (12) 3947-5953 sandro@ita.br Cleverson

More information

Powered Space Flight

Powered Space Flight Powered Space Flight KOIZUMI Hiroyuki ( 小泉宏之 ) Graduate School of Frontier Sciences, Department of Advanced Energy & Department of Aeronautics and Astronautics ( 基盤科学研究系先端エネルギー工学専攻, 工学系航空宇宙工学専攻兼担 ) Scope

More information

is a revolution relative to a fixed celestial position. is the instant of transit of mean equinox relative to a fixed meridian position.

is a revolution relative to a fixed celestial position. is the instant of transit of mean equinox relative to a fixed meridian position. PERIODICITY FORMULAS: Sidereal Orbit Tropical Year Eclipse Year Anomalistic Year Sidereal Lunar Orbit Lunar Mean Daily Sidereal Motion Lunar Synodical Period Centenial General Precession Longitude (365.25636042

More information

PATHFINDING AND V-INFININTY LEVERAGING FOR PLANETARY MOON TOUR MISSIONS

PATHFINDING AND V-INFININTY LEVERAGING FOR PLANETARY MOON TOUR MISSIONS AAS 09-222 PATHFINDING AND V-INFININTY LEVERAGING FOR PLANETARY MOON TOUR MISSIONS Adam T. Brinckerhoff * and Ryan P. Russell The well established technique of V-infinity leveraging is applied to the phasefixed

More information

Earth-Mars Halo to Halo Low Thrust

Earth-Mars Halo to Halo Low Thrust Earth-Mars Halo to Halo Low Thrust Manifold Transfers P. Pergola, C. Casaregola, K. Geurts, M. Andrenucci New Trends in Astrodynamics and Applications V 3 June / -2 July, 28 Milan, Italy Outline o Introduction

More information

Radial Acceleration. recall, the direction of the instantaneous velocity vector is tangential to the trajectory

Radial Acceleration. recall, the direction of the instantaneous velocity vector is tangential to the trajectory Radial Acceleration recall, the direction of the instantaneous velocity vector is tangential to the trajectory 1 Radial Acceleration recall, the direction of the instantaneous velocity vector is tangential

More information

AP Physics C Textbook Problems

AP Physics C Textbook Problems AP Physics C Textbook Problems Chapter 13 Pages 412 416 HW-16: 03. A 200-kg object and a 500-kg object are separated by 0.400 m. Find the net gravitational force exerted by these objects on a 50.0-kg object

More information

MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30.

MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30. MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 4 Due Thursday, July 30. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details that

More information

orbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG

orbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits

More information

18. Kepler as a young man became the assistant to A) Nicolaus Copernicus. B) Ptolemy. C) Tycho Brahe. D) Sir Isaac Newton.

18. Kepler as a young man became the assistant to A) Nicolaus Copernicus. B) Ptolemy. C) Tycho Brahe. D) Sir Isaac Newton. Name: Date: 1. The word planet is derived from a Greek term meaning A) bright nighttime object. B) astrological sign. C) wanderer. D) nontwinkling star. 2. The planets that were known before the telescope

More information

GRAVITATION. F = GmM R 2

GRAVITATION. F = GmM R 2 GRAVITATION Name: Partner: Section: Date: PURPOSE: To explore the gravitational force and Kepler s Laws of Planetary motion. INTRODUCTION: Newton s law of Universal Gravitation tells us that the gravitational

More information

CHAPTER 8 PLANETARY MOTIONS

CHAPTER 8 PLANETARY MOTIONS 1 CHAPTER 8 PLANETARY MOTIONS 8.1 Introduction The word planet means wanderer (πλάνητες αστέρες wandering stars); in contrast to the fixed stars, the planets wander around on the celestial sphere, sometimes

More information

t S 18. Determining Planetary Co-ordinates

t S 18. Determining Planetary Co-ordinates 8. Determining Planetary Co-ordinates θ θ 0 ω R In the heliocentric reference frame which rotates with the Earth s orbital motion, suppose that initially a planet s unplanet line makes an angle θ 0 with

More information