Visual walkthrough — Orbital perturbations — J2 effect (oblateness), derivation of nodal precession

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

Yeh page parent note ka headline result scratch se rebuild karta hai. Koi bhi symbol use nahi hoga jab tak uska plain-word meaning aur picture dono na ho. Hum ek equation chase kar rahe hain, aur end tak aap uska har piece dekh lenge:

$\dot{\Omega }=-\frac{3}{2}n{J}_2{\left(\frac{R_E}{p}\right)}^2\cos i$

Chalte chalte destination ko pehle words mein decode karte hain. $\Omega$ (Greek capital "Omega") ke upar dot ka matlab hai "har second $\Omega$ kitni tezi se badal raha hai" — ek rate. $\Omega$ khud ek angle hai jo batata hai ki tilted orbit ring kis taraf point kar rahi hai. Toh poori equation ek hi sawaal hai: orbit ring kitni tezi se swivel karti hai? Neeche har ek step right-hand side ke har symbol ko earn karta hai, ek step at a time.

---

Step 1 — Orbit kya hai, aur hum kya change kar sakte hain?

KYA. Ek satellite Earth ke around ek ellipse par loop karta hai. Us loop ko space mein pin karne ke liye kuch numbers chahiye. Unme se chaar size, shape aur tilt describe karte hain; do loop kis taraf point karta hai yeh batate hain.

KYUN. Isse pehle ki hum puchen "yeh kitni tezi se swivel karta hai," hume pata hona chahiye ki kya swivel karta hai. Bulge, nikalta hai, size, shape aur tilt freeze kar deta hai lekin do pointing angles rotate karta hai. Toh pehle hume saare choh se milna hoga.

PICTURE. Figure mein, flat gray disk Earth ka equatorial plane hai — Earth ke mote belly ka plane. Blue ellipse orbit hai, us plane se upar ki taraf tilted hai.

Jis line par tilted orbit equatorial plane ko slice karti hai woh line of nodes hai. Jis point par satellite us line ko northward jaate hue cross karta hai woh ascending node hai — orange dot. Jab $\Omega$ change hota hai, nodes ki poori line rotate hoti hai. Woh rotation nodal precession hai. Two-body problem & Keplerian orbits se sab kuch ye chhe numbers deta hai; yahan hum poochte hain ki bulge inke saath kya karta hai.

---

Step 2 — Satellite kahan baitha hai? Geometry se Latitude

KYA. Hum $\varphi$, geocentric latitude introduce karte hain — Earth ke center se dekhe gaye satellite ka equatorial plane ke upar ya neeche ka angle. Hum dikhayenge ki $\sin \varphi =\sin i\sin u,u=\omega +\nu .$

KYUN. Bulge sirf yeh care karta hai ki satellite equator ke upar hai ya neeche — exactly yahi $\varphi$ measure karta hai. Hume $\varphi$ ko un chezon ke terms mein likhna hai jo hamare paas already hain ($i$, aur hum loop mein kitni door hain). Jab tak hum yeh nahi jaante ki $\varphi$ kya depend karta hai, hum kisi formula mein $\varphi$ use nahi kar sakte.

PICTURE. Figure celestial sphere par baitha spherical triangle dikhata hai: ek corner ascending node par, ek satellite par, ek equator par us point par jo satellite ke seedha "neeche" hai.

Yahan $u=\omega +\nu$ argument of latitude hai — ek single angle jo measure karta hai "ascending node se loop ke around kitni door." Figure padh ke:

$\underset{\text{height above equator}}{\underbrace{\sin \varphi }}=\underset{\text{how steep the ramp is}}{\underbrace{\sin i}}\cdot \underset{\text{how far up the ramp we’ve walked}}{\underbrace{\sin u}}$

Yeh pure spherical trigonometry hai — wohi triangle jo aap Two-body problem & Keplerian orbits mein milte ho.

---

Step 3 — Bulge gravity mein ek extra term add karta hai

