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

3.2.33 · Physics › Orbital Mechanics & Astrodynamics

Earth perfect sphere nahi hai. Yeh equator par thoda bahar nikla hua hai (bulge). Woh chhoti si bulge quietly satellite orbits ko drag karti hai, poore orbital plane ko rotate karti rehti hai. Yeh note exactly derive karta hai ki yeh kitni tezi se hota hai.

The Big Picture (80/20)

20% jo tumhe poori tarah samajhna chahiye: equatorial bulge ek chhota extra potential term $J_2$ add karta hai. Iska time-averaged torque orbit ki size ($a$), shape ($e$), aur tilt ($i$) ko chhoda nahi karta, lekin do orientation angles ko dheere dheere rotate karta hai:

• right ascension of ascending node $\Omega$ (orbit plane swivel karta hai) → nodal precession
• argument of perigee $\omega$ (ellipse apne hi plane mein rotate karta hai)

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

[!intuition] WHY does a bulge cause precession?

Ek perfect sphere bilkul apne center ki taraf pull karta hai → gravity central hoti hai → orbit plane fixed rehta hai (angular momentum $\vec{L}$ conserved rehta hai).

Equatorial bulge equator ke around extra mass add karta hai. Jab satellite equatorial plane ke upar ya neeche hota hai, toh bulge use thoda back toward the equator kheenchta hai — force ka ek component jo center ki taraf point nahi karta. Yeh off-axis pull orbit par ek torque exert karta hai. Kyunki satellite continuously sweep karta rehta hai, torque average hokar ek steady drift ban jaata hai: orbital plane randomly tumble nahi kar sakta, isliye woh ek spinning top ki tarah gravity ke under precesses karta hai.

[!definition] The J2 coefficient

$J_2$ Earth ke gravitational potential ke multipole (spherical-harmonic) expansion ka sabse bada coefficient hai. Yeh oblateness (equatorial bulge) measure karta hai. Earth ke liye $J_2\approx 1.0826\times 10^{-3}$ — chhota, lekin low orbits ke liye yeh baaki saari perturbations par dominate karta hai.

$U(r,\varphi )=-\frac{\mu }{r}\left[1-J_2{\left(\frac{R_E}{r}\right)}^2\frac{3{\sin }^2\varphi -1}{2}-\dots \right]$

• $\mu =G{M}_E$, $R_E$ = equatorial radius, $\varphi$ = geocentric latitude.
• $P_2(\sin \varphi )=\frac{1}{2}(3{\sin }^2\varphi -1)$ degree 2 ka Legendre polynomial hai.
• Leading $-\mu /r$ point-mass (Keplerian) term hai.

[!formula] Derivation of nodal precession $\dot{\Omega }$

Step 1 — The disturbing function $R$

Potential ko Kepler + perturbation mein split karo. Disturbing function woh extra piece hai: $R=\frac{\mu J_2{R}_E^2}{2r^3}\left(3{\sin }^2\varphi -1\right)$ Yeh step kyun? Lagrange's planetary equations $R$ par feed karti hain = point-mass baseline ke upar potential energy per mass.

Step 2 — Express latitude via orbital geometry

Inclination $i$ ke orbit mein argument of latitude $u=\omega +\nu$ par ek satellite ke liye: $\sin \varphi =\sin i\sin u$ Yeh step kyun? Satellite ki latitude sirf iss baat par depend karti hai ki woh orbit mein kitni dur hai ($u$) aur orbit kitna tilted hai ($i$). Yeh orbital triangle ki seedhi spherical trigonometry hai.

Toh: $R=\frac{\mu J_2{R}_E^2}{2r^3}\left(3{\sin }^{2}i{\sin }^{2}u-1\right)$

Step 3 — Average over one orbit

Revolution per perturbations tiny hain; sirf secular (non-repeating) drift matter karta hai. Ek period par average karo. Do facts: $⟨{\sin }^{2}u⟩=\frac{1}{2},⟨\frac{1}{r^3}⟩=\frac{1}{a^3(1-e^2{)}^{3/2}}=\frac{1}{p^3/...}\Rightarrow \frac{1}{a^3(1-e^2{)}^{3/2}}$ Yeh step kyun? Short-period wiggles ek full loop mein cancel ho jaate hain; bacha hua constant part woh hai jo orbit ke baad orbit accumulate hota hai.

