42 min read

> "The Earth is the cradle of humanity, but one cannot live in the cradle forever."

Prerequisites

  • 6
  • 10

Learning Objectives

  • Apply the patched-conic approximation and compute a body's sphere of influence, and say when the approximation is and isn't valid.
  • Extend the Hohmann transfer to a heliocentric Earth-to-Mars trajectory and compute its departure and arrival hyperbolic excess velocities, transfer time, and total delta-v.
  • Define the synodic period and compute the launch window and phase angle for an interplanetary transfer.
  • Define characteristic energy C3 and hyperbolic excess velocity, and size the departure burn from a low parking orbit.
  • Describe the arrival options — orbit insertion, aerobraking, and aerocapture — and compute a Mars orbit-insertion burn.
  • Explain how a gravity assist changes a spacecraft's heliocentric energy by re-pointing its velocity in a moving planet's frame.

Chapter 11: Interplanetary Trajectories

"The Earth is the cradle of humanity, but one cannot live in the cradle forever." — Konstantin Tsiolkovsky

Overview

For ten chapters we have orbited one planet. Now we leave.

This is the chapter where everything we have built — the rocket equation, orbital energy, the vis-viva equation, the Hohmann transfer — is turned loose on the oldest ambition in the book: going to another world. And here is the wonderful thing. We do not need any new physics to do it. The same inverse-square gravity that holds a satellite in low Earth orbit holds Earth in orbit around the Sun, and holds Mars, and holds a probe drifting between them. Interplanetary flight is not a different subject from orbital mechanics; it is orbital mechanics played on a bigger table, with the Sun at the center instead of the Earth. The one genuinely new idea is a bookkeeping trick — patched conics — for handling the fact that a spacecraft on its way to Mars falls, at different times, under the sway of three different gravitational masters: Earth as it departs, the Sun for the long cruise, and Mars as it arrives.

We are going to cash in a promise made back in Chapter 6. There we computed the delta-v to escape Earth and hinted at the "trans-Mars injection" that appears on the delta-v map; we even quoted a hyperbolic excess speed of about $2.94\ \text{km/s}$ and promised to justify it here. This is where we make good. By the end of the chapter you will have computed the entire Earth-to-Mars Hohmann trajectory from first principles: how fast you must be going when you leave, how long the coast takes (spoiler: most of a year), why you can only leave every twenty-six months, what energy the launch vehicle must deliver, and how hard you must brake when you arrive. This is the climax of the Hohmann-to-Mars thread that has run through the book since Chapter 6.

Then we do something that still feels like cheating the first time you see it. We let a spacecraft steal energy from a planet — flying close past Jupiter and coming away faster than it arrived, for free — the gravity assist that flung the two Voyager probes out of the solar system. It is the single most elegant trick in all of orbital mechanics, and it is nothing but Newton, applied with nerve.

In this chapter, you will learn to:

  • Stitch a trajectory together out of two-body conics using the sphere of influence — patched conics.
  • Compute a complete Earth-to-Mars Hohmann transfer: the departure and arrival speeds, the ~259-day cruise, and the ~5.6 km/s of heliocentric delta-v it demands.
  • Find the launch window from the synodic period, and see why Mars opens its door only every 26 months.
  • Use $C_3$ and hyperbolic excess velocity to size the real departure burn from low Earth orbit.
  • Understand arrival — orbit insertion, aerobraking, aerocapture — and the physics of the gravity assist.

Learning Paths

🚀 Space Enthusiast: Read 11.1 for the big picture (spheres of influence), then spend your time on 11.2 (the trip to Mars), 11.3 (why you must wait for a window), and 11.6 (the gravity assist — the most beautiful idea in the chapter). You can take the departure and arrival delta-v numbers in 11.4–11.5 as given.

📐 Engineering Student: Read everything and reproduce the Mars transfer in 11.2 with a calculator before reading our numbers. Sections 11.4 and 11.5 are the ones you will use in real trajectory work; the ⭐⭐/⭐⭐⭐ exercises build the porkchop-plot intuition that professionals live by.

🎮 KSP Player: You have flown all of this. Focus on 11.1 (KSP's "sphere of influence" is exactly this chapter's idea, made literal), 11.3 (the phase-angle indicator on your transfer window), and 11.6 (every Jool slingshot you have ever done). The numbers here are the real solar system's version of your transfer-window planner.

🛰️ Industry Prep: The vocabulary of this chapter — $C_3$, $v_\infty$, launch window, DLA — is the language of a launch-services contract and a mission's trajectory design. Section 11.4 explains the number a launch provider actually sells you. Your Mission Design Checkpoint adds the interplanetary leg.


11.1 Patched conics and spheres of influence

Here is the problem, stated honestly. A spacecraft flying from Earth to Mars is pulled on, at every instant, by the Sun, by Earth, by Mars, and by every other body in the solar system. Even keeping only the three that matter — Sun, Earth, spacecraft — we have the notorious three-body problem, which, unlike the clean two-body orbits of Chapter 8, has no closed-form solution. You cannot write down an equation for the spacecraft's position as a function of time. (We will look this beautiful monster in the eye in Chapter 15.) If interplanetary flight really required solving the full many-body problem before you could plan a trajectory, no probe would ever have left Earth.

It does not, because of a saving observation about the sizes of things. Although the Sun's gravity reaches everywhere, near a planet the planet's gravity dominates. A spacecraft in low Earth orbit feels Earth pulling it roughly a thousand times harder than the Sun does; out in deep space between the planets, the Sun wins by a mile. There is, for each planet, a fuzzy boundary — a bubble around it — inside which you may pretend the planet is the only gravitating body, and outside which you may pretend the Sun is. We give that bubble a name.

Definition (sphere of influence). The sphere of influence (SOI) of a planet is the region around it within which the planet's gravity, not the Sun's, is treated as the dominant force on a spacecraft. Its radius is well approximated by $$r_{\text{SOI}} \approx a_{\text{planet}}\left(\frac{m_{\text{planet}}}{m_{\text{Sun}}}\right)^{2/5},$$ where $a_{\text{planet}}$ is the planet's distance from the Sun and $m_{\text{planet}}/m_{\text{Sun}}$ is its mass ratio to the Sun. Inside the SOI we model the trajectory as a two-body orbit about the planet; outside it, as a two-body orbit about the Sun.

The exponent of $2/5$ comes from comparing the planet's gravitational pull to the Sun's tidal disturbance, and you should treat the formula as a well-earned rule of thumb rather than a sharp edge — the "sphere" is not really a sphere and the boundary is not really a wall. But it gives excellent numbers, and it licenses the trick that makes this whole chapter possible.

Definition (patched conics). The patched-conic approximation models an interplanetary trajectory as a sequence of two-body conic-section orbits — one per gravitational body — joined ("patched") at the boundaries of the spheres of influence. A Mars mission becomes three conics: a departure hyperbola relative to Earth, a heliocentric transfer ellipse about the Sun, and an arrival hyperbola relative to Mars. Each piece is a clean two-body problem we already know how to solve; we hand the spacecraft off from one to the next at each SOI crossing.