KYA. Ek perfect sphere $-\mu /r$ potential se khinchta hai, jahan $\mu =G{M}_E$ (Earth ki mass times gravitational constant) aur $r$ center se distance hai. Bulge ek correction add karta hai. Poora: $U(r,\varphi )=-\frac{\mu }{r}\left[1-J_2{\left(\frac{R_E}{r}\right)}^2\frac{3{\sin }^2\varphi -1}{2}\right].$

KYUN. Ek lumpy Earth se gravity correction shells add karke build hoti hai — Spherical harmonics of gravity fields expansion. $J_2$ sabse badi correction hai: yeh equatorial bulge encode karta hai. Hum sirf yahi ek rakhte hain kyunki low orbits ke liye yeh baaki sab par dominate karta hai.

PICTURE. Figure round Earth (perfectly circular pull, gray) aur true squashed Earth (equator ke paas extra outward pull, orange arrows) ko contrast karta hai.

Term by term: $U=-\frac{\mu }{r}\left[1\underset{\begin{array}{c}\text{bulge size,}\\ \approx 1.08\times 10^{-3}\end{array}}{\underbrace{-J_2}}\underset{\begin{array}{c}\text{fades fast}\\ \text{with distance}\end{array}}{\underbrace{{\left(\frac{R_E}{r}\right)}^2}}\underset{\begin{array}{c}\text{latitude shape}\\ =P_2(\sin \varphi )\end{array}}{\underbrace{\frac{3{\sin }^2\varphi -1}{2}}}\right]$

• $R_E$ = Earth ka equatorial radius. Ratio $(R_E/r{)}^2$ kehta hai ki bulge tab sabse zyada matter karta hai jab aap close skim karte ho.
• $\frac{1}{2}(3{\sin }^2\varphi -1)$ Legendre polynomial $P_2$ hai (dekho Legendre polynomials). Yeh poles ke paas positive aur equator ke paas negative hai — "fat belly, flat poles" ka mathematical fingerprint.

---

Step 4 — Geometry ko disturbing function mein daalo

KYA. Step 2 ka $\sin \varphi =\sin i\sin u$ ko $R$ mein substitute karo: $R=\frac{\mu J_2{R}_E^2}{2r^3}\left(3{\sin }^{2}i{\sin }^{2}u-1\right).$

KYUN. $R$ latitude $\varphi$ use karke likha gaya tha, lekin latitude hamare chhe orbit numbers mein se ek nahi hai — $i$ aur $u$ hain. Lagrange planetary equations ki machinery use karne ke liye hume $R$ purely orbital quantities mein express karna hoga. Yeh substitution woh translation hai.

PICTURE. Figure ek loop mein $R$ plot karta hai (jaise $u$ $0\to 360°$ run karta hai): yeh per orbit do baar upar neeche oscillate karta hai, kyunki ${\sin }^{2}u$ ke do humps hain.

$R=\frac{\mu J_2{R}_E^2}{2r^3}(3\underset{\text{tilt, fixed}}{\underbrace{{\sin }^{2}i}}\underset{\text{position, wiggling}}{\underbrace{{\sin }^{2}u}}-1)$

Structure notice karo: ek part fixed tilt $i$ par depend karta hai, doosra rapidly changing position $u$ par. Woh split agale step ki key hai.

---

Step 5 — Ek loop par average karo: sirf woh rakho jo accumulate hota hai

KYA. Wiggling quantities ko ek full orbit par unke averages se replace karo. Do facts kaam karte hain: $⟨{\sin }^{2}u⟩=\frac{1}{2},⟨\frac{1}{r^3}⟩=\frac{1}{a^3(1-e^2{)}^{3/2}}.$ Result hai mean disturbing function $\bar{R}$, $(3{\sin }^{2}i-2)$ ke proportional.

KYUN. Ek loop mein wiggles orbit ko ek taraf push karte hain, phir wapas pull karte hain — woh cancel ho jaate hain. Jo bach jaata hai, woh non-cancelling constant, woh hai jo loop ke baad loop ke baad quietly build up hota hai. Woh leftover secular drift hai. Hum sirf drift chase karte hain, toh hum average karte hain.

PICTURE. Figure Step 4 ki wiggly $R(u)$ curve dikhata hai jisme uske mean se hoke ek flat dashed line hai — dashed line ke upar aur neeche ke wiggles equal area rakhte hain, toh sirf dashed level survive karta hai.