Averaging (careful $\nu$-weighting ke saath, $r^{2}d\nu =\sqrt{\mu p}dt$ use karke) mean disturbing function deta hai: $\bar{R}=\frac{n^2{J}_2{R}_E^2}{4(1-e^2{)}^{3/2}}\left(3{\sin }^{2}i-2\right)\cdot (\text{const})$ jahan $n=\sqrt{\mu /a^3}$ mean motion hai aur $p=a(1-e^2)$.

Step 4 — Lagrange's planetary equation for $\Omega$

$\frac{d\Omega }{dt}=\frac{1}{n{a}^2\sqrt{1-e^2}\sin i}\frac{\partial \bar{R}}{\partial i}$ Yeh step kyun? Node motion is baat se driven hoti hai ki jab tum orbit ko tilt karo toh averaged energy kaise change hoti hai — yahi $\partial \bar{R}/\partial i$ ka matlab hai.

Kyunki $\bar{R}\propto (3{\sin }^{2}i-2)$ hai, humein $\partial \bar{R}/\partial i\propto 6\sin i\cos i$ milta hai. $\sin i$ cancel ho jaata hai, aur $\cos i$ bacha rehta hai:

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

jahan $p=a(1-e^2)$, $n=\sqrt{\mu /a^3}$.

Companion result (same machinery, $\partial \bar{R}/\partial e$... actually $\partial /\partial i$ for $\omega$):

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

[!intuition] Reading the two results

Nodal precession $\dot{\Omega }\propto -\cos i$:

• $i=0°$ (equatorial): $\cos i=1$, maximum westward ($\dot{\Omega }<0$).
• $i=90°$ (polar): $\cos i=0$, koi precession nahi.
• $i>90°$ (retrograde): $\cos i<0$, toh $\dot{\Omega }>0$ — plane eastward drift karta hai. Yahi trick hai Sun-synchronous orbits ke peeche: $i\approx 98°$ choose karo taaki $\dot{\Omega }=+0.9856°/\text{day}$ ho, exactly Earth ki Sun ke around orbit → orbit plane Sun ke saath fixed angle rakhta hai.

Apsidal rotation $\dot{\omega }\propto (5{\cos }^{2}i-1)$:

• Zero hota hai jab ${\cos }^{2}i=1/5\Rightarrow i=63.4°$ (ya $116.6°$). Yeh critical inclination perigee ko freeze karta hai — Molniya orbits isse use karti hain taaki apogee northern hemisphere ke upar tikti rahe.

[!example] Worked: precession of the ISS

Given: $a\approx 6.79\times 10^6$ m, $e\approx 0$, $i=51.6°$, $\mu =3.986\times 10^{14}$, $R_E=6.378\times 10^6$, $J_2=1.0826\times 10^{-3}$.

1. Mean motion: $n=\sqrt{\mu /a^3}=\sqrt{3.986\times 10^{14}/(6.79\times 10^6{)}^3}=1.13\times 10^{-3}$ rad/s. Kyun? Base timescale set karta hai; precession $n$ ke saath scale karta hai.
2. $p=a(1-e^2)=a$ (nearly circular). Toh $(R_E/p{)}^2=(6.378/6.79{)}^2=0.882$. Kyun? Yeh "bulge field mein kitna deep" wala factor hai — bade orbits J2 kam feel karte hain.
3. $\dot{\Omega }=-\frac{3}{2}(1.13\times 10^{-3})(1.0826\times 10^{-3})(0.882)\cos 51.6°$ $=-\frac{3}{2}(1.13\times 10^{-3})(1.0826\times 10^{-3})(0.882)(0.622)$ $\approx -1.01\times 10^{-6}$ rad/s.
4. Convert: $\times 86400\times \frac{180}{\pi }\approx -5.0°/\text{day}$. ✔ (ISS node ~5°/day regress karta hai, reality se match karta hai.)

[!example] Worked: designing a Sun-synchronous orbit

Chahiye $\dot{\Omega }=+1.991\times 10^{-7}$ rad/s ($=0.9856°/$day). $a=7.2\times 10^6$ m, $e=0$ ke liye:

