3.2.12 · Physics › Orbital Mechanics & Astrodynamics
Har orbit mein ek conserved quantity hoti hai jise specific angular momentum h kehte hain (angular momentum per unit mass ). Yeh encode karta hai ki orbit kitni "wide" hai. Relation h = GM p us conserved number ko orbit ki geometric size se semi-latus rectum p ke through link karta hai. Toh: shape measure karo (p ), aur turant spin pata chal jaayega (h ) — aur vice versa.
Definition Specific angular momentum
Ek body of mass m jo bahut badi mass M ke around orbit kar rahi hai, uska angular momentum hai L = r × m v . Hum m se divide kar dete hain taaki ek per-mass quantity mile jo orbiting object ki mass par depend na kare:
h ≡ r × v , h = ∣ h ∣
Units: m 2 /s (note: not kg m 2 /s — mass pehle hi divide ho chuki hai).
m se divide kyun karte hain? Kyunki two-body problem mein chhoti mass equation of motion se cancel ho jaati hai (r ¨ = − r 2 GM r ^ ). Dynamics sirf r , v par depend karti hai, isliye natural conserved quantities bhi per-mass hoti hain.
Intuition Koi torque nahi ⟹
h constant
Gravity seedha r ke along point karti hai (central force). Torque hai τ = r × F , aur r × ( kuch bhi ∥ r ) = 0 . Koi torque nahi ⟹ angular momentum mein koi change nahi.
Scratch se derivation.
d t d h = d t d ( r × v ) = = 0 v × v + r × a
Yeh step kyun? Cross products ke liye product rule; v × v = 0 kyunki ek vector apne aap ke parallel hota hai.
Acceleration purely radial hai: a = − r 2 GM r ^ . Toh
r × a = − r 2 GM ( r × r ^ ) = 0
Yeh step kyun? r = r r ^ , aur r ^ × r ^ = 0 . Isliye
d t d h = 0 ⇒ h = const.
Constant h ⟹ orbit ek fixed plane mein rehti hai (motion h ke perpendicular rehti hai) aur Kepler's 2nd law hold karta hai (d t d A = 2 h ).
Humein orbit equation chahiye. Two-body problem ke conic-section solution se shuru karte hain (neeche ek baar derive kiya gaya hai).
Polar coordinates use karte hain. Radial equation of motion hai
r ¨ − r θ ˙ 2 = − r 2 GM
Kyun? Polar form mein standard acceleration; RHS gravity hai.
Angular momentum deta hai h = r 2 θ ˙ (r × v ka z -component). Toh θ ˙ = h / r 2 .
Trick: u = 1/ r substitute karo aur variable θ mein change karo.
r ˙ = d θ d r θ ˙ = − u 2 1 d θ d u ⋅ h u 2 = − h d θ d u
Yeh step kyun? r = 1/ u ⇒ d θ d r = − u − 2 d θ d u , aur θ ˙ = h u 2 .
Dobara differentiate karo:
r ¨ = − h d θ 2 d 2 u θ ˙ = − h 2 u 2 d θ 2 d 2 u
Radial equation mein substitute karo (jahan r θ ˙ 2 = u 1 ( h u 2 ) 2 = h 2 u 3 ):
− h 2 u 2 u ′′ − h 2 u 3 = − GM u 2
− h 2 u 2 se divide karo:
u ′′ + u = h 2 GM
Yeh step kyun? Yeh ek clean linear ODE hai — Binet's equation.
Iska solution: particular u p = GM / h 2 , homogeneous A cos ( θ − θ 0 ) :
u = h 2 GM ( 1 + e cos θ ) ⇒ r = 1 + e c o s θ h 2 / GM
r = 1 + e cos θ p se compare karo:
p = GM h 2 ⟹ h = GM p
Intuition Semi-latus rectum HOTA kya hai?
Orbit equation mein θ = 9 0 ∘ set karo: cos 9 0 ∘ = 0 toh r = p . Toh p woh orbital radius hai jo major axis ke perpendicular, focus se measure kiya jaata hai — ellipse ki "half-width" focus par. Bada p ⟹ motii orbit ⟹ bada h .
Key checkpoints (sab r = 1 + e c o s θ p se):
Periapsis (θ = 0 ): r p = 1 + e p
Apoapsis (θ = 18 0 ∘ ): r a = 1 − e p
Harmonic-mean relation: p 2 = r p 1 + r a 1 .
Worked example 1 — Circular Low-Earth orbit
G M ⊕ = 3.986 × 1 0 14 m 3 / s 2 , r = 7000 km = 7 × 1 0 6 m , circular hai toh p = r .
h = GM p = 3.986 × 1 0 14 × 7 × 1 0 6 = 2.79 × 1 0 21 ≈ 5.28 × 1 0 10 m 2 /s
Yeh step kyun? Circle ⟹ e = 0 ⟹ p = a = r , toh bas radius plug in karo.
Check: v = GM / r = 7547 m/s , aur h = v r = 7547 × 7 × 1 0 6 = 5.28 × 1 0 10 ✓ (match karta hai).
Worked example 2 — Periapsis speed aur radius se
Ek probe periapsis par hai jahan r p = 8000 km , v p = 9.0 km/s . Periapsis par velocity purely tangential hoti hai, isliye h = r p v p .
h = 8 × 1 0 6 × 9000 = 7.2 × 1 0 10 m 2 /s
Phir p = GM h 2 = 3.986 × 1 0 14 ( 7.2 × 1 0 10 ) 2 = 1.30 × 1 0 7 m = 13 , 000 km .
Yeh step kyun? Periapsis par v ⊥ r hota hai matlab ∣ r × v ∣ = r v (koi sin θ loss nahi).
Worked example 3 — Eccentricity nikalo
Same probe: r p = p / ( 1 + e ) se, e = r p p − 1 = 8000 13000 − 1 = 0.625 .
Yeh step kyun? Orbit equation ka periapsis form use karo, ab jab p pata hai.
h ki units kg·m²/s hain."
Kyun sahi lagta hai: ordinary angular momentum L = m v r mein mass hoti hai. Fix: h specific hai (per unit mass), h = L / m , units m 2 /s . Poora point hi mass-independent rehna hai.
Common mistake "Orbit mein har jagah
h = r v hota hai."
Kyun sahi lagta hai: peri/apoapsis par sach hai toh general lagta hai. Fix: generally h = r v sin ϕ jahan ϕ flight-path angle hai r aur v ke beech. Sirf apsides par ϕ = 9 0 ∘ hota hai, isliye sin ϕ = 1 .
Common mistake "Ellipse ke liye
p = a use karo."
Kyun sahi lagta hai: circles ke liye sach hai, aur a famous "size" parameter hai. Fix: p = a ( 1 − e 2 ) . a tab hi use karo jab e = 0 .
h 2 = GM p toh h = GM p /2 ya kuch aisa."
Kyun sahi lagta hai: sloppy algebra. Fix: h = GM p — square root, half nahi.
Recall Feynman: ek 12-saal ke bachche ko samjhao
Ek ball ko dhage se bandho aur ghuma do. Ball ka "spin amount" wahi rehta hai jab tak tum use sideways push na karo — gravity ek aise dhage ki tarah hai jo sirf centre ki taraf kheenchta hai, kabhi sideways nahi, isliye planet ka spin amount hamesha same rehta hai. Woh spin amount hi h hai. Bada loop matlab bada spin amount. Formula bas yeh kehta hai: batao loop kitna wide hai (woh hai p ) aur central pull kitna strong hai (woh hai GM ), aur main spin bataa deta hoon: unhe multiply karo aur square root lo.
"He Got My Pull" → h = GM p . He(h) Got(√) My(GM) Pull(p).
Aur geometry ke liye: "p is the perpendicular peek" — p woh radius hai jab tum 9 0 ∘ par seedha bahar dekhte ho.
Kisi bhi central force mein h constant kyun rehta hai?
h ki units kya hain, aur kg kyun nahi?
Orbit ke kin points par exactly h = r v hota hai?
Ek line of reasoning mein p = h 2 / GM derive karo.
Specific angular momentum h kya hai? h = r × v , angular momentum per unit mass, units
m 2 /s .
Orbital motion mein h conserved kyun hai? Gravity central hai (
r ke parallel), isliye torque
r × F = 0 , hence
d h / d t = 0 .
p ke terms mein h batao.h = GM p , jahan
p semi-latus rectum hai.
p ko a aur e se express karo.p = a ( 1 − e 2 ) , toh
h = GM a ( 1 − e 2 ) .
Semi-latus rectum geometrically kya hai? True anomaly 9 0 ∘ par orbital radius: r = p jab cos θ = 0 .
h = r v exactly kab hold karta hai?Periapsis aur apoapsis par, jahan
v ⊥ r (flight-path angle
= 9 0 ∘ ).
h , r , v ke beech general relation?h = r v sin ϕ jahan
ϕ ,
r aur
v ke beech angle hai.
Radius r ki circular orbit mein h ? h = GM r = v c i r c r , kyunki
p = a = r .
Orbit mein u = 1/ r kaunsa ODE satisfy karta hai? Binet's equation u ′′ + u = GM / h 2 .
h ke zariye Kepler's 2nd law?d A / d t = h /2 = const (equal times mein equal areas).
Two-body problem — jahan se orbit equation r = p / ( 1 + e cos θ ) aati hai.
Kepler's Second Law — d A / d t = h /2 constant h ka direct consequence hai.
Specific orbital energy — h ke saath pair karta hai orbit ko fully fix karne ke liye (a energy se, e dono se).
Eccentricity vector — e = GM v × h − r ^ directly h use karta hai.
Flight-path angle — h = r v sin ϕ mein sin ϕ explain karta hai.
Vis-viva equation — v 2 = GM ( 2/ r − 1/ a ) orbit ke along speeds dene ke liye h se combine hota hai.
Angular momentum L = r x m v
Specific ang. mom. h = r x v
Kepler 2nd law dA/dt = h/2
Orbit equation r = p / 1+e cos theta