Orbital mechanics
How spacecraft move, from Kepler's laws to Lambert targeting, gravity assists and a full mission. Short sections, worked examples and graded practice, with every number checked by a script.
19 chapters · 34 worked examples · 57 practice problems · 2257 script checks
How to use this pack
- 1Read 'The idea in one minute' and the 'You will be able to' list.
- 2Work through the short sections and the worked examples with a pencil.
- 3Try each practice problem before opening the hint; open the full solution only to check.
- 4Do the three-question check and the spaced recap before moving on.
Chapters
- 1Two-body problem and Kepler's laws30 min
- 2Conic sections and orbital elements35 min
- 3Energy, vis-viva and angular momentum35 min
- 4Time of flight and Kepler's equation40 min
- 5From state vectors to elements and back40 min
- 6Perturbations: J₂, drag, third body and pressure from sunlight40 min
- 7Designer orbits: sun-synchronous, frozen and Molniya35 min
- 8Ground tracks30 min
- 9Coplanar manoeuvres: Hohmann, bi-elliptic and phasing40 min
- 10Plane changes and combined manoeuvres35 min
- 11Rendezvous: Clohessy-Wiltshire40 min
- 12Lambert's problem45 min
- 13Porkchop plots35 min
- 14Patched conics, sphere of influence and the Oberth effect40 min
- 15Gravity assists35 min
- 16Lagrange points35 min
- 17GEO and station-keeping35 min
- 18Launch windows and launch azimuth30 min
- 19Mission case end to end45 min
Your ticks are saved in this browser only. Original work, and a study aid, not assignment answers.
Two-body problem and Kepler's laws
Two bodies pull on each other along the line between them. Subtract one equation of motion from the other and you get a single equation for the satellite relative to Earth, with one constant, μ. Its solutions are conic sections. Because the force is central, angular momentum is conserved and the satellite sweeps equal areas in equal times. Because the force falls as 1/r², the period depends only on the size of the orbit.
You will be able to
- Reduce two gravitating bodies to one relative equation of motion
- State Kepler's three laws and say which conservation law each one comes from
- Compute circular speed and period, and explain why Earth-fixed geometry uses the sidereal day
From two bodies to one
Newton's law gives each body an acceleration towards the other. Their difference is the relative acceleration r̈ = −μ r/r³ with μ = G(M + m). For a satellite m is negligible, so μ is simply Earth's GM = 398600.4418 km³/s². Measured μ is known far better than G or M separately, which is why we always work with μ.[R1]
Kepler's three laws
First: the orbit is an ellipse with Earth's centre at one focus. Second: the line from Earth to the satellite sweeps equal areas in equal times, which is conservation of angular momentum. Third: the square of the period is proportional to the cube of the semi-major axis, T = 2π√(a³/μ).
Circular orbits
For a circle, gravity supplies exactly the centripetal acceleration: v²/r = μ/r², so v = √(μ/r). Higher orbits are slower and take longer. Near the surface a lap takes about an hour and a half; at 42,164 km it takes one sidereal day.
Which day?
Earth turns once relative to the stars in a sidereal day of 86164.09 s. The solar day is longer because Earth has moved along its own orbit and must turn a little more to face the Sun again. Orbit planes are fixed relative to the stars, so Earth-fixed geometry uses the sidereal day.[R2, R5]
Which clock?
Orbit work uses three clocks. UT1 follows Earth's actual rotation, so it wanders slightly. TAI is the steady atomic count. UTC, the clock on every ephemeris file, ticks with TAI but is held within 0.9 s of UT1 by inserting leap seconds; it now runs 37 s behind TAI. Use UTC to stamp an epoch and UT1 (UTC plus the published DUT1 correction) whenever Earth's orientation matters, such as a ground track.[R26]
Key equations
where
where
where
where
for a circle a = r
where
Worked examples
How long is a lap of the Earth?[R1]
Given
A space station flies a circular orbit 415 km above the equatorial radius RE = 6378.137 km. Find its speed and period. Then check Kepler's third law on an eccentric orbit (a = 9500 km, e = 0.25) by integrating the equations of motion for one predicted period.
new symbols
Working
- r = RE + h = 6378.137 + 4156793.137 km
- v = √(μ/r) = √(398600.4418 / 6793.137)7.6601 km/s
- T = 2π √(r³/μ) = 5572.1 s92.87 min
- Eccentric orbit: T = 2π √(9500³ / 398600.4418) = 9215.0 s2.560 h
- Numerical integration (DOP853) for 9215.0 s ends 0.00019 m from the starting point
Answer. v = 7.660 km/s and T = 92.9 min (15.51 laps per day). The eccentric orbit closes after exactly 2.560 h.
Sanity check. The period depends only on a, not on e: the integrated eccentric orbit came back to within 0.00019 m of where it started, after the time Kepler's third law predicts.
r = 6793.137 km, v = 7.6601 km/s, T = 5572.1 s = 92.87 min laps per day = 15.51 eccentric orbit: T = 9215.0 s = 2.560 h after one Kepler period the integrated orbit is 0.0000 m from its start
Why the solar day is four minutes longer[R2, R5]
Given
Earth spins once relative to the stars in one sidereal day, 86164.0905 s. It also moves round the Sun once in 365.2422 days. How long until the Sun is overhead again?
Working
- Spin rate relative to the stars: 1/86164.0905 rev/s; the Sun's apparent drift: 1/31556926 rev/s in the same direction
- Rate relative to the Sun = 1/86164.0905 − 1/31556926 rev/s
- Solar day = 1 / (that rate)86400.0 s
- Difference = 235.9 s3.93 min
Answer. Solar day = 86400 s, 3.93 min longer than the sidereal day.
Sanity check. Orbits are fixed relative to the stars, so anything tied to Earth's rotation (ground tracks, GEO) uses the sidereal day.
solar day = 86400.0 s difference = 235.9 s = 3.93 min
Practice
An Earth-observation orbit
Given
A satellite circles at 550 km altitude. Find its speed, its period, and how many orbits it completes per day.
Add RE first: every formula uses distance from Earth's centre, not altitude.
Answer. v = 7.585 km/s, T = 95.6 min, 15.05 orbits per day.
Working
- r = 6378.137 + 5506928.137 km
- v = √(398600.4418 / 6928.137)7.5851 km/s
- T = 2π √(6928.137³ / 398600.4418) = 5739.0 s95.65 min
- Orbits per day = 86400 / 5739.015.05
Answer. v = 7.585 km/s, T = 95.6 min, 15.05 orbits per day.
Sanity check. Slightly slower and longer than the station at 415 km (7.660 km/s, 92.9 min): higher orbits are slower.
r = 6928.137 km, v = 7.5851 km/s T = 5739.0 s = 95.65 min orbits per day = 15.055
Choose the altitude from the period[R2]
Given
A mission wants a circular orbit with a period of exactly 100 min. What altitude is that? What altitude gives exactly 14 orbits per sidereal day?
Invert Kepler's third law: a = (μ (T/2π)²)^(1/3).
Answer. h = 758 km for a 100-min period; h = 881 km for 14 orbits per sidereal day.
Working
- a = (398600.4418 × (6000/2π)²)^(1/3)7136.64 km
- h = a − RE = 7136.64 − 6378.137758.5 km
- 14 per sidereal day: T = 86164.0905/146154.58 s
- a = (398600.4418 × (6154.58/2π)²)^(1/3) = 7258.69 km, h880.6 km
Answer. h = 758 km for a 100-min period; h = 881 km for 14 orbits per sidereal day.
Sanity check. Period grows as a^(3/2): a 4 % longer period needs only about 2.7 % more radius, so altitudes are sensitive to the period you ask for.
100-min period: a = 7136.64 km, h = 758.5 km 14 per sidereal day: T = 6154.58 s, a = 7258.69 km, h = 880.6 km
Kepler's third law from the second
Given
An orbit has perigee altitude 1000 km and apogee altitude 40000 km. The perigee speed is √(μ(2/rp − 1/a)). Use the area law (areal rate h/2 is constant) and the ellipse area πab to find the period, and show it equals 2π√(a³/μ).
h = rp vp at perigee (velocity is perpendicular to the radius there). b = a√(1 − e²).
Answer. T = 43854 s = 12.182 h by either route.
Working
- rp = 7378.137 km, ra = 46378.137 km; a = 26878.137 km, e0.72550
- b = a√(1 − e²) = 18498.223 km; vp9.65501 km/s
- h = rp vp = 71235.97 km²/s; areal rate h/235617.98 km²/s
- T = πab / (h/2)43854.05 s
- Kepler III: 2π√(a³/μ)43854.05 s(equal)
Answer. T = 43854 s = 12.182 h by either route.
Sanity check. Using h = √(μ a (1 − e²)) the two expressions are algebraically identical; matching to many digits confirms the arithmetic.
a = 26878.137 km, e = 0.72550, b = 18498.223 km, vp = 9.65501 km/s, h = 71235.97 km^2/s area law: T = 43854.05 s = 12.1817 h Kepler III: T = 43854.05 s
Common mistakes
- Using altitude instead of radius in √(μ/r). Always add RE = 6378.137 km.
- Mixing metres and kilometres: μ = 398600.4418 is in km³/s², so r must be in km and v comes out in km/s.
- Assuming a higher orbit is faster. It is slower, even though it took more energy to reach.
- Using 86400 s for anything tied to Earth's rotation relative to the stars; use the sidereal day 86164.09 s.
Check yourself
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R5Williams, D.R., "Earth Fact Sheet", NASA NSSDCA, Goddard Space Flight Center, updated 15 Nov 2024.
- R26NIST Time and Frequency Division, "Leap second and UT1-UTC information" (created 6 Aug 2009, updated 28 Oct 2026).
Where in each source, and links: References.
Conic sections and orbital elements
Every two-body orbit is a conic: circle, ellipse, parabola or hyperbola. One equation describes them all, r = p/(1 + e cos ν), and the eccentricity e picks the shape. To place the orbit in space you need three more angles, and to place the satellite on it, one more. Those six numbers are the classical orbital elements.
You will be able to
- Use the orbit equation r = p/(1 + e cos ν) on any conic
- Classify an orbit by its eccentricity
- Name the six classical elements and move a position from the perifocal to the ECI frame
The orbit equation
Measured from Earth's centre, the radius depends only on the true anomaly ν, the angle from perigee. p, the semi-latus rectum, is the radius at ν = 90°. For an ellipse p = a(1 − e²), perigee is a(1 − e) and apogee a(1 + e).
The conic family
e = 0 is a circle and 0 < e < 1 an ellipse: both closed. e = 1 is a parabola, the boundary case that just escapes. e > 1 is a hyperbola: the satellite leaves with speed to spare and its path bends by a finite angle.
Six numbers place an orbit
a and e give size and shape. Inclination i tilts the plane from the equator. Right ascension of the ascending node Ω turns the line where the plane crosses the equator northbound. Argument of perigee ω points perigee within the plane. True anomaly ν places the satellite.
Two frames
The perifocal frame lies in the orbit plane with x towards perigee, so position is simply (r cos ν, r sin ν, 0). The Earth-centred inertial frame has z along Earth's spin axis and x towards the vernal equinox, fixed relative to the stars. Three rotations (ω, i, Ω) carry one into the other.
Where it sits in the sky
Telescopes and ground stations use two more frames. The horizon frame gives altitude above the horizon and azimuth clockwise from north. The equatorial frame projects Earth's equator onto the sky: declination is the angle north or south of it, right ascension the angle east from the vernal equinox, counted in hours (1 h = 15°). The link between them is sidereal time. Local sidereal time is the right ascension now on your meridian, and a star's hour angle is HA = LST − RA. Greenwich mean sidereal time at J2000.0 (2000 Jan 1, 12:00 UT1) is 18.6974 h (18 h 41.8 min); add your east longitude divided by 15 to get LST. At London (0.1278° W) that is 18.6889 h, so a star at RA 6 h then has HA = 12.6889 h.[R27]
Frames that drift
The equinox is not fixed. The Sun and Moon pull on the equatorial bulge and the spin axis precesses: the general precession runs at about 50.3 arcsec per year, a full turn in about 25,770 years. A smaller wobble, nutation, rides on top, its largest term following the Moon's 18.6-year node cycle. So every catalogue of right ascensions and every inertial frame names an epoch; the standard one is J2000.0.[R28]
| Element | Symbol | Sets | Range |
|---|---|---|---|
| Semi-major axis | a | size (and period) | a > 0 closed, a < 0 hyperbola |
| Eccentricity | e | shape | 0 ≤ e |
| Inclination | i | tilt of the plane | 0° to 180° |
| Right ascension of ascending node | Ω | where the plane crosses the equator | 0° to 360° |
| Argument of perigee | ω | where perigee sits in the plane | 0° to 360° |
| True anomaly | ν | where the satellite is now | 0° to 360° |
i above 90° means retrograde: the satellite moves against Earth's rotation.
Key equations
where
where
where
where
where
only for e > 1
Worked examples
Shape of a transfer ellipse
Given
A transfer orbit has perigee altitude 250 km and apogee at geostationary radius, 42164.2 km. Find a, e, p and the radius a quarter-turn after perigee (ν = 90°).
new symbols
Working
- rp = 6378.137 + 250 = 6628.137 km; ra42164.170 km
- a = (rp + ra)/224396.153 km
- e = (ra − rp)/(ra + rp)0.72831
- p = a(1 − e²)11455.490 km
- r(90°) = p / (1 + e cos 90°) = p11455.490 km
Answer. a = 24396 km, e = 0.7283, p = 11455 km, r(90°) = 11455 km (altitude 5077 km).
Sanity check. p sits between rp and ra, much nearer rp: an eccentric orbit stays low for only a short arc around perigee.
rp = 6628.137 km, ra = 42164.170 km a = 24396.153 km, e = 0.72831, p = 11455.490 km r(90 deg) = 11455.490 km, altitude 5077.4 km
Where is it, in Earth-centred coordinates?[R1]
Given
A satellite has a = 7200 km, e = 0.01, i = 51.6°, Ω = 120°, ω = 80°, ν = 30°. Find its position in the Earth-centred inertial (ECI) frame.
Working
- p = a(1 − e²) = 7199.280 km; r = p/(1 + e cos ν)7137.468 km
- In the perifocal frame (x to perigee): r = (6181.23, 3568.73, 0) km
- Rotate by ω about z, by i about the new x (the node line), then by Ω about z
- ECI: r = (-2387.33, -4197.13, 5256.25) km
Answer. r = (-2387.3, -4197.1, 5256.3) km, |r| = 7137.5 km.
Sanity check. z = r sin(ω + ν) sin i = 5256.25 km matches the third component: the height above the equator depends only on the argument of latitude and i.
p = 7199.280 km, r = 7137.468 km perifocal r = (6181.23, 3568.73, 0) km ECI r = (-2387.33, -4197.13, 5256.25) km, |r| = 7137.47 km
Practice
From altitudes to a and e
Given
An orbit has perigee altitude 500 km and apogee altitude 1500 km. Find a and e.
Convert altitudes to radii first, then a = (rp + ra)/2 and e = (ra − rp)/(ra + rp).
Answer. a = 7378.1 km, e = 0.0678 (dimensionless).
Working
- rp = 6878.137 km, ra7878.137 km
- a = (6878.137 + 7878.137)/27378.137 km
- e = (7878.137 − 6878.137) / (7878.137 + 6878.137)0.06777
Answer. a = 7378.1 km, e = 0.0678 (dimensionless).
Sanity check. A 1000 km altitude difference is a small eccentricity because it is compared with radii of about 7,000 km.
rp = 6878.137 km, ra = 7878.137 km a = 7378.137 km, e = 0.06777
A highly elliptical orbit
Given
An orbit has a = 26600 km and e = 0.72. Find the radius at ν = 135°, and the true anomalies at which r = 20000 km.
Rearrange the orbit equation: cos ν = (p/r − 1)/e. There are two answers, ±ν.
Answer. r(135°) = 26097 km; r = 20000 km at ν = ±120.0°.
Working
- p = a(1 − e²) = 26600 × (1 − 0.72²)12810.56 km
- r(135°) = 12810.56 / (1 + 0.72 cos 135°)26097.0 km
- cos ν = (12810.56/20000 − 1)/0.72-0.49927
- ν = ±119.95°
Answer. r(135°) = 26097 km; r = 20000 km at ν = ±120.0°.
Sanity check. rp = 7448 km and ra = 45752 km bracket both radii, as they must.
p = 12810.56 km, r(135 deg) = 26097.0 km cos(nu) = -0.49927, so r = 20000 km at nu = 119.951 deg and at minus that
An escaping hyperbola
Given
A probe leaves Earth on a hyperbola with perigee altitude 300 km and e = 1.3. Find p, the true anomaly of the asymptote ν∞, the angle the path turns through, and the radius at ν = 100°.
The radius goes to infinity where 1 + e cos ν = 0. The turn angle is 2ν∞ − 180°.
Answer. p = 15360 km, ν∞ = 140.3°, turn 100.6°, r(100°) = 19838 km.
Working
- rp = 6678.137 km; p = rp(1 + e)15359.72 km
- ν∞ = acos(−1/e) = acos(-0.76923)140.28°
- Turn angle δ = 2ν∞ − 180° = 100.57° (equivalently 2 asin(1/e))
- r(100°) = 15359.72 / (1 + 1.3 cos 100°)19838.0 km
Answer. p = 15360 km, ν∞ = 140.3°, turn 100.6°, r(100°) = 19838 km.
Sanity check. Only true anomalies inside ±ν∞ exist on a hyperbola; 100° is inside, and its radius is large but finite.
p = 15359.72 km nu_inf = 140.28 deg turn angle = 100.57 deg (check 2 asin(1/e) = 100.57 deg) r(100 deg) = 19838.0 km
Common mistakes
- Using altitudes in e = (ra − rp)/(ra + rp). Radii are needed; altitudes give the wrong e.
- Forgetting the second solution when inverting the orbit equation: r is the same at +ν and −ν.
- Measuring ν from the ascending node. ν is measured from perigee; the angle from the node is the argument of latitude u = ω + ν.
- Thinking i = 100° is a 100° tilt the 'normal' way. Above 90° the orbit is retrograde.
Check yourself
Recap from earlier chapters
What is the period of a circular orbit at 800 km altitude?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R27U.S. Naval Observatory, Astronomical Applications Dept., "Computing Approximate Sidereal Time".
- R28Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS TN 36, Ch. 5 (10 Aug 2012 update).
Where in each source, and links: References.
Energy, vis-viva and angular momentum
Two numbers never change on a two-body orbit: the energy per kilogram and the angular momentum per kilogram. Energy fixes the size of the orbit, angular momentum fixes its shape. From one position and velocity you can therefore find the whole orbit, and the speed at any other radius follows from one line: the vis-viva equation.
You will be able to
- Compute specific energy and use vis-viva to get speed anywhere on an orbit
- Find escape speed and the hyperbolic excess speed
- Use h = r v cos γ for the flight-path angle, and check a propagator by its conserved quantities
Specific energy
Kinetic plus potential energy per unit mass, ε = v²/2 − μ/r, is constant. It depends only on the semi-major axis: ε = −μ/(2a). Negative means bound (ellipse), zero is a parabola, positive is a hyperbola.
Vis-viva
Put the two energy forms together and you get v² = μ(2/r − 1/a). It gives the speed at any radius once you know a. It is the single most useful line in this pack: every manoeuvre chapter uses it.
Angular momentum and flight-path angle
h = |r × v| = r v cos γ, where γ is the angle of the velocity above the local horizontal. γ is zero at perigee and apogee, positive while climbing. h fixes the parameter p = h²/μ, and with a this gives e.
Escape
Setting ε = 0 gives the escape speed √(2μ/r), √2 times the circular speed. With more speed than that, the probe arrives at infinity with a leftover speed v∞, the hyperbolic excess speed: ε = v∞²/2.
Key equations
where
where
where
where
where
where
Worked examples
Speed and climb angle partway up
Given
A satellite is on an orbit with perigee altitude 500 km and apogee altitude 2000 km. As it climbs through 1200 km altitude, find its energy, speed and flight-path angle.
new symbols
Working
- rp = 6878.137 km, ra = 8378.137 km, a = 7628.137 km, e0.09832
- ε = −μ/(2a) = −398600.4418/(2 × 7628.137)-26.1270 km²/s²
- Vis-viva at r = 7578.137 km: v = √(μ(2/r − 1/a))7.2762 km/s
- h = √(μp) = √(398600.4418 × 7554.397)54874.27 km²/s
- cos γ = h/(r v) = 0.995176, so γ5.630°(climbing)
Answer. ε = -26.127 km²/s², v = 7.276 km/s, γ = 5.63°.
Sanity check. Check with the anomaly form: ν = 91.83° and tan γ = e sin ν/(1 + e cos ν) gives the same γ. A small e gives a small climb angle.
a = 7628.137 km, e = 0.09832 energy = -26.1270 km^2/s^2 v = 7.2762 km/s, h = 54874.27 km^2/s, gamma = 5.630 deg
Watch a computer conserve energy
Given
Integrate the equation r̈ = −μr/r³ numerically for 5 revolutions of an orbit with a = 12000 km, e = 0.5, starting at perigee. Track energy and angular momentum.
Working
- Start: rp = 6000 km, vp = 9.9825 km/s; period 13082.3 s
- Integrate with an 8th-order Runge-Kutta scheme (DOP853, tolerance 10⁻¹²) for 65411 s
- Energy ε = v²/2 − μ/r at every output: start -16.608352 km²/s², largest relative change 2.6e-11
- |h| = |r × v|: start 59894.941 km²/s, largest relative change 2.4e-11
- Speed swings between 3.3275 and 9.9825 km/s
Answer. ε = -16.6084 km²/s² and h = 59894.9 km²/s stay constant to better than 1 part in 10⁹, while the speed changes by a factor 3.00.
Sanity check. Both drifts (3e-11 and 2e-11) are below 10⁻⁹. Checking conserved quantities is the standard way to trust a propagator.
rp = 6000 km, vp = 9.9825 km/s, period 13082.3 s energy = -16.608352 km^2/s^2, largest relative change 1.1e-12 h = 59894.941 km^2/s, largest relative change 1.9e-14 both constant to better than 1 part in 10^9: True speed from 3.3275 to 9.9825 km/s, a factor 3.00
Practice
Escape from an 800 km orbit[R1]
Given
A spacecraft is on a circular orbit at 800 km altitude. Find its speed, the local escape speed and the Δv needed to escape with a single prograde burn.
Escape speed makes ε = 0: vesc = √(2μ/r).
Answer. vc = 7.452 km/s, vesc = 10.538 km/s, Δv = 3.087 km/s.
Working
- r = 7178.137 km; vc = √(μ/r)7.4518 km/s
- vesc = √(2μ/r)10.5385 km/s
- Δv = 10.5385 − 7.45183.0866 km/s
Answer. vc = 7.452 km/s, vesc = 10.538 km/s, Δv = 3.087 km/s.
Sanity check. vesc/vc = 1.4142 = √2 at every radius.
v_c = 7.4518 km/s, v_esc = 10.5385 km/s escape burn dv = 3.0866 km/s
Orbit from one radar fix
Given
Tracking finds a satellite at r = 8000 km moving at 7.6 km/s, 5° above the local horizontal. Find a, e, perigee and apogee altitudes.
Energy gives a. Angular momentum h = r v cos γ gives p. Then e = √(1 − p/a).
Answer. a = 9515 km, e = 0.1810, perigee altitude 1415 km, apogee altitude 4860 km.
Working
- ε = v²/2 − μ/r = 7.6²/2 − 398600.4418/8000-20.9451 km²/s²
- a = −μ/(2ε)9515.38 km
- h = r v cos γ = 8000 × 7.6 × cos 5° = 60568.64 km²/s; p = h²/μ9203.60 km
- e = √(1 − p/a)0.18101
- rp = a(1 − e) = 7793.0 km (altitude 1415 km); ra11237.8 km(altitude 4860 km)
Answer. a = 9515 km, e = 0.1810, perigee altitude 1415 km, apogee altitude 4860 km.
Sanity check. r = 8000 km lies between rp and ra, and the climb angle is positive, so the satellite is on its way to apogee.
energy = -20.9451 km^2/s^2, a = 9515.38 km h = 60568.64 km^2/s, p = 9203.60 km, e = 0.18101 perigee altitude 1414.8 km, apogee altitude 4859.7 km
A probe on its way out
Given
After its departure burn a probe is at perigee, 400 km altitude, moving at 11.6 km/s. Find its energy, its hyperbolic excess speed v∞, its eccentricity, and its speed when it is 1,000,000 km from Earth's centre.
ε = v²/2 − μ/r is conserved. At infinite distance ε = v∞²/2. At perigee e = 1 + rp v∞²/μ.
Answer. ε = 8.47 km²/s², v∞ = 4.117 km/s, e = 1.2882, v = 4.212 km/s at a million km.
Working
- ε = 11.6²/2 − 398600.4418/6778.1378.4732 km²/s²(positive: escaping)
- v∞ = √(2ε)4.1166 km/s
- e = 1 + rp v∞²/μ = 1 + 6778.137 × 4.1166² / 398600.44181.28817
- At r = 1e+06 km: v = √(2(ε + μ/r))4.2123 km/s
Answer. ε = 8.47 km²/s², v∞ = 4.117 km/s, e = 1.2882, v = 4.212 km/s at a million km.
Sanity check. At a million km the speed is already within 2 % of v∞: almost all the slowing happens close to Earth.
energy = 8.4732 km^2/s^2 (positive: escaping) v_inf = 4.1166 km/s, e = 1.28817 speed at 1000000 km: 4.2123 km/s
Common mistakes
- Using v cos γ when you need the radial speed: v sin γ is radial, v cos γ is horizontal.
- Forgetting that a is negative for a hyperbola. Then −μ/(2a) is positive, as it should be.
- Using vis-viva with an altitude in place of r.
- Treating energy drift in a numerical propagation as physics. In pure two-body motion any drift is integration error.
Check yourself
Recap from earlier chapters
Kepler's second law is a statement of which conservation law?
What is p for an ellipse with a = 10000 km and e = 0.2?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
Where in each source, and links: References.
Time of flight and Kepler's equation
Position as a function of angle is easy: the orbit equation. Position as a function of time is harder, because the satellite speeds up and slows down. Kepler's trick is an auxiliary angle, the eccentric anomaly E, that turns time into a simple equation, M = E − e sin E. Going from angle to time is direct; going from time to angle needs a few steps of Newton's method.
You will be able to
- Move between true, eccentric and mean anomaly
- Solve Kepler's equation with Newton's method and know when it has converged
- Find the time between two points, and the position at a given time, on an ellipse or a hyperbola
Three anomalies
True anomaly ν is the real angle from perigee. Mean anomaly M = n t grows uniformly with time, as if the orbit were a circle. Eccentric anomaly E is the geometric bridge between them, measured on the circle that just encloses the ellipse.
Kepler's equation
M = E − e sin E. Given ν, compute E then M then t: no iteration. Given t, you must invert it. Newton's method converges in a handful of steps from the guess E₀ = M + e sin M, and from E₀ = π for very eccentric orbits.
The hyperbolic version
For e > 1 the same idea uses a hyperbolic angle F: M = e sinh F − F, with M = √(μ/(−a)³) t. A Newton solve works the same way.
Checking your answer
Two checks catch nearly every error. Put the answer back into Kepler's equation and see the residual vanish. Or integrate the motion numerically for the same time and compare positions.
| Step | E (rad) | Residual E − e sin E − M (rad) |
|---|---|---|
| 0 | 1.1480830276 | -1.88e-01 |
| 1 | 1.4124278609 | 2.31e-02 |
| 2 | 1.3865107052 | 2.32e-04 |
| 3 | 1.3862447747 | 2.43e-08 |
| 4 | 1.3862447468 | 3.33e-16 |
Newton iterations for e = 0.7, M = 40°. The residual roughly squares at every step: quadratic convergence.
Key equations
where
where
where
where
where
where
a < 0 for a hyperbola
Worked examples
Solve Kepler's equation by Newton's method
Given
An orbit has e = 0.7. Find E and ν when the mean anomaly is M = 40°.
new symbols
Working
- M = 40° = 0.698132 rad. Start with E₀ = M + e sin M1.148083 rad
- Update E ← E − (E − e sin E − M)/(1 − e cos E)
- Iteration 1: E = 1.4124278609 rad, residual 2.31e-02 rad
- Iteration 2: E = 1.3865107052 rad, residual 2.32e-04 rad
- Iteration 3: E = 1.3862447747 rad, residual 2.43e-08 rad
- Iteration 4: E = 1.3862447468 rad, residual 3.33e-16 rad
- E = 79.4260°; ν = 2 atan(√((1+e)/(1−e)) tan(E/2))126.3430°
Answer. E = 79.426°, ν = 126.343°, in 5 Newton steps.
Sanity check. ν > E > M, as always between perigee and apogee: the satellite runs ahead of the uniform mean motion near perigee.
step 1: E = 1.4124278609 rad, residual 2.31e-02 rad step 2: E = 1.3865107052 rad, residual 2.32e-04 rad step 3: E = 1.3862447747 rad, residual 2.43e-08 rad step 4: E = 1.3862447468 rad, residual 3.33e-16 rad step 5: E = 1.3862447468 rad, residual 1.11e-16 rad E = 79.4260 deg, nu = 126.3430 deg, in 5 Newton steps
How long to climb to ν = 150°?
Given
A highly elliptical orbit has perigee altitude 600 km and apogee altitude 39000 km. How long after perigee does the satellite reach ν = 150°? Check by integrating the motion.
Working
- a = 26178.137 km, e = 0.73344, n = √(μ/a³) = 1.490599e-04 rad/s, T11.709 h
- E = 2 atan(√((1−e)/(1+e)) tan(ν/2))1.942745 rad
- M = E − e sin E1.259460 rad
- t = M/n = 8449.4 s2.347 h
- Integrating from perigee for 8449.4 s lands 2.7e-05 m from the Kepler position
Answer. t = 2.35 h after perigee, 20.0 % of the period for 42 % of the angle.
Sanity check. The satellite covers the first 150° quickly and spends most of the period near apogee, which is the whole point of such an orbit.
a = 26178.137 km, e = 0.73344, T = 11.709 h E = 1.942745 rad, M = 1.259460 rad t = 8449.4 s = 2.347 h after perigee that is 20.0 % of the period for 42 % of the angle
Practice
A quarter-turn on a transfer orbit
Given
On a transfer orbit from 250 km altitude to GEO radius (42164.2 km), how long after perigee is the satellite at ν = 90°?
ν → E → M → t. Use the half-angle formula for E.
Answer. t = 25.7 min after perigee.
Working
- a = 24396.153 km, e0.72831
- E = 2 atan(√((1−e)/(1+e)) tan 45°)0.754941 rad
- M = E − e sin E0.255868 rad
- t = M √(a³/μ) = 1544.3 s25.7 min
Answer. t = 25.7 min after perigee.
Sanity check. The whole climb to apogee takes 316 min; reaching 90° takes only 8 % of that.
a = 24396.153 km, e = 0.72831, E = 0.754941 rad, M = 0.255868 rad t = 1544.3 s = 25.74 min after perigee
Where will it be in 40 minutes?
Given
A satellite passes perigee of an orbit with a = 9000 km, e = 0.3. Where is it (ν and r) 40 min later?
M = n t, then solve M = E − e sin E for E by Newton, then convert E to ν.
Answer. ν = 131.57°, r = 10226 km (altitude 3847 km).
Working
- n = √(μ/a³) = 7.394437e-04 rad/s; M = n t1.774665 rad
- Newton on M = E − e sin E: E = 2.041975 rad after 4 steps
- ν = 2 atan(√((1+e)/(1−e)) tan(E/2))131.573°
- r = a(1 − e cos E)10225.63 km
Answer. ν = 131.57°, r = 10226 km (altitude 3847 km).
Sanity check. 28 % of the period has passed, and ν = 132° is well over 102°: near perigee the satellite moves faster than average.
n = 7.394437e-04 rad/s, M = 1.774665 rad, E = 2.041975 rad nu = 131.573 deg, r = 10225.63 km (altitude 3847.5 km)
Leaving on a hyperbola
Given
A probe passes perigee at 300 km altitude on a hyperbola with e = 1.8. How long until it is 100,000 km from Earth's centre?
a = rp/(1 − e) is negative. r = a(1 − e cosh F) gives F; then M = e sinh F − F and t = M √(−a³/μ).
Answer. t = 3.42 h after perigee, at ν = 116.9°.
Working
- rp = 6678.137 km; a = rp/(1 − e)-8347.67 km
- cosh F = (1 − r/a)/e = 7.210772, so F2.663880
- M = e sinh F − F10.190090
- t = M √((−a)³/μ) = 12310.0 s3.419 h
- True anomaly there: ν = 2 atan(√((e+1)/(e−1)) tanh(F/2))116.85°(asymptote at 123.75°)
Answer. t = 3.42 h after perigee, at ν = 116.9°.
Sanity check. A universal-variable propagation of the perigee state for the same time lands at the same radius, and ν is below the asymptote angle as it must be.
a = -8347.67 km, F = 2.663880, M = 10.190090 t = 12310.0 s = 3.419 h nu = 116.85 deg (asymptote at 123.75 deg)
Common mistakes
- Working in degrees inside Kepler's equation. E − e sin E only makes sense with E in radians.
- Getting the wrong quadrant from atan: use the half-angle formula (or atan2), which keeps E and ν in the same half of the orbit.
- Stopping Newton after a fixed number of steps without checking the residual.
- Using M = n t with t measured from the wrong event. t is time since perigee passage.
Check yourself
Recap from earlier chapters
What is the vis-viva equation?
What is the perigee radius of an ellipse in terms of a and e?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
Where in each source, and links: References.
From state vectors to elements and back
Tracking gives you a position and a velocity at one instant: six numbers. The classical elements are six other numbers describing the same orbit in a more meaningful way. Going from r, v to elements is a short recipe built on three vectors (h, n, e) and a few arc cosines; going back is a rotation from the orbit's own frame into the Earth-centred one.
You will be able to
- Compute the six classical elements from a position and velocity
- Apply the three quadrant checks that every acos needs
- Build r and v from elements, and handle the circular and equatorial cases
Three helper vectors
h = r × v is perpendicular to the orbit plane and gives i. The node vector n = k × h points to the ascending node and gives Ω. The eccentricity vector e points to perigee, its length is e, and it gives ω and ν.[R11]
Quadrant checks
acos returns 0° to 180° only. Three sign tests fix this: if ny < 0 then Ω = 360° − Ω; if ez < 0 then ω = 360° − ω; if r·v < 0 (falling towards perigee) then ν = 360° − ν.
Back to r and v
In the perifocal frame the state is simple: r = (r cos ν, r sin ν, 0) and v = √(μ/p)(−sin ν, e + cos ν, 0). Rotate by ω, then i, then Ω and you are in ECI.
When elements break
On a circular orbit perigee does not exist, so ω and ν are undefined; use the argument of latitude u = ω + ν. On an equatorial orbit the node does not exist; use the longitude of perigee or true longitude. Numerically, near-circular orbits make ω and ν very sensitive to noise.
Key equations
where
where
where
where
where
where
Worked examples
Six elements from one state vector
Given
A tracking station gives r = (4506.847, 6754.624, 4127.819) km and v = (-1.220029, 3.493813, -5.021162) km/s in the ECI frame. Find the classical elements.
new symbols
Working
- |r| = 9109.089 km, |v| = 6.237569 km/s, r·v-2625.539 km²/s(negative: falling towards perigee)
- h = r × v = (-48337.9, 17593.6, 23986.9) km²/s, |h|56757.88 km²/s
- i = acos(hz/|h|)65.000°
- n = k × h = (-17593.6, -48337.9, 0); Ω = acos(nx/|n|) = 110.000°, but ny < 0 so Ω = 360° − that250.000°
- e = ((v² − μ/r) r − (r·v) v)/μ: |e|0.120000
- ω = acos(n·e/(|n||e|)) = 50.000°, but ez < 0 so ω310.000°
- ν = acos(e·r/(|e||r|)) = 160.000°, but r·v < 0 so ν200.000°
- a = −μ/(2ε) with ε = v²/2 − μ/r = -24.30491 km²/s²: a8200.00 km
Answer. a = 8200.0 km, e = 0.1200, i = 65.00°, Ω = 250.00°, ω = 310.00°, ν = 200.00°.
Sanity check. Every acos returns 0° to 180°; each of the three sign tests moved an angle into the other half. Feeding the elements back through coe2rv returns the same state.
|r| = 9109.089 km, |v| = 6.237569 km/s, r.v = -2625.539, |h| = 56757.88 km^2/s a = 8200.00 km, e = 0.12000, i = 65.000 deg RAAN = 250.000 deg, argp = 310.000 deg, nu = 200.000 deg
Elements to a state vector
Given
An Earth-observation satellite has a = 7000 km, e = 0.0012, i = 97.8°, Ω = 30°, ω = 90°, ν = 45°. Find r and v in ECI.
Working
- p = a(1 − e²) = 6999.990 km; r = p/(1 + e cos ν)6994.055 km
- Perifocal: r = (4945.544, 4945.544, 0) km, v = √(μ/p)(−sin ν, e + cos ν, 0) = (-5.33587, 5.34492, 0) km/s
- Multiply by Q = R₃(−Ω) R₁(−i) R₃(−ω)
- ECI: r = (-3947.373, -3054.037, 4899.787) km
- ECI: v = (-4.99092, -2.04532, -5.28650) km/s
Answer. r = (-3947.4, -3054.0, 4899.8) km, v = (-4.991, -2.045, -5.287) km/s.
Sanity check. A rotation never changes length: |r| and |v| are the same in both frames. With u = ω + ν = 135° the satellite has passed its most northerly point and z = 4900 km is still north of the equator.
p = 6999.990 km, r = 6994.055 km ECI r = (-3947.373, -3054.037, 4899.787) km ECI v = (-4.99092, -2.04532, -5.28650) km/s
Practice
Inclination from a state
Given
A satellite is at r = (6800, 0, 0) km with v = (0, 5.2, 5.6) km/s. Find h and the inclination.
h = r × v; i = acos(hz/|h|).
Answer. |h| = 51966 km²/s, i = 47.12°.
Working
- h = r × v = (0·5.6 − 0·5.2, 0·0 − 6800·5.6, 6800·5.2 − 0) = (0, -38080, 35360) km²/s
- |h|51965.53 km²/s
- i = acos(35360/51965.53)47.121°
Answer. |h| = 51966 km²/s, i = 47.12°.
Sanity check. The satellite is on the x axis, in the equatorial plane, so it is at a node and i is simply the angle of v above the equator: atan(5.6/5.2) = 47.12°.
h = (0, -38080, 35360) km^2/s |h| = 51965.53 km^2/s, i = 47.121 deg
All six elements
Given
Find all six classical elements for r = (-7902.47, -296.38, 5413.27) km, v = (-1.29711, -6.84390, 0.06230) km/s.
Follow the recipe: h, i, n, Ω, e, ω, ν, a. Do the three sign checks (ny, ez, r·v) at the end.
Answer. a = 11500 km, e = 0.250, i = 35.00°, Ω = 80.00°, ω = 40.00°, ν = 60.00°.
Working
- |r| = 9583.338 km, |v| = 6.96601 km/s, r·v12616.01 km²/s
- h = (37029.4, -6529.3, 53699.3); i35.000°
- n = (6529.3, 37029.4, 0); ny ≥ 0 so Ω80.000°
- e vector = (-0.09638, 0.21146, 0.09217); e = 0.25000; ez ≥ 0 so ω40.000°
- r·v > 0 so ν60.000°
- ε = -17.33040 km²/s², a = −μ/(2ε)11500.04 km
Answer. a = 11500 km, e = 0.250, i = 35.00°, Ω = 80.00°, ω = 40.00°, ν = 60.00°.
Sanity check. Check: rp = a(1 − e) must be above Earth's surface; here it is 2247 km altitude. Converting back with coe2rv reproduces the state to rounding.
energy = -17.33040 km^2/s^2 a = 11500.04 km, e = 0.25000, i = 35.000 deg RAAN = 80.000 deg, argp = 40.000 deg, nu = 60.000 deg
When perigee is almost undefined
Given
A near-circular orbit has a = 6930 km, e = 2.0e-04, i = 97.5°, Ω = 300°, ω = 120°, ν = 250°. Convert to r, v and back. Which recovered quantity is reliable, and why? Find it directly from the node vector.
With e tiny, the direction of perigee is set by a vector only 10⁻⁴ long. The angle from the node to the satellite, u = ω + ν, does not need perigee at all.
Answer. u = 10.000° is reliable; ω and ν separately are ill-conditioned when e ≈ 0, though in exact arithmetic the round trip still closes (error about 10⁻¹² km, floating-point level).
Working
- coe2rv: r = (3276.554, -5989.325, 1193.168) km, |r|6930.474 km
- rv2coe gives ω = 120.0000° and ν250.0000°
- u = ω + ν10.000000°(given 10.000000°)
- Directly: u = acos(n·r/(|n||r|)) = 10.000000°, with z ≥ 0 deciding the half
- Round trip back to r, v: position error 1.5e-12 km
Answer. u = 10.000° is reliable; ω and ν separately are ill-conditioned when e ≈ 0, though in exact arithmetic the round trip still closes (error about 10⁻¹² km, floating-point level).
Sanity check. Real tracking noise of a few metres changes the e vector's direction by tens of degrees at e = 10⁻⁴, but barely moves u. That is why near-circular orbits are described with u (or equinoctial elements).
r = (3276.554, -5989.325, 1193.168) km recovered argp = 120.0000 deg, nu = 250.0000 deg, e = 2.0e-04 u = argp + nu = 10.000000 deg; directly from the node vector u = 10.000000 deg round-trip position error 6.4e-13 km
Common mistakes
- Skipping the quadrant checks. Half of all random states give a wrong Ω, ω or ν without them.
- Using the r·v test for ω or the ez test for ν: each angle has its own test.
- Trusting ω and ν on a near-circular orbit. Report u = ω + ν instead.
- Applying the rotations in the wrong order. The perifocal-to-ECI matrix is R₃(−Ω) R₁(−i) R₃(−ω).
Check yourself
Recap from earlier chapters
Write the orbit equation.
How many Newton steps does Kepler's equation usually need at e = 0.7?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R11Vallado, D.A., companion software to Fundamentals of Astrodynamics and Applications (CelesTrak), GitHub repo CelesTrak/fundamentals-of-astrodynamics, commit 49df047 (AGPL-3.0; values read, no code copied).
Where in each source, and links: References.
Perturbations: J₂, drag, third body and pressure from sunlight
A real orbit is a Kepler ellipse that slowly bends. Earth's equatorial bulge (J₂) turns the orbit plane and the ellipse inside it, steadily, by degrees per day. Air drag, thin as it is, takes energy and lowers the orbit. The Moon and the Sun pull a little differently on the satellite than on Earth, and sunlight itself pushes. Each is tiny compared with Earth's central pull, but they add up over days and years, and each one wins at a different height.
You will be able to
- Compute the J₂ drift of the node and the perigee and say which way they turn
- Estimate drag decay per orbit and a rough orbital lifetime
- Rank J₂, drag, the Moon, the Sun and sunlight pressure for any altitude
Earth's bulge: J₂
Earth is wider at the equator by about 21 km. The extra mass pulls an inclined orbit towards the equator, and the orbit responds like a gyroscope: the plane precesses. The node moves west for a prograde orbit and east for a retrograde one. Inside the plane the ellipse rotates too, except at 63.4°, the critical inclination.[R1]
Air drag
Drag removes energy, so the orbit shrinks and, oddly, the satellite speeds up. Density falls roughly exponentially with height, so the decay rate depends hugely on altitude and on the Sun's activity, which heats and swells the upper atmosphere.[R14, R15]
The Moon, the Sun and sunlight
A third body matters through the difference in its pull on the satellite and on Earth: a tidal acceleration that grows with distance from Earth. Sunlight carries momentum and pushes on every lit surface, by an amount set by the area-to-mass ratio, not by altitude.[R2, R8]
What a tracker sees: the period shrinks
Radar measures the period more easily than the height. From T = 2π√(a³/μ), dT/da = 3T/(2a), and with the drag decay per revolution this gives dT/dt = −3π CD (A/m) ρ a (SI units, circular orbit, density held fixed). For the small satellite of worked example 6.2 at 400 km (ρ = 2.8e-12 kg/m³) the period shortens by 0.34 s per day. Run the same relation backwards and a measured period decay gives the density along the orbit.[R14]
Key equations
where
where
zero at the critical inclination 63.4°
where
circular orbit, SI units
where
where
valid when r ≪ d
where
Worked examples
How fast does the orbit plane turn?[R1]
Given
A circular orbit has altitude 800 km and inclination 70°. Find the J₂ node rate and the perigee rate in degrees per day.
Working
- a = RE + h = 6378.137 + 800 = 7178.137 km; n = √(μ/a³)1.038129e-03 rad/s
- Common factor n J₂ (RE/p)² = 1.038129e-03 × 0.00108262668 × (6378.137/7178.137)²8.8735e-07 rad/s
- Ω̇ = −1.5 × that × cos 70° = -4.5524e-07 rad/s-2.254 °/day
- ω̇ = 0.75 × that × (5 cos² 70° − 1) = -2.7626e-07 rad/s-1.368 °/day
Answer. Ω̇ = -2.25 °/day (westward), ω̇ = -1.37 °/day.
Sanity check. The node moves about 823° in a year: J₂ is not a small effect over a mission. A negative ω̇ is expected above the critical inclination of 63.4°.
n = 1.038129e-03 rad/s, factor = 8.8735e-07 rad/s node rate = -4.5524e-07 rad/s = -2.254 deg/day perigee rate = -2.7626e-07 rad/s = -1.368 deg/day
How fast does drag pull a small satellite down?[R14, R15]
Given
A small satellite (CD = 2.2, A/m = 0.01 m²/kg) is in a circular orbit at 400 km, where the density is about 2.8e-12 kg/m³. Find the height lost per orbit and per day, then estimate the time to fall to 200 km.
new symbols
Working
- a = 6778.137 km6.7781e+06 m
- Δa per rev = −2π CD (A/m) ρ a² = −2π × 2.2 × 0.01 × 2.80e-12 × (6.7781e+06)²-17.78 m
- Revolutions per day = 86,400 / 5553.6 s = 15.56, so Δa per day-277 m
- Local scale height from the density table: H = 54.6 km. Density grows quickly as the orbit drops, so the decay speeds up.
- Integrating da/dt = −√(μa) ρ(h) CD (A/m) with the table down to 200 km gives 185 days
Answer. About 18 m per orbit (0.28 km per day) at first; roughly 185 days (6 months) to reach 200 km.
Sanity check. The decay is not linear: at the start the satellite would need 723 days at the initial rate, yet it falls in 185. Halve the density (solar minimum) and the lifetime roughly doubles.
height lost: 17.78 m per orbit, 15.56 orbits/day, so 0.277 km per day scale height near 400 km: 54.6 km time to fall from 400 to 200 km: 184.9 days (about 6 months)
Practice
Node drift of a mid-inclination orbit[R1]
Given
A circular orbit at 600 km has inclination 45°. How many degrees per day does its ascending node move, and in which direction?
Use Ω̇ = −1.5 n J₂ (RE/p)² cos i with p = a for a circle, then convert rad/s to °/day.
Answer. Ω̇ = -5.14 °/day: the node drifts west by about 5.1 deg every day.
Working
- a = 6978.137 km, n1.083078e-03 rad/s
- Ω̇ = −1.5 × 1.083078e-03 × 0.00108262668 × (6378.137/6978.137)² × cos 45°-1.0390e-06 rad/s
- × (180/π) × 86,400-5.144 °/day
Answer. Ω̇ = -5.14 °/day: the node drifts west by about 5.1 deg every day.
Sanity check. Prograde orbit (i < 90°) gives a negative rate, as it must. 154° in a month is why constellation planes are spaced with J₂ in mind.
n = 1.083078e-03 rad/s, node rate = -1.0390e-06 rad/s = -5.144 deg/day: the node drifts west by 5.14 deg per day
Drag or sunlight: which pushes harder at 700 km?[R8, R14]
Given
A small satellite has CD = 2.2, CR = 1.3 and A/m = 0.01 m²/kg. At 700 km the density is about 3.1e-14 kg/m³. Compare the drag acceleration with the solar-radiation-pressure acceleration (solar flux 1361 W/m²).
Assumption: Density at 700 km from the U.S. Standard Atmosphere 1976 table (moderate solar activity).
Radiation pressure is P = S/c in N/m². Drag is ½ρv²CD(A/m); use the circular speed in m/s.
Answer. adrag = 1.9e-08 m/s² and aSRP = 5.9e-08 m/s²: at 700 km they are within a factor of 3.1 of each other.
Working
- P = S/c = 1361 / 2997924584.540e-06 N/m²
- aSRP = P CR (A/m) = 4.540e-06 × 1.3 × 0.015.90e-08 m/s²
- v = √(μ/r)7.504 km/s
- adrag = ½ × 3.10e-14 × (7504)² × 2.2 × 0.011.92e-08 m/s²
- Ratio drag / SRP0.33
Answer. adrag = 1.9e-08 m/s² and aSRP = 5.9e-08 m/s²: at 700 km they are within a factor of 3.1 of each other.
Sanity check. Both are below 10⁻⁷ m/s², about 6e-09 g: tiny, but they act for years. The density can change by tens of times with solar activity, so the winner at this height is not fixed.
P = 4.540e-06 N/m^2, a_SRP = 5.90e-08 m/s^2 v = 7.504 km/s, a_drag = 1.92e-08 m/s^2 drag/SRP = 0.33: within a factor of 3.1
Reboosting a space station[R14, R21]
Given
A large station (CD = 2.2, A/m = 0.0032 m²/kg) flies a circular orbit at 410 km. Assume the density stays at 2.36e-12 kg/m³. How much height does it lose in 30 days, and what total Δv does a two-burn reboost need to restore it?
Assumption: Density held constant over the month (it really rises as the orbit drops and changes with solar activity).
Find Δa per revolution, multiply by the number of revolutions, then do a small Hohmann transfer back up.
Answer. About 2.2 km lost in 30 days; the reboost costs about 1.3 m/s.
Working
- Δa per rev = −2π × 2.2 × 0.0032 × 2.36e-12 × (6.7881e+06 m)²-4.8 m
- Revolutions in 30 days = 30 × 86,400 / 5565.9465.7
- Height lost = 4.8 m × 465.72.24 km
- Hohmann from 407.76 km back to 410 km: Δv₁ = 0.63 m/s, Δv₂0.63 m/s
- Total = 1.27 m/s; small-change shortcut v Δa / (2a)1.27 m/s
Answer. About 2.2 km lost in 30 days; the reboost costs about 1.3 m/s.
Sanity check. A few metres per second per month is a modest budget, which is why stations reboost every few weeks rather than carrying huge tanks. The two methods agree to within 1 %. NASA's own ISS ephemeris pairs this area with CD = 1.1, so CD = 2.2 makes this an upper-end estimate.
da per rev = 4.81 m, 465.7 revs, height lost = 2.24 km reboost: 0.63 + 0.63 = 1.26 m/s shortcut v da / (2a) = 1.26 m/s
Common mistakes
- Using degrees inside n or mixing rad/s with °/day. Compute in rad/s, convert at the end.
- Forgetting the sign: prograde orbits regress west (Ω̇ < 0), retrograde ones east.
- Treating drag decay as linear in time. Density rises as the orbit drops, so the end comes fast.
- Putting kilometres into the drag formula. Use metres and kg/m³ together.
Check yourself
Recap from earlier chapters
(Chapter 1) What is the period of a circular orbit at 400 km?
(Chapter 3) What is the specific energy of that orbit?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R8Mamajek, E.E., Prša, A., Torres, G. et al., "IAU 2015 Resolution B3 on Recommended Nominal Conversion Constants for Selected Solar and Planetary Properties", arXiv:1510.07674 (2015).
- R14NOAA, NASA and USAF, U.S. Standard Atmosphere, 1976, NASA-TM-X-74335 (NTRS 19770009539), Washington DC, Oct 1976.
- R15Boring, J.W. and Humphris, R.R., "Drag coefficients for spheres in free molecular flow in O at satellite velocities", NASA CR-2233, 1973.
- R21NASA JSC Flight Operations Directorate, TOPO office, ISS public ephemeris ISS.OEM_J2K_EPH.txt (CCSDS OEM, created 2026-11-02).
Where in each source, and links: References.
Designer orbits: sun-synchronous, frozen and Molniya
Chapter 6 treated J₂ as a nuisance. Designers turn it into a tool. Choose the inclination so the node turns exactly once a year and the orbit keeps the same angle to the Sun all year: sun-synchronous. Choose a small eccentricity and put perigee over the pole, and J₂ and J₃ cancel so the ellipse stops wandering: a frozen orbit. Choose 63.4° and the perigee stays fixed, so a long, slow apogee can hang over the far north: Molniya.
You will be able to
- Find the sun-synchronous inclination for any altitude
- Compute the frozen eccentricity that keeps altitude repeatable over each latitude
- Design a Molniya orbit and explain why it flies at 63.4°
Sun-synchronous
The node must move east by 360° per year, about one degree per day. J₂ gives an eastward drift only for retrograde orbits, so the inclination ends up just past polar, between about 96° and 102° for low orbits. The reward is a fixed local solar time at every pass: the same lighting for every image.[R5, R22]
Frozen
Odd zonal harmonics (mostly J₃) make the eccentricity and perigee oscillate. With perigee at 90° and one particular small eccentricity, the oscillation is zero, and the altitude over each latitude repeats orbit after orbit.
Molniya
A 12-hour orbit with a high apogee spends most of its time near apogee. Flown at the critical inclination, J₂ does not move the perigee, so the apogee stays over the northern hemisphere for years.
Key equations
where
where
where
where
Worked examples
A sun-synchronous orbit at 600 km[R5]
Given
Find the inclination that makes a circular 600 km orbit sun-synchronous, and the node rate that it gives.
new symbols
Working
- Sun-synchronous means the node turns once a year, eastward, to follow the Sun: Ω̇ = 2π / (365.2422 days)
- Required Ω̇ = 1.9911e-07 rad/s0.9856 °/day
- n = √(μ/a³) = 1.083078e-03 rad/s, a6978.137 km
- cos i = −Ω̇ / (1.5 n J₂ (RE/a)²)-0.13550
- i97.788°
Answer. i = 97.79°; the node advances 0.9856 °/day, exactly the Sun's mean motion.
Sanity check. Slightly retrograde, as every sun-synchronous orbit must be: only a retrograde orbit has an eastward node drift.
required node rate = 1.9911e-07 rad/s = 0.9856 deg/day n = 1.083078e-03 rad/s, cos i = -0.13550, i = 97.788 deg
A frozen orbit
Given
For a sun-synchronous orbit with a = 7078.137 km (i = 98.19°), find the frozen eccentricity with the perigee held at ω = 90°.
new symbols
Working
- J₂ and J₃ together try to change e and ω. They cancel when ω = 90° and e = −J₃ RE sin i / (2 J₂ a)
- e = −(-2.53215e-06) × 6378.137 × sin 98.19° / (2 × 0.00108262668 × 7078.137)0.00104
- Perigee altitude a(1 − e) − RE = 692.6 km; apogee altitude a(1 + e) − RE707.4 km
Answer. e = 0.0010 with ω = 90°: the altitude over each latitude repeats every orbit, varying between 693 and 707 km.
Sanity check. The orbit is still nearly circular (a 15 km swing), but now the swing is fixed relative to latitude, which keeps image scale and drag repeatable.
i = 98.19 deg, frozen e = 0.00104 altitude varies from 692.6 to 707.4 km
Practice
Sun-synchronous at 900 km[R5]
Given
What inclination makes a circular orbit at 900 km sun-synchronous?
Set −1.5 n J₂ (RE/a)² cos i equal to 2π per year and solve for cos i.
Answer. i = 99.03°.
Working
- a = 7278.137 km, n1.016807e-03 rad/s
- Ω̇ needed1.9911e-07 rad/s
- cos i = −1.9911e-07 / (1.5 × 1.016807e-03 × 0.00108262668 × (6378.137/7278.137)²)-0.15701
- i99.033°
Answer. i = 99.03°.
Sanity check. Higher than the 97.79° needed at 600 km, as the curve in the figure shows.
n = 1.016807e-03 rad/s, cos i = -0.15701 i = 99.033 deg
When the inclination is a little off
Given
A satellite meant to be sun-synchronous at 700 km is injected 0.1° too high in inclination. How far does the local time of its ascending node drift in a year?
Compute Ω̇ at the wrong inclination and subtract the Sun's rate. One degree of node offset is 4 minutes of local time.
Answer. The node runs ahead of the Sun by 4.37° per year, a local-time drift of about 17 min per year.
Working
- Correct i = 98.188°; actual i98.288°
- Ω̇(actual) = 2.015212e-07 rad/s; Sun's rate1.991064e-07 rad/s
- Difference × one year = 4.366° per year
- Local time: 4.366° × 4 min/° = 17.5 min per year
Answer. The node runs ahead of the Sun by 4.37° per year, a local-time drift of about 17 min per year.
Sanity check. A higher inclination makes cos i more negative, so the eastward drift grows and the node runs ahead of the Sun. A tenth of a degree matters: injection accuracy is part of the mission design.
correct i = 98.188 deg, actual i = 98.288 deg node rate 2.015212e-07 rad/s vs Sun 1.991064e-07 rad/s drift = 4.366 deg/yr = 17.5 min of local time per year
Design a Molniya orbit
Given
A Molniya orbit repeats twice per sidereal day (86164.0905 s) and has a perigee altitude of 1000 km. Find its apogee altitude and eccentricity, the fraction of each orbit spent in the half nearest apogee (eccentric anomaly from 90° to 270°), and the perigee drift if it were flown at 60° instead of the critical inclination.
a comes from the period. For the dwell, use Kepler's equation M = E − e sin E at E = 90° and 270°.
Answer. Apogee altitude 39367 km, e = 0.722; 73.0 % of each orbit is spent in the apogee half; at 60° the perigee would drift 0.037 deg/day.
Working
- a = (μ (T/2π)²)^(1/3) with T = 43082.05 s: a26561.8 km
- ra = 2a − rp = 2 × 26561.8 − 7378.137 = 45745.4 km, apogee altitude 39367 km
- e = (ra − rp)/(ra + rp)0.7222
- M(90°) = π/2 − e, M(270°) = 3π/2 + e, so the fraction is (π + 2e)/(2π) = ½ + e/π0.7299
- At i = 60°: ω̇ = 0.75 n J₂ (RE/p)² (5 cos² i − 1)0.0369 °/day
Answer. Apogee altitude 39367 km, e = 0.722; 73.0 % of each orbit is spent in the apogee half; at 60° the perigee would drift 0.037 deg/day.
Sanity check. Over a year that drift is 13°, which would move the slow apogee away from the high northern latitudes it is meant to serve. That is why Molniya orbits fly at 63.4°.
a = 26561.8 km, apogee radius 45745.4 km, apogee altitude 39367.3 km e = 0.7222, apogee-half dwell = 0.7299 (73.0 %) perigee drift at 60 deg = 0.0369 deg/day
Common mistakes
- Expecting a prograde sun-synchronous orbit. The node must move east, which needs cos i < 0.
- Using the solar day (86,400 s) for the Molniya period. It repeats on the sidereal day, so T = 43,082 s.
- Forgetting that the frozen condition also fixes ω = 90°: the eccentricity alone is not enough.
- Assuming sun-synchronous orbits stay put forever: drag lowers a, which changes the needed inclination.
Check yourself
Recap from earlier chapters
(Chapter 6) Which way does the node of a 45° orbit drift, and why?
(Chapter 2) A Molniya orbit has e = 0.722. How many times further out is apogee than perigee?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R5Williams, D.R., "Earth Fact Sheet", NASA NSSDCA, Goddard Space Flight Center, updated 15 Nov 2024.
- R22ESA, "Types of orbits", 30 Mar 2020.
Where in each source, and links: References.
Ground tracks
The ground track is the path of the point directly under the satellite. The orbit is (nearly) fixed in space while Earth turns beneath it, so each pass lands further west. Choose the period so that after a whole number of orbits Earth has turned a whole number of times, and the track repeats exactly: the basis of every mapping mission.
You will be able to
- Compute the westward shift between successive passes, with J₂
- Explain the sidereal day and why it is 4 minutes shorter than the solar day
- Find the altitude of a repeating ground track by root finding
Two kinds of day
Earth turns once relative to the stars in 86164.0905 s, the sidereal day. Because Earth also moves round the Sun, it must turn a little more to face the Sun again: the solar day of 86,400 s. Orbits care about the sidereal day; clocks follow the solar one.[R2]
The westward shift
During one nodal period Earth turns east by ω_E TN, while J₂ swings the node by Ω̇ TN. The equator crossing moves west by their difference. The track reaches a highest latitude equal to the inclination (or 180° minus it for a retrograde orbit).
Repeating tracks
If k orbits take exactly as long as d turns of Earth relative to the orbit plane, the track closes after d days. Because J₂ changes both the nodal period and the node rate, the altitude that repeats is found numerically.
Key equations
where
where
where
where
Worked examples
How far west does each pass land?
Given
A station orbits at 420 km, inclination 51.6°. Find the nodal period (equator to equator) and how far west each new ascending pass crosses the equator.
new symbols
Working
- Kepler period 2π√(a³/μ)5578.2 s
- With J₂, the argument of latitude advances at n + Ṁ_J₂ + ω̇ = 1.127252e-03 rad/s, so TN5573.9 s
- Node rate Ω̇ = -1.0001e-06 rad/s; Earth turns at ω_E7.292115e-05 rad/s
- Shift per orbit = (ω_E − Ω̇) TN = (7.292115e-05 − (-1.0001e-06)) × 5573.9 = 0.41203 rad23.608°
Answer. TN = 92.90 min; each pass crosses the equator 23.61° further west (2628 km at the equator).
Sanity check. Without J₂ the shift would be ω_E T = 23.31°. The node's own westward regression adds a little to Earth's eastward turn.
Kepler period 5578.2 s; nodal period 5573.9 s = 92.90 min node rate -1.0001e-06 rad/s shift per orbit = 0.41203 rad = 23.608 deg = 2628 km at the equator
Practice
Spacing between passes
Given
For a circular orbit at 700 km and 98° inclination, how far apart in longitude are successive equator crossings? Compare with the answer that ignores J₂.
Without J₂ the shift is just Earth's rotation during one period. With J₂ use (ω_E − Ω̇) TN.
Answer. 24.72° with J₂ against 24.76° without: the passes are about 24.7° apart.
Working
- Kepler period T = 5926.4 s; ω_E T24.761°
- With J₂: TN = 5933.6 s, Ω̇1.9457e-07 rad/s(eastward, retrograde orbit)
- (ω_E − Ω̇) TN24.725°
Answer. 24.72° with J₂ against 24.76° without: the passes are about 24.7° apart.
Sanity check. For a retrograde orbit the node follows Earth's spin, so the shift is slightly smaller than Earth's rotation alone.
without J2: T = 5926.4 s, spacing 24.761 deg with J2: T_N = 5933.6 s, node rate 1.9457e-07 rad/s, spacing 24.725 deg
The figure-8 of an inclined geosynchronous orbit
Given
A geosynchronous satellite has drifted to 3° inclination (circular, period one sidereal day). Its ground track is a small figure-8. Find its north-south and east-west half-sizes.
Latitude swings by ±i. In longitude, along-track motion is slightly faster than Earth at the nodes and slower near the extremes of latitude; for small i the half-width is i²/4 (radians).
Answer. ±3° in latitude but only ±0.039° in longitude: a tall, thin figure-8.
Working
- North-south: the track reaches latitude ±i = ±3°
- East-west, small-i result: half-width = i²/4 = (0.05236)²/4 = 6.854e-04 rad0.0393°
- Numerical ground track over one sidereal day: longitude range ±0.0393°
Answer. ±3° in latitude but only ±0.039° in longitude: a tall, thin figure-8.
Sanity check. 0.039° at GEO radius is 29 km: a fixed ground antenna with a wide beam still sees it; a narrow beam must track north-south.
north-south half-size: 3.000 deg east-west half-size: i^2/4 = 6.854e-04 rad = 0.0393 deg numerical ground track: longitude range +/- 0.0393 deg
A sun-synchronous orbit that repeats every two days
Given
Find the altitude of a circular sun-synchronous orbit that makes exactly 29 revolutions in 2 days and then retraces its ground track. Include J₂ in both the nodal period and the node rate.
Repeat condition: k TN = d × 2π/(ω_E − Ω̇). For a sun-synchronous orbit ω_E − Ω̇ is exactly 2π per solar day. Solve numerically: a and i depend on each other.
Answer. h = 720.0 km at i = 98.27°: 29 orbits in 2 days, then the track repeats.
Working
- For SSO, ω_E − Ω̇ = 2π/86,400 s, so the condition is TN = 2 × 86,400 / 295958.62 s
- For a trial altitude h: i = SSO inclination(h), then TN(h, i) from n + Ṁ_J₂ + ω̇
- Root of 29 TN(h) − 2 × 86,400 = 0 by Brent's method: h719.95 km
- Check: i = 98.270°, TN = 5958.621 s, 29 × TN = 172800.0 s2 days
Answer. h = 720.0 km at i = 98.27°: 29 orbits in 2 days, then the track repeats.
Sanity check. 14.5 revolutions per day sits between the 14 and 15 rev/day families, as expected for an orbit near 700 km.
h = 719.95 km, i = 98.270 deg, T_N = 5958.621 s check: 29 x T_N = 172800.0 s = 2.0000 days
Common mistakes
- Using 86,400 s for Earth's rotation. Earth turns 360° in a sidereal day.
- Forgetting the node rate: for a near-polar orbit it changes the shift only a little, but it decides whether a track repeats.
- Expecting a track above the inclination: the maximum latitude equals i (or 180° − i).
Check yourself
Recap from earlier chapters
(Chapter 7) Why does a sun-synchronous orbit's repeat condition use the solar day?
(Chapter 1) What is the ISS-like orbit's Kepler period at 420 km?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
Where in each source, and links: References.
Coplanar manoeuvres: Hohmann, bi-elliptic and phasing
Every in-plane manoeuvre is a change of speed at a point, which changes the energy and so the semi-major axis. The Hohmann transfer uses two tangential burns at opposite ends of an ellipse and is the cheapest two-burn transfer between circles. For very large changes, going out much further first and coming back is cheaper still. And to move along your own orbit, you change your lap time: slow down to catch up.
You will be able to
- Size a Hohmann transfer and the propellant it needs
- Say when a bi-elliptic transfer beats Hohmann, and why 11.94 is the number
- Design a phasing manoeuvre to move along the same orbit
Hohmann
Burn prograde to stretch the orbit into an ellipse that just touches the target circle, coast half an orbit, burn prograde again to round it off. Both burns are along the velocity, which wastes nothing.[R1]
Bi-elliptic
Three burns: out to a far apoapsis, raise the periapsis there where the speed is tiny, then come in and circularise. When the radius ratio is above 11.94 this can beat Hohmann; above 15.58 every bi-elliptic does. It always takes much longer.
Phasing
To meet a point ahead of you in the same orbit, enter a slightly smaller ellipse with a shorter period, lap once and return. To fall behind, use a larger one. The rocket equation converts every Δv into propellant mass.
One-tangent burn
A Hohmann transfer is the cheapest two-burn route but the slowest. Aim the transfer ellipse past the target and the second burn happens where the ellipse cuts the target orbit at an angle. From a 300 km circle to GEO with the transfer apogee at 1.5 GEO radii: 5.280 km/s in 3.01 h, against Hohmann's 3.893 km/s in 5.28 h. Time costs Δv.[R1]
Staging and gravity loss
The rocket equation applies stage by stage and the Δv values add. Dropping an empty stage stops you accelerating dead mass. With exhaust speed 3.0 km/s, 100 t at lift-off, 2 t of payload and structure 10 % of the stage mass, one stage gives 6.41 km/s; the same masses split into two equal stages give 7.75 km/s. Burning against gravity costs extra: for a vertical burn the ideal Δv is reduced by g₀ tb, 1.18 km/s for a 120 s burn, which is why rockets pitch over early.[R29]
Key equations
where
where
where
where
Worked examples
Raising a small satellite with a Hohmann transfer[R1, R9]
Given
A 150 kg satellite moves from a circular 400 km orbit to a circular 1200 km orbit with two burns of a hydrazine thruster (Isp = 220 s). Find Δv, time of flight and propellant mass.
new symbols
Working
- r₁ = 6778.137 km, r₂ = 7578.137 km, at = (r₁ + r₂)/27178.137 km
- Burn 1: √(μ(2/r₁ − 1/at)) − √(μ/r₁) = 7.8793 − 7.6686210.8 m/s
- Burn 2: √(μ/r₂) − √(μ(2/r₂ − 1/at)) = 7.2525 − 7.0475205.0 m/s
- Total Δv = 415.7 m/s; time of flight = half the transfer period50.4 min
- Rocket equation: m₀/mf = exp(Δv/(Isp g₀)) = exp(415.7/(220 × 9.80665))1.2125
- Propellant = m₀(1 − mf/m₀) = 150 × (1 − 1/1.2125)26.3 kg
Answer. Δv = 416 m/s in 50 min, using 26.3 kg of hydrazine.
Sanity check. That is 18 % of the satellite's mass for an 800 km climb: small satellites rarely have this much, which is why they are launched close to their final orbit.
burn 1 = 210.8 m/s, burn 2 = 205.0 m/s, total 415.7 m/s time of flight = 50.4 min mass ratio 1.2125, propellant 26.29 kg
Slow down to catch up: a phasing manoeuvre[R16]
Given
A satellite in a circular 800 km orbit must move 20° ahead of its current slot in one revolution of a phasing orbit. Find the phasing orbit and the Δv.
new symbols
Working
- Circular period T = 6052.4 s. To arrive 20° early, the phasing orbit must be shorter: Tph = T(1 − 20/360)5716.2 s
- aph = (μ(Tph/2π)²)^(1/3) = 6909.756 km; the burn point becomes apogee, so perigee rp = 2aph − r6641.374 km(altitude 263 km)
- Burn 1 (retrograde): Δv = √(μ/r) − √(μ(2/r − 1/aph)) = 7.4518 − 7.3057146.2 m/s
- After one phasing orbit, burn 2 (prograde) restores the circle: another 146.2 m/s
- Total 292.3 m/s
Answer. Drop perigee to 263 km with a 146 m/s braking burn, wait one orbit, then speed up again: 292 m/s in total.
Sanity check. Braking to get ahead feels wrong, but a lower orbit is a shorter, faster lap. Spreading the move over several phasing orbits cuts the Δv almost in proportion.
T = 6052.4 s, T_ph = 5716.2 s, a_ph = 6909.756 km phasing perigee altitude 263.2 km braking burn 146.2 m/s, same again to restore: total 292.3 m/s
Practice
A short Hohmann climb[R1]
Given
Find the two burns and the total Δv to go from a circular 300 km orbit to a circular 900 km orbit.
Transfer semi-major axis is the mean of the two radii; use vis-viva at each end.
Answer. Δv₁ = 164 m/s, Δv₂ = 161 m/s, total 325 m/s.
Working
- at = (6678.137 + 7278.137)/26978.137 km
- Δv₁ = 7.8901 − 7.7258164.3 m/s
- Δv₂ = 7.4005 − 7.2396160.8 m/s
- Total325.1 m/s
Answer. Δv₁ = 164 m/s, Δv₂ = 161 m/s, total 325 m/s.
Sanity check. The first burn is slightly bigger, as always for a Hohmann climb. Sanity: the difference of circular speeds is 325 m/s, close to the total for a small change.
a_t = 6978.137 km dv1 = 164.3 m/s, dv2 = 160.8 m/s, total 325.1 m/s
Falling back in a navigation-satellite orbit[R1]
Given
A satellite in a circular 20200 km orbit must drop 30° behind its slot, spread over 3 revolutions of a phasing orbit. Find the phasing orbit's semi-major axis and the total Δv.
To fall behind you need a longer lap: burn prograde. Each phasing revolution must be longer by (Δθ/360°)/k of a period.
Answer. aph = 27068 km and a total Δv of 69.8 m/s over 3 revolutions.
Working
- T = 43121.9 s; Tph = T(1 + 30/(360 × 3))44319.7 s
- aph = 27068.07 km; the burn point is perigee and apogee radius27558.0 km
- Burn 1 = √(μ(2/r − 1/aph)) − √(μ/r) = 34.89 m/s; burn 2 is the same, retrograde
- Total69.78 m/s
Answer. aph = 27068 km and a total Δv of 69.8 m/s over 3 revolutions.
Sanity check. The total time is 36.9 h. Doing it in one revolution would cost about 3 times more.
T = 43121.9 s, T_ph = 44319.7 s a_ph = 27068.07 km, apogee radius 27558.0 km each burn 34.89 m/s, total 69.78 m/s
When going further out is cheaper[R1]
Given
A spacecraft goes from a circular orbit of radius 6978.1 km to one 15 times larger (104672 km). Compare a Hohmann transfer with a bi-elliptic transfer whose intermediate apoapsis is 314016 km (three times the target). Give Δv and time of flight for both.
Bi-elliptic: burn to an ellipse with apoapsis rb, at rb raise periapsis to r₂, at r₂ circularise. Times of flight are half-periods of the two ellipses.
Answer. Bi-elliptic 4014 m/s against Hohmann 4053 m/s, a saving of 39 m/s, but 221 h instead of 18 h.
Working
- Hohmann: Δv = 2791.2 + 1261.5 = 4052.7 m/s; t18.23 h
- Bi-elliptic burn 1 (to rb): 3013.8 m/s
- Burn 2 at rb (periapsis up to r₂): 561.7 m/s
- Burn 3 at r₂ (circularise, braking): 438.6 m/s
- Bi-elliptic total 4014.1 m/s; t221.3 h
Answer. Bi-elliptic 4014 m/s against Hohmann 4053 m/s, a saving of 39 m/s, but 221 h instead of 18 h.
Sanity check. R = 15 is above the 11.94 crossover, so a bi-elliptic can win; the saving is small (1.0 %) and the trip is 12 times longer.
Hohmann: 4052.7 m/s, 18.23 h bi-elliptic: 3013.8 + 561.7 + 438.6 = 4014.1 m/s, 221.3 h saving 38.6 m/s
Common mistakes
- Subtracting speeds at different radii. Each Δv compares two speeds at the same point.
- Speeding up to catch a target ahead. That raises the orbit and lengthens the lap: you fall further behind.
- Using kilometres per second with Isp g₀ in metres per second inside the rocket equation.
- Choosing a bi-elliptic for a mission with a deadline: the Δv saving costs days of flight.
Check yourself
Recap from earlier chapters
(Chapter 3) What is the speed at the perigee of the 400-to-1200 km transfer ellipse?
(Chapter 4) How long does the Hohmann coast take?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R9NIST, CODATA 2022 recommended values, "standard acceleration of gravity".
- R16NASA Small Spacecraft Systems Virtual Institute, State-of-the-Art of Small Spacecraft Technology, Ch. 4 "In-Space Propulsion", 2026 edition.
- R29NASA Glenn Research Center, "Ideal Rocket Equation", Beginner's Guide to Aeronautics, updated 21 Nov 2023.
Where in each source, and links: References.
Plane changes and combined manoeuvres
Turning an orbit plane means turning the velocity vector, and that costs in proportion to the speed you turn. At low-orbit speed a few degrees is very expensive; far out, where you move slowly, it is cheap. Better still, fold the turn into a burn you already need: velocity changes add as vectors, so the extra cost of a small turn inside a large burn is tiny.
You will be able to
- Compute the cost of a pure plane change and say where to do it
- Combine a speed change and a plane change in one burn
- Optimise how a plane change is split between two burns
A pure turn
To rotate the plane by Δi without changing speed, the burn must cancel part of the old velocity and add part of the new: Δv = 2v sin(Δi/2). It is done at a node, where the old and new planes intersect.
Turn where you are slow
The cost scales with v, so turns belong at apogee or far from the planet. For a geostationary mission the turn goes with the apogee burn into GEO.
One burn, two jobs
When the speed also changes, one burn along the vector difference is cheapest: the law of cosines. Moving a few degrees of the turn into the perigee burn reduces the total a little more.
Key equations
where
where
where
Worked examples
The same turn at two speeds
Given
Find the Δv to turn an orbit plane by 10° in a circular 500 km orbit and in GEO.
new symbols
Working
- Pure plane change: Δv = 2v sin(Δi/2)
- LEO: v = 7.6126 km/s, Δv = 2 × 7.6126 × sin 5.0°1.3270 km/s
- GEO: v = 3.0747 km/s, Δv = 2 × 3.0747 × sin 5.0°0.5359 km/s
Answer. 1327 m/s in LEO, 536 m/s in GEO.
Sanity check. The ratio 2.48 is just the ratio of speeds. A 10° turn in LEO costs more than a Hohmann climb from 500 km to 1300 km (407 m/s).
LEO: v = 7.6126 km/s, dv = 1.3270 km/s = 1327.0 m/s GEO: v = 3.0747 km/s, dv = 0.5359 km/s = 535.9 m/s
One burn at apogee: circularise and turn together
Given
A transfer orbit from a 28.5° latitude site has perigee altitude 250 km, apogee at GEO radius and inclination 28.5°. Compare one combined burn at apogee with circularising first and then turning the plane.
new symbols
Working
- a = (6628.137 + 42164.17)/2 = 24396.15 km; apogee speed va = √(μ(2/ra − 1/a))1.6026 km/s
- GEO speed vGEO3.0747 km/s
- Combined: Δv = √(v_a² + v_GEO² − 2 va vGEO cos Δi) = √(1.6026² + 3.0747² − 2 × 1.6026 × 3.0747 × cos 28.5°)1.8333 km/s
- Separate: (3.0747 − 1.6026) + 2 × 3.0747 × sin 14.25°2.9857 km/s
Answer. Combined 1833 m/s against 2986 m/s for two separate burns: a saving of 1152 m/s.
Sanity check. Velocity changes are vectors. One burn along the vector difference is always shorter than two burns along two sides of the triangle.
apogee speed 1.6026 km/s, GEO speed 3.0747 km/s combined burn 1833.3 m/s, separate burns 2985.7 m/s saving 1152.4 m/s
Practice
A small plane trim
Given
A satellite in a circular 700 km orbit must change its inclination by 5°. What Δv does that cost?
Δv = 2v sin(Δi/2) with the circular speed.
Answer. Δv = 655 m/s.
Working
- v = √(μ/r) = √(398600.4418/7078.137)7.5043 km/s
- Δv = 2 × 7.5043 × sin 2.5°0.6547 km/s
Answer. Δv = 655 m/s.
Sanity check. About 8.7 % of the orbital speed for a 5° turn; for small angles Δv ≈ v Δi = 655 m/s.
v = 7.5043 km/s, dv = 0.6547 km/s = 654.7 m/s
Why launch sites near the equator win for GEO[R25]
Given
Using the transfer orbit of worked example 10.2 (apogee speed 1.6026 km/s), find the combined apogee burn into GEO for transfer orbits inclined at 5.2° and at 45° (the site latitudes).
Same law of cosines, only the angle changes.
Answer. 1486 m/s from 5.2° and 2248 m/s from 45°, a difference of 762 m/s.
Working
- Δv = √(v_a² + v_GEO² − 2 va vGEO cos Δi), va = 1.6026, vGEO3.0747 km/s
- Δi = 5.2°: Δv1.4857 km/s
- Δi = 45°: Δv2.2480 km/s
Answer. 1486 m/s from 5.2° and 2248 m/s from 45°, a difference of 762 m/s.
Sanity check. The 5.2° case is close to the pure circularisation cost of 1472 m/s, as it should be for a small angle.
di = 5.2 deg: combined apogee burn 1485.7 m/s di = 45.0 deg: combined apogee burn 2248.0 m/s difference 762.2 m/s
Where to do the plane change: an optimisation
Given
A spacecraft goes from a circular 300 km orbit at 40° inclination to an equatorial GEO with a Hohmann transfer. A share s of the plane change is done in the perigee burn and the rest in the apogee burn. Find the best s and the saving compared with doing it all at apogee.
Write the total as two law-of-cosines burns, one at each end, as a function of s. Minimise numerically (bounded scalar search).
Answer. Do about 2.6° of the turn at perigee (s = 0.066); total 4504 m/s, saving 34 m/s.
Working
- Perigee: vc1 = 7.7258, vp = 10.1515 km/s. Apogee: va = 1.6078, vGEO3.0747 km/s
- Δv(s) = √(v_c1² + v_p² − 2v_c1vp cos(sΔi)) + √(v_a² + v_GEO² − 2v_avGEO cos((1 − s)Δi))
- All at apogee, s = 0: Δv4.5387 km/s
- Bounded minimisation: s = 0.0659 (2.64° at perigee), Δv4.5043 km/s
Answer. Do about 2.6° of the turn at perigee (s = 0.066); total 4504 m/s, saving 34 m/s.
Sanity check. A small turn added to a big burn costs almost nothing (the vector barely gets longer), so the optimum always shifts a few degrees to perigee.
vc1 = 7.7258, vp = 10.1515, va = 1.6078, vGEO = 3.0747 km/s all at apogee: 4538.7 m/s best s = 0.0659 (2.64 deg at perigee): 4504.3 m/s saving 34.4 m/s
Common mistakes
- Using degrees inside the small-angle form. Δv ≈ v Δi needs Δi in radians.
- Doing the plane change in LEO because it is the first burn. It is the most expensive place.
- Adding the speed change and the plane change as scalars: they are vectors.
- Burning away from a node: a plane change has to happen where the two planes cross.
Check yourself
Recap from earlier chapters
(Chapter 9) What is the circularisation Δv at GEO for the 250 km-perigee transfer, with no plane change?
(Chapter 2) Which element does a pure plane change at the node alter?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R25ESA, "Europe's Spaceport".
Where in each source, and links: References.
Rendezvous: Clohessy-Wiltshire
Rendezvous is orbital mechanics seen from the target. Put the target at the origin and ride along with it: radial up, along-track forward, cross-track out of the plane. Close to the target the equations become linear and have a neat closed-form solution, the Clohessy-Wiltshire (CW) equations. They contain the great paradox of orbits: to catch up with something ahead of you, you slow down. A slower craft drops to a lower orbit, a lower orbit is faster, and it overtakes.
You will be able to
- Describe relative motion in the target's LVLH frame (radial, along-track, cross-track)
- Explain why a retrograde burn makes a chaser move ahead
- Plan a two-burn approach with the Clohessy-Wiltshire solution and check it against the full orbit
The LVLH frame
Local vertical, local horizontal: x points radially out from Earth through the target, y points along the target's velocity, z completes the right-handed set along the orbit normal. The line along y through the target is the V-bar; the line along x is the R-bar. Approaches to a space station come in along one of them.
The CW equations
For a target on a circular orbit and a chaser a few kilometres away, gravity differences and the rotating frame give three linear equations. The radial and along-track motions are coupled through Coriolis terms 2n; the cross-track motion is a simple oscillation at the orbit rate. Because they are linear, the motion after any burn is a matrix times the state: the state transition matrix.
Slow down to catch up
A burn along +y raises the opposite side of the orbit; the craft climbs, its period grows and it falls behind. A burn along −y does the reverse: it drops, speeds up relative to the target and drifts ahead by 6π Δv/n per orbit. The relative path is a looping cycloid, as in the figure.
Checking the linear model
CW assumes a circular target orbit, small separations and no perturbations. Measure x as the difference in radius and y as the distance along the curved orbit, and for a 10 km separation at 400 km altitude the linear prediction is within a few metres of the full two-body answer; at hundreds of kilometres it breaks down, and you go back to Lambert's problem (Chapter 12).
Key equations
where
valid for a circular target orbit and separations much smaller than the orbit radius
where
where
where
where
Worked examples
Slow down to catch up
Given
A target is on a circular orbit at 400 km. A chaser sits 10 km behind it on the same orbit (on the V-bar). It wants to arrive at the target after exactly one orbit with two along-track burns. Find the burn, the total Δv and how low the chaser dips.
new symbols
Working
- Mean motion n = √(μ/a³) = √(398600.4418 / 6778.137³) = 1.131367e-03 rad/s; period T = 2π/n92.56 min
- An along-track burn Δvy moves the chaser by Δy = −6π Δvy / n per orbit (CW solution at t = T)
- We need Δy = +10 km, so Δvy = −10 × 1.131367e-03 / (6π)-0.600 m/s(a retrograde burn)
- After one orbit the CW solution gives position (0.0e+00, 1.8e-15) km and along-track rate -0.600 m/s, so a second burn of +0.600 m/s stops it
- Total Δv = 2 × 0.600 = 1.200 m/s; the lowest point is x = 4Δvy/n = -2122 m, half an orbit after the burn
- Full two-body check: the chaser arrives 3.1 m from the target, 0.03 % of the gap
Answer. Burn 0.60 m/s retrograde, wait one orbit (92.6 min), burn 0.60 m/s prograde: 1.20 m/s in total. The chaser dips 2122 m below the target's orbit.
Sanity check. Slowing down lowers the orbit, a lower orbit has a shorter period, and so the chaser laps ahead. The tiny dip (hundreds of metres for a 10 km gap) is why phasing near a station is so cheap.
n = 1.131367e-03 rad/s, T = 92.56 min burn 1 = -0.600 m/s (retrograde), i.e. 0.60 m/s after one orbit: x = -0.0e+00 km, y = 1.8e-15 km, ydot = -0.600 m/s burn 2 = 0.60 m/s prograde; total 1.20 m/s the chaser dips 2122 m below the target's orbit
A two-burn approach in half an orbit
Given
The chaser is 1 km below and 8 km behind the target (same 400 km orbit for the target). Being on a lower circular orbit it already drifts forward at ẏ = −(3/2)n x = 1.70 m/s. Plan a two-burn CW transfer that reaches the target in half an orbit.
new symbols
Working
- Transfer time t = T/2 = 46.28 min, so nt = π
- Required velocity just after burn 1: v₀⁺ = −Φ_rv⁻¹ Φ_rr r₀ = (-1.596, 1.980, 0) m/s
- Burn 1: Δv₁ = v₀⁺ − v₀ = (-1.596, 0.283, 0) m/s, size 1.621 m/s
- Arrival velocity from the CW solution; burn 2 cancels it: |Δv₂|1.621 m/s
- Total = 3.242 m/s. Full two-body propagation of the same burn misses by 0.8 m
Answer. Δv₁ = 1.62 m/s, Δv₂ = 1.62 m/s, total 3.24 m/s, arriving 46.3 min later.
Sanity check. Metres per second, not kilometres per second: relative motion near a target is a game of very small burns. The linear CW model is good to a few metres here because the separation is tiny compared with the orbit radius.
drift rate before burn 1: 1.70 m/s velocity after burn 1: (-1.596, 1.980) m/s dv1 = 1.621 m/s, dv2 = 1.621 m/s, total 3.242 m/s arrival 46.28 min after burn 1
Practice
How far does a tiny burn move you?
Given
A satellite on a circular orbit at 500 km makes an along-track burn of -0.2 m/s (retrograde). How far, and in which direction, has it moved along-track relative to where it would have been, after one orbit?
Use Δy = −6π Δvy / n with n in rad/s and Δv in km/s.
Answer. Δy = +3.41 km: about 3.4 km ahead after one orbit.
Working
- a = 6378.137 + 500 = 6878.137 km; n = √(μ/a³)1.106783e-03 rad/s
- Δy = −6π × (-2.0e-04) / 1.106783e-033.406 km
- Positive y is ahead of the reference point
Answer. Δy = +3.41 km: about 3.4 km ahead after one orbit.
Sanity check. 0.2 m/s is about 1/38,063 of the orbital speed, yet the effect grows every orbit: drift is how operators move satellites along a constellation.
n = 1.106783e-03 rad/s dy = +3.406 km (ahead)
Drifting in from below
Given
A chaser is on a circular orbit 2 km below a target at 450 km, 20 km behind it. How fast does it close the gap, and how long until it reaches the target's along-track position?
A slightly lower circular orbit drifts forward at ẏ = −(3/2) n x.
Answer. ẏ = 3.36 m/s, so the gap closes in 1.7 h.
Working
- n = √(μ/(6378.137 + 450)³)1.118963e-03 rad/s
- ẏ = −1.5 × 1.118963e-03 × (-2)3.357 m/s
- Time = 20 km / 3.357 m/s = 5958 s1.65 h
- Exact check: (n_chaser − n_target) × a3.358 m/s
Answer. ẏ = 3.36 m/s, so the gap closes in 1.7 h.
Sanity check. The linear rate and the exact difference of angular rates agree to 0.04 %, as they should for a 2 km offset.
n = 1.118963e-03 rad/s, drift 3.357 m/s (exact 3.358 m/s) time to close 20 km: 5958 s = 1.655 h
Coming down the R-bar
Given
A chaser is released 3 km directly above a target (on the R-bar) at 400 km, with zero velocity relative to the LVLH frame. Plan a two-burn CW transfer that reaches the target in a quarter of an orbit. Give both burns and the total.
Use the CW state transition matrix at nt = π/2: v₀⁺ = −Φ_rv⁻¹ Φ_rr r₀. The second burn cancels the arrival velocity.
Answer. Δv₁ = 6.27 m/s, Δv₂ = 2.19 m/s, total 8.47 m/s.
Working
- n = 1.131367e-03 rad/s, t = T/4 = 23.14 min, nt = π/2
- v₀⁺ = (-4.130, -4.723, 0) m/s, so Δv₁6.274 m/s
- Arrival velocity cancelled by Δv₂2.192 m/s
- Total8.466 m/s
- Two-body check: the chaser arrives 6.0 m from the target
Answer. Δv₁ = 6.27 m/s, Δv₂ = 2.19 m/s, total 8.47 m/s.
Sanity check. The first burn has a large along-track part, not just a radial push down: in orbit, radial and along-track motion are coupled by the Coriolis terms.
transfer time 23.14 min, velocity after burn 1 (-4.130, -4.723) m/s dv1 = 6.274 m/s, dv2 = 2.192 m/s, total 8.466 m/s
Common mistakes
- Burning towards the target to catch it. A forward burn raises the orbit and the chaser falls behind.
- Using CW for large separations or an elliptical target orbit. Check the separation is a small fraction of the radius.
- Mixing the inertial and rotating-frame velocities. CW velocities are measured in the rotating LVLH frame.
- Forgetting the second burn. A two-impulse transfer arrives with a relative velocity that must be cancelled.
Check yourself
Recap from earlier chapters
(Chapter 3) What is the circular speed at 400 km altitude?
(Chapter 1) What is the orbital period at 400 km?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
Where in each source, and links: References.
Lambert's problem
Most orbit questions are 'given where I am and how fast I'm going, where will I be?'. Mission design asks the opposite: I am here, I must be there at that time; what velocity do I need? That is Lambert's problem. Two position vectors and a flight time define exactly one conic (for less than one revolution, going the short or the long way). Every launch window, every interplanetary porkchop plot and every long-range rendezvous is a grid of Lambert solves.
You will be able to
- Say why Lambert's problem sits behind launch windows, targeting and rendezvous
- Solve it with the universal-variable method and Stumpff functions
- Verify a solver against a Hohmann case and by propagating the answer
Why it matters
Targeting a planet, choosing a launch date, intercepting a satellite, or picking up a sample at a given time are all boundary-value problems: the end points and the time are fixed, the orbit is unknown. Lambert's problem answers them, and its output (v₁ and v₂) gives the burns at both ends by simple subtraction.[R13]
The universal-variable method
Write the unknown orbit with the universal variable z, which is positive for ellipses, zero for parabolas and negative for hyperbolas. The flight time becomes a single monotonic function of z built from the Stumpff functions C(z) and S(z). Solve F(z) = 0 with a bracketed root finder, then the Lagrange coefficients f, g and ġ give both velocities.[R13]
Short way, long way, and 180°
For a given pair of points you can go the short way (Δθ < 180°) or the long way (Δθ > 180°); the direction of motion decides which. At exactly 180° the two vectors do not define a plane, so the problem is singular; real tools handle it with a small plane offset. Allowing extra full revolutions adds further solutions, each pair with its own minimum time.
Trust, but verify
A Lambert solver should pass three tests: it rebuilds a known ellipse from two of its points, its answer propagated forward lands on r₂, and the orbit's energy is the same at both ends. This chapter's solver passes all three, and reproduces a standard textbook-style test case to four decimals.[R12, R13]
Key equations
where
where
A = 0 at Δθ = 180°: the singular case
where
where
where
where
Worked examples
Lambert finds the Hohmann ellipse[R1]
Given
Take the Hohmann transfer from a 300 km circular orbit to GEO (r = 42164.2 km). Pick the departure point (perigee) and the point 170° further round the ellipse, and the flight time between them. Lambert's solver, given only the two positions and the time, must rebuild the same ellipse.
Working
- Transfer ellipse: at = (6678.137 + 42164.2)/2 = 24421.2 km, e0.72654
- Time from ν = 0 to ν = 170° by Kepler's equation (Chapter 4): Δt = 14532.2 s4.037 h
- Lambert(r₁, r₂, Δt) returns |v₁|10.151493 km/s
- Vis-viva at perigee: √(μ(2/rp − 1/at))10.151493 km/s
- Orbit from Lambert's v₁: a = 24421.2 km, e0.72654
Answer. Lambert gives 10.1515 km/s, exactly the Hohmann perigee speed, and the same ellipse (a = 24421 km).
Sanity check. Two positions plus a time define one conic. If you feed in points from a known ellipse, you must get that ellipse back: this is the first test any Lambert solver should pass.
transfer ellipse a = 24421.2 km, e = 0.72654, flight time 14532.2 s = 4.037 h Lambert |v1| = 10.151493 km/s, vis-viva perigee speed 10.151493 km/s orbit from Lambert: a = 24421.2 km, e = 0.72654
Aim for a point, arrive on time
Given
A spacecraft on a circular orbit of radius 7000 km wants to reach the point at radius 12000 km that lies 110° ahead in the same plane, in exactly 60 min. Find the departure velocity, the transfer orbit and the departure burn.
Working
- Inputs: r₁ = 7000 km, r₂ = 12000 km, Δθ = 110°, Δt3600 s(prograde, short way)
- A = sin Δθ √(r₁r₂/(1 − cos Δθ))7434.4 km
- Solve √μ Δt = (y/C)3/2 S + A√y for z: z5.57899(z > 0, so an ellipse)
- Lagrange coefficients give v₁ = (2.7041, 7.8177, 0) km/s, |v₁|8.2721 km/s
- Transfer orbit from (r₁, v₁): a = 8768.6 km, e0.3784
- Circular speed at r₁ is 7.5461 km/s along +y, so the burn is |v₁ − vc|2.7177 km/s
- Check: numerical propagation of v₁ for 3600 s misses r₂ by 7.0e-10 km
Answer. v₁ = 8.272 km/s on an ellipse with a = 8769 km and e = 0.378; departure burn 2.718 km/s.
Sanity check. The burn has a radial part as well as an along-track part: to arrive at a fixed time you generally leave at a flight-path angle, not tangentially.
A = 7434.4 km, z = 5.57899 v1 = (2.7041, 7.8177, 0) km/s, |v1| = 8.2721 km/s transfer orbit a = 8768.6 km, e = 0.3784 departure burn 2.7177 km/s
Practice
The cheapest-energy ellipse
Given
Two points at r₁ = 8000 km and r₂ = 10000 km are 90° apart. Find the chord, the semi-perimeter and the semi-major axis of the minimum-energy transfer ellipse between them.
Chord by the law of cosines; semi-perimeter s = (r₁ + r₂ + c)/2; the minimum-energy ellipse has a = s/2.
Answer. c = 12806 km, s = 15403 km, amin = 7702 km.
Working
- c = √(r₁² + r₂² − 2r₁r₂ cos Δθ)12806.2 km
- s = (8000 + 10000 + 12806.2)/215403.1 km
- amin = s/27701.6 km
- Its flight time (Lambert's theorem) is 3259 s; a Lambert solve at that time returns a7701.6 km
Answer. c = 12806 km, s = 15403 km, amin = 7702 km.
Sanity check. amin must exceed half the larger radius (5000 km) and is less than the mean radius of a Hohmann-like path: 7702 km fits.
chord c = 12806.2 km, semi-perimeter s = 15403.1 km minimum-energy a = 7701.6 km
Solve a Lambert problem[R13]
Given
A craft on a circular orbit at r₁ = 7000 km must reach r₂ = 14000 km, 140° ahead in the same plane, in 1.5 h. Find |v₁|, the transfer a and e, and the departure burn.
Use the universal-variable method: compute A, solve F(z) = 0, then f, g and v₁ = (r₂ − f r₁)/g.
Answer. |v₁| = 8.657 km/s, a = 10235 km, e = 0.401, departure burn 2.389 km/s.
Working
- A = sin Δθ √(r₁r₂/(1 − cos Δθ))4788.3 km
- Root of F(z): z8.33762
- v₁ = (2.2458, 8.3605, 0) km/s, |v₁|8.6569 km/s
- Transfer orbit: a = 10235.3 km, e0.4006
- Departure burn |v₁ − vc| with vc = 7.5461 km/s: 2.3890 km/s
Answer. |v₁| = 8.657 km/s, a = 10235 km, e = 0.401, departure burn 2.389 km/s.
Sanity check. r₂ = 14000 km must lie between the transfer perigee (6135 km) and apogee (14336 km), and it does.
A = 4788.3 km, z = 8.33762 v1 = (2.2458, 8.3605, 0) km/s, |v1| = 8.6569 km/s a = 10235.3 km, e = 0.4006, departure burn 2.3890 km/s
Intercept: choose the flight time
Given
A chaser on a circular orbit at 7000 km wants to meet a target on a coplanar circular orbit at 9000 km that is currently 30° ahead. For flight times of 30, 45 and 60 min, find the total Δv (departure plus matching the target's velocity on arrival) and pick the cheapest.
The target moves while you fly: its angle at arrival is the lead angle plus n_target × Δt. Solve one Lambert problem per flight time.
Answer. The 60 min transfer is cheapest at 0.89 km/s total.
Working
- Target mean motion n7.394437e-04 rad/s
- Δt = 30 min: target at 106.3°, Δvdep = 1.006, Δvarr = 1.797, total 2.803 km/s
- Δt = 45 min: target at 144.4°, Δvdep = 0.523, Δvarr = 0.739, total 1.262 km/s
- Δt = 60 min: target at 182.5°, Δvdep = 0.459, Δvarr = 0.430, total 0.889 km/s
- Cheapest: 60 min with 0.889 km/s
Answer. The 60 min transfer is cheapest at 0.89 km/s total.
Sanity check. For comparison a Hohmann transfer between the two radii costs 0.888 km/s but needs the right phasing; paying for a fixed arrival time always costs at least that much.
30 min: target at 106.3 deg, dep 1.006 + arr 1.797 = 2.803 km/s 45 min: target at 144.4 deg, dep 0.523 + arr 0.739 = 1.262 km/s 60 min: target at 182.5 deg, dep 0.459 + arr 0.430 = 0.889 km/s cheapest: 60 min with 0.889 km/s
Common mistakes
- Letting the solver pick the wrong direction. Decide prograde or retrograde first; it sets whether Δθ is the short or long angle.
- Feeding a transfer angle of exactly 180°. The plane is undefined; offset one point slightly or treat the Hohmann case separately.
- Forgetting that the target moves. In an intercept the arrival position is where the target will be, not where it is now.
- Trusting a solver without a test. Propagate v₁ forward and check you land on r₂.
Check yourself
Recap from earlier chapters
(Chapter 4) What does Kepler's equation relate?
(Chapter 9) How long does a Hohmann transfer from 300 km to GEO take?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R12poliastro project, tests/test_iod.py::test_curtis52 (Curtis, Orbital Mechanics for Engineering Students, Example 5.2), commit 21fd771. Secondary record of a textbook case.
- R13Izzo, D., "Revisiting Lambert's problem", Celestial Mechanics and Dynamical Astronomy 121 (2015) 1-15, doi:10.1007/s10569-014-9587-y; open preprint arXiv:1403.2705.
Where in each source, and links: References.
Porkchop plots
A porkchop plot is a map of launch dates against flight times (or arrival dates), coloured by how much energy the trip needs. Each point is one Lambert problem: where Earth is on the departure day, where the target planet is on the arrival day, and the flight time between. The cheap region forms a closed valley that, drawn as contours, looks like a pork chop. Everything here uses a SIMPLIFIED model: circular, coplanar planetary orbits and an arbitrary model day 0 when both planets are lined up. It shows the idea, not real dates.
You will be able to
- Explain what a porkchop plot shows and how it is built from Lambert solves
- Compute synodic periods and Hohmann phase angles
- Read the cost of a late launch from the plot
What is plotted
On each grid point solve Lambert's problem around the Sun. Subtract Earth's velocity from v₁ to get the hyperbolic excess speed at departure; C3 = v∞² is the launch-vehicle energy, in km²/s². Subtracting the planet's velocity from v₂ gives the arrival v∞, which sets the capture burn.
Why windows recur
The geometry repeats every synodic period, the time for Earth to lap the target. For Mars that is about 26 months; for Venus about 19. Between windows the planets are badly placed and C3 is huge.
A simplified model, clearly labelled
Real porkchops use planetary ephemerides with eccentric, inclined orbits. Inclination turns the 180° transfer line into an expensive ridge and splits the valley in two (type I and type II transfers). The circular, coplanar model keeps the physics and removes the detail.[R4]
Key equations
where
where
where
where
Worked examples
The Hohmann benchmark to Mars[R2, R4]
Given
In the simplified model (circular, coplanar orbits at 1.000000 and 1.523679 AU), find the Hohmann flight time, the hyperbolic excess speed and C3 at Earth, and where Mars must be at departure.
new symbols
Working
- at = (rE + rM)/2 = 1.8877e+08 km; tH = π√(a_t³/μ_Sun)258.9 days
- v∞ at Earth = √(μ_Sun(2/rE − 1/at)) − √(μ_Sun/rE) = 2.945 km/s, so C3 = v∞²8.67 km²/s²
- v∞ at Mars2.649 km/s
- Mars moves nM tH = 135.7° during the flight, so at departure it must lead Earth by φ = 180° − that44.3°
Answer. tH = 259 days, C3 = 8.67 km²/s², v∞ = 2.94 km/s at Earth and 2.65 km/s at Mars, Mars 44.3° ahead at departure.
Sanity check. This lead angle recurs once per synodic period, which is why Mars windows come about every 26 months.
a_t = 1.8877e+08 km, t_H = 258.9 days v_inf at Earth = 2.945 km/s, C3 = 8.671 km^2/s^2 v_inf at Mars = 2.649 km/s Mars moves 135.7 deg, so it must lead Earth by 44.34 deg
Reading the porkchop
Given
Solve Lambert's problem on a grid of departure days (every 5 days around the window) and flight times (120 to 420 days), compute C3 at each point and find the minimum.
Working
- 2914 grid points solved (points within 2° of a 180° transfer are skipped)
- Grid minimum (5-day steps): day 684, 255 days, C38.68 km²/s²
- Refined by a simplex search: day 683.8, flight 255.0 days, C38.68 km²/s²
- Hohmann benchmark: C3 = 8.67 km²/s² at 259 days
Answer. Best departure is model day 684 with a 255 day flight and C3 = 8.68 km²/s², 0.1 % above Hohmann.
Sanity check. In this coplanar model the true optimum is the Hohmann transfer itself (a 180° arc), which the grid must skip; the best usable point sits just beside it. With real inclined orbits the 180° line becomes an expensive ridge and the optimum splits into a short-way and a long-way valley.
2914 grid points; grid minimum day 684, 255 days, C3 = 8.68 km^2/s^2 refined: day 683.8, flight 255.0 days, C3 = 8.676 km^2/s^2 Hohmann: C3 = 8.671 km^2/s^2 at 259 days; the optimum is 0.1 % above it
Practice
How often do windows open?[R4]
Given
Using the simplified circular orbits, find the orbital periods of Earth (1.000000 AU) and Mars (1.523679 AU) and the Earth-Mars synodic period.
T = 2π√(a³/μ_Sun); 1/Tsyn = |1/TE − 1/TM|.
Answer. Tsyn = 780 days, about 26 months between Mars windows.
Working
- TE = 365.26 days, TM686.97 days
- 1/Tsyn = |1/365.26 − 1/686.97|0.001282 per day
- Tsyn = 779.9 days25.6 months
Answer. Tsyn = 780 days, about 26 months between Mars windows.
Sanity check. Close to two Earth years, because Mars takes nearly two years to go round: Earth must lap it.
T_E = 365.26 days, T_M = 686.97 days synodic period = 779.9 days = 25.6 months
Lining up for Venus[R4]
Given
For a Hohmann transfer from Earth inward to Venus (0.723332 AU) in the simplified model, find the flight time, where Venus must be relative to Earth at departure, and how often that geometry recurs.
Same logic as Mars: φ = π − nV tH. A negative angle means Venus must be behind Earth.
Answer. tH = 146 days, Venus 54.0° behind Earth at departure, windows every 584 days.
Working
- at = (rE + rV)/2 = 1.2890e+08 km; tH = π√(a_t³/μ_Sun)146.1 days
- nV tH = 234.0°, so φ = 180° − 234.0°-54.0°
- Synodic period Earth-Venus584 days
Answer. tH = 146 days, Venus 54.0° behind Earth at departure, windows every 584 days.
Sanity check. Venus is faster than Earth, so it must start behind and catch up to the arrival point: a negative lead angle is the expected sign.
t_H = 146.1 days, Venus moves 234.0 deg phi = -54.03 deg: Venus 54.03 deg behind Earth at departure windows every 584 days
The cost of a late launch
Given
Take the optimum found in this chapter (model day 683.8, 255 day flight, C3 = 8.68 km²/s²). The launch slips 30 days and the flight time is kept. Solve the new Lambert problem and find the new C3 and the percentage increase.
Departure day d + 30 fixes Earth's position; arrival day d + 30 + t fixes Mars's. Solve Lambert with μ_Sun and subtract Earth's velocity.
Answer. C3 rises to 14.4 km²/s², 66 % more than the optimum.
Working
- Earth at model day 714, Mars at model day 969; transfer angle 164.2°
- Lambert gives v∞ = 3.790 km/s, so C314.37 km²/s²
- Increase = 14.37/8.68 − 166 %
Answer. C3 rises to 14.4 km²/s², 66 % more than the optimum.
Sanity check. A month's slip moves you off the valley floor of the porkchop. Launch vehicles quote payload against C3, so this directly costs payload mass; it is why planetary windows are only a few weeks wide.
refined optimum: day 683.8, flight 255.0 days optimum C3 = 8.68 km^2/s^2 late: Earth day 714, Mars day 969, transfer angle 164.2 deg v_inf = 3.791 km/s, C3 = 14.37 km^2/s^2, 65.7 % more
Common mistakes
- Reading the simplified plot as real dates. Model day 0 is an arbitrary reference; real windows need ephemerides.
- Quoting C3 in km/s. C3 is a speed squared: km²/s².
- Forgetting arrival. The cheapest launch can bring a large arrival v∞ and an expensive capture.
- Solving exactly on the 180° line. It is singular for Lambert; skip it or offset slightly.
Check yourself
Recap from earlier chapters
(Chapter 12) What are the inputs and outputs of Lambert's problem?
(Chapter 9) Why is a Hohmann transfer only a benchmark?
Sources for this chapter
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R4Standish, E.M. and Williams, J.G., "Approximate Positions of the Planets", NASA JPL SSD.
- R7IAU 2012 Resolution B2, "Re-definition of the astronomical unit of length", XXVIII General Assembly, Beijing, Aug 2012 (text hosted by SYRTE, Observatoire de Paris).
Where in each source, and links: References.
Patched conics, sphere of influence and the Oberth effect
An interplanetary trip has three bodies pulling at once, but almost everywhere one of them dominates. Patched conics uses that: near Earth, the craft follows a hyperbola around Earth; far from planets, an ellipse around the Sun; near Mars, a hyperbola around Mars. The pieces are glued at the sphere of influence by matching velocities. The link between them is the hyperbolic excess speed v∞: what is left of the speed after climbing out of a planet's gravity well. Burning deep in the well, where you are fast, buys more v∞ than burning far away: the Oberth effect.
You will be able to
- Split an interplanetary trip into planet-centred and Sun-centred conics
- Size a departure or capture burn from a hyperbolic excess speed
- Explain and compute the Oberth effect
Sphere of influence
Inside a radius r = a(m/M)2/5 the planet's pull dominates the Sun's disturbance of the motion around it. For Earth that is under a million kilometres; outside it, the Sun takes over. It is an approximation for patching, not a physical wall.[R2, R4]
Hyperbolic excess speed
A craft leaving a planet on a hyperbola keeps a speed v∞ at large distance; its energy v²/2 − μ/r equals v∞²/2. In the Sun's frame its velocity is the planet's velocity plus v∞. C3 = v∞² is the number launch vehicles quote.
Departure and capture
The burn happens at periapsis of the planet-centred hyperbola. The required speed is √(v∞² + 2μ/rp); subtract the parking-orbit speed for departure, or the target-orbit speed for capture.
The Oberth effect
Energy gained by a burn is v·Δv (to first order), so the same Δv adds more energy where v is large. Planetary departures burn at the lowest point of the parking orbit for exactly this reason.
Key equations
where
where
where
where
where
where
Worked examples
Leaving Earth for Mars[R2, R3]
Given
A spacecraft in a 300 km circular parking orbit must leave Earth with the Hohmann excess speed to Mars, v∞ = 2.945 km/s (Chapter 13). Find C3, the speed needed at periapsis of the departure hyperbola and the burn.
new symbols
Working
- C3 = v∞²8.671 km²/s²
- Periapsis radius rp = 6378.137 + 3006678.137 km
- vp = √(v∞² + 2μ/rp) = √(8.671 + 119.375)11.3157 km/s
- Circular speed vc = √(μ/rp)7.7258 km/s
- Δv = vp − vc = 3.5900 km/s; hyperbola eccentricity e = 1 + rp v∞²/μ1.1453
Answer. C3 = 8.67 km²/s², vp = 11.316 km/s, burn Δv = 3.590 km/s.
Sanity check. Escape alone from 300 km costs 3.200 km/s; adding a 2.94 km/s excess costs only 0.390 km/s more. That bargain is the Oberth effect.
v_inf = 2.945 km/s, C3 = 8.671 km^2/s^2 v_p = 11.3157 km/s, v_c = 7.7258 km/s, burn = 3.5900 km/s, e = 1.1453
The Oberth effect in numbers
Given
A spacecraft is on a parabolic (just-escaping) trajectory. It has 1.0 km/s of Δv to spend. Compare burning at periapsis, 300 km up, with burning far from Earth where its speed has dropped to almost zero.
Working
- At periapsis the speed is the escape speed √(2μ/r)10.9259 km/s
- After the burn: 11.9259 km/s; far away v∞ = √(v² − v_esc²)4.7804 km/s
- Far away the speed is about 0, so the same burn gives v∞1.0000 km/s
- Gain factor4.78
Answer. Burning at periapsis gives v∞ = 4.78 km/s; burning far away gives 1.00 km/s: 4.8 times as much.
Sanity check. Kinetic energy grows with v², so a fixed Δv adds the most energy where you are already fastest: deep in the gravity well, at periapsis.
periapsis speed 10.9259 km/s, after the burn 11.9259 km/s burn at periapsis: v_inf = 4.7804 km/s; far away: v_inf = 1.0000 km/s gain factor 4.78
Practice
How big is Mars's neighbourhood?[R2, R4]
Given
Find the radius of the sphere of influence of Mars (orbit radius 1.523679 AU, μ = 42828 km³/s²) and of Earth, using r = a(μ_p/μ_Sun)2/5.
Convert AU to km first (1 AU = 149,597,870.7 km).
Answer. Mars SOI ≈ 577,225 km (170 Mars radii); Earth SOI ≈ 924,647 km.
Working
- Mars: a = 2.2794e+08 km; (μ_M/μ_Sun)0.4 = 2.5324e-03; r5.7723e+05 km
- Earth: a = 1.4960e+08 km; r9.2465e+05 km
- Mars SOI in planet radii: 170 RM
Answer. Mars SOI ≈ 577,225 km (170 Mars radii); Earth SOI ≈ 924,647 km.
Sanity check. Earth's SOI is about 2.4 times the Moon's distance, which is why the Moon is firmly Earth's.
Mars SOI = 577225.3 km = 170 Mars radii Earth SOI = 924646.8 km
Capture at Mars[R2, R3]
Given
A spacecraft arrives at Mars with v∞ = 2.6489 km/s (Hohmann arrival, from Chapter 13) and periapsis 400 km. Find the capture burn into (a) a 400 km circular orbit and (b) an ellipse 400 × 20,000 km. Mars radius 3396.2 km, μ = 42828 km³/s².
Speed at periapsis of the arrival hyperbola is √(v∞² + 2μ/rp). Subtract the speed of the orbit you want at that point.
Answer. (a) 2.08 km/s into the circular orbit; (b) 1.03 km/s into the ellipse.
Working
- rp = 3796.2 km; vp = √(7.017 + 22.564)5.4388 km/s
- (a) circular speed √(μ/rp) = 3.3588 km/s, Δv2.0799 km/s
- (b) ellipse a = (3796.2 + 23396.2)/2 = 13596.2 km; periapsis speed 4.4061 km/s; Δv1.0327 km/s
Answer. (a) 2.08 km/s into the circular orbit; (b) 1.03 km/s into the ellipse.
Sanity check. Capturing into a long ellipse saves 1.05 km/s; many missions capture loosely and lower the orbit later, often by aerobraking.
v_inf = 2.649 km/s, periapsis speed 5.4388 km/s (a) circular capture 2.0799 km/s; (b) ellipse capture 1.0327 km/s
The best height to capture[R2, R3]
Given
For capture into a circular orbit at Mars with v∞ = 2.6489 km/s (Hohmann arrival, from Chapter 13), the burn Δv(r) = √(v∞² + 2μ/r) − √(μ/r) depends on the capture radius r. Find the radius that minimises it, the altitude, the minimum Δv, and compare with the 400 km case.
Differentiate Δv(r) with respect to r and set it to zero; the algebra gives r = 2μ/v∞².
Answer. ropt = 12,208 km (altitude 8,811 km), Δvmin = 1.873 km/s, 207 m/s less than at 400 km.
Working
- dΔv/dr = 0 gives ropt = 2μ/v∞² = 2 × 42828 / 7.016612208 km
- Altitude = 12208 − 3396.28811 km
- Δvmin = v∞/√21.8731 km/s
- At 400 km it was 2.0799 km/s; a numerical scan confirms the minimum at 12208 km
Answer. ropt = 12,208 km (altitude 8,811 km), Δvmin = 1.873 km/s, 207 m/s less than at 400 km.
Sanity check. Here the Oberth effect and the cost of a fast circular orbit pull in opposite directions; the optimum is a balance, and it is shallow, which is why mission designers rarely capture that high.
r_opt = 12207.5 km (scan: 12208 km), altitude 8811.3 km dv_min = 1.8731 km/s (= v_inf/sqrt 2 = 1.8731) at 400 km: 2.0799 km/s, so 207 m/s more
Common mistakes
- Adding v∞ to the escape speed. Speeds add in squares: vp = √(v∞² + v_esc²).
- Using the Sun's μ inside the sphere of influence, or the planet's outside it.
- Treating the sphere of influence as a hard boundary. It is a patching approximation.
- Burning far from the planet to save propellant. The Oberth effect says the opposite.
Check yourself
Recap from earlier chapters
(Chapter 3) What is the escape speed at 300 km?
(Chapter 13) What does C3 measure?
Sources for this chapter
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R3NASA JPL Solar System Dynamics, "Planetary Physical Parameters" (radii from Archinal et al. 2018, CMDA 130:22; change log 2019-12-12).
- R4Standish, E.M. and Williams, J.G., "Approximate Positions of the Planets", NASA JPL SSD.
Where in each source, and links: References.
Gravity assists
A gravity assist looks like free energy, and in the Sun's frame it is. Relative to the planet, the craft comes in and goes out at the same speed; the planet only turns its direction. But the planet is itself moving fast around the Sun. Add the planet's velocity to the turned v∞ and the craft's speed around the Sun has changed, up or down depending on which side it passed. The energy comes from the planet, which slows by an immeasurably small amount.
You will be able to
- Compute the turn angle of a flyby from v∞ and the periapsis
- Build the velocity triangle and find the speed change around the Sun
- Explain why passing behind a planet adds energy and in front removes it
In the planet's frame
The flyby is a hyperbola with eccentricity e = 1 + rp v∞²/μ. The incoming and outgoing asymptotes are separated by the turn angle δ = 2 asin(1/e). A closer pass or a slower approach turns more.
In the Sun's frame
Velocities add: before, v = Vp + v∞,in; after, v = Vp + v∞,out. Since both v∞ vectors have the same length, the outgoing tip lies on a circle centred on the tip of Vp. Where on the circle it lands is set by δ and the side of the pass.
Leading and trailing passes
Passing behind the planet (its trailing side) turns v∞ towards the planet's direction of motion: the craft speeds up. Passing in front (leading side) does the reverse: the craft slows, which is how missions fall inwards to Mercury or the Sun.
Key equations
where
where
where
where
where
Worked examples
A Jupiter flyby, from both sides[R3, R8]
Given
A probe on a solar orbit with perihelion 1 AU and aphelion 6 AU meets Jupiter (circular orbit, 5.2026 AU) on its way out. It passes at 20 Jupiter radii. Find v∞, the turn angle, and the probe's speed around the Sun after passing behind Jupiter and in front of it.
new symbols
Working
- Probe at r = 7.7830e+08 km: v = 9.358 km/s, split into radial 5.602 and along-track 7.496 km/s; Jupiter moves at 13.058 km/s
- v∞,in = v − VJ = (5.602, -5.562) km/s, |v∞|7.894 km/s
- e = 1 + rp v∞²/μ = 1 + 1429840 × 62.322 / 1.2669e+08 = 1.7034; δ = 2 asin(1/e)71.9°
- Behind Jupiter (trailing side) v∞ turns towards Jupiter's motion: speed after18.077 km/s
- In front (leading side) it turns away: speed after6.974 km/s
- Check: |v∞,out| = 7.894 km/s, unchanged
Answer. v∞ = 7.89 km/s turned by 72°: the Sun-frame speed goes from 9.36 to 18.08 km/s behind Jupiter, or down to 6.97 km/s in front.
Sanity check. Relative to Jupiter nothing is gained: the hyperbola is symmetric. The gain comes from adding a rotated v∞ to Jupiter's own large velocity; Jupiter loses a negligible amount of orbital energy in return.
behind Jupiter: Sun-frame speed after = 18.077 km/s in front of Jupiter: Sun-frame speed after = 6.974 km/s before: v = 9.358 km/s (radial 5.602, along-track 7.496), Jupiter 13.058 km/s v_inf = 7.894 km/s, e = 1.7034, turn = 71.9 deg
Practice
How far does Venus bend you?[R2, R3]
Given
A probe passes Venus (R = 6051.8 km, μ = 324859 km³/s²) with v∞ = 5.0 km/s at a periapsis altitude of 500 km. Find the eccentricity of the flyby hyperbola and the turn angle.
e = 1 + rp v∞²/μ with rp measured from the planet's centre; δ = 2 asin(1/e).
Answer. e = 1.504, turn angle δ = 83.3°.
Working
- rp = 6051.8 + 5006551.8 km
- e = 1 + 6551.8 × 25.00 / 3248591.5042
- δ = 2 asin(1/1.5042)83.33°
Answer. e = 1.504, turn angle δ = 83.3°.
Sanity check. e > 1 confirms a hyperbola; slower or closer passes turn more.
r_p = 6551.8 km, e = 1.5042, turn angle = 83.33 deg
The free Δv of a flyby
Given
For the pass of O15.1, find the size of the velocity change the flyby delivers, and the largest possible value at Venus for the same v∞ (a grazing pass, rp = R).
The velocity change is the vector difference of v∞,out and v∞,in, both of length v∞: Δv = 2 v∞ sin(δ/2).
Answer. Δv = 6.65 km/s at 500 km, up to 6.82 km/s for a grazing pass.
Working
- Δv = 2 × 5.0 × sin(41.67°)6.648 km/s
- Grazing: e = 1 + 6051.8 × 25.00 / 324859 = 1.4657, δ86.04°
- Δvmax = 2 × 5.0 × sin(43.02°)6.823 km/s
Answer. Δv = 6.65 km/s at 500 km, up to 6.82 km/s for a grazing pass.
Sanity check. Several km/s for free, compared with chemical stages that cost tonnes of propellant per km/s: the reason Cassini, Galileo and Messenger all used Venus.
at 500 km: e = 1.5042, turn 83.33 deg, dv = 6.648 km/s grazing: e = 1.4657, turn 86.04 deg, dv = 6.823 km/s
Thrown out of the Solar System?[R3, R8]
Given
The probe of the worked example instead passes behind Jupiter at 4 Jupiter radii. Find the turn angle, its speed around the Sun afterwards, and whether it now escapes the Sun.
Rotate v∞,in by δ towards Jupiter's direction of motion, add Jupiter's velocity, then compare with the solar escape speed √(2μ_Sun/r) at Jupiter's distance.
Answer. δ = 122°, speed 20.84 km/s against an escape speed of 18.47 km/s: the probe escapes the Sun.
Working
- |v∞| = 7.894 km/s; rp = 285968 km; e = 1.1407; δ122.5°
- v∞,out = (1.683, 7.713) km/s; vout = VJ + v∞,out = (1.683, 20.771) km/s, |v|20.839 km/s
- Solar escape speed at 5.2026 AU: √(2μ_Sun/r)18.467 km/s
- Heliocentric energy afterwards: v²/2 − μ_Sun/r46.62 km²/s²
Answer. δ = 122°, speed 20.84 km/s against an escape speed of 18.47 km/s: the probe escapes the Sun.
Sanity check. A probe that arrived on a bound ellipse can leave on a hyperbola: Jupiter supplied the energy. This is how the outer-planet missions reached escape.
v_inf = 7.894 km/s, e = 1.1407, turn = 122.5 deg speed after = 20.839 km/s, solar escape speed = 18.467 km/s energy = 46.62 km^2/s^2: escapes the Sun
Common mistakes
- Expecting a speed change relative to the planet. |v∞| is the same before and after.
- Measuring rp from the surface. The flyby formulas use the distance from the planet's centre.
- Forgetting the side of the pass. The same δ can add or remove energy.
- Adding speeds instead of vectors. Build the triangle Vp + v∞.
Check yourself
Recap from earlier chapters
(Chapter 14) What is v∞?
(Chapter 2) Which eccentricities give a hyperbola?
Sources for this chapter
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R3NASA JPL Solar System Dynamics, "Planetary Physical Parameters" (radii from Archinal et al. 2018, CMDA 130:22; change log 2019-12-12).
- R4Standish, E.M. and Williams, J.G., "Approximate Positions of the Planets", NASA JPL SSD.
- R8Mamajek, E.E., Prša, A., Torres, G. et al., "IAU 2015 Resolution B3 on Recommended Nominal Conversion Constants for Selected Solar and Planetary Properties", arXiv:1510.07674 (2015).
Where in each source, and links: References.
Lagrange points
Put a third, tiny body near two big ones that circle each other, and look from a frame that turns with them. There are five places where gravity from both and the outward pull of the turning frame cancel. A spacecraft parked there keeps its place relative to the pair. Three points sit on the line through the two bodies, two lead and trail the smaller body by 60 degrees.
You will be able to
- Explain what a Lagrange point is in the rotating frame
- Find the L1 and L2 distances exactly and with the Hill approximation
- State which points are stable and why L1/L2 missions need station-keeping
The rotating frame
In the circular restricted three-body problem the two primaries move on circles about their common centre of mass, and the third body is too light to disturb them. In a frame rotating at their orbital rate the primaries stand still, and an equilibrium is a point where gravity and the centrifugal term balance.
Where the points are
L1 lies between the bodies, L2 just beyond the smaller one, L3 on the far side of the larger one. For a small mass ratio, L1 and L2 sit at about a(μ*/3)^(1/3) from the smaller body. L4 and L5 form equilateral triangles with the two primaries.[R22, R24]
Stable and unstable
The three collinear points are unstable: a small error grows within weeks. L4 and L5 are stable when the mass ratio is below about 0.0385, which holds for Sun-Earth, Earth-Moon and Sun-Jupiter (the Trojan asteroids).[R23]
Why missions go there
Sun-Earth L2 keeps the Sun, Earth and Moon behind the spacecraft and gives a steady thermal environment, ideal for infrared telescopes. Real missions fly large halo or Lissajous orbits around the point and trim them with small burns.
Living near L1 and L2
Nothing parks exactly on L1 or L2. A spacecraft nudged off the point drifts away on a timescale of about 23 days, so missions fly large loops around the point instead, called halo or Lissajous orbits, and trim them with small periodic burns. Webb's loop around Sun-Earth L2 is comparable in size to the Moon's orbit and takes about 6 months. Solar observatories fly similar loops around L1.[R23, R24]
Key equations
where
where
three roots: L1, L2, L3
where
good when μ* ≪ 1
where
L4 and L5 are linearly stable below it
Worked examples
How far away is Sun-Earth L2?[R22, R24]
Given
Treat the Sun and Earth as a circular restricted three-body system with separation a = 149,597,871 km and mass ratio μ* = μ_Earth / (μ_Sun + μ_Earth). Find the distances of L1 and L2 from Earth, exactly (root of the collinear balance) and with the Hill approximation.
new symbols
Working
- μ* = 398600.4418 / (1.327124e+11 + 398600.4418)3.0035e-06
- Hill: r ≈ a (μ*/3)^(1/3) = 149,597,871 × (3.0035e-06/3)^(1/3)1,496,557 km
- Exact root of the collinear balance: L1 is 0.009970 a = 1,491,551 km from Earth
- L2 is 0.010037 a = 1,501,532 km from Earth (on the night side)
- Hill sits between them: L1 is 0.33 % closer, L2 0.33 % farther
Answer. L1 ≈ 1.492 million km and L2 ≈ 1.502 million km from Earth; the Hill estimate 1.497 million km is within 0.3 %.
Sanity check. Both are about 3.9 times the Moon's distance, still deep inside the Sun's pull, which is why the points exist at all: there the combined pull matches the rotating frame.
mu* = 3.0035e-06 L1: 0.009970 a = 1,491,551 km = 1.492 million km L2: 0.010037 a = 1,501,532 km = 1.502 million km Hill: 1,496,557 km = 1.497 million km, within 0.33 %
Earth-Moon L1, where Hill starts to slip[R2, R6]
Given
Repeat for the Earth-Moon system: a = 384,400 km, μ_Moon = 4902.8 km³/s². Where is L1, measured from the Moon and from Earth, and how good is the Hill estimate here?
new symbols
Working
- μ* = 4902.8 / (398600.4418 + 4902.8)0.01215
- Exact L1: 0.15093 a = 58,019 km from the Moon, so 326,381 km from Earth
- Hill: 384,400 × (0.01215/3)^(1/3)61,274 km
- Error of Hill5.6 %(the mass ratio is no longer tiny)
Answer. Earth-Moon L1 is 58,019 km from the Moon (326,381 km from Earth); Hill overestimates it by 5.6 %.
Sanity check. μ* = 0.012 is about 4,046 times the Sun-Earth value, so the first-order Hill formula degrades from 0.3 % to several per cent.
mu* = 0.01215 L1 is 58,019 km from the Moon, 326,381 km from Earth Hill: 61,274 km, an overestimate of 5.6 %
Practice
Sun-Mars L1 with one line
Given
Mars orbits the Sun at a = 1.523679 AU (1 AU = 149,597,870.7 km) with μ_Mars = 42,828 km³/s². Estimate the distance of Sun-Mars L1 from Mars with the Hill approximation.
Compute the mass ratio first; it is tiny, so r ≈ a(μ*/3)^(1/3) is accurate.
Answer. r ≈ 1,084,053 km (the exact root gives 1,082,332 km).
Working
- a = 1.523679 × 149,597,870.7227,939,134 km
- μ* = 42,828 / 1.327125e+113.2271e-07
- r ≈ 227,939,134 × (3.2271e-07/3)^(1/3)1,084,053 km
Answer. r ≈ 1,084,053 km (the exact root gives 1,082,332 km).
Sanity check. Mars is about a tenth of Earth's mass but farther out, so its L1 (1.08 million km) is a bit under Earth's 1.49 million km.
a = 227,939,134 km, mass ratio = 3.2271e-07 Hill: r = 1,084,053 km (exact root 1,082,332 km)
Earth-Moon L2 for a relay[R2, R6]
Given
A relay for the lunar far side would sit near Earth-Moon L2. Using a = 384,400 km and μ* = 0.01215, find L2's distance beyond the Moon (exact) and the percentage error of the Hill estimate 61,274 km.
L2 lies beyond the smaller body; solve the collinear balance between x = 1 − μ* and 1.5 (the library does this by bracketing).
Answer. L2 is 64,515 km beyond the Moon; Hill is off by -5.0 %.
Working
- Exact root: xL2 = 1.155682, Moon at 1 − μ*0.987849
- d = (1.155682 − 0.987849) × 384,40064,515 km
- Hill error = (61,274 − 64,515) / 64,515-5.0 %
Answer. L2 is 64,515 km beyond the Moon; Hill is off by -5.0 %.
Sanity check. L2 is farther from the Moon than L1 (58,019 km): on the far side the Moon's pull adds to Earth's, so balance needs a larger distance.
x_L2 = 1.155682, Moon at 0.987849 L2 is 64,515 km beyond the Moon Hill 61,274 km is off by -5.02 %
Why L2 needs station-keeping[R23]
Given
Linearise the rotating-frame motion about Sun-Earth L2. With c₂ = μ*/|x − 1 + μ*|³ + (1 − μ*)/|x + μ*|³, the in-plane unstable root is λ² = [c₂ − 2 + √(9c₂² − 8c₂)]/2, in units of the orbital rate n = 2π per year. Find c₂, λ and the e-folding time of a small error, in days.
Use xL2 from the exact collinear root; the time unit is 1/n with n = 2π / (one year in seconds).
Answer. λ = 2.484 n, so a small error grows by a factor e every 23 days.
Working
- xL2 = 1.0100341, μ*3.0035e-06
- c₂3.9408
- λ² = (3.9408 − 2 + √(9 × 3.9408² − 8 × 3.9408))/2 = 6.1723, λ2.4844
- n = 2π / 31556926 s1.9911e-07 rad/s
- e-folding time = 1/(λ n) = 2.022e+06 s23.4 days
Answer. λ = 2.484 n, so a small error grows by a factor e every 23 days.
Sanity check. In 70 days an error grows about 20 times: missions at L2 trim their orbit every few weeks with small burns, a few m/s per year.
x_L2 = 1.0100341, c2 = 3.9408, lambda = 2.4844 n e-folding time = 2.022e+06 s = 23.4 days
Common mistakes
- Measuring L1/L2 from the larger body: the Hill formula gives the distance from the smaller one.
- Thinking a spacecraft sits exactly at L2. It orbits around it, often hundreds of thousands of km wide.
- Assuming every Lagrange point is stable. Only L4 and L5 are, and only below the Routh mass ratio.
- Using the Hill formula for Earth-Moon and expecting sub-per-cent accuracy; the mass ratio is too large.
Check yourself
Recap from earlier chapters
What is the circular speed at Earth's distance from the Sun?
Why does a period depend only on a, not on e?
Sources for this chapter
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R6Williams, D.R., "Moon Fact Sheet", NASA NSSDCA, updated 11 Jan 2024.
- R22ESA, "Types of orbits", 30 Mar 2020.
- R23Cornish, N.J., "What is a Lagrange Point?", NASA Science, 27 Mar 2018, updated 1 Nov 2024.
- R24NASA Science, "Webb Orbit".
Where in each source, and links: References.
GEO and station-keeping
A geostationary satellite circles Earth exactly as fast as Earth turns, so it hangs over one longitude. That only works on one circle, 42,164 km from Earth's centre, in the equator plane. Nature keeps pushing it off: the Sun and Moon tilt the orbit, Earth's slightly oval equator drags it in longitude. Station-keeping is the small, regular set of burns that pushes it back.
You will be able to
- Derive the GEO radius from the sidereal day
- Turn a semi-major axis error into a longitude drift rate
- Budget north-south, east-west and end-of-life Δv for a GEO mission
The one orbit
Match the period to the sidereal day, 23 h 56 min 4 s, not the 24 h solar day. Zero eccentricity and zero inclination complete the conditions.[R1, R2, R22]
Drift in longitude
A slightly higher orbit has a longer period, so the satellite falls behind Earth and drifts west; lower drifts east. The drift rate is −1.5 n Δa/a, about 0.013 deg/day per km.
North-south: the big one
The Sun's and Moon's pull tilt the orbit by about a degree per year. Removing it costs about 50 m/s per year, by far the largest station-keeping item.[R18, R19]
East-west and end of life
Earth's equator is slightly elliptical, so most longitudes feel a small along-track force; correcting it costs one or two m/s per year. At end of life the satellite is raised a few hundred km into a graveyard orbit to free the slot.[R18, R20]
Key equations
where
where
where
where
from Δa = 2aΔv/v
where
Worked examples
Why GEO is 42,164 km and not 42,241 km[R1, R2, R11]
Given
A geostationary satellite must turn once per rotation of Earth relative to the stars: the sidereal day, 86164.0905 s. Find the GEO radius and speed, and the radius you would get by mistake from the 86,400 s solar day.
new symbols
Working
- a = (μ (T/2π)²)^(1/3) = (398600.4418 × (86164.0905/2π)²)^(1/3)42,164.17 km
- Altitude = 42,164.17 − 6378.137 = 35,786.0 km; speed v = √(μ/a)3.0747 km/s
- With T = 86400 s: a = 42,241.10 km, 76.9 km too high
- That error drifts at -0.988 deg/day
Answer. RGEO = 42,164.2 km (altitude 35,786 km), v = 3.075 km/s; using the solar day would put it 77 km too high.
Sanity check. The error drifts the satellite about one degree per day westward, roughly a full circle per year: the solar day already contains one extra turn per year.
R_GEO = 42,164.17 km, altitude 35,786.0 km, v = 3.0747 km/s with the solar day: a = 42,241.10 km, 76.9 km too high that error drifts at -0.988 deg/day
Fifteen years of north-south station-keeping[R17, R18]
Given
Without control, the Sun and Moon tilt a GEO orbit by about 0.85 deg per year. Find the yearly Δv to remove it, the 15-year total including 1.5 m/s/yr east-west, and the propellant for a 3000 kg satellite with Isp = 315 s.
new symbols
Working
- ΔvNS per year = 2 v sin(Δi/2) = 2 × 3074.7 × sin(0.425 deg)45.61 m/s
- Per year with east-west: 45.61 + 1.547.11 m/s
- 15 years: 707 m/s
- Mass ratio = exp(Δv / (Isp g₀)) = exp(707 / (315 × 9.80665))1.2571
- Propellant = m₀ (1 − 1/MR) = 3000 × (1 − 1/1.2571)613 kg
Answer. About 46 m/s per year north-south, 707 m/s over 15 years, needing 613 kg of propellant.
Sanity check. North-south dominates (it is 30 times the east-west budget): that is why electric propulsion, with Isp ten times higher, is now common for station-keeping.
north-south 45.61 m/s per year, 47.11 with east-west 15 years: 706.7 m/s, mass ratio 1.2571 propellant 613.5 kg
Practice
A satellite 2.5 km too high
Given
After a manoeuvre a GEO satellite has a semi-major axis 2.5 km above RGEO = 42,164.2 km. Find its longitude drift rate in deg/day and how long it takes to leave a ±0.05 deg box if it starts in the centre.
λ̇ = −1.5 n Δa / a with n = 2π / sidereal day; positive Δa means slower, so it drifts west.
Answer. λ̇ = -0.0321 deg/day westward; it leaves the box after 1.6 days.
Working
- n = 2π / 86164.09057.292116e-05 rad/s
- λ̇ = −1.5 × 7.292116e-05 × 2.5 / 42,164.2-6.4855e-09 rad/s
- = -0.0321 deg/day (west)
- Time to move 0.05 deg = 0.05 / 0.03211.56 days
Answer. λ̇ = -0.0321 deg/day westward; it leaves the box after 1.6 days.
Sanity check. A few km of error leaves a tight box in days, which is why GEO operators trim the orbit every one to two weeks.
drift = -6.4855e-09 rad/s = -0.0321 deg/day (west) leaves the box after 1.56 days
Retiring to the graveyard[R20]
Given
At end of life a 1900 kg GEO satellite is raised 300 km with a two-burn Hohmann transfer. Find each burn, the total Δv and the propellant (Isp = 315 s).
The two radii are almost equal, so both burns are small and nearly equal; use the Hohmann formulas exactly anyway.
Answer. Δv ≈ 10.9 m/s in two near-equal burns, about 6.7 kg of propellant.
Working
- r₁ = 42,164.2 km, r₂ = 42,464.2 km, at42,314.2 km
- Δv₁ = 5.445 m/s, Δv₂5.435 m/s
- Total10.88 m/s
- Propellant = 1900 × (1 − exp(−10.88 / (315 × 9.80665)))6.7 kg
Answer. Δv ≈ 10.9 m/s in two near-equal burns, about 6.7 kg of propellant.
Sanity check. Check: for a small raise Δv ≈ v Δr / (2r) = 10.94 m/s, matching. Cheap, but it must be reserved: an empty tank cannot retire.
dv1 = 5.445 m/s, dv2 = 5.435 m/s, total 10.88 m/s propellant 6.68 kg
Moving a satellite 12 degrees east
Given
An operator moves a GEO satellite 12 deg east in 20 days. A small tangential burn Δv changes the semi-major axis by Δa = 2aΔv/v, so λ̇ = −3nΔv/v. Find the burn to start the drift, the semi-major axis of the drift orbit, and the total Δv (start plus stop).
East means faster, so the burn is retrograde (lower orbit). Required λ̇ = Δλ / t in rad/s.
Answer. Burn 1.70 m/s retrograde, drift orbit 46.7 km below GEO, total 3.4 m/s.
Working
- λ̇ = 12 deg / 20 days = 0.60 deg/day1.2120e-07 rad/s
- Δv = λ̇ v / (3n) = 1.2120e-07 × 3074.7 / (3 × 7.292116e-05)1.703 m/s(retrograde)
- Vis-viva after the burn: a = 1/(2/r − (v − Δv)²/μ) = 42,117.51 km, 46.66 km below GEO
- Stop burn equal and opposite: total = 2 × 1.7033.41 m/s
Answer. Burn 1.70 m/s retrograde, drift orbit 46.7 km below GEO, total 3.4 m/s.
Sanity check. Faster relocation costs proportionally more: half the time doubles the Δv. Check the sign: a lower orbit has a shorter period, so it moves ahead (east).
drift rate 1.2120e-07 rad/s, burn 1.703 m/s retrograde drift orbit a = 42,117.51 km, 46.66 km below GEO total with the stop burn: 3.41 m/s
Common mistakes
- Using 86,400 s for GEO: the radius comes out about 77 km high and the satellite drifts a degree a day.
- Getting the drift sign backwards. Higher and slower means it drifts west.
- Forgetting the end-of-life reserve; it is small but must exist when the tank is almost empty.
- Treating north-south and east-west as similar. North-south is roughly thirty times larger.
Check yourself
Recap from earlier chapters
Hohmann total for a small raise Δr at radius r is about?
Why change plane where the speed is lowest?
Write the rocket equation.
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R2NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).
- R11Vallado, D.A., companion software to Fundamentals of Astrodynamics and Applications (CelesTrak), GitHub repo CelesTrak/fundamentals-of-astrodynamics, commit 49df047 (AGPL-3.0; values read, no code copied).
- R17ArianeGroup GmbH, "400N bipropellant apogee motors" product page (S400-12, S400-15), © 2002-2024.
- R18Lovell, R.R. and O'Malley, T.A., Station Keeping of High Power Communication Satellites, NASA TM X-2136, Lewis Research Center, Nov 1970.
- R19Decoto, J. and Loerch, P. (Orbital ATK), "Technique for GEO RSO Station Keeping Characterization and Maneuver Detection", AMOS Technical Conference, 2015.
- R20Inter-Agency Space Debris Coordination Committee, IADC Space Debris Mitigation Guidelines, IADC-02-01 Rev. 2, Mar 2020.
- R22ESA, "Types of orbits", 30 Mar 2020.
Where in each source, and links: References.
Launch windows and launch azimuth
A rocket can only put a satellite into a plane that passes over the launch site. Launching due east gives the lowest inclination, equal to the site's latitude, and borrows Earth's rotation for free. Any other inclination means steering north- or south-east, and any particular plane only sweeps over the site twice a day: those are the launch windows.
You will be able to
- Find the launch azimuth for a target inclination, inertial and relative to the rotating ground
- Explain why the lowest reachable inclination equals the site latitude
- Time the two daily launch windows into a given orbit plane
Inclination and latitude
The orbit plane goes through Earth's centre and the site. Its inclination can never be less than the site latitude; it equals it when the rocket flies due east.
Azimuth
Spherical trigonometry gives sin β = cos i / cos φ. Close to i = φ the azimuth is 90 deg; for a polar orbit it is 0 or 180 deg; retrograde orbits need a westward component.
The rotating ground
The pad already moves east at ω_E RE cos φ, about 465 m/s at the equator. The rocket only adds the difference, which shifts the ground-relative azimuth a little and saves or costs a few hundred m/s.[R1, R25]
Launch windows
A fixed orbit plane turns relative to the Earth once per sidereal day, so the site passes under it twice: once on the ascending side and once on the descending side. Rendezvous or a sun-synchronous local time pins down which plane, so the window opens at a fixed time each day.
Key equations
where
needs i ≥ φ
where
where
where
second window at 180° − Δλ
Worked examples
Launching to a 51.6 deg station orbit
Given
A site at latitude 28.5 deg launches to a circular 200 km orbit inclined 51.6 deg. Find the inertial azimuth, the azimuth and speed needed relative to the rotating ground, and the time between the two daily windows.
new symbols
Working
- sin β = cos i / cos φ = cos 51.6 / cos 28.5 = 0.7068, β44.98 deg(north-east)
- Orbit speed 7784 m/s; ground speed eastward = ω_E RE cos φ408.7 m/s
- Relative to the ground: east 5502 − 409, north 5507: azimuth 42.77 deg, speed 7501 m/s
- Window geometry: sin Δλ = tan φ / tan i = 0.4303, Δλ25.49 deg
- Gap between the ascending and descending windows = (180 − 2Δλ)/360 × sidereal day8.58 h
Answer. Inertial azimuth 45.0 deg, ground-relative 42.8 deg; Earth's rotation saves 283 m/s; windows 8.6 h and 15.4 h apart.
Sanity check. The rotating-frame azimuth is a little north of the inertial one, because the ground already carries the rocket east.
inertial azimuth 44.98 deg; orbit speed 7784 m/s; ground speed 408.7 m/s relative to the ground: azimuth 42.77 deg, speed 7501 m/s, saving 283.3 m/s windows 8.58 h and 15.36 h apart
Practice
Free speed from Earth's rotation[R1, R25]
Given
Find the eastward speed a launch site carries at latitudes 5.2, 28.5 and 45.0 deg (ω_E = 7.292115e-05 rad/s, RE = 6378.137 km), and the difference between the first and the last.
v = ω_E RE cos φ; keep RE in metres for m/s.
Answer. 463, 409 and 329 m/s; the near-equatorial site gains 134 m/s for an eastward launch.
Working
- Equator: ω_E RE = 7.292115e-05 × 6378137465.1 m/s
- φ = 5.2 deg: 465.1 × cos 5.2463.2 m/s
- φ = 28.5 deg: 465.1 × cos 28.5408.7 m/s
- φ = 45.0 deg: 465.1 × cos 45.0328.9 m/s
- Difference134.3 m/s
Answer. 463, 409 and 329 m/s; the near-equatorial site gains 134 m/s for an eastward launch.
Sanity check. The gain is only a few per cent of the 7784 m/s orbit speed, but the bigger prize of a low latitude is reaching low inclination without a plane change.
equator: 465.1 m/s latitude 5.2 deg: 463.2 m/s latitude 28.5 deg: 408.7 m/s latitude 45.0 deg: 328.9 m/s difference 5.2 vs 45.0 deg: 134.3 m/s
The price of latitude[R25]
Given
A satellite must end up in an equatorial 200 km orbit. Each site launches due east (i = latitude) and removes the inclination with one plane change in LEO. Compare sites at 5.2 and 28.5 deg.
Launching due east gives the minimum inclination, equal to the latitude. Then Δv = 2v sin(Δi/2).
Answer. 706 m/s versus 3832 m/s: the higher latitude costs 3126 m/s more in LEO.
Working
- vLEO = √(μ/r)7.7843 km/s
- 5.2 deg: Δv = 2 × 7.7843 × sin(2.6 deg)0.706 km/s
- 28.5 deg: Δv = 2 × 7.7843 × sin(14.25 deg)3.832 km/s
- Penalty of the higher site3.126 km/s
Answer. 706 m/s versus 3832 m/s: the higher latitude costs 3126 m/s more in LEO.
Sanity check. That is so expensive that nobody does it in LEO; the change is made at GTO apogee where the speed is low (Chapter 19).
v_LEO = 7.7843 km/s 5.2 deg site: 706.2 m/s 28.5 deg site: 3832.2 m/s penalty of the higher site: 3126.0 m/s
A sun-synchronous launch from 57 deg north
Given
A site at 57 deg latitude launches to a 550 km sun-synchronous orbit. Find the inclination, the northbound and southbound inertial azimuths, and how much extra speed relative to the ground the retrograde orbit needs because Earth turns the wrong way.
Use the J2 sun-synchronous condition from Chapter 7. With i > 90 deg the formula gives a negative azimuth (north-west); the southbound option is 180 − β.
Answer. i = 97.59 deg; azimuths 346.0 deg (north) or 194.0 deg (south); the retrograde orbit costs 65 m/s extra.
Working
- iSSO97.593 deg
- sin β = cos 97.593 / cos 57 = -0.2426, β = -14.04 deg → northbound 345.96 deg
- Southbound azimuth = 180 − (-14.04)194.04 deg
- Ground speed east = 253.3 m/s; orbit speed 7585 m/s
- Speed needed relative to the ground = 7650.5 m/s, penalty 65.4 m/s
Answer. i = 97.59 deg; azimuths 346.0 deg (north) or 194.0 deg (south); the retrograde orbit costs 65 m/s extra.
Sanity check. Small penalty because the orbit is only 7.6 deg past polar and the site's ground speed is only 253 m/s at 57 deg.
i = 97.593 deg azimuth northbound 345.96 deg, southbound 194.04 deg orbit speed 7585 m/s, ground speed 253.3 m/s, needed 7650.5 m/s penalty 65.4 m/s
Common mistakes
- Expecting to reach an inclination below the site latitude without a plane change.
- Mixing the inertial azimuth with the azimuth relative to the ground.
- Treating every launch time as equal: a target plane passes over the site only twice a day.
- Forgetting that retrograde orbits (like sun-synchronous) pay for Earth's rotation instead of gaining.
Check yourself
Recap from earlier chapters
Sun-synchronous inclination at 550 km?
Δv for a 10 deg plane change at 7.7 km/s?
Sources for this chapter
- R1Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.
- R25ESA, "Europe's Spaceport".
Where in each source, and links: References.
Mission case end to end
Every chapter so far solved one piece. A real mission is all of them in a row: the launch site fixes the starting inclination, a transfer orbit lifts the satellite, the plane change goes where it is cheapest, station-keeping spends a little every year, and a final burn clears the slot. Add the pieces, apply the rocket equation, and you have the propellant the satellite must carry.
You will be able to
- Chain launch, transfer, plane change, station-keeping and disposal into one Δv budget
- Find the best split of a plane change between two burns numerically
- Turn a Δv budget into propellant mass and check it with an independent propagation
The plan
A communications satellite launches due east from 28.5 deg into a 200 km parking orbit. The upper stage fires at perigee to enter a geostationary transfer orbit. The satellite fires at apogee to circularise and remove the inclination.
Where to change the plane
Most of the plane change belongs at apogee, where the speed is low. A few degrees at perigee, folded into the large perigee burn, cost almost nothing extra and make the apogee burn cheaper. A one-variable minimisation finds the best split.
Check it twice
The coast from perigee to apogee is half the transfer period by Kepler's equation; a numerical propagation of the same orbit must reach the GEO radius with constant energy.
From Δv to kilograms
Station-keeping and disposal are added to the apogee burn, and the rocket equation turns the total into propellant. About half the separated mass ends up as propellant.
| Phase | Who burns | Δv (m/s) |
|---|---|---|
| Parking orbit → GTO perigee burn | upper stage | 2478 |
| Apogee burn: circularise + plane change | spacecraft | 1789 |
| North-south station-keeping, 15 years | spacecraft | 684 |
| East-west station-keeping, 15 years | spacecraft | 22 |
| Graveyard raise | spacecraft | 11 |
| Spacecraft total | 2507 |
Site latitude 28.5 deg, 200 km parking orbit, 2.17 deg of plane change at perigee. Isp 315 s and separation mass 5000 kg are assumptions.
Key equations
where
where
where
where
Worked examples
A GEO satellite from launch to retirement[R17, R18, R20]
Given
Launch due east from a 28.5 deg site into a 200 km circular parking orbit, transfer to GEO by a Hohmann-type GTO, take out the 28.5 deg of inclination with the best split between perigee and apogee, then keep station for 15 years and retire 300 km above GEO. Spacecraft: 5000 kg at separation, Isp 315 s.
new symbols
Working
- Launch azimuth: sin β = cos 28.5 / cos 28.5 = 1, β = 90 deg; the pad adds 409 m/s
- Parking orbit: r = 6,578.137 km, v = 7.7843 km/s, specific energy −μ/2r-30.297 km²/s²
- GTO: a = 24,371.2 km, e = 0.7301; perigee 10.2388 km/s, apogee 1.5974 km/s; coast 5.26 h (half the period, Kepler)
- Best split (numerical minimum): 2.169 deg at perigee, 26.331 deg at apogee
- Perigee burn 2478 m/s; apogee burn √(v_a² + v_GEO² − 2 va vGEO cos Δi)1789 m/s
- Saving against doing all the plane change at apogee: 24.1 m/s
- Station-keeping: 684 m/s N-S + 22 m/s E-W; graveyard 10.9 m/s
- Spacecraft total 2507 m/s; mass ratio exp(2507 / (315 × 9.80665))2.2513
- Propellant = 5000 − 5000/2.2513 = 2779 kg; 2221 kg left at end of life
Answer. Spacecraft Δv 2507 m/s, of which the apogee burn is 1789 m/s; propellant 2779 kg of the 5000 kg separated mass.
Sanity check. The numerical propagation of the GTO reaches 42,164.2 km at apogee with energy constant to 1e-11: the coast is consistent with Kepler. About 56 % of the satellite is propellant, typical of chemical GEO platforms.
pad speed 409 m/s; GTO a = 24,371.2 km, perigee 10.2388, apogee 1.5974 km/s best split 2.169 deg at perigee: perigee burn 2478 m/s, apogee burn 1789.3 m/s, saving 24.1 m/s station-keeping 684 + 22 m/s, graveyard 10.9 m/s spacecraft total 2506.8 m/s, mass ratio 2.2513 propellant 2779.1 kg of 5000 kg; 2221 kg left at end of life
Practice
The same satellite from 5.2 deg[R25]
Given
Launch due east from a 5.2 deg site into the same 200 km parking orbit and GTO. Doing the whole plane change at apogee, find the apogee burn and compare it with 1836 m/s from 28.5 deg.
Only the plane change differs; use the combined-burn law of cosines with va and vGEO from the worked example.
Answer. Apogee burn 1491 m/s, 346 m/s less than from 28.5 deg.
Working
- va = 1.5974 km/s, vGEO3.0747 km/s
- Δv = √(1.5974² + 3.0747² − 2 × 1.5974 × 3.0747 × cos 5.2 deg)1491 m/s
- Saving = 1836 − 1491346 m/s
Answer. Apogee burn 1491 m/s, 346 m/s less than from 28.5 deg.
Sanity check. Even with zero plane change the burn would be vGEO − va = 1477 m/s, so 5.2 deg costs only 14 m/s extra.
v_a = 1.5974 km/s, v_GEO = 3.0747 km/s apogee burn from 5.2 deg: 1490.9 m/s; from 28.5 deg: 1836.5 m/s saving 345.6 m/s
Ten seconds of Isp[R17]
Given
Keep the spacecraft Δv budget of 2507 m/s and separation mass 5000 kg. How much propellant is saved if the engine's Isp rises from 315 s to 325 s?
Propellant = m₀(1 − exp(−Δv/(Isp g₀))) for each Isp.
Answer. 56 kg less propellant, mass that can become payload.
Working
- Isp 315 s: mp = 5000 × (1 − exp(−2507/(315 × 9.80665)))2779.1 kg
- Isp 325 s: mp = 5000 × (1 − exp(−2507/(325 × 9.80665)))2722.9 kg
- Saving56.2 kg
Answer. 56 kg less propellant, mass that can become payload.
Sanity check. A 3.2 % rise in Isp saves 2.0 % of the propellant: the exponential makes Isp valuable on big budgets.
Isp 315 s: propellant 2779.2 kg Isp 325 s: propellant 2723.0 kg saving 56.2 kg
An Earth-observation variant[R14, R15]
Given
A small satellite goes to a 550 km sun-synchronous orbit from a 57 deg site. Find the inclination, check the J2 node rate is one turn per year, then estimate the yearly decay Δa = 2π CD (A/m) ρ a² per revolution × revolutions per year (CD = 2.2, A/m = 0.01 m²/kg, ρ = 2e-13 kg/m³) and the 5-year reboost Δv.
Assumption: ρ = 2e-13 kg/m³ at 550 km (moderate solar activity)
Work Δa in metres (a in m). A small reboost costs Δv ≈ v Δa / (2a). Southbound azimuth = 180 − β.
Answer. i = 97.59 deg; decay about 8.8 km per year; 24.0 m/s keeps the altitude for 5 years.
Working
- i = 97.593 deg; southbound azimuth from 57 deg: 194.04 deg
- Ω̇ = −1.5 n J₂ (RE/a)² cos i = 1.9911e-07 rad/s → × one year360.0000 deg
- Δa per rev = 2π × 2.2 × 0.01 × 2e-13 × (6.9281e+06)²1.592 m
- Revolutions per year = 5499 → Δa8.76 km/yr
- Reboost Δv = 7585.1 × 8.76 / (2 × 6,928.1) = 4.79 m/s per year → 5 years: 24.0 m/s
Answer. i = 97.59 deg; decay about 8.8 km per year; 24.0 m/s keeps the altitude for 5 years.
Sanity check. The density is the weak link: at solar maximum it can be tens of times higher, so real budgets carry a large margin.
i = 97.593 deg, southbound azimuth 194.04 deg, node turns 360.0000 deg/yr decay 1.592 m/rev x 5499 revs = 8.76 km/yr reboost 4.79 m/s per year, 5 years: 24.0 m/s
Common mistakes
- Doing the whole plane change in LEO: it costs kilometres per second instead of hundreds of m/s.
- Adding Δv from different phases but applying the rocket equation with the wrong starting mass.
- Forgetting the end-of-life reserve or the station-keeping years in the budget.
- Trusting one method alone. Cross-check analytic coasts with a numerical propagation.
Check yourself
Recap from earlier chapters
Specific energy of a 200 km circular orbit?
Hohmann LEO-to-GEO total from 200 km, coplanar?
North-south station-keeping per year?
Sources for this chapter
- R14NOAA, NASA and USAF, U.S. Standard Atmosphere, 1976, NASA-TM-X-74335 (NTRS 19770009539), Washington DC, Oct 1976.
- R15Boring, J.W. and Humphris, R.R., "Drag coefficients for spheres in free molecular flow in O at satellite velocities", NASA CR-2233, 1973.
- R17ArianeGroup GmbH, "400N bipropellant apogee motors" product page (S400-12, S400-15), © 2002-2024.
- R18Lovell, R.R. and O'Malley, T.A., Station Keeping of High Power Communication Satellites, NASA TM X-2136, Lewis Research Center, Nov 1970.
- R20Inter-Agency Space Debris Coordination Committee, IADC Space Debris Mitigation Guidelines, IADC-02-01 Rev. 2, Mar 2020.
- R25ESA, "Europe's Spaceport".
Where in each source, and links: References.
Formula sheet
Two-body basics
Time and state
Perturbations
Manoeuvres
Rendezvous
Lambert
Interplanetary
Flybys
Lagrange points
GEO
Launch
Budgets
Glossary
- Apsidal rotation
- The turning of the ellipse's major axis within the orbit plane, also caused by J₂.
- Argument of latitude (u)
- Angle from the ascending node to the satellite, u = ω + ν; defined even for circular orbits.
- Bi-elliptic transfer
- Three-burn transfer via a distant apoapsis; cheaper than Hohmann for radius ratios above 11.94.
- C3
- Characteristic energy, the square of the hyperbolic excess speed, in km²/s².
- Circular restricted three-body problem
- Two massive bodies on circular orbits about their barycentre and a third body too light to affect them.
- Clohessy-Wiltshire equations
- Linearised equations of relative motion about a target on a circular orbit.
- Critical inclination
- 63.4° (or 116.6°): the inclination at which J₂ leaves the argument of perigee fixed.
- Eccentric anomaly (E)
- Auxiliary angle that turns Kepler's equation into M = E − e sin E.
- Eccentricity (e)
- Shape of a conic: 0 circle, below 1 ellipse, 1 parabola, above 1 hyperbola.
- Flight-path angle (γ)
- Angle of the velocity above the local horizontal.
- Frozen orbit
- An orbit whose small eccentricity and perigee at 90° make the J₂/J₃ drift of e and ω vanish.
- Geostationary orbit (GEO)
- Circular, equatorial orbit with a period of one sidereal day, radius about 42,164 km.
- Graveyard orbit
- Disposal orbit a few hundred km above GEO used at end of life.
- Gravitational parameter (μ)
- G times the mass of the central body; for Earth 398600.4418 km³/s².
- Gravity assist
- A planetary flyby that turns the craft's velocity relative to the planet and so changes its speed around the Sun.
- Ground track
- The path traced on Earth's surface by the point directly beneath a satellite.
- Halo orbit
- A periodic three-dimensional orbit around a collinear Lagrange point.
- Hill sphere
- Region around a smaller body where its gravity dominates the tidal pull of the larger one; radius about a(μ*/3)^(1/3).
- Hohmann transfer
- Two-burn, tangential transfer between coplanar circles along a half-ellipse; cheapest two-burn option.
- Hyperbolic excess speed (v∞)
- Speed left over far from the planet on an escape hyperbola.
- J₂
- Coefficient describing Earth's equatorial bulge; the largest perturbation for most Earth orbits.
- Lagrange point
- An equilibrium of the restricted three-body problem in the rotating frame; there are five, L1 to L5.
- Lambert's problem
- Find the orbit (and velocities) joining two positions in a given flight time.
- Launch azimuth
- Direction of flight at launch, measured clockwise from north.
- Launch window
- The time interval in which a launch reaches the required orbit plane or target.
- LVLH frame
- Local vertical, local horizontal: axes radial, along-track and cross-track, riding with the target.
- Mean anomaly (M)
- Angle that grows uniformly with time, M = n t.
- Molniya orbit
- A 12-hour, highly elliptical orbit at 63.4° whose long apogee dwell serves high latitudes.
- Nodal regression
- The steady turning of the orbit plane's ascending node caused by J₂; westward for prograde orbits.
- Oberth effect
- A burn adds more orbital energy when made at high speed, deep in a gravity well.
- Perifocal frame
- Orbit-plane frame with x towards perigee and z along h.
- Phasing manoeuvre
- Temporary change of period that moves a satellite along its own orbit.
- Plane change
- A burn that rotates the orbit plane; costs 2v sin(Δi/2) at speed v.
- Porkchop plot
- Contour map of launch energy (C3) or arrival speed over departure date and flight time.
- Repeat ground track
- An orbit that makes k revolutions in exactly d days and then retraces its path over the ground.
- Semi-latus rectum (p)
- Radius of an orbit at 90° from perigee; p = h²/μ = a(1 − e²).
- Sidereal day
- One rotation of Earth relative to the stars, 86164.09 s.
- Sphere of influence
- Region around a planet inside which its gravity is treated as the only one, for patched conics.
- Station-keeping
- Small regular burns that hold a satellite in its assigned orbit or longitude box.
- Sun-synchronous orbit
- A near-polar retrograde orbit whose node turns 360° per year, keeping a fixed local solar time.
- Synodic period
- Time between repeats of the same relative geometry of two orbiting bodies.
- True anomaly (ν)
- Angle at Earth's centre from perigee to the satellite.
- Universal variable
- A single variable z that describes ellipses, parabolas and hyperbolas in one set of equations.
- V-bar / R-bar
- The target's along-track line and radial line; standard approach corridors for rendezvous.
- Vis-viva equation
- v² = μ(2/r − 1/a): speed anywhere on an orbit from its size.
Answers
| Problem | Answer |
|---|---|
| O1.1 | v = 7.585 km/s, T = 95.6 min, 15.05 orbits per day. |
| O1.2 | h = 758 km for a 100-min period; h = 881 km for 14 orbits per sidereal day. |
| O1.3 | T = 43854 s = 12.182 h by either route. |
| O2.1 | a = 7378.1 km, e = 0.0678 (dimensionless). |
| O2.2 | r(135°) = 26097 km; r = 20000 km at ν = ±120.0°. |
| O2.3 | p = 15360 km, ν∞ = 140.3°, turn 100.6°, r(100°) = 19838 km. |
| O3.1 | vc = 7.452 km/s, vesc = 10.538 km/s, Δv = 3.087 km/s. |
| O3.2 | a = 9515 km, e = 0.1810, perigee altitude 1415 km, apogee altitude 4860 km. |
| O3.3 | ε = 8.47 km²/s², v∞ = 4.117 km/s, e = 1.2882, v = 4.212 km/s at a million km. |
| O4.1 | t = 25.7 min after perigee. |
| O4.2 | ν = 131.57°, r = 10226 km (altitude 3847 km). |
| O4.3 | t = 3.42 h after perigee, at ν = 116.9°. |
| O5.1 | |h| = 51966 km²/s, i = 47.12°. |
| O5.2 | a = 11500 km, e = 0.250, i = 35.00°, Ω = 80.00°, ω = 40.00°, ν = 60.00°. |
| O5.3 | u = 10.000° is reliable; ω and ν separately are ill-conditioned when e ≈ 0, though in exact arithmetic the round trip still closes (error about 10⁻¹² km, floating-point level). |
| O6.1 | Ω̇ = -5.14 °/day: the node drifts west by about 5.1 deg every day. |
| O6.2 | adrag = 1.9e-08 m/s² and aSRP = 5.9e-08 m/s²: at 700 km they are within a factor of 3.1 of each other. |
| O6.3 | About 2.2 km lost in 30 days; the reboost costs about 1.3 m/s. |
| O7.1 | i = 99.03°. |
| O7.2 | The node runs ahead of the Sun by 4.37° per year, a local-time drift of about 17 min per year. |
| O7.3 | Apogee altitude 39367 km, e = 0.722; 73.0 % of each orbit is spent in the apogee half; at 60° the perigee would drift 0.037 deg/day. |
| O8.1 | 24.72° with J₂ against 24.76° without: the passes are about 24.7° apart. |
| O8.2 | ±3° in latitude but only ±0.039° in longitude: a tall, thin figure-8. |
| O8.3 | h = 720.0 km at i = 98.27°: 29 orbits in 2 days, then the track repeats. |
| O9.1 | Δv₁ = 164 m/s, Δv₂ = 161 m/s, total 325 m/s. |
| O9.2 | aph = 27068 km and a total Δv of 69.8 m/s over 3 revolutions. |
| O9.3 | Bi-elliptic 4014 m/s against Hohmann 4053 m/s, a saving of 39 m/s, but 221 h instead of 18 h. |
| O10.1 | Δv = 655 m/s. |
| O10.2 | 1486 m/s from 5.2° and 2248 m/s from 45°, a difference of 762 m/s. |
| O10.3 | Do about 2.6° of the turn at perigee (s = 0.066); total 4504 m/s, saving 34 m/s. |
| O11.1 | Δy = +3.41 km: about 3.4 km ahead after one orbit. |
| O11.2 | ẏ = 3.36 m/s, so the gap closes in 1.7 h. |
| O11.3 | Δv₁ = 6.27 m/s, Δv₂ = 2.19 m/s, total 8.47 m/s. |
| O12.1 | c = 12806 km, s = 15403 km, amin = 7702 km. |
| O12.2 | |v₁| = 8.657 km/s, a = 10235 km, e = 0.401, departure burn 2.389 km/s. |
| O12.3 | The 60 min transfer is cheapest at 0.89 km/s total. |
| O13.1 | Tsyn = 780 days, about 26 months between Mars windows. |
| O13.2 | tH = 146 days, Venus 54.0° behind Earth at departure, windows every 584 days. |
| O13.3 | C3 rises to 14.4 km²/s², 66 % more than the optimum. |
| O14.1 | Mars SOI ≈ 577,225 km (170 Mars radii); Earth SOI ≈ 924,647 km. |
| O14.2 | (a) 2.08 km/s into the circular orbit; (b) 1.03 km/s into the ellipse. |
| O14.3 | ropt = 12,208 km (altitude 8,811 km), Δvmin = 1.873 km/s, 207 m/s less than at 400 km. |
| O15.1 | e = 1.504, turn angle δ = 83.3°. |
| O15.2 | Δv = 6.65 km/s at 500 km, up to 6.82 km/s for a grazing pass. |
| O15.3 | δ = 122°, speed 20.84 km/s against an escape speed of 18.47 km/s: the probe escapes the Sun. |
| O16.1 | r ≈ 1,084,053 km (the exact root gives 1,082,332 km). |
| O16.2 | L2 is 64,515 km beyond the Moon; Hill is off by -5.0 %. |
| O16.3 | λ = 2.484 n, so a small error grows by a factor e every 23 days. |
| O17.1 | λ̇ = -0.0321 deg/day westward; it leaves the box after 1.6 days. |
| O17.2 | Δv ≈ 10.9 m/s in two near-equal burns, about 6.7 kg of propellant. |
| O17.3 | Burn 1.70 m/s retrograde, drift orbit 46.7 km below GEO, total 3.4 m/s. |
| O18.1 | 463, 409 and 329 m/s; the near-equatorial site gains 134 m/s for an eastward launch. |
| O18.2 | 706 m/s versus 3832 m/s: the higher latitude costs 3126 m/s more in LEO. |
| O18.3 | i = 97.59 deg; azimuths 346.0 deg (north) or 194.0 deg (south); the retrograde orbit costs 65 m/s extra. |
| O19.1 | Apogee burn 1491 m/s, 346 m/s less than from 28.5 deg. |
| O19.2 | 56 kg less propellant, mass that can become payload. |
| O19.3 | i = 97.59 deg; decay about 8.8 km per year; 24.0 m/s keeps the altitude for 5 years. |
Assumptions and sources
Each input names the source it was checked against; R-ids link to the References page. Verdict: CONFIRMED means the source gives the same value; ADJUSTED means the value was changed to match the source; UNVERIFIED means no readable open source was found (the reason is given), so check it before using it in a real design. This pack: 19 confirmed, 4 adjusted, 4 unverified.
| Assumption | Source | Verdict |
|---|---|---|
| Sun gravitational parameter μ_Sun = 1.32712440018 × 10¹¹ km³/s² | R2: JPL SSD Astrodynamic Parameters (DE440) (1.32712440041e11); R8: IAU 2015 B3 | CONFIRMEDagrees to 2e-12 relative; ours is the earlier DE405 value |
| Moon gravitational parameter μ_Moon = 4902.800 km³/s² | R2: JPL SSD Astrodynamic Parameters (DE440) | CONFIRMED |
| Mean Earth-Moon distance 384,400 km | R6: NASA NSSDCA Moon Fact Sheet (2024) | CONFIRMED |
| Venus orbit radius 0.723332 AU (circular model) | R4: Standish and Williams, Approximate Positions of the Planets, Table 1 | CONFIRMED5e-6 relative |
| Venus μ = 324,859 km³/s² | R2: JPL SSD Astrodynamic Parameters (DE440) | CONFIRMED |
| Venus radius 6051.8 km | R3: JPL SSD Planetary Physical Parameters (Archinal et al. 2018) | CONFIRMED |
| Earth orbit radius 1.000000 AU (circular model) | R4: Standish and Williams, Approximate Positions of the Planets, Table 1 | CONFIRMED3e-6 |
| Mars orbit radius 1.523679 AU (circular model) | R4: Standish and Williams, Approximate Positions of the Planets, Table 1 | CONFIRMED2e-5 |
| Mars μ = 42,828 km³/s² | R2: JPL SSD Astrodynamic Parameters (DE440) | CONFIRMED |
| Mars equatorial radius 3396.2 km | R3: JPL SSD Planetary Physical Parameters (Archinal et al. 2018) | CONFIRMED |
| Jupiter orbit radius 5.2026 AU (circular model) | R4: Standish and Williams, Approximate Positions of the Planets, Table 1 | CONFIRMED6e-5; rounded value |
| Jupiter μ = 1.26686 × 10⁸ km³/s² | [R2] R8: IAU 2015 B3 nominal (GM)J = 1.2668653e8 km³/s² | CONFIRMED |
| Jupiter equatorial radius 71,492 km | R3: JPL SSD Planetary Physical Parameters (Archinal et al. 2018); R8: IAU 2015 B3 nominal | ADJUSTEDto 71,492 km |
| Drag coefficient CD = 2.2 for a compact satellite in free-molecular flow (representative; measured sphere values are just above 2) | R15: Boring and Humphris, NASA CR-2233 (1973), abstract | CONFIRMEDas representative (2.2 sits just above the measured sphere value; real values depend on shape and surface, typically 2-2.5). Note: ISS operational fit uses CD = 1.10 with its own reference area (R21), so CD only has meaning together with the area it is paired with |
| Radiation-pressure coefficient CR = 1.3 (1 = fully absorbing, 2 = perfect mirror) | UNVERIFIED: the 1-2 bounds are definitional; no readable source for 1.3 as a typical value | UNVERIFIEDthe bounds are definitional; the typical value 1.3 has no readable primary. Kept, labelled representative |
| Area-to-mass ratio A/m = 0.01 m²/kg for a typical small satellite | UNVERIFIED: representative design input, not a measured value | UNVERIFIEDdesign input for a generic small satellite, not a physical constant. Kept, labelled representative |
| Area-to-mass ratio A/m = 0.0032 m²/kg for a large crewed station (ISS drag area over mass) | R21: NASA JSC TOPO, ISS public ephemeris header (2 Nov 2026): 1511.8 m², 473,406.5 kg | ADJUSTEDto 0.0032 m²/kg. With the pack's CD = 2.2 the product CD A/m is still about 2x the TOPO fit (0.0070 vs 0.0035), so the station decay example is an upper-end estimate; say so in the text |
| Atmospheric density table, 200-1000 km, U.S. Standard Atmosphere 1976 (moderate solar activity; real density varies by tens of times over the solar cycle at these heights) | R14: U.S. Standard Atmosphere 1976 (NASA-TM-X-74335), Table I, pp. 69-73 | ADJUSTEDrows 250 → 6.1e-11, 900 → 5.8e-15, 1000 → 3.6e-15 (other nine rows already equal USSA76 to 2 s.f.). Cross-check: Vallado's exponential model (R11 ATMOSEXP.DAT) is 10-25 % higher in 300-600 km, consistent with the "varies about 10x over the cycle" caveat. Re-cite to R14 (the SMAD citation was never read) |
| Monopropellant hydrazine thruster specific impulse Isp = 220 s | R16: NASA Small Spacecraft State of the Art 2026, ch. 4, §4.6.1 (200-235 s) | CONFIRMED |
| Lunisolar inclination drift of an uncontrolled GEO satellite about 0.85 deg/yr | R18: Lovell and O'Malley, NASA TM X-2136 (1970), Summary; R19: Decoto and Loerch, AMOS 2015, §2 | CONFIRMEDfor 0.85. The 0.75-0.95 range over the 18.6-yr lunar node cycle is UNVERIFIED (no readable source); Phase 2 drops the range from the label rather than keep an unsourced number |
| Typical GEO east-west station-keeping budget 1.5 m/s per year (sourced bound: under 1.75 m/s/yr for the triaxiality drift; about 2 m/s/yr often quoted) | R18: NASA TM X-2136 (1970), Summary and Introduction | CONFIRMED1.5 lies inside the sourced bound; range 1-2 consistent with R18 |
| End-of-life re-orbit about 300 km above GEO (IADC minimum 235 km + 1000 CR A/m; 300 km covers A/m up to 0.05 m²/kg at CR = 1.3) | R20: IADC Space Debris Mitigation Guidelines, IADC-02-01 Rev. 2 (2020), §5.3.1, p. 11 | CONFIRMEDwith CR = 1.3, 300 km covers A/m up to 0.05 m²/kg (formula gives 261 km at 0.02, 300 km at 0.05) |
| Spacecraft bipropellant apogee engine and thrusters, Isp = 315 s (slightly below 400 N apogee motors at 318-321 s) | R17: ArianeGroup 400 N apogee motors S400-12 / S400-15 | CONFIRMED315 s is slightly conservative for the apogee engine; small biprop thrusters are lower |
| Mission case: spacecraft mass at separation in GTO 5000 kg | UNVERIFIED: representative mid-size GEO platform (design input) | UNVERIFIEDmission-case design input, not a fact; labelled "representative mid-size GEO platform" |
| Mission case: air density at 550 km about 2.4 × 10⁻¹³ kg/m³ (U.S. Standard Atmosphere 1976, moderate solar activity; varies by about 10x over the cycle) | R14: U.S. Standard Atmosphere 1976 (NASA-TM-X-74335), Table I, p. 72 | ADJUSTEDto 2.4e-13 kg/m³ (consistent with the corrected table; our 5e-13 was 2.1x the moderate-activity reference). Note: the ADCS v2 pack uses the same 5e-13 (not edited here) |
| Drag coefficient CD = 2.2 for a compact satellite | R15: Boring and Humphris, NASA CR-2233 (1973), abstract | CONFIRMEDas representative |
| Mission case: area-to-mass ratio A/m = 0.01 m²/kg | UNVERIFIED: representative small satellite (design input) | UNVERIFIEDdesign input |
Fixed constants: μ = 398600.4418 km³/s² (IERS Conventions 2010), Rₑ = 6378.137 km (GRS80), J₂ = 0.00108262668 (EGM96), J₃ = -2.53215306e-06 (EGM2008 to 4 s.f.), Earth rotation 7.292115e-05 rad/s (IERS), sidereal day 86164.0905 s (JPL SSD); AU = 149597870.7 km (IAU 2012 Resolution B2). Planet, Sun and Moon values are listed as assumptions with their sources; full references are listed below.
Built by engineers, for engineers. Every worked number is computed and asserted by the build script.
References
Every reference the pack cites. "Where" names the section, table or page that holds the cited value; an entry says so when the source was not read in full or is not openly available.
| Id | Reference |
|---|---|
| R1 | Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS Technical Note 36, Ch. 1 "General definitions and numerical standards" (update of 16 Nov 2017). IERS.Where: Table 1.1 (GM⊕, c, au, ERA rate), Table 1.2 (GRS80: a_E, J2, ω), pp. 3-4https://iers-conventions.obspm.fr/content/chapter1/icc1.pdfCited in: Chapters 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 17, 18 |
| R2 | NASA JPL Solar System Dynamics, "Astrodynamic Parameters" (planetary GMs from DE440; Park et al. 2021, AJ 161:105).Where: GM Sun, GM Moon, planet system GMs, au, mean sidereal dayhttps://ssd.jpl.nasa.gov/astro_par.htmlCited in: Chapters 1, 6, 8, 13, 14, 15, 16, 17; assumptions |
| R3 | NASA JPL Solar System Dynamics, "Planetary Physical Parameters" (radii from Archinal et al. 2018, CMDA 130:22; change log 2019-12-12).Where: equatorial radii of Venus, Mars, Jupiterhttps://ssd.jpl.nasa.gov/planets/phys_par.htmlCited in: Chapters 14, 15; assumptions |
| R4 | Standish, E.M. and Williams, J.G., "Approximate Positions of the Planets", NASA JPL SSD.Where: Table 1 (1800-2050 AD), semi-major axeshttps://ssd.jpl.nasa.gov/planets/approx_pos.htmlCited in: Chapters 13, 14, 15; assumptions |
| R5 | Williams, D.R., "Earth Fact Sheet", NASA NSSDCA, Goddard Space Flight Center, updated 15 Nov 2024.Where: GM, R_E, J2, sidereal rotation, tropical year, e, obliquityhttps://nssdc.gsfc.nasa.gov/planetary/factsheet/earthfact.htmlCited in: Chapters 1, 7 |
| R6 | Williams, D.R., "Moon Fact Sheet", NASA NSSDCA, updated 11 Jan 2024.Where: semi-major axis 0.3844 × 10⁶ kmhttps://nssdc.gsfc.nasa.gov/planetary/factsheet/moonfact.htmlCited in: Chapter 16; assumptions |
| R7 | IAU 2012 Resolution B2, "Re-definition of the astronomical unit of length", XXVIII General Assembly, Beijing, Aug 2012 (text hosted by SYRTE, Observatoire de Paris).Where: au = 149 597 870 700 m exactlyhttps://syrte.obspm.fr/IAU_resolutions/IAUResol_2012_0.htmlCited in: Chapter 13 |
| R8 | Mamajek, E.E., Prša, A., Torres, G. et al., "IAU 2015 Resolution B3 on Recommended Nominal Conversion Constants for Selected Solar and Planetary Properties", arXiv:1510.07674 (2015).Where: p. 2: S = 1361 W/m², (GM)☉ = 1.3271244e20, R_eJ = 7.1492e7 m, (GM)J = 1.2668653e17 m³/s²https://arxiv.org/abs/1510.07674Cited in: Chapters 6, 15; assumptions |
| R9 | NIST, CODATA 2022 recommended values, "standard acceleration of gravity".Where: g_n = 9.80665 m/s², exacthttps://physics.nist.gov/cgi-bin/cuu/Value?gnCited in: Chapter 9 |
| R10 | MathWorks, Aerospace Toolbox documentation, gravityzonal (Earth zonal coefficients standardised to EGM2008, Pavlis et al. 2012). Secondary.Where: Algorithms table: J2 = 1.08262617e-3, J3 = −2.53241e-6, J4 = −1.61990e-6https://www.mathworks.com/help/aerotbx/ug/gravityzonal.html |
| R11 | Vallado, D.A., companion software to Fundamentals of Astrodynamics and Applications (CelesTrak), GitHub repo CelesTrak/fundamentals-of-astrodynamics, commit 49df047 (AGPL-3.0; values read, no code copied).Where: software/python/tests/astro/twobody/test_frame_conversions.py (Ex. 2-5 rv2coe values); software/python/src/valladopy/constants.py (EGM-08 constants, sidereal day); software/misc/pascal/ATMOSEXP.DAT (exponential atmosphere table)https://github.com/CelesTrak/fundamentals-of-astrodynamicsCited in: Chapters 5, 17 |
| R12 | poliastro project, tests/test_iod.py::test_curtis52 (Curtis, Orbital Mechanics for Engineering Students, Example 5.2), commit 21fd771. Secondary record of a textbook case.Where: r₁, r₂, 1 h, expected v₁, v₂https://github.com/poliastro/poliastro/blob/main/tests/test_iod.pyCited in: Chapter 12 |
| R13 | Izzo, D., "Revisiting Lambert's problem", Celestial Mechanics and Dynamical Astronomy 121 (2015) 1-15, doi:10.1007/s10569-014-9587-y; open preprint arXiv:1403.2705.Where: method and test philosophy for Lambert solvershttps://arxiv.org/abs/1403.2705Cited in: Chapter 12 |
| R14 | NOAA, NASA and USAF, U.S. Standard Atmosphere, 1976, NASA-TM-X-74335 (NTRS 19770009539), Washington DC, Oct 1976.Where: Table I, geometric altitude, metric units, density column: pp. 69-73 (pdf pp. 84-88)https://ntrs.nasa.gov/citations/19770009539Cited in: Chapters 6, 19; assumptions |
| R15 | Boring, J.W. and Humphris, R.R., "Drag coefficients for spheres in free molecular flow in O at satellite velocities", NASA CR-2233, 1973.Where: abstract: measured C_D "slightly greater than 2" at orbital energieshttps://ntrs.nasa.gov/citations/19730010970Cited in: Chapters 6, 19; assumptions |
| R16 | NASA Small Spacecraft Systems Virtual Institute, State-of-the-Art of Small Spacecraft Technology, Ch. 4 "In-Space Propulsion", 2026 edition.Where: §4.6.1: hydrazine thrusters typically 200-235 s (1 N class and up); Table 4-1https://www.nasa.gov/smallsat-institute/sst-soa/in-space_propulsion/Cited in: Chapter 9; assumptions |
| R17 | ArianeGroup GmbH, "400N bipropellant apogee motors" product page (S400-12, S400-15), © 2002-2024.Where: Isp at nominal point 318 s (S400-12), 321 s (S400-15), MMH/MONhttps://space-propulsion.com/spacecraft-propulsion/apogee-motors/index.htmlCited in: Chapters 17, 19; assumptions |
| R18 | Lovell, R.R. and O'Malley, T.A., Station Keeping of High Power Communication Satellites, NASA TM X-2136, Lewis Research Center, Nov 1970.Where: Summary (pdf p. 4): inclination grows 0.85 deg/yr; triaxiality E-W Δv < 1.75 m/s/yr; Introduction (pdf p. 5): 46 m/s N-S and 2 m/s E-W per year "documented in a number of reports"https://ntrs.nasa.gov/citations/19710001770Cited in: Chapters 17, 19; assumptions |
| R19 | Decoto, J. and Loerch, P. (Orbital ATK), "Technique for GEO RSO Station Keeping Characterization and Maneuver Detection", AMOS Technical Conference, 2015.Where: §2: about 0.85 deg/yr for a beginning-of-life GEO satellite (citing Soop 1994), about 46 m/s/yr N-Shttps://amostech.com/TechnicalPapers/2015/SSA/Decoto.pdfCited in: Chapter 17; assumptions |
| R20 | Inter-Agency Space Debris Coordination Committee, IADC Space Debris Mitigation Guidelines, IADC-02-01 Rev. 2, Mar 2020.Where: §5.3.1, p. 11: ΔH_min = 235 km + 1000 C_R A/m, e ≤ 0.003https://orbitaldebris.jsc.nasa.gov/library/iadc-space-debris-guidelines-revision-2.pdfCited in: Chapters 17, 19; assumptions |
| R21 | NASA JSC Flight Operations Directorate, TOPO office, ISS public ephemeris ISS.OEM_J2K_EPH.txt (CCSDS OEM, created 2026-11-02).Where: header: MASS = 473406.50 kg, DRAG_AREA = 1511.80 m², DRAG_COEFF = 1.10https://nasa-public-data.s3.amazonaws.com/iss-coords/current/ISS_OEM/ISS.OEM_J2K_EPH.txtCited in: Chapter 6; assumptions |
| R22 | ESA, "Types of orbits", 30 Mar 2020.Where: GEO altitude 35,786 km; LEO below 2000 km; SSO 600-800 km; L1/L2 about 1.5 million kmhttps://www.esa.int/Enabling_Support/Space_Transportation/Types_of_orbitsCited in: Chapters 7, 16, 17 |
| R23 | Cornish, N.J., "What is a Lagrange Point?", NASA Science, 27 Mar 2018, updated 1 Nov 2024.Where: L1-L3 unstable (about 23 days), L4/L5 stable if mass ratio > 24.96https://science.nasa.gov/resource/what-is-a-lagrange-point/Cited in: Chapter 16 |
| R24 | NASA Science, "Webb Orbit".Where: L2 about 1.5 million km from Earth; orbit about L2 takes about 6 months; small periodic thrusthttps://science.nasa.gov/mission/webb/orbit/Cited in: Chapter 16 |
| R25 | ESA, "Europe's Spaceport".Where: Kourou at latitude 5°3′ N; Earth rotation adds 460 m/shttps://www.esa.int/Enabling_Support/Space_Transportation/Europe_s_Spaceport/Europe_s_Spaceport2Cited in: Chapters 10, 18, 19 |
| R26 | NIST Time and Frequency Division, "Leap second and UT1-UTC information" (created 6 Aug 2009, updated 28 Oct 2026).Where: TAI − UTC = 37 s; leap seconds keep UT1 − UTC within ±0.9 s; DUT1 broadcasthttps://www.nist.gov/pml/time-and-frequency-division/time-realization/leap-secondsCited in: Chapter 1 |
| R27 | U.S. Naval Observatory, Astronomical Applications Dept., "Computing Approximate Sidereal Time".Where: GMST formula in hours from days since J2000.0; LST = GMST + east longitude/15; sidereal time = hour angle of the vernal equinoxhttps://aa.usno.navy.mil/faq/GASTCited in: Chapter 2 |
| R28 | Petit, G. and Luzum, B. (eds.), IERS Conventions (2010), IERS TN 36, Ch. 5 (10 Aug 2012 update).Where: eq. (5.43), argument F14: general precession in longitude p_A = 0.02438175 t rad (t in Julian centuries) → 50.29″/yr, about 25,770 yr per turnhttps://iers-conventions.obspm.fr/content/chapter5/icc5.pdfCited in: Chapter 2 |
| R29 | NASA Glenn Research Center, "Ideal Rocket Equation", Beginner's Guide to Aeronautics, updated 21 Nov 2023.Where: Δu = V_eq ln(m_f/m_e) = Isp g₀ ln MR; with gravity Δu = V_eq ln MR − g₀ t_bhttps://www1.grc.nasa.gov/beginners-guide-to-aeronautics/ideal-rocket-equation/Cited in: Chapter 9 |