$⟨{\sin }^{2}u⟩=\frac{1}{2}\Rightarrow (3{\sin }^{2}i\cdot \frac{1}{2}-1)=\frac{1}{2}(3{\sin }^{2}i-2).$

• $⟨\cdot ⟩$ matlab "ek orbit par average."
• $1/r^3$ average $p=a(1-e^2)$, semi-latus rectum laata hai — in purposes ke liye orbit ka "effective radius."
• $n=\sqrt{\mu /a^3}$ mean motion hai: average angular speed, yaani $2\pi$ divided by orbital period. Yeh isliye appear karta hai kyunki time par averaging positions ko weight karta hai by kitni der satellite wahan ruka.

$\bar{R}\propto \frac{n^2{J}_2{R}_E^2}{(1-e^2{)}^{3/2}}(3{\sin }^{2}i-2).$

---

Step 6 — Woh lever jo tilt ko swivel mein convert karta hai

KYA. Node ke liye Lagrange's planetary equation hai $\frac{d\Omega }{dt}=\frac{1}{n{a}^2\sqrt{1-e^2}\sin i}\frac{\partial \bar{R}}{\partial i}.$

KYUN. Yeh angular-momentum bookkeeping disguise mein hai. Node tab swivel karta hai jab averaged energy $\bar{R}$ change hoti agar aap orbit ko thoda tilt karte. Woh "change per tilt" partial derivative $\partial \bar{R}/\partial i$ hai — yeh answer deta hai "bulge energy inclination ke liye kitni sensitive hai?" Derivative sahi tool hai kyunki hum chahte hain ek rate of change with respect to a small nudge, jo exactly wahi hai jo derivative measure karta hai.

PICTURE. Figure orbit ring ko ek gyroscope ki tarah dikhata hai: bulge ek torque (red arrow) apply karta hai dono spin axis aur equator ke perpendicular, ring ko sideways precess karne par force karta hai tumble karne ki jagah.

Ab differentiate karo. Kyunki $\bar{R}\propto (3{\sin }^{2}i-2)$: $\frac{\partial }{\partial i}(3{\sin }^{2}i-2)=6\sin i\cos i.$

Wapas slot karo aur cancellation dekho: $\dot{\Omega }\propto \frac{1}{\sin i}\cdot \underset{\text{from }{\partial }_i\bar{R}}{\underbrace{6\sin i\cos i}}=6\cos i.$

Denominator mein $\sin i$ cancel ho jaata hai derivative ke $\sin i$ se, aur sirf $\cos i$ bachta hai. Har constant carefully collect karke boxed result milta hai:

$\boxed{\dot{\Omega }=-\frac{3}{2}n{J}_2{\left(\frac{R_E}{p}\right)}^2\cos i}$

$\dot{\Omega }=\underset{\text{sign + size}}{\underbrace{-\frac{3}{2}}}\underset{\text{orbit speed}}{\underbrace{n}}\underset{\text{bulge}}{\underbrace{J_2}}\underset{\text{bulge field mein kitne deep}}{\underbrace{{\left(\frac{R_E}{p}\right)}^2}}\underset{\text{tilt factor}}{\underbrace{\cos i}}$

Right side par har symbol ab dekha ja chuka hai: $n$ loop timing se (Step 5), $J_2$ aur $R_E/p$ bulge potential se (Steps 3, 5), aur sab-se-important $\cos i$ lever cancellation se (yeh step).

---

Step 7 — Har case, koi gap nahi: $\cos i$ padhna

KYA. Har possible tilt se guzro aur $\dot{\Omega }$ ka sign aur size check karo.

KYUN. Contract: reader ko kabhi aisa scenario nahi milna chahiye jo humne nahi dikhaya. $\cos i$ sign change karta hai, toh physics flip ho jaati hai — hume sab cover karna hai.

PICTURE. Figure $\dot{\Omega }$ versus inclination $i$ ko $0°$ se $180°$ tak plot karta hai, har regime mark karta hai.

