Worked examples — Orbital perturbations — J2 effect (oblateness), derivation of nodal precession

The parent J2 & nodal precession note gave us two master formulas. This page stress-tests them: every inclination sign, every degenerate input, every limiting case, plus real word problems and an exam twist. By the end you will have seen every cell of the scenario grid solved.

The two formulas we will keep reaching for (built in the parent, so we may use them here):

Everything else you need — what $n$, $p$, $i$, $\Omega$, $\omega$ mean — is built in Two-body problem & Keplerian orbits and the parent. The engine that produces these rates is Lagrange planetary equations; the $J_2$ term itself comes from Spherical harmonics of gravity fields and Legendre polynomials.

Constants used everywhere on this page (Earth): $\mu =3.986\times 10^{14}\frac{{\text{m}}^3}{{\text{s}}^2},R_E=6.378\times 10^6\text{m},J_2=1.0826\times 10^{-3}.$

---

The scenario matrix

Every problem J2 can throw at you lands in exactly one of these cells. The worked examples below are labelled with the cell they cover.

Cell	What makes it special	Example
A. Prograde, $0<i<90°$	$\cos i>0$ → $\dot{\Omega }<0$ (westward)	Ex 1 (ISS)
B. Equatorial limit $i=0°$	$\cos i=1$ → maximum precession	Ex 2
C. Polar limit $i=90°$	$\cos i=0$ → $\dot{\Omega }=0$ exactly	Ex 3
D. Retrograde $i>90°$	$\cos i<0$ → $\dot{\Omega }>0$ (eastward) → Sun-sync	Ex 4 (design)
E. Eccentric orbit, $e\ne 0$	must use $p=a(1-e^2)$, not $a$	Ex 5 (Molniya)
F. Critical inclination for $\dot{\omega }$	$5{\cos }^{2}i-1=0$ → apsides freeze	Ex 6
G. Size scaling / limiting	how $\dot{\Omega }$ dies as $a\to \infty$	Ex 7
H. Real-world word problem	"how many days to drift 90°?"	Ex 8
I. Exam twist	given $\dot{\Omega }$, back-solve $i$ with sign trap	Ex 9

---

[!example] Ex 1 — Cell A: prograde, the ISS

Statement. $a=6.79\times 10^6$ m, $e\approx 0$, $i=51.6°$. Find $\dot{\Omega }$ in degrees/day.

Steps.

1. Mean motion $n=\sqrt{\mu /a^3}=\sqrt{3.986\times 10^{14}/(6.79\times 10^6{)}^3}=1.131\times 10^{-3}$ rad/s. Why this step? $\dot{\Omega }$ scales directly with $n$; it sets the base clock speed of the orbit.
2. $p=a(1-e^2)=a$ (circular), so ${\left(R_E/p\right)}^2=(6.378/6.79{)}^2=0.8825$. Why this step? This factor is "how deep in the bulge field" — bigger orbits feel a weaker bulge.
3. $\cos 51.6°=0.6216$.
4. $\dot{\Omega }=-\frac{3}{2}(1.131\times 10^{-3})(1.0826\times 10^{-3})(0.8825)(0.6216)=-1.008\times 10^{-6}$ rad/s. Why this step? Just assembling the master formula.
5. Convert: $\times 86400\times \frac{180}{\pi }=-4.99°/$day.

---

[!example] Ex 2 — Cell B: the equatorial limit

Statement. Same ISS orbit size ($a=6.79\times 10^6$ m, $e=0$) but $i=0°$ (equatorial). Find $\dot{\Omega }$ and compare to Ex 1.

Steps.

1. $n$ and $(R_E/p{)}^2$ are unchanged: $1.131\times 10^{-3}$ rad/s and $0.8825$. Why this step? Only $i$ changed; the size-dependent factors are identical.
2. $\cos 0°=1$.
3. $\dot{\Omega }=-\frac{3}{2}(1.131\times 10^{-3})(1.0826\times 10^{-3})(0.8825)(1)=-1.621\times 10^{-6}$ rad/s.
4. Convert: $=-8.03°/$day.

---

[!example] Ex 3 — Cell C: the polar limit (degenerate zero)

Statement. Same size again, but $i=90°$ (polar). Find $\dot{\Omega }$.

Steps.

1. $\cos 90°=0$. Why this step? Nodal precession is the plane rotating about the polar axis; that needs the orbit's angular momentum to have a component along the polar axis — a polar orbit's angular momentum lies in the equatorial plane, so zero component, zero torque about that axis.
2. $\dot{\Omega }=-\frac{3}{2}n{J}_2(R_E/p{)}^2\cdot 0=0$ exactly.

---

[!example] Ex 4 — Cell D: retrograde design (Sun-synchronous)