1. $n=\sqrt{\mu /a^3}=1.033\times 10^{-3}$ rad/s. Kyun? Formula mein chahiye.
2. Solve karo $\cos i=\frac{\dot{\Omega }}{-\frac{3}{2}n{J}_2(R_E/p{)}^2}$.
3. $(R_E/p{)}^2=(6.378/7.2{)}^2=0.785$. Denominator $=-\frac{3}{2}(1.033\times 10^{-3})(1.0826\times 10^{-3})(0.785)=-1.317\times 10^{-6}$.
4. $\cos i=1.991\times 10^{-7}/(-1.317\times 10^{-6})=-0.151\Rightarrow i=98.7°$. ✔ Negative kyun? Prograde-plus effect ke liye retrograde inclination chahiye.

[!mistake] Steel-manning the classic errors

"J2 orbit ki size aur energy change karta hai." Kyun sahi lagta hai: yeh ek extra force hai, aur forces usually energy/speed change karte hain. Fix: $\dot{a}=\dot{e}=0$ ke secular (averaged) parts; J2 average par orbital energy conserve karta hai aur sirf orientation angles $\Omega ,\omega$ rotate karta hai (plus ek chhota $\dot{M}$). Size aur shape sirf short periods mein wiggle karte hain jo har orbit mein cancel ho jaate hain.

"Polar orbits sabse tezi se precess karte hain." Kyun sahi lagta hai: polar orbits equatorial bulge ko sabse steeply cross karti hain. Fix: polar axis ke baare mein torque ke liye orbit ko us axis ke along angular momentum chahiye — yahi $\cos i$ hai. $i=90°$ par, $\cos i=0$: koi nodal precession nahi bilkul bhi.

"Sign matter nahi karta." Kyun sahi lagta hai: precession toh precession hai. Fix: $\dot{\Omega }$ ka sign ($-\cos i$ se) exactly wahi cheez hai jo retrograde Sun-synchronous orbits possible banata hai. Sign galat karo aur tumhara satellite design fail ho jaata hai.

"$p$ ki jagah $a$ har jagah use karo." Kyun sahi lagta hai: circular orbits ke liye dono equal hain. Fix: eccentric orbits ke liye $p=a(1-e^2)$ use karo; andar chhupa hua $(1-e^2{)}^2$ factor eccentric orbits ke liye precession significantly speed up karta hai.

[!recall]- Feynman: explain to a 12-year-old

Socho Earth thoda squished ball hai — apne belly ke around thoda fatla. Ek satellite iske around looping marble ki tarah hai. Kyunki belly bulge marble ko thoda sideways kheenchta hai jab bhi woh equator ke upar ya neeche hota hai, poora loop ek hi direction mein nahi reh sakta — woh dheere dheere murtaa hai, jaise ek spinning top seedha khade rehne ki jagah dheere dheere circle karta hai. Yeh kitna murtaa hai woh iss baat par depend karta hai ki loop kitna tilted hai: poora seedha upar-neeche (poles ke upar) waala loop bilkul nahi murta, jabki equator ke paas flat loop sabse zyada murta hai. Space engineers ise purpose se use karte hain — woh exact tilt choose karte hain taaki loop itna ghume ki hamesha Sun ki taraf face karta rahe.

[!mnemonic]

"NOde precession ek COsine hai; PErigee Five-COs-squared hai."

• $\dot{\Omega }\sim \cos i$ → "Node—Cosine"
• $\dot{\omega }\sim 5{\cos }^{2}i-1$ → "Perigee—Five"
• Dono minus/plus aur $\frac{3}{2}n{J}_2(R_E/p{)}^2$ front factor carry karte hain. Node $-\frac{3}{2}$ hai, perigee $+\frac{3}{4}$ hai.

Flashcards

Connections

• Two-body problem & Keplerian orbits — woh unperturbed baseline jo yeh modify karta hai.
• Lagrange planetary equations — woh machinery jo $\bar{R}$ ko element rates mein convert karti hai.
• Spherical harmonics of gravity fields — jahan se $J_2,J_3,C_{nm},S_{nm}$ aate hain.
• Sun-synchronous orbits — $\dot{\Omega }$ ka direct engineering use.
• Molniya & critical inclination orbits — $\dot{\omega }=0$ ka application.
• Angular momentum & torque — woh deep reason kyun precession hota hai.
• Legendre polynomials — potential mein $P_2$.

Concept Map