Inclination $i$	$\cos i$	$\dot{\Omega }$	Matlab
$0°$ (equatorial, prograde)	$+1$	sabse zyada negative	ring westward sabse fast swivel karti hai
$0°$ aur $90°$ ke beech	positive	negative	westward, $i$ badhne par slow hota hai
$90°$ (polar)	$0$	zero	ring bilkul swivel nahi karti
$90°$ aur $180°$ ke beech	negative	positive	ring eastward swivel karti hai
$180°$ (equatorial, retrograde)	$-1$	sabse zyada positive	eastward, fastest

---

Step 8 — Saathi: perigee rotation aur critical angle

KYA. Wohi machinery, lekin $\Omega$ ki jagah $\omega$ ke liye differentiate karke deta hai $\dot{\omega }=\frac{3}{4}n{J}_2{\left(\frac{R_E}{p}\right)}^2\left(5{\cos }^{2}i-1\right).$

KYUN. Hum ise include karte hain kyunki yeh ek doosra special inclination chupaata hai, aur kyunki yeh dikhata hai ki method general hai — ek one-trick derivation nahi.

PICTURE. Figure tilt factor $(5{\cos }^{2}i-1)$ plot karta hai; yeh $i=63.4°$ aur $116.6°$ par zero cross karta hai.

$5{\cos }^{2}i-1=0\Rightarrow {\cos }^{2}i=\frac{1}{5}\Rightarrow i=63.4°\text{ or }116.6°.$

Is critical inclination par perigee rotate karna band kar deta hai — yeh "freeze" ho jaata hai. Molniya orbits exactly isi ko exploit karte hain taaki unka apogee ghanton tak northern hemisphere ke upar hang kare. Dekho Molniya & critical inclination orbits.

---

Ek-picture summary

Yeh single figure poori derivation compress karta hai: bulge → off-center pull (Step 3) → tilted ring par torque (Step 6) → ek loop par average (Step 5) → $\cos i$ se scale kiya gaya node ka steady swivel (Steps 6–7).

Recall Poore walkthrough ki Feynman retelling

Ek thoda squashed Earth imagine karo — beech mein mota. Ek satellite ek tilted ring par uske around ride karta hai. Jab bhi satellite us mote belly ke upar ya neeche drift karta hai, extra belly-mass use thoda sideways khinchti hai, equator ki taraf wapas (Step 3). Woh sideways tug Earth ke center ki taraf seedha nahi hai, toh yeh ring ko twist karta hai — ek torque (Step 6). Ek full loop mein, zyaadatar twisting cancel ho jaati hai; sirf ek tiny steady leftover bachta hai (Step 5). Woh leftover slowly ring ki pointing direction turn karta hai — node angle $\Omega$ drift karta hai. Kitni tezi se? Yeh orbit ki speed $n$, bulge $J_2$ ki size, aap kitne close skim karte ho ($R_E/p$), aur sabse badhkar $\cos i$ ke through tilt par depend karta hai (Step 6). Seedhe poles ke upar ($i=90°$): bilkul turn nahi. Equator ke paas flat: sabse fast turn, westward. Vertical se aage lean karo ($i>90°$): turn eastward flip ho jaata hai — aur agar aap us eastward turn ko Earth ki Sun ke around yearly loop se match karo, aapka orbit Sun ke saath hamesha same face rakhta hai (Step 7). Wohi idea, ring ki jagah ellipse ke perigee par apply karo, toh ek magic tilt $63.4°$ milta hai jahan ellipse rotate karna band kar deta hai (Step 8). Ek bulge, do beautiful design tricks.

Recall Quick self-test

Kaunsa orbit tilt zero nodal precession deta hai? ::: $i=90°$, polar orbit, kyunki $\cos 90°=0$. Derivative se $\sin i$ $\dot{\Omega }$ mein kyun nahi appear karta? ::: Yeh Lagrange's equation ke denominator mein $\sin i$ ko cancel kar deta hai, pure $\cos i$ chhodta hai. Retrograde ($i>90°$) orbits Sun-synchronicity ke liye special kyun hain? ::: $\cos i<0$ $\dot{\Omega }$ ko positive (eastward) flip karta hai, jise Earth ki Sun ke around yearly motion se match kiya ja sakta hai.