Statement. Design a Sun-synchronous orbit at $a=7.2\times 10^6$ m, $e=0$. Required $\dot{\Omega }=+1.991\times 10^{-7}$ rad/s ($=+0.9856°/$day, matching Earth's march around the Sun). Find $i$.

Steps.

1. $n=\sqrt{\mu /a^3}=\sqrt{3.986\times 10^{14}/(7.2\times 10^6{)}^3}=1.0334\times 10^{-3}$ rad/s. Why this step? Appears in the denominator when we back-solve for $\cos i$.
2. $(R_E/p{)}^2=(6.378/7.2{)}^2=0.7849$. Why this step? Bigger orbit than the ISS → weaker bulge factor (0.78 vs 0.88).
3. Prefactor $D=-\frac{3}{2}n{J}_2(R_E/p{)}^2=-\frac{3}{2}(1.0334\times 10^{-3})(1.0826\times 10^{-3})(0.7849)=-1.317\times 10^{-6}$ rad/s.
4. $\cos i=\frac{\dot{\Omega }}{D}=\frac{+1.991\times 10^{-7}}{-1.317\times 10^{-6}}=-0.1512$. Why this step? Inverting the master formula for the one unknown, $i$.
5. $i=\arccos (-0.1512)=98.7°$.

---

[!example] Ex 5 — Cell E: eccentric orbit, the $p$ vs $a$ trap (Molniya)

Statement. A Molniya orbit: $a=2.656\times 10^7$ m, $e=0.74$, $i=63.4°$. Find $\dot{\Omega }$ in degrees/day. Show what you'd get wrongly if you used $a$ instead of $p$.

Steps.

1. $n=\sqrt{\mu /a^3}=\sqrt{3.986\times 10^{14}/(2.656\times 10^7{)}^3}=1.4589\times 10^{-4}$ rad/s. Why this step? Large orbit → slow mean motion.
2. $p=a(1-e^2)=2.656\times 10^7(1-0.74^2)=2.656\times 10^7(0.4524)=1.2016\times 10^7$ m. Why this step? $p$, not $a$, is the correct length scale in the formula (parent's fourth "mistake").
3. $(R_E/p{)}^2=(6.378\times 10^6/1.2016\times 10^7{)}^2=(0.5308{)}^2=0.2818$.
4. $\cos 63.4°=0.4478$.
5. $\dot{\Omega }=-\frac{3}{2}(1.4589\times 10^{-4})(1.0826\times 10^{-3})(0.2818)(0.4478)=-2.991\times 10^{-8}$ rad/s.
6. Convert: $=-0.148°/$day.
7. Wrong version (using $a$): $(R_E/a{)}^2=(6.378/26.56{)}^2=0.05766$, giving ${\dot{\Omega }}_{\text{wrong}}=-6.12\times 10^{-9}$ rad/s $=-0.0303°/$day.

---

[!example] Ex 6 — Cell F: critical inclination (apsides freeze)

Statement. At what inclination does $\dot{\omega }=0$? Verify it kills apsidal drift for the Molniya orbit of Ex 5.

Steps.

1. Set $5{\cos }^{2}i-1=0\Rightarrow {\cos }^{2}i=1/5$. Why this step? $\dot{\omega }=0$ is exactly the bracket vanishing; the prefactor is never zero for a real orbit.
2. $\cos i=\pm \sqrt{1/5}=\pm 0.4472\Rightarrow i=63.43°$ or $i=116.57°$.
3. Check the Molniya value $i=63.4°$: $5{\cos }^2(63.4°)-1=5(0.4478{)}^2-1=5(0.2005)-1=0.0026\approx 0$.

---

[!example] Ex 7 — Cell G: size scaling, the far-orbit limit

Statement. Show how $\dot{\Omega }$ dies as $a\to \infty$. Compute the ratio of $\dot{\Omega }$ at $a=4R_E$ vs $a=2R_E$ (both circular, same $i$).

Steps.

1. Write the $a$-dependence: $n=\sqrt{\mu /a^3}\propto a^{-3/2}$, and $(R_E/p{)}^2=(R_E/a{)}^2\propto a^{-2}$ (circular). Why this step? Isolate every factor's power of $a$ so the ratio is clean.
2. Product $\propto a^{-3/2-2}=a^{-7/2}$.
3. Ratio $\frac{\dot{\Omega }(2R_E)}{\dot{\Omega }(4R_E)}={\left(\frac{4R_E}{2R_E}\right)}^{7/2}=2^{7/2}=11.31$.

---

[!example] Ex 8 — Cell H: real-world word problem

Statement. How many days does the ISS take to drift its ascending node a full $90°$ around the equator?

Steps.

1. Use $|\dot{\Omega }|=4.99°/$day from Ex 1. Why this step? Same orbit; we already computed its rate.
2. Time $=\frac{90°}{4.99°/\text{day}}=18.04$ days. Why this step? Constant secular rate → linear accumulation, so just divide angle by rate.

---

[!example] Ex 9 — Cell I: exam twist (given rate, sign trap)

Statement. A satellite at $a=7.0\times 10^6$ m, $e=0$ is observed with $\dot{\Omega }=-6.5°/$day. Find its inclination. Then state which of the two mathematical solutions is physical and why the sign matters.

Steps.

1. Convert the target: $\dot{\Omega }=-6.5°/\text{day}\times \frac{\pi }{180}\times \frac{1}{86400}=-1.313\times 10^{-6}$ rad/s. Why this step? The formula lives in SI rad/s.
2. $n=\sqrt{\mu /a^3}=\sqrt{3.986\times 10^{14}/(7.0\times 10^6{)}^3}=1.0796\times 10^{-3}$ rad/s.
3. $(R_E/p{)}^2=(6.378/7.0{)}^2=0.8302$.
4. Prefactor $D=-\frac{3}{2}(1.0796\times 10^{-3})(1.0826\times 10^{-3})(0.8302)=-1.4557\times 10^{-6}$ rad/s.
5. $\cos i=\dot{\Omega }/D=(-1.313\times 10^{-6})/(-1.4557\times 10^{-6})=0.9020$. Why this step? Inverting the master formula; both negatives cancel → positive cosine → prograde, as forecast.
6. $i=\arccos (0.9020)=25.6°$.

---

[!recall]- One-line summary of the whole matrix