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

3.2.33 · D3 · Physics › Orbital Mechanics & Astrodynamics › Orbital perturbations — J2 effect (oblateness), derivation o

Parent J2 & nodal precession note ne humein do master formulas di hain. Yeh page unhe stress-test karta hai: har inclination sign, har degenerate input, har limiting case, saath mein real word problems aur ek exam twist bhi. Is page ke end tak tumne scenario grid ka har cell solve hua dekh liya hoga.

Do formulas jo hum baar baar use karenge (parent mein derive ki gayi hain, toh hum yahan seedha use kar sakte hain):

Baaki jo kuch bhi chahiye — $n$, $p$, $i$, $\Omega$, $\omega$ ka matlab — woh Two-body problem & Keplerian orbits aur parent note mein diya hua hai. Yeh rates jo engine produce karta hai woh hai Lagrange planetary equations; $J_2$ term khud Spherical harmonics of gravity fields aur Legendre polynomials se aata hai.

Is poore page par use hone wale constants (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}.$

---

Scenario matrix

J2 ka har problem in cells mein se exactly ek mein fit hota hai. Neeche ke worked examples us cell ke naam se label kiye gaye hain jo woh cover karta hai.

Cell	Kya cheez special hai	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$	$p=a(1-e^2)$ use karna zaroori hai, $a$ nahi	Ex 5 (Molniya)
F. Critical inclination for $\dot{\omega }$	$5{\cos }^{2}i-1=0$ → apsides freeze	Ex 6
G. Size scaling / limiting	$\dot{\Omega }$ kaise $a\to \infty$ par khatam hota hai	Ex 7
H. Real-world word problem	"90° drift hone mein kitne din lagenge?"	Ex 8
I. Exam twist	diya gaya $\dot{\Omega }$, back-solve $i$ sign trap ke saath	Ex 9

---

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

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

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. Yeh step kyun? $\dot{\Omega }$ directly $n$ ke saath scale karta hai; yeh orbit ki base clock speed set karta hai.
2. $p=a(1-e^2)=a$ (circular), toh ${\left(R_E/p\right)}^2=(6.378/6.79{)}^2=0.8825$. Yeh step kyun? Yeh factor batata hai "bulge field mein kitna andar" — bade orbits kamzor bulge feel karte hain.
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. Yeh step kyun? Bas master formula assemble kar rahe hain.
5. Convert: $\times 86400\times \frac{180}{\pi }=-4.99°/$day.

---

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

Statement. Same ISS orbit size ($a=6.79\times 10^6$ m, $e=0$) lekin $i=0°$ (equatorial). $\dot{\Omega }$ nikalo aur Ex 1 se compare karo.

Steps.

1. $n$ aur $(R_E/p{)}^2$ unchanged hain: $1.131\times 10^{-3}$ rad/s aur $0.8825$. Yeh step kyun? Sirf $i$ badla hai; size-dependent factors same hain.
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: polar limit (degenerate zero)

Statement. Same size phir se, lekin $i=90°$ (polar). $\dot{\Omega }$ nikalo.

Steps.

1. $\cos 90°=0$. Yeh step kyun? Nodal precession matlab plane ka polar axis ke around rotate karna; iske liye orbit ke angular momentum ka polar axis ke saath ek component hona chahiye — polar orbit ka angular momentum equatorial plane mein hota hai, toh woh component zero, toh us axis ke around torque bhi zero.
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. $a=7.2\times 10^6$ m, $e=0$ par ek Sun-synchronous orbit design karo. Required $\dot{\Omega }=+1.991\times 10^{-7}$ rad/s ($=+0.9856°/$day, jo Earth ke Sun ke around march se match karta hai). $i$ nikalo.

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. Yeh step kyun? Jab hum $\cos i$ ke liye back-solve karte hain toh denominator mein aata hai.
2. $(R_E/p{)}^2=(6.378/7.2{)}^2=0.7849$. Yeh step kyun? ISS se bada orbit → kamzor 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$. Yeh step kyun? Ek unknown $i$ ke liye master formula ko invert kar rahe hain.
5. $i=\arccos (-0.1512)=98.7°$.

---

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

Statement. Ek Molniya orbit: $a=2.656\times 10^7$ m, $e=0.74$, $i=63.4°$. $\dot{\Omega }$ degrees/day mein nikalo. Dikhao ki agar $p$ ki jagah $a$ use karte toh kya galat answer aata.

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. Yeh step kyun? Bada 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. Yeh step kyun? Formula mein correct length scale $p$ hai, $a$ nahi (parent ki chauthi "galti").
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. Galat version ($a$ use karke): $(R_E/a{)}^2=(6.378/26.56{)}^2=0.05766$, jisse ${\dot{\Omega }}_{\text{wrong}}=-6.12\times 10^{-9}$ rad/s $=-0.0303°/$day milta.

---

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

Statement. Kis inclination par $\dot{\omega }=0$ hota hai? Ex 5 ke Molniya orbit ke liye verify karo ki yeh apsidal drift khatam kar deta hai.

Steps.

1. Set $5{\cos }^{2}i-1=0\Rightarrow {\cos }^{2}i=1/5$. Yeh step kyun? $\dot{\omega }=0$ exactly bracket ke vanish hone se hota hai; prefactor kisi real orbit ke liye kabhi zero nahi hota.
2. $\cos i=\pm \sqrt{1/5}=\pm 0.4472\Rightarrow i=63.43°$ ya $i=116.57°$.
3. Molniya value $i=63.4°$ check karo: $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, far-orbit limit

Statement. Dikhao ki $a\to \infty$ par $\dot{\Omega }$ kaise khatam hota hai. $a=4R_E$ vs $a=2R_E$ par $\dot{\Omega }$ ka ratio nikalo (dono circular, same $i$).

Steps.

1. $a$-dependence likho: $n=\sqrt{\mu /a^3}\propto a^{-3/2}$, aur $(R_E/p{)}^2=(R_E/a{)}^2\propto a^{-2}$ (circular). Yeh step kyun? Har factor ka $a$ ka power isolate karo taaki ratio clean ho.
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. ISS ko apna ascending node equator ke around poore $90°$ drift karne mein kitne din lagte hain?

Steps.

1. Ex 1 se $|\dot{\Omega }|=4.99°/$day use karo. Yeh step kyun? Same orbit hai; humne already uski rate compute ki hai.
2. Time $=\frac{90°}{4.99°/\text{day}}=18.04$ din. Yeh step kyun? Constant secular rate → linear accumulation, toh bas angle ko rate se divide karo.

---

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

Statement. $a=7.0\times 10^6$ m, $e=0$ wale satellite par $\dot{\Omega }=-6.5°/$day observe kiya gaya hai. Uski inclination nikalo. Phir batao ki do mathematical solutions mein se kaun sa physical hai aur sign kyun matter karta hai.

Steps.

1. Target convert karo: $\dot{\Omega }=-6.5°/\text{day}\times \frac{\pi }{180}\times \frac{1}{86400}=-1.313\times 10^{-6}$ rad/s. Yeh step kyun? Formula SI rad/s mein hai.
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$. Yeh step kyun? Master formula invert kar rahe hain; dono negatives cancel → positive cosine → prograde, jaise forecast kiya tha.
6. $i=\arccos (0.9020)=25.6°$.

---

[!recall]- Poore matrix ka one-line summary