Chapter 12 — Key Takeaways (Perturbations)
A one-page reference. Reread this before an exam, or before you size any orbit's station-keeping budget.
The perturbed equation of motion
$$\ddot{\mathbf{r}} = -\frac{\mu}{r^3}\mathbf{r} + \mathbf{a}_p, \qquad \mathbf{a}_p = \mathbf{a}_{J2} + \mathbf{a}_{\text{drag}} + \mathbf{a}_{\text{3-body}} + \mathbf{a}_{\text{SRP}} + \cdots$$
The two-body term is huge; every perturbation is tiny. What matters is not magnitude but whether the effect accumulates (secular) or averages away (periodic).
Perturbation hierarchy (order of magnitude, per unit mass)
| Effect | At LEO (400 km) | At GEO | Character |
|---|---|---|---|
| Two-body gravity | $8.7\ \text{m/s}^2$ | $0.22\ \text{m/s}^2$ | the orbit |
| J2 oblateness | $\sim 1\times10^{-2}$ | $\sim 8\times10^{-6}$ | secular ($\Omega, \omega$); $\propto 1/r^4$ |
| Drag | $\sim 10^{-6}$ (±10×) | $\to 0$ | secular in energy; only low orbits |
| Third-body (Moon+Sun) | $\sim 10^{-6}$ | $\sim 10^{-5}$ | grows as $r$; dominates beyond GEO |
| SRP | $\sim 10^{-7}$ | $\sim 10^{-7}$ | depends on $A/m$ |
Ranking flips with altitude: J2 rules LEO; luni-solar $\gtrsim$ J2 at GEO; drag matters only low down.
The two J2 secular rates (memorize the structure)
| Effect | Rate | What it does | Sign / special value |
|---|---|---|---|
| Nodal regression | $\dot\Omega = -\dfrac{3}{2}\dfrac{nJ_2R_\oplus^2}{(1-e^2)^2 a^2}\cos i$ | swivels the plane ($\Omega$) | west if prograde; 0 at $i=90^\circ$ |
| Apsidal precession | $\dot\omega = \dfrac{3}{4}\dfrac{nJ_2R_\oplus^2}{(1-e^2)^2 a^2}(5\cos^2 i - 1)$ | rotates the ellipse ($\omega$) | 0 at $i=63.4^\circ$ (critical) |
with $n = \sqrt{\mu/a^3}$. J2 does NOT change $a$ or $e$ secularly — it only reorients the orbit. Both rates scale as $a^{-7/2}$, so they are ~600× stronger at LEO than GEO.
The three J2 payoffs (Chapter 9's promises, delivered)
| Design | Condition | Result |
|---|---|---|
| Sun-synchronous orbit | tune $i$ so $\dot\Omega = +0.9856^\circ/\text{day}$ (eastward → $i>90^\circ$) | constant lighting; $i \approx 98^\circ$ at 700–800 km |
| Molniya frozen apogee | fly at critical inclination $i = 63.4^\circ$ | $\dot\omega = 0$; apogee stays over the north |
| Polar orbit | $i = 90^\circ$ | $\dot\Omega = 0$; plane fixed in inertial space |
Sun-synchronous inclination: $\cos i = -\dot\Omega_{\text{req}}/\left[\frac{3}{2}nJ_2(R_\oplus/a)^2\right]$, $\dot\Omega_{\text{req}} = 1.991\times10^{-7}\ \text{rad/s}$.
Drag and orbital decay
| Quantity | Relation | Notes |
|---|---|---|
| Drag deceleration | $a_D = \frac{1}{2}\rho v^2/\beta$ | $\beta = m/(C_D A)$ = ballistic coefficient (🔗 Ch. 5/7) |
| Decay rate (circular) | $\dot a = -\rho v a/\beta$ | $\propto 1/\beta$: light/broad decays fast |
| Per revolution | $\Delta a_{\text{rev}} = -2\pi\rho a^2/\beta$ | — |
Drag is strongest at perigee → lowers apogee first → circularizes, then spirals in. Density at 400 km varies ~10× over the solar cycle — lifetime predictions carry wide error bars. Rough lifetimes: 200 km → days; 400 km → months (reboost); 800 km → decades; ≥1000 km → centuries (self-cleaning fails).
Solar radiation pressure & third-body
- SRP: $P = S/c = 4.54\ \mu\text{Pa}$ (absorbed; $\times 2$ for a mirror); $a_{\text{SRP}} = (S/c)(1+r)\,A/m$. Matters for high area-to-mass and at high altitude (no drag).
- Third-body (tidal): $a_3 \sim 2\mu_3 r/d^3$. Moon $\approx 2\times$ Sun (proximity beats mass). Drives GEO inclination up $\sim 0.85^\circ/\text{yr}$ → the expensive north–south station-keeping.
Station-keeping delta-v (per year)
| Orbit | Fights | $\Delta v$ |
|---|---|---|
| GEO north–south | luni-solar inclination drift | $\sim 45$–$50\ \text{m/s/yr}$ (~90% of budget) |
| GEO east–west | Earth triaxiality ($J_{22}$) | $\sim 2$–$4\ \text{m/s/yr}$ |
| LEO reboost | drag | a few to tens of m/s/yr (altitude-dependent) |
| GEO disposal | one-time | $\sim 11\ \text{m/s}$ to graveyard |
15-yr GEO total $\approx 0.78\ \text{km/s}$ → propellant via $m_p = m_0(1-e^{-\Delta v/v_e})$. Electric ($I_{sp}\sim1500$) vs chemical ($\sim300$) saves hundreds of kg → often decides mission lifetime (theme 4).
Numbers worth memorizing
- $J_2 = 1.0826\times10^{-3}$; ISS nodal regression $\approx -5^\circ/\text{day}$ (plane loops in ~72 days).
- Sun-synchronous $i \approx 98^\circ$; critical inclination $63.4^\circ$; polar $\dot\Omega = 0$.
- GEO north–south station-keeping $\approx 50\ \text{m/s/yr}$ ($\approx 0.75\ \text{km/s}$ over 15 yr).
- SRP $\approx 4.5\ \mu\text{Pa}$; Moon perturbs $\approx 2\times$ the Sun.
Common pitfalls
| Pitfall | Reality |
|---|---|
| "J2 lowers/circularizes the orbit." | J2 rotates $\Omega, \omega$; it does not change $a$ or $e$. Drag does that. |
| "Polar = max nodal regression." | Polar ($90^\circ$) gives zero regression; max is near equatorial. |
| "Sun-synchronous = polar." | SSO is near $98^\circ$ (retrograde) and about timing; polar is about reach. |
| "Station-keeping is a GEO thing." | LEO fights drag, GEO fights luni-solar — different perturbation, same idea. |
| "You can thrust a plane sun-synchronous." | It would cost ~47 km/s/yr; J2 does it for free at the right $i$. |
Mission / astrotools additions this chapter
- MDR: added a station-keeping budget line (dominant perturbation × annual $\Delta v$ × life → propellant → lifetime limit).
orbits.py:j2_nodal_rate(a, e, i),sun_sync_inclination(a), a station-keeping estimator, and a numerical-propagation (Cowell's method) note.