💡 Intuition: Think of interplanetary travel like a series of connecting flights, each operated by a different gravitational "airline." Earth's gravity flies you out to the edge of Earth's SOI and hands you off to the Sun. The Sun flies you across the vast heliocentric distance to the edge of Mars's SOI and hands you off to Mars. Mars's gravity flies the final approach. At each handoff you simply match up the position and velocity — the "connecting gate" is the SOI boundary — and switch which body you call "the center." No single airline flies the whole route, but the trip still gets planned, one clean leg at a time.

Let us compute a sphere of influence, because the numbers reveal why patching works so well.

Worked Example: Earth's sphere of influence. Earth orbits at $a_\oplus = 1.496\times10^8\ \text{km}$ (1 AU), and its mass is $m_\oplus/m_\odot = 3.00\times10^{-6}$ of the Sun's. Then $$r_{\text{SOI}} \approx 1.496\times10^8 \times \left(3.00\times10^{-6}\right)^{2/5}\ \text{km}.$$ Take the two-fifths power: $\left(3.00\times10^{-6}\right)^{2/5}$. The fifth root of $3.00\times10^{-6}$ is about $0.0786$, and squaring gives $\approx 6.18\times10^{-3}$. So $$r_{\text{SOI}} \approx 1.496\times10^8 \times 6.18\times10^{-3} \approx 9.25\times10^5\ \text{km} > \approx 924{,}000\ \text{km}.$$ Earth's SOI reaches only about $924{,}000\ \text{km}$ — roughly $2.4$ times the distance to the Moon. Sanity-check the size: that is $9.24\times10^5 / 1.496\times10^8 \approx 0.006\ \text{AU}$, a mere 0.6% of the way to the Sun. The spacecraft spends the overwhelming majority of its journey out in the Sun's domain; the Earth and Mars hyperbolas are brief episodes at the very start and end. That enormous separation of scales is exactly why treating the pieces independently works so beautifully.

For reference, the same formula gives Mars an SOI of about $577{,}000\ \text{km}$ and the Moon an SOI of about $66{,}000\ \text{km}$ (relative to Earth). Each is tiny compared with the distance the spacecraft travels between them — which is the geometric fact patched conics leans on.

🔧 Engineering Reality: Patched conics is an approximation, and mission designers know exactly where it lies. It ignores the fact that inside Earth's SOI the Sun is still tugging (a little), and that near the SOI boundary the two pulls are comparable. For preliminary design — sizing the delta-v budget, picking a launch window, choosing a launch vehicle — patched conics is accurate to a percent or two and is used universally. For the actual flown trajectory, navigators throw it away and integrate the full equations of motion numerically, including every planet, the Sun's finite size, solar radiation pressure, and more (the perturbations of Chapter 12). The rule: patched conics to design the mission; full numerical integration to fly it. This chapter designs.

🔄 Check Your Understanding 1. Why can't we just solve the exact equations of motion for a spacecraft, the Sun, Earth, and Mars all at once? 2. Earth's SOI radius is about $924{,}000\ \text{km}$, only $0.6\%$ of an AU. Why is that smallness good news for the patched-conic method?

Answers

  1. Because three or more mutually gravitating bodies form the three-body (or $n$-body) problem, which has no closed-form solution — there is no formula for position versus time. We can only integrate it numerically. Patched conics sidesteps this by never having more than one gravitating body active at a time. 2. Because it means the "planet-dominated" legs are geometrically tiny compared with the Sun-dominated cruise. The scales are so separated that the small error made by pretending the boundary is sharp — and that only one body acts inside each region — barely affects the answer. If the SOIs were large and overlapping, patching would fail.

11.2 The interplanetary Hohmann

We now compute the trip to Mars. In Chapter 10 you derived the Hohmann transfer — the fuel-cheapest way to move between two circular, coplanar orbits: fire once to enter an ellipse tangent to the inner orbit at one end and the outer orbit at the other, coast halfway around, and fire again to circularize. Every word of that carries over to interplanetary flight. The only change is the cast: the central body is no longer Earth but the Sun, and the two "circular orbits" are no longer LEO and GEO but the orbits of Earth and Mars around the Sun.

Strategy first. Ignore Earth's and Mars's gravity for the moment (patched conics lets us — we are now out in the Sun's domain, between the SOIs). Treat Earth's orbit as a circle of radius $r_1 = 1\ \text{AU}$ and Mars's orbit as a circle of radius $r_2 = 1.524\ \text{AU}$. Build the Hohmann transfer ellipse between them exactly as in Chapter 10, but with the Sun's gravitational parameter. Use vis-viva (Chapter 6) to get the spacecraft's speed at each end, and compare it with the planet's own orbital speed. The difference between them is the speed the spacecraft needs relative to the planet — which is exactly what the departure and arrival hyperbolas must supply. That is the whole plan.

Our central-body constant is the Sun's gravitational parameter, $\mu_\odot = 1.327\times10^{11}\ \text{km}^3/\text{s}^2$ (Appendix B). Distances, converting AU to km with $1\ \text{AU} = 1.496\times10^8\ \text{km}$:

$$r_1 = 1.496\times10^8\ \text{km} \quad(\text{Earth}), \qquad r_2 = 1.524\times1.496\times10^8 = 2.280\times10^8\ \text{km} \quad(\text{Mars}).$$

Step 1 — the planets' orbital speeds. Each planet is (to good approximation) on a circular heliocentric orbit, so its speed is the circular velocity $\sqrt{\mu_\odot/r}$:

$$v_1 = \sqrt{\frac{\mu_\odot}{r_1}} = \sqrt{\frac{1.327\times10^{11}}{1.496\times10^8}} = \sqrt{887.0} = 29.78\ \text{km/s} \quad(\text{Earth}),$$ $$v_2 = \sqrt{\frac{\mu_\odot}{r_2}} = \sqrt{\frac{1.327\times10^{11}}{2.280\times10^8}} = \sqrt{582.0} = 24.13\ \text{km/s} \quad(\text{Mars}).$$

Sanity check: Earth's mean orbital speed is famously about $29.8\ \text{km/s}$, and Mars's about $24\ \text{km/s}$ — both land exactly. Already a good sign.

Step 2 — the transfer ellipse. Its semi-major axis is half the sum of the two orbital radii:

$$a_t = \frac{r_1 + r_2}{2} = \frac{1.496\times10^8 + 2.280\times10^8}{2} = 1.888\times10^8\ \text{km} \quad(= 1.262\ \text{AU}).$$

Step 3 — the spacecraft's speeds at the ends of the transfer. From vis-viva, $v = \sqrt{\mu_\odot\left(\tfrac{2}{r} - \tfrac{1}{a_t}\right)}$. At perihelion (the low end, at Earth's orbit — this is departure):

$$v_p = \sqrt{1.327\times10^{11}\left(\frac{2}{1.496\times10^8} - \frac{1}{1.888\times10^8}\right)} = \sqrt{1.327\times10^{11}\times 8.072\times10^{-9}} = \sqrt{1071} = 32.73\ \text{km/s}.$$

At aphelion (the high end, at Mars's orbit — this is arrival):

$$v_a = \sqrt{1.327\times10^{11}\left(\frac{2}{2.280\times10^8} - \frac{1}{1.888\times10^8}\right)} = \sqrt{1.327\times10^{11}\times 3.476\times10^{-9}} = \sqrt{461.2} = 21.48\ \text{km/s}.$$

Here already is the whole story of the transfer, sitting in four numbers. At Earth's orbit the spacecraft must move at $32.73\ \text{km/s}$, but Earth is only moving at $29.78$; so the spacecraft must be going faster than Earth by

$$v_{\infty,\oplus} = v_p - v_1 = 32.73 - 29.78 = 2.95\ \text{km/s}.$$

And at Mars's orbit the spacecraft is moving at only $21.48\ \text{km/s}$ while Mars sails past at $24.13$; so the spacecraft arrives slower than Mars by

$$v_{\infty,\text{Mars}} = v_2 - v_a = 24.13 - 21.48 = 2.65\ \text{km/s}.$$

Those two leftover speeds — $2.95\ \text{km/s}$ at Earth and $2.65\ \text{km/s}$ at Mars — are the hyperbolic excess velocities, the speeds the spacecraft has relative to each planet at the edge of its sphere of influence. They are the hinge between the heliocentric transfer and the planet-centered hyperbolas, and we will spend Sections 11.4 and 11.5 turning them into actual rocket burns. Notice the first one matches the $\approx 2.94\ \text{km/s}$ that Chapter 6 quoted on faith — promise kept.

   Earth at departure                                   Mars at arrival
   (perihelion of transfer)                          (aphelion of transfer)
        ●  v_p = 32.73 km/s                              ★  v_a = 21.48 km/s
         \  Earth moves 29.78  ──► faster by             /  Mars moves 24.13 ──► slower by
          \       v∞ = 2.95 km/s (C3 ≈ 8.7)             /        v∞ = 2.65 km/s
           \                                           /
            \            transfer ellipse             /
             \        a_t = 1.262 AU, half-orbit     /
              \                                      /
               \______________  ☉  _______________/
                r1 = 1 AU      Sun      r2 = 1.524 AU
              (the probe coasts 180° around the Sun, engine off, for ~259 days)

Step 4 — the total heliocentric delta-v. In the idealized picture where the spacecraft is the planet's orbit (we will refine this in 11.4–11.5), the delta-v is just the sum of the two speed gaps:

$$\Delta v_{\text{helio}} = v_{\infty,\oplus} + v_{\infty,\text{Mars}} = 2.95 + 2.65 = 5.60\ \text{km/s}.$$

That $\approx 5.6\ \text{km/s}$ is the canonical "heliocentric cost of a Mars Hohmann," and you will see it quoted throughout mission design.

Step 5 — the transfer time. The spacecraft coasts exactly half of the transfer ellipse, so the trip takes half the ellipse's period. Kepler's third law (Chapter 8) gives the period as $T = 2\pi\sqrt{a_t^3/\mu_\odot}$, so the transfer time is half of that:

$$t_{\text{trans}} = \pi\sqrt{\frac{a_t^3}{\mu_\odot}} = \pi\sqrt{\frac{(1.888\times10^8)^3}{1.327\times10^{11}}} = \pi\sqrt{5.071\times10^{13}} = \pi\times7.121\times10^6 = 2.237\times10^7\ \text{s}.$$

Convert to days: $2.237\times10^7\ \text{s} \div 86{,}400\ \text{s/day} = 259\ \text{days}$ — about 8.5 months. Sanity-check against reality: real Mars missions cruise for roughly six to nine months, and a minimum-energy Hohmann sits right at the long end. The number is right.

🔗 Connection: the counterintuitive arrival. Look again at what happens at Mars. The spacecraft slows down on the way out — it is climbing the Sun's gravity well from 1 to 1.524 AU, trading kinetic energy for potential exactly as Chapter 6 taught — so by the time it reaches Mars's orbit it is dawdling along at $21.5\ \text{km/s}$, well below Mars's $24.1$. Mars runs into the spacecraft from behind. This is the same "higher is slower" law that governs a LEO-to-GEO transfer, now stretched across the solar system. The spacecraft does not chase Mars down; it climbs to Mars's altitude, arrives moving slowly, and lets Mars catch up to it.

📜 From History: Hohmann's paper trip. Walter Hohmann worked out the minimum-energy transfer in his 1925 book Die Erreichbarkeit der Himmelskörper ("The Attainability of Celestial Bodies"), decades before any rocket could fly it. Like Tsiolkovsky's rocket equation, the interplanetary Hohmann transfer was pure theory that sat waiting for the hardware to catch up. When Mariner 4 flew past Mars in 1965 — humanity's first successful Mars encounter — it flew, in essence, the ellipse Hohmann had drawn on paper forty years earlier. The mathematics of this section is not a modern computer's invention; it is a century old, and it predicted the space age before the space age existed.

🔄 Check Your Understanding 1. At Mars's orbit the spacecraft moves at $21.5\ \text{km/s}$ but Mars moves at $24.1\ \text{km/s}$. Does the spacecraft need to speed up or slow down (relative to the Sun) to match Mars's circular orbit? What does that tell you about the direction of the arrival burn? 2. Why is the transfer time fixed at 259 days once you have chosen a Hohmann transfer — what would you have to give up to make the trip faster?

Answers

  1. It must speed up relative to the Sun (from $21.5$ to $24.1\ \text{km/s}$) to join Mars's faster circular orbit. So, in the Sun's frame, the arrival maneuver adds prograde speed. (In Mars's frame, though, the spacecraft is overtaken by the planet, so it arrives with $v_\infty$ and must brake into a bound orbit — Section 11.5.) 2. The Hohmann ellipse is fixed the moment you fix its two ends ($r_1$ and $r_2$); its semi-major axis, and therefore its period and half-period, are then determined. To arrive sooner you must fly a faster, higher-energy trajectory — a non-tangential transfer with a larger semi-major axis or a different shape — which costs more delta-v at both ends. Speed costs fuel; the Hohmann is slow precisely because it is cheap.

11.3 Launch windows and the synodic period

Suppose your rocket is fueled and your Mars probe is bolted on. Can you launch today? Almost certainly not — and the reason is pure geometry. Your spacecraft will spend 259 days coasting to the point in space where Mars's orbit crosses the far end of your transfer ellipse. For the mission to work, Mars must actually be at that point when the spacecraft arrives. Launch on the wrong day and the spacecraft sails through Mars's orbit into empty space, months ahead of or behind the planet, and coasts on forever having accomplished nothing.

So Earth and Mars must be in the right relative positions at departure. Since the spacecraft sweeps $180^\circ$ around the Sun while Mars moves along its own orbit, Mars must start out ahead of Earth by just the right lead angle so that planet and probe reach the rendezvous point together. That lead angle — the phase angle at departure — works out to about $44^\circ$ for Earth-to-Mars. (Mars, moving at its own angular rate, covers about $136^\circ$ during the 259-day cruise; since the probe covers $180^\circ$, Mars must begin $180^\circ - 136^\circ = 44^\circ$ ahead.) The narrow span of dates around which this alignment holds is the launch window.

Definition (launch window). A launch window is the interval of time during which a spacecraft may depart and still reach its target with the planned trajectory and delta-v budget. For an interplanetary Hohmann transfer, it is set by the requirement that the departure and destination planets be correctly phased — the destination leading by the right angle — so the spacecraft and the planet arrive at the transfer point together. Miss the window and you must wait for the geometry to recur.

How long until the geometry recurs? Both planets keep orbiting, at different rates — Earth once a year, Mars once every $687$ days — so their alignment is not fixed; it drifts and then, eventually, repeats. The time between successive identical alignments is the synodic period.

Definition (synodic period). The synodic period $T_{\text{syn}}$ of two bodies orbiting the same primary is the time between successive identical relative alignments (for example, from one launch window to the next). It is set by the difference of their angular rates: $$\frac{1}{T_{\text{syn}}} = \left|\frac{1}{T_1} - \frac{1}{T_2}\right|,$$ where $T_1$ and $T_2$ are the two orbital periods. The faster inner body laps the slower outer one once per synodic period.

💡 Intuition: Two runners on a circular track, one faster than the other, line up side by side only every so often — not once per lap, but once per extra lap the faster runner gains on the slower. The synodic period is exactly that "time to lap." Earth is the faster inner runner; it must gain a full $360^\circ$ on Mars to restore the same starting geometry. Because their per-day angular speeds differ only a little, gaining a whole lap takes far longer than either runner's own lap time — which is why the Earth-Mars synodic period is longer than either planet's year.

Worked Example: the Earth-Mars launch window. Earth's orbital period is $T_1 = 365.25\ \text{days}$; Mars's is $T_2 = 686.98\ \text{days}$. Then $$\frac{1}{T_{\text{syn}}} = \frac{1}{365.25} - \frac{1}{686.98} > = (2.7379 - 1.4557)\times10^{-3} = 1.2822\times10^{-3}\ \text{day}^{-1},$$ so $$T_{\text{syn}} = \frac{1}{1.2822\times10^{-3}} = 780\ \text{days} \approx 2.14\ \text{years} > \approx 26\ \text{months}.$$ Mars launch opportunities open only about every 26 months. Sanity-check against the historical record: Mars missions do indeed cluster in windows roughly two years apart — 2018, 2020, 2022, 2024, and so on. (Not every window is used, but the opportunities recur on this 26-month drumbeat.) If you miss one, you wait more than two years for the next. This single number shapes the entire economics and politics of Mars exploration: a slipped schedule does not cost weeks, it costs a two-year cycle.

🔧 Engineering Reality: Real launch windows are days-to-weeks wide, not instantaneous, because a mission accepts a little extra delta-v to leave a few days early or late. Trajectory designers capture this with a porkchop plot — a contour map of required $C_3$ (and arrival $v_\infty$) as a function of launch date and arrival date, whose nested contours look like a pork chop. The bottom of the "chop" is the cheapest departure; the mission picks a launch period around it that trades a bit of performance for schedule robustness. We will not draw one here, but every real interplanetary mission is planned on one, and the single point at its minimum is the idealized Hohmann we just computed.

⚠️ Common Misconception: "Why not just launch whenever and adjust with the engine?" Because the adjustment is ruinously expensive. To leave off-window and still hit Mars, you would have to fly a steeply non-optimal trajectory — arriving where Mars will be rather than where the cheap ellipse goes — and the extra delta-v can dwarf the entire mission budget. Recall the tyranny of the rocket equation (Chapter 3): a few extra km/s is not a few extra percent of fuel, it is often a different, impossible rocket. It is almost always cheaper to wait for the geometry than to pay to fight it. Patience is a propellant.


11.4 $C_3$ and hyperbolic departure

We have the heliocentric picture. Now we zoom back in to Earth and ask the practical question: what must the launch vehicle actually do to put the spacecraft on that transfer ellipse? This is where the departure hyperbola — the first patched conic — earns its keep.

The spacecraft does not begin in interplanetary space moving at Earth's $29.78\ \text{km/s}$. It begins in a low parking orbit around Earth, and it must leave Earth's sphere of influence with a leftover speed of $v_{\infty,\oplus} = 2.95\ \text{km/s}$ relative to Earth, pointed the right way, so that once it is out in the Sun's domain it adds to Earth's orbital motion and rides the transfer ellipse. Inside Earth's SOI, a trajectory that reaches the boundary still moving is — by the energy classification of Chapter 6 — a hyperbola. The speed it retains at the far edge is what we now name formally.

Definition (hyperbolic excess velocity). The hyperbolic excess velocity $v_\infty$ is the speed a spacecraft retains relative to a body after it has climbed entirely out of that body's gravity well — its speed "at infinity," i.e. at the edge of the sphere of influence. It is set by the trajectory's specific energy through $\varepsilon = v_\infty^2/2$: a bound orbit has $\varepsilon<0$ and no $v_\infty$ (it never escapes); a parabola has $\varepsilon = 0$ and $v_\infty = 0$ (it just barely escapes); a hyperbola has $\varepsilon > 0$ and $v_\infty > 0$. For an interplanetary departure, $v_\infty$ is exactly the excess speed the heliocentric transfer demanded — here, $2.95\ \text{km/s}$.

Launch providers, however, do not quote you a $v_\infty$. They quote a $C_3$ — and the two are the same idea, squared.

Definition (characteristic energy, $C_3$). The characteristic energy $C_3$ is twice the specific orbital energy of a departure trajectory — equivalently, the square of the hyperbolic excess velocity: $$C_3 \equiv v_\infty^2 = 2\varepsilon = -\frac{\mu}{a}.$$ It has units of $\text{km}^2/\text{s}^2$ (energy per unit mass). A larger $C_3$ means a more energetic departure — a faster escape, reaching a more distant destination or arriving sooner. $C_3 = 0$ is a parabolic (just-barely-escape) trajectory. $C_3$ is the number a launch vehicle's performance is rated against: every launcher has a published curve of deliverable payload mass versus $C_3$.

For our Mars departure,

$$C_3 = v_{\infty,\oplus}^2 = (2.95\ \text{km/s})^2 \approx 8.7\ \text{km}^2/\text{s}^2.$$

That single number — "we need a $C_3$ of about $8.7\ \text{km}^2/\text{s}^2$" — is what a mission tells a launch provider, and the provider's payload-vs-$C_3$ chart says how much mass they can send. (Real Earth-to-Mars departures run $C_3 \approx 8$–$17\ \text{km}^2/\text{s}^2$ depending on the year's geometry; our $8.7$ is a favorable, near-minimum window — an illustrative value, since the true figure is date-dependent and needs an ephemeris.)

Now the burn itself. The spacecraft sits in a low circular parking orbit — take $300\ \text{km}$ altitude, so $r_{\text{park}} = 6{,}371 + 300 = 6{,}671\ \text{km}$, using Earth's mean radius from Appendix B. Its circular speed there is

$$v_{\text{circ}} = \sqrt{\frac{\mu_\oplus}{r_{\text{park}}}} = \sqrt{\frac{3.986\times10^5}{6{,}671}} = \sqrt{59.75} = 7.73\ \text{km/s}.$$

To be on the departure hyperbola, the spacecraft needs, at that same perigee radius, a speed set by energy conservation. Its hyperbola has specific energy $\varepsilon = v_\infty^2/2$, and $\varepsilon = v_{\text{peri}}^2/2 - \mu_\oplus/r_{\text{park}}$, so

$$v_{\text{peri}} = \sqrt{v_\infty^2 + \frac{2\mu_\oplus}{r_{\text{park}}}} = \sqrt{v_\infty^2 + v_{\text{esc}}^2},$$

where $v_{\text{esc}} = \sqrt{2\mu_\oplus/r_{\text{park}}} = \sqrt{119.5} = 10.93\ \text{km/s}$ is the escape speed from the parking orbit. Plugging in:

$$v_{\text{peri}} = \sqrt{(2.95)^2 + (10.93)^2} = \sqrt{8.68 + 119.5} = \sqrt{128.2} = 11.32\ \text{km/s}.$$

The trans-Mars injection (TMI) burn is the jump from the parking-orbit speed to this hyperbola speed:

$$\Delta v_{\text{TMI}} = v_{\text{peri}} - v_{\text{circ}} = 11.32 - 7.73 = 3.59\ \text{km/s}.$$

Sanity check: the delta-v map of Chapter 3 listed "LEO → trans-Mars injection $\approx 3.6\ \text{km/s}$." We just derived it from scratch — $3.59\ \text{km/s}$. Another promise kept.

Now stare at a small miracle in those numbers. We needed $2.95\ \text{km/s}$ of excess speed at infinity, yet the burn to produce it was only $3.59\ \text{km/s}$ — not the $2.95$ plus the full $10.93$ escape speed you might have feared. Where did the discount come from? From burning deep in Earth's gravity well, while already moving fast.

🚪 Threshold Concept: the Oberth effect makes deep burns cheap. A burn adds kinetic energy at a rate that depends on how fast you are already going: the energy gained is $\Delta\varepsilon = v\,\Delta v + \tfrac{1}{2}\Delta v^2$, and the first term rewards a large existing speed $v$. In the parking orbit you are screaming along at $7.73\ \text{km/s}$, so every km/s you add buys a lot of energy. If instead you coasted slowly out to the edge of the SOI and burned there, moving at a crawl, the same $\Delta v$ would buy far less energy, and you would need much more of it to reach the same $v_\infty$. This is the Oberth effect, met first in Chapter 6: burn low, burn fast. It is why interplanetary missions inject from a low parking orbit rather than a high one, why a powered flyby is so potent, and why "get deep, then burn hard" is a mantra of trajectory design. Leaving a planet is not about overpowering gravity — it is about spending your delta-v where it is worth the most.

Worked Example: reading a launcher's $C_3$ chart (illustrative). Suppose a launch vehicle's data sheet says it can send $4{,}000\ \text{kg}$ to $C_3 = 0$ (Earth escape) but only $2{,}500\ \text{kg}$ to $C_3 = 15\ \text{km}^2/\text{s}^2$. Your Mars mission needs $C_3 = 8.7$; interpolating (very roughly) between those points, the deliverable mass is around $4{,}000 - \tfrac{8.7}{15}(4{,}000-2{,}500) \approx 4{,}000 - 870 = 3{,}130\ \text{kg}$. If your fueled spacecraft masses $3{,}000\ \text{kg}$, you close — with thin margin. This is the actual conversation between a mission and a rocket: not "how much delta-v" but "what $C_3$, and how much mass at that $C_3$." The numbers here are illustrative (Tier 3), but the logic is exactly how a launch is bought. (The real payload-vs-$C_3$ curve is convex, not linear — this straight-line interpolation is only a back-of-envelope sketch.)

🔄 Check Your Understanding 1. A mission's $C_3$ is $8.7\ \text{km}^2/\text{s}^2$. What is its hyperbolic excess velocity, and what would $C_3 = 0$ physically mean? 2. Why is the trans-Mars injection burn ($3.59\ \text{km/s}$) so much less than $v_\infty + v_{\text{esc}} = 2.95 + 10.93 = 13.9\ \text{km/s}$?

Answers

  1. $v_\infty = \sqrt{C_3} = \sqrt{8.7} = 2.95\ \text{km/s}$. $C_3 = 0$ means a parabolic trajectory — just barely escaping Earth with zero speed left over, the boundary between bound and unbound. 2. Because the spacecraft is already moving at $7.73\ \text{km/s}$ in the parking orbit; you only pay for the difference between the hyperbola's perigee speed and the circular speed, not for the whole escape from a standstill. And the speeds add in quadrature ($v_{\text{peri}} = \sqrt{v_\infty^2 + v_{\text{esc}}^2}$), not linearly, because they are rooted in energy, which goes as speed squared — the Oberth discount.

11.5 Arrival: orbit insertion and aerocapture

The spacecraft coasts for 259 days and crosses into Mars's sphere of influence. Now the patched-conic handoff runs in reverse: relative to Mars, the spacecraft arrives with the hyperbolic excess velocity we computed, $v_{\infty,\text{Mars}} = 2.65\ \text{km/s}$. Left alone, it would whip once around Mars on an arrival hyperbola and fly straight back out into solar orbit — a flyby, which is exactly what the early Mariners did. To stay — to become a Mars orbiter or to land — the spacecraft must shed enough energy to drop from that hyperbola ($\varepsilon > 0$) onto a bound ellipse ($\varepsilon < 0$). There are two ways to pay: with the engine, or with the atmosphere.

Option 1 — propulsive orbit insertion. Fire the engine at closest approach (periapsis) to slow down into a captured orbit. As at departure, burning at periapsis is cheapest — the Oberth effect again. Let us size the burn to capture into a low circular orbit at $400\ \text{km}$ altitude, so $r_{\text{arr}} = 3{,}390 + 400 = 3{,}790\ \text{km}$ (Mars's mean radius is $3{,}390\ \text{km}$, and $\mu_{\text{Mars}} = 4.283\times10^4\ \text{km}^3/\text{s}^2$, both from Appendix B).

The speed at periapsis on the arrival hyperbola:

$$v_{\text{peri}} = \sqrt{v_{\infty,\text{Mars}}^2 + \frac{2\mu_{\text{Mars}}}{r_{\text{arr}}}} = \sqrt{(2.65)^2 + \frac{2\times4.283\times10^4}{3{,}790}} = \sqrt{7.02 + 22.60} = \sqrt{29.62} = 5.44\ \text{km/s}.$$

The circular speed we want to end up at:

$$v_{\text{circ}} = \sqrt{\frac{\mu_{\text{Mars}}}{r_{\text{arr}}}} = \sqrt{\frac{4.283\times10^4}{3{,}790}} = \sqrt{11.30} = 3.36\ \text{km/s}.$$

So the Mars orbit-insertion (MOI) burn to circularize is

$$\Delta v_{\text{MOI}} = v_{\text{peri}} - v_{\text{circ}} = 5.44 - 3.36 = 2.08\ \text{km/s}.$$

That is expensive — comparable to the entire trans-Mars injection. Which is why almost no real mission captures directly into a low circular orbit. Instead they capture into a highly elliptical orbit first, braking only enough to become barely bound. The minimum capture burn just nudges the periapsis speed below the local escape speed, $v_{\text{esc}}(r_{\text{arr}}) = \sqrt{22.60} = 4.75\ \text{km/s}$, so the smallest possible insertion is only $5.44 - 4.75 = 0.69\ \text{km/s}$ — a third of the cost. Real Mars orbiters (Mars Reconnaissance Orbiter, MAVEN) do exactly this, spending roughly $1$–$1.5\ \text{km/s}$ to drop onto a long ellipse, then lowering it over weeks with the atmosphere.

Option 2 — let the atmosphere do the braking. Mars has an atmosphere, thin but real, and a spacecraft can trade the engine for aerodynamic drag — dumping orbital energy as heat exactly as a returning capsule does at Earth (Chapter 7). Two flavors:

  • Aerobraking: after a cheap propulsive capture into a high ellipse, dip the periapsis into the upper atmosphere on each pass. Each grazing pass sheds a little energy to drag, slowly shrinking the orbit over weeks or months. It is slow but nearly free of propellant — the method that circularized MRO and Mars Odyssey.
  • Aerocapture: the aggressive version — fly directly into the atmosphere on arrival, in one pass, shedding enough energy to go straight from the arrival hyperbola to a bound orbit without a big capture burn at all. It could save that entire $\sim 2\ \text{km/s}$ of MOI, but it demands a heat shield and pinpoint guidance through a razor-thin corridor; it has been used at Earth and studied intensively for Mars but not yet flown there.

🔗 Connection: the same energy, paid two ways. Whether you brake with a rocket or with an atmosphere, you are solving the identical problem — removing the positive energy that makes your orbit a hyperbola. The engine converts that energy back into (unburned) propellant's chemical energy you chose not to spend; the atmosphere converts it into heat. This is the arrival mirror of the departure Oberth trick, and it is the direct setup for the Mars mission's "seven minutes of terror" — entry, descent, and landing — that we dissect in Chapter 34.

🐛 Find the Error. A student computes the Mars arrival speed by taking the spacecraft's heliocentric aphelion speed, $21.48\ \text{km/s}$, and using that as the $v_\infty$ at Mars: "the spacecraft hits Mars's atmosphere at $21.5\ \text{km/s}$, so the heat shield must handle $21.5\ \text{km/s}$." Where is the error, and what is the right approach speed?

Answer

The $21.48\ \text{km/s}$ is the spacecraft's speed relative to the Sun, not relative to Mars. What matters for arrival is the speed relative to Mars, which is the hyperbolic excess $v_\infty = v_2 - v_a = 24.13 - 21.48 = 2.65\ \text{km/s}$ — the two heliocentric velocities nearly cancel because Mars and the spacecraft are moving the same direction at similar speeds. The spacecraft then falls into Mars's well, so at the top of the atmosphere ($r \approx 3{,}490\ \text{km}$) it is doing $\sqrt{v_\infty^2 + 2\mu_{\text{Mars}}/r} = \sqrt{2.65^2 + 2(4.283\times10^4)/3{,}490} \approx \sqrt{7.0 + 24.5} = 5.6\ \text{km/s}$, not $21.5$. Confusing the heliocentric speed with the planet-relative speed is the single most common patched-conic mistake — always ask "relative to which body?"


11.6 Gravity assists: stealing energy from a planet

We come to the most elegant trick in orbital mechanics — the one that looks, the first time you meet it, like a violation of energy conservation. A spacecraft flies close past a planet, fires no engine, and comes away moving faster. Voyager 2 did it four times, riding flybys of Jupiter, Saturn, Uranus, and Neptune to reach the edge of the solar system on a trickle of propellant. Where does the extra speed come from? Nowhere magical — it comes, quite literally, out of the planet's orbit. Let us see how.

Definition (gravity assist). A gravity assist (or gravitational slingshot) is a maneuver in which a spacecraft flies through a planet's sphere of influence and exploits the planet's motion around the Sun to change its own heliocentric speed and direction — gaining or losing orbital energy with no propellant. The energy is transferred to or from the planet's orbital motion.

The key is to watch the flyby from two frames at once.

In the planet's frame, the encounter is a plain hyperbola — the arrival conic of Section 11.5, but one we do not brake out of. The spacecraft falls in from the edge of the SOI at speed $v_\infty$, whips around the planet, and climbs back out to the edge. Because the planet's frame is (locally) inertial and no propellant is burned, energy is conserved in this frame: the spacecraft leaves with exactly the same speed $v_\infty$ it arrived with. All the flyby does, in the planet's frame, is rotate the velocity vector by some turning angle $\delta$ — like a ball on a string swinging around a post. Nothing is gained; the speed relative to the planet is untouched.

In the Sun's frame, the planet is moving — Jupiter at $13\ \text{km/s}$, Mars at $24$, and so on. The spacecraft's heliocentric velocity is the vector sum of its planet-relative velocity and the planet's own velocity: $\mathbf{v}_{\text{helio}} = \mathbf{v}_\infty + \mathbf{v}_{\text{planet}}$. The flyby did not change the length of $\mathbf{v}_\infty$, but it changed its direction — and when you add a rotated $\mathbf{v}_\infty$ to the planet's velocity, the resulting heliocentric speed can be very different. Turn $\mathbf{v}_\infty$ to point more along the planet's motion, and the spacecraft speeds up in the Sun's frame. Turn it against, and the spacecraft slows down. The planet's velocity is the free ingredient the spacecraft borrows.

💡 Intuition: an elastic bounce off a moving wall. Throw a tennis ball at $10\ \text{m/s}$ toward a wall that is rushing toward you at $13\ \text{m/s}$. In the wall's frame the ball comes in at $23\ \text{m/s}$ and (elastically) bounces back at $23\ \text{m/s}$. Back in your frame, the wall's $13\ \text{m/s}$ adds on again: the ball returns at $23 + 13 = 36\ \text{m/s}$ — it gained $2\times13 = 26\ \text{m/s}$, twice the wall's speed, without any engine. A gravity assist is the gravitational, frictionless, curved version of that bounce: the planet is the moving wall, gravity is the "spring" that turns the ball around, and the speed the spacecraft steals comes from the wall's motion. The ball slowed the wall by an utterly imperceptible amount — and so does the spacecraft slow the planet.

Let us make it concrete with a Jupiter flyby, the workhorse of the outer solar system.

Worked Example: a Jupiter gravity assist (illustrative). Let Jupiter move in the $+x$ direction at its orbital speed $v_J = 13\ \text{km/s}$. A spacecraft approaches with a hyperbolic excess of $v_\infty = 10\ \text{km/s}$ relative to Jupiter — say, coming in "downward," $\mathbf{v}_\infty^{\text{in}} = (0,\,-10)\ \text{km/s}$. Its heliocentric velocity before the flyby is $$\mathbf{v}_{\text{before}} = \mathbf{v}_\infty^{\text{in}} + \mathbf{v}_J = (0,-10) + (13,0) = (13,\,-10)\ \text{km/s},$$ with speed $\sqrt{13^2 + 10^2} = \sqrt{269} = 16.40\ \text{km/s}$. Now suppose the flyby geometry turns the planet-relative velocity by $90^\circ$, so it leaves pointing along Jupiter's motion: $\mathbf{v}_\infty^{\text{out}} = (10,\,0)\ \text{km/s}$ — same length $10$, new direction. The heliocentric velocity after is $$\mathbf{v}_{\text{after}} = \mathbf{v}_\infty^{\text{out}} + \mathbf{v}_J = (10,0) + (13,0) = (23,\,0)\ \text{km/s},$$ with speed $23\ \text{km/s}$. The spacecraft's heliocentric speed jumped from $16.40$ to $23\ \text{km/s}$ — a gain of $6.6\ \text{km/s}$, and its specific orbital energy about the Sun rose from $\tfrac{1}{2}(16.40)^2 = 134\ \text{km}^2/\text{s}^2$ to $\tfrac{1}{2}(23)^2 = 265\ \text{km}^2/\text{s}^2$, nearly doubling — all without burning a drop of propellant. In the planet's frame the speed never changed ($v_\infty = 10\ \text{km/s}$ in and out); in the Sun's frame the probe robbed Jupiter of a whisper of orbital energy and ran off with it.

How big a turn can the planet provide? The turning angle depends on how close and how fast you pass:

$$\sin\!\left(\frac{\delta}{2}\right) = \frac{1}{e}, \qquad e = 1 + \frac{r_p\,v_\infty^2}{\mu_{\text{planet}}},$$

where $e$ is the eccentricity of the flyby hyperbola and $r_p$ is the periapsis distance. A closer pass (smaller $r_p$) or a slower approach (smaller $v_\infty$) bends the path more.

Worked Example: how hard can Jupiter turn you? Take our $v_\infty = 10\ \text{km/s}$ and a close pass at $r_p = 350{,}000\ \text{km}$ (about 5 Jupiter radii, roughly where Voyager 1 crossed), with $\mu_{\text{Jupiter}} = 1.267\times10^8\ \text{km}^3/\text{s}^2$. The eccentricity is $$e = 1 + \frac{(350{,}000)(10^2)}{1.267\times10^8} = 1 + \frac{3.50\times10^7}{1.267\times10^8} > = 1 + 0.276 = 1.276,$$ so $\sin(\delta/2) = 1/1.276 = 0.784$, giving $\delta/2 = 51.6^\circ$ and a turn of $\delta \approx 103^\circ$. Jupiter can swing a fast spacecraft's velocity through more than a right angle in a single pass — which is exactly why it, the most massive planet, is the great enabler of outer-solar-system missions. Our $90^\circ$ turn in the previous example is comfortably within Jupiter's reach. (Values are illustrative — Tier 3 — but of the right size for a real Jovian flyby.)

🚪 Threshold Concept: the assist is free to you, paid by the planet. Energy is not created in a gravity assist; it is transferred. By Newton's third law, the same gravitational tug that speeds the spacecraft up slows the planet down — the spacecraft gains exactly the orbital energy the planet loses. But the planet outmasses the spacecraft by a factor of order $10^{24}$, so its share of the exchange is unmeasurably small: a one-tonne probe that gains $6.6\ \text{km/s}$ slows Jupiter's orbital motion by about $(10^3/1.9\times10^{27})\times 6.6\ \text{km/s} \approx 3\times10^{-21}\ \text{km/s}$ — Jupiter's orbit shrinks by less than the width of a proton over the age of the universe. The books balance perfectly; the planet simply never notices. Once you see that the "free" energy is really a withdrawal from the planet's colossal orbital bank account, the paradox dissolves — and a whole new class of missions opens up. This is the idea we will generalize in Chapter 15 into the "interplanetary superhighway" of low-energy transfers.

📜 From History: the Grand Tour. In the mid-1960s a graduate student, Michael Minovitch, worked out that a single spacecraft could chain gravity assists to tour the outer planets, and a rare alignment of Jupiter, Saturn, Uranus, and Neptune in the late 1970s — one that recurs only every 176 years — made it possible to visit all four with one launch. NASA flew it as Voyager. Voyager 2 used Jupiter's assist to reach Saturn, Saturn's to reach Uranus, and Uranus's to reach Neptune, each flyby both bending the path to the next target and adding the energy to get there. Without gravity assists, the propellant to visit all four would have been impossible; with them, two modest probes launched in 1977 are now in interstellar space. We will return to the Grand Tour as inspiration in Chapter 39. It remains the most elegant single demonstration of everything in this chapter.

🔄 Check Your Understanding 1. In the planet's frame, a gravity assist does not change the spacecraft's speed at all. So how does the spacecraft end up faster in the Sun's frame? 2. Could you use a gravity assist to slow down — say, to help a probe fall toward the inner solar system? How?

Answers

  1. The flyby rotates the spacecraft's planet-relative velocity vector $\mathbf{v}_\infty$ (same length, new direction). Heliocentric velocity is $\mathbf{v}_\infty + \mathbf{v}_{\text{planet}}$; adding the rotated $\mathbf{v}_\infty$ to the planet's velocity yields a different (here, larger) heliocentric speed. The planet's own motion is the borrowed ingredient. 2. Yes — turn $\mathbf{v}_\infty$ to point against the planet's motion, and the heliocentric speed drops. Missions to Mercury (MESSENGER, BepiColombo) use exactly this: repeated Venus and Mercury flybys bleed off heliocentric energy so the probe can spiral inward and be captured, because arriving at Mercury with too much speed would make orbit insertion impossibly expensive. A slingshot can brake as well as boost.

Mission Design Checkpoint: your interplanetary leg and interplanetary.py

If your mission leaves Earth's neighborhood — a Mars orbiter (Track C) or an asteroid rendezvous (Track D) — this is the chapter where your Mission Design Review (MDR) gains its interplanetary leg. (Tracks A and B stay in Earth-Moon space; you can still build the module and run it on the Mars numbers for practice, and you will need the launch-window idea for any deep-space extension.)

The design. Add an "Interplanetary transfer" block to your MDR with: the departure and arrival hyperbolic excess velocities $v_\infty$; the departure $C_3 = v_\infty^2$; the transfer time; the synodic period (your launch cadence); and the arrival delta-v for orbit insertion (or a note that you will aerocapture). For a Track-C Mars orbiter, that block reads $v_{\infty,\oplus} = 2.95\ \text{km/s}$, $C_3 = 8.7\ \text{km}^2/\text{s}^2$, cruise $\approx 259$ days, window every $\approx 780$ days, $v_{\infty,\text{Mars}} = 2.65\ \text{km/s}$, MOI $\approx 0.7$–$2.1\ \text{km/s}$ depending on the capture orbit. Add the TMI ($3.6\ \text{km/s}$) and MOI to the delta-v budget you started in Chapter 3. (Track D: rerun the transfer with your asteroid's semi-major axis in place of Mars's $1.524\ \text{AU}$.)

The code. Create astrotools/interplanetary.py. It extends the geocentric maneuvers.py from Chapter 10 to heliocentric transfers and adds the timing and $C_3$ tools of this chapter. Keep these signatures stable — the capstone's mission.py will call them.

import math

MU_SUN = 1.327e11    # km^3/s^2, Sun's gravitational parameter
AU = 1.496e8         # km per astronomical unit

def hohmann_transfer(mu, r1, r2):
    """Heliocentric Hohmann r1 -> r2 (km). Returns (v_inf_depart, v_inf_arrive,
    tof_seconds): hyperbolic-excess speeds (km/s) at each planet and transfer time (s)."""
    v1 = math.sqrt(mu / r1)                  # departure-planet circular speed
    v2 = math.sqrt(mu / r2)                  # arrival-planet circular speed
    a  = (r1 + r2) / 2                       # transfer-ellipse semi-major axis
    vp = math.sqrt(mu * (2 / r1 - 1 / a))    # perihelion (departure) speed
    va = math.sqrt(mu * (2 / r2 - 1 / a))    # aphelion (arrival) speed
    return abs(vp - v1), abs(v2 - va), math.pi * math.sqrt(a**3 / mu)

def synodic_period(T1, T2):
    """Time between launch windows for two orbits of periods T1, T2 (same units)."""
    return 1.0 / abs(1.0 / T1 - 1.0 / T2)

def c3_required(v_inf):
    """Characteristic energy C3 (km^2/s^2) from hyperbolic excess speed v_inf (km/s)."""
    return v_inf ** 2

if __name__ == "__main__":
    v_dep, v_arr, tof = hohmann_transfer(MU_SUN, 1.000 * AU, 1.524 * AU)
    print("v_inf depart (Earth):", round(v_dep, 2), "km/s")
    print("v_inf arrive (Mars): ", round(v_arr, 2), "km/s")
    print("departure C3:        ", round(c3_required(v_dep), 2), "km^2/s^2")
    print("transfer time:       ", round(tof / 86400), "days")
    print("synodic period:      ", round(synodic_period(365.25, 686.98)), "days")
    # Expected output:
    # v_inf depart (Earth): 2.95 km/s
    # v_inf arrive (Mars):  2.65 km/s
    # departure C3:         8.68 km^2/s^2
    # transfer time:        259 days
    # synodic period:       780 days

Every number this module prints, we derived by hand in this chapter. By the capstone, you will point hohmann_transfer at your destination, feed its $v_\infty$ into the rocket equation to size the injection stage, and feed its transfer time into the mission timeline. The road to another planet is now a function call.


Summary

Interplanetary flight is orbital mechanics with the Sun at the center, stitched together by patched conics. Carry these forward:

Idea The essential fact
Patched conics Model the trajectory as two-body conics — Earth hyperbola, Sun ellipse, Mars hyperbola — joined at the spheres of influence. Design tool, not flight tool.
Sphere of influence $r_{\text{SOI}} \approx a_{\text{pl}}(m_{\text{pl}}/m_\odot)^{2/5}$. Earth: $\approx 924{,}000\ \text{km}$ ($0.6\%$ of an AU) — tiny, which is why patching works.
Interplanetary Hohmann Heliocentric transfer ellipse tangent to both planets' orbits. Earth→Mars: $v_{\infty,\oplus}=2.95$, $v_{\infty,\text{Mars}}=2.65\ \text{km/s}$, $\Delta v_{\text{helio}}\approx5.6\ \text{km/s}$, cruise $\approx 259$ days.
Launch window / synodic period Depart only when the planets are phased ($\approx44^\circ$ lead for Mars). $1/T_{\text{syn}}=|1/T_1-1/T_2|$; Earth-Mars $\approx 780$ days ($26$ months).
$C_3$ and $v_\infty$ $v_\infty$ = speed left at the SOI edge; $C_3 = v_\infty^2 = 2\varepsilon = -\mu/a$. Mars departure $C_3\approx 8.7\ \text{km}^2/\text{s}^2$. What a launcher is rated against.
Hyperbolic departure TMI from a $300\ \text{km}$ parking orbit: $v_{\text{peri}}=\sqrt{v_\infty^2+v_{\text{esc}}^2}=11.3\ \text{km/s}$; burn $=11.3-7.7=3.6\ \text{km/s}$. Burn low, burn fast (Oberth).
Arrival Brake the arrival hyperbola onto a bound orbit: MOI $\approx 2.1\ \text{km/s}$ to low circular, $\approx 0.7\ \text{km/s}$ to a loose ellipse; or aerobrake / aerocapture.
Gravity assist Flyby preserves $v_\infty$ in the planet's frame but rotates it; adding to the planet's motion changes heliocentric speed. Free to the probe, paid (imperceptibly) by the planet.

Key numbers worth memorizing: Earth heliocentric speed $\approx 29.8\ \text{km/s}$; Mars $\approx 24.1\ \text{km/s}$; Earth→Mars Hohmann cruise $\approx 259$ days; synodic period $\approx 26$ months; Mars departure $C_3 \approx 8.7\ \text{km}^2/\text{s}^2$; TMI from LEO $\approx 3.6\ \text{km/s}$.


Spaced Review

Retrieval strengthens memory. Answer from memory before checking, then look back at the cited section. This chapter revisits Chapter 6 and Chapter 10.

  1. (§11.4, Ch. 6) Chapter 6 classified trajectories by the sign of their specific orbital energy $\varepsilon$. A departure trajectory has $C_3 = 8.7\ \text{km}^2/\text{s}^2$. What is the sign of $\varepsilon$, what shape is the trajectory, and how are $C_3$ and $\varepsilon$ related?
  2. (§11.2, Ch. 6) Using vis-viva with $\mu_\odot$, why is the spacecraft moving slower than Mars when it reaches Mars's orbit, even though it started out moving faster than Earth?
  3. (§11.2, Ch. 10) The interplanetary Hohmann is the same maneuver you derived in Chapter 10 for LEO→GEO. Name the one thing that changes and two things that stay the same.
  4. (§11.4, Ch. 10) In Chapter 10 a Hohmann needed two burns. Where are the "two burns" of the interplanetary Hohmann, and why does each now happen deep inside a planet's gravity well rather than out in the Sun's?

Answers

  1. $\varepsilon > 0$ (positive), so the trajectory is a hyperbola — it escapes with speed to spare. They are related by $C_3 = 2\varepsilon = v_\infty^2$; here $\varepsilon = 4.35\ \text{km}^2/\text{s}^2$ and $v_\infty = \sqrt{8.7} = 2.95\ \text{km/s}$. 2. Because moving out from 1 AU to 1.524 AU means climbing the Sun's gravity well: the spacecraft trades kinetic energy for potential energy (vis-viva with fixed transfer-ellipse energy), so it slows from $32.7\ \text{km/s}$ at perihelion to $21.5\ \text{km/s}$ at aphelion — below Mars's $24.1$. Higher is slower, solar-system-scale. 3. Changes: the central body — the Sun instead of Earth (so $\mu_\odot$ replaces $\mu_\oplus$, and the "orbits" are the planets'). Stay the same: the tangential-burn Hohmann geometry (ellipse touching both circles), and the tools — vis-viva for the speeds, Kepler's third law for the time. 4. The two burns are the trans-Mars injection (departure, §11.4) and the Mars orbit insertion (arrival, §11.5). Each happens at periapsis deep in a planet's well because that is where the spacecraft is moving fastest, and the Oberth effect makes a burn there buy the most energy per unit delta-v — burn low, burn fast.

What's Next

We have built the ideal trajectory to Mars — clean two-body conics, patched at sharp spheres of influence, on circular coplanar planetary orbits. It is a magnificent approximation, and it is wrong in every detail. Real planetary orbits are elliptical and tilted relative to one another; Earth is not a point mass but an oblate, lumpy body whose bulge tugs orbits around; the atmosphere a spacecraft skims never fully lets go; sunlight itself pushes on a spacecraft hard enough, over months, to matter. Every one of these is a perturbation — a small, persistent departure from the perfect Kepler orbit — and small though each is, over an interplanetary distance or a multi-year mission they accumulate into the difference between hitting Mars and missing it. In Chapter 12 we let go of the fiction of the perfect two-body orbit and ask what really happens when gravity is not quite an inverse square and the spacecraft is not quite alone — the physics that turns a designed trajectory into a flown one.