3.2.10 · D3 · Physics › Orbital Mechanics & Astrodynamics › Vis-viva equation v² = GM(2 - r − 1 - a) — derivation
Intuition Yeh page kis liye hai
Vis-viva equation v 2 = GM ( r 2 − a 1 ) ek formula jaisi lagti hai, lekin orbit ki shape (a positive, infinite, ya negative) isse bilkul alag tarike se behave karati hai. Yeh page har tarah ki orbit aur har woh trap dhundti hai jo formula laga sakta hai, aur har ek ko scratch se solve karti hai. End tak aap koi bhi vis-viva scenario nahin dekhenge jo aapne pehle na dekha ho.
Shuru karne se pehle, teen simple-language reminders taaki koi bhi symbol anjana na lage:
Definition Teen symbols jo hum baar baar use karte hain
r = bade body M ke center se aapki current distance (har moment badle).
a = semi-major axis : poori orbit ki sabse lambi width ka aadha — orbit ke liye ek fixed number.
GM = standard gravitational parameter μ : pull kitna strong hai. Earth ke liye G M ⊕ = 398600 km 3 / s 2 .
Hum hamesha km aur seconds mein kaam karenge, isliye GM ko km 3 / s 2 mein rakhne se har answer km/s mein aata hai. Metres aur km mix karna #1 silent error hai — hum isse ek baar km choose karke aur kabhi switch na karke bachate hain.
Har vis-viva problem in cells mein se kisi ek mein hoti hai. Niche ke examples tagged hain us cell se jo woh cover karte hain.
Cell
Orbit type
Sign / value of a
Kya special hai
Example
C1
Circular
a = r (finite > 0 )
r 2 − a 1 = r 1 ; speed constant
Ex 1
C2
Elliptical, perigee par
a > 0 , r = r p < a
sabse fast point
Ex 2
C3
Elliptical, apogee par
a > 0 , r = r a > a
sabse slow point
Ex 2
C4
Elliptical, general point
a > 0 , r p < r < r a
beech ki speed
Ex 3
C5
Parabolic (escape boundary)
a → ∞ , a 1 → 0
ε = 0 exactly
Ex 4
C6
Hyperbolic (unbound)
a < 0
v 2 > 2 GM / r , excess speed
Ex 5
C7
Ulta solve karein a ke liye
given v , r → orbit size dhundho
orbit classify karo
Ex 6
C8
Real-world word problem
mission Δ v
Hohmann transfer
Ex 7
C9
Exam twist / degenerate
r → r at a → r limit, sign trap
galti pakdo
Ex 8
Worked example Ex 1 — Circular orbit (Cell C1)
Ek GPS satellite Earth ke around r = 26560 km radius par circle kar rahi hai. Woh kitni fast ja rahi hai?
Forecast: Pehle guess karo — kya yeh low satellite ki ∼ 7.7 km/s speed se fast hogi ya slow? (Zyada upar orbit → zyada slow.)
Step 1. Pehchano ki circular hai ⇒ a = r .
Yeh step kyun? Ek circle ek aisi orbit hai jis ki "sabse lambi width" uska diameter hi hai, toh uska aadha, a , radius r ke barabar hota hai. Yeh sirf yahi case hai jahan r aur a same number hote hain.
Step 2. a = r ko vis-viva mein substitute karo:
v 2 = GM ( r 2 − r 1 ) = r GM .
Yeh step kyun? Dono fractions collapse ho jaate hain kyunki unka same denominator hai — yeh "circular matlab constant speed" ka algebraic fingerprint hai.
Step 3. Numbers plug karo:
v = 26560 398600 = 15.008 ≈ 3.874 km/s .
Verify: Units: km 3 / s 2 ÷ km = km 2 / s 2 = km/s ✓. Aur 3.87 < 7.7 km/s — LEO se slow, exactly jaisa ek badi orbit mein hona chahiye. ✓
Worked example Ex 2 — Ellipse dono apsides par (Cells C2 & C3)
Ek elliptical orbit ka perigee r p = 7000 km aur apogee r a = 42000 km hai. Perigee (sabse paas) aur apogee (sabse door) par speed dhundho.
Forecast: Kaun zyada hoga, v p ya v a ? (Paas mein whip-fast, door mein drift-slow — toh v p > v a .)
Step 1. Apsides se a nikalo.
Yeh step kyun? Vis-viva ko hamesha a chahiye, aur dono apsides long axis ke opposite ends par hote hain, toh unki distances ka sum poori length 2 a hota hai:
a = 2 r p + r a = 2 7000 + 42000 = 24500 km .
Step 2. Perigee par speed, r = r p = 7000 km:
v p 2 = 398600 ( 7000 2 − 24500 1 ) = 398600 ( 2.8571 × 1 0 − 4 − 4.0816 × 1 0 − 5 ) .
Yeh step kyun? Hum sirf r = r p aur fixed a formula mein daalte hain — kuch aur nahin badalta.
v p 2 = 398600 × 2.4490 × 1 0 − 4 = 97.62 ⇒ v p ≈ 9.88 km/s .
Step 3. Apogee par speed, r = r a = 42000 km, same a :
v a 2 = 398600 ( 42000 2 − 24500 1 ) = 398600 ( 4.7619 × 1 0 − 5 − 4.0816 × 1 0 − 5 ) .
v a 2 = 398600 × 6.803 × 1 0 − 6 = 2.712 ⇒ v a ≈ 1.647 km/s .
Verify: Angular momentum match karna chahiye: r p v p = 7000 × 9.88 = 6.92 × 1 0 4 aur r a v a = 42000 × 1.647 = 6.92 × 1 0 4 . Equal ✓ (yeh Conservation of Angular Momentum action mein hai). Aur v p > v a jaise forecast tha. ✓
Worked example Ex 3 — Ellipse ek general beech ke point par (Cell C4)
Ex 2 wali same orbit (a = 24500 km). Jab satellite r = 20000 km par hai — perigee aur apogee ke beech kahin — tab speed kya hai?
Forecast: 20000 km, 7000 aur 42000 ke beech hai, toh speed 9.88 aur 1.65 km/s ke beech expect karo.
Step 1. Same a = 24500 km reuse karo.
Yeh step kyun? a poori orbit describe karta hai; yeh sirf isliye nahin badalta ki hum move ho gaye. Sirf r badalta hai.
Step 2. r = 20000 substitute karo:
v 2 = 398600 ( 20000 2 − 24500 1 ) = 398600 ( 1.0000 × 1 0 − 4 − 4.0816 × 1 0 − 5 ) .
v 2 = 398600 × 5.918 × 1 0 − 5 = 23.59 ⇒ v ≈ 4.857 km/s .
Verify: 1.65 < 4.86 < 9.88 — perigee/apogee range ke andar jaise forecast tha ✓. Units km/s ✓.
Worked example Ex 4 — Parabolic escape,
a → ∞ limit (Cell C5)
Ek probe r = 7000 km par hai. Isse Earth se just barely escape karne ke liye (infinity tak pahunchne ke liye zero speed ke saath) kitni speed chahiye?
Forecast: Escape speed famously circular speed ki 2 times hoti hai — toh is radius par circular speed ka lagbhag 1.41 × guess karo.
Step 1. "Just barely escape" matlab orbit parabolic hai, yaani a → ∞ .
Yeh step kyun? Ek escaping-but-not-overshooting path ka total specific energy ε = 0 hota hai (bound aur free ke beech ka exact boundary). Kyunki ε = − GM / ( 2 a ) , ε = 0 set karne se a → ∞ force hota hai, toh a 1 → 0 .
Step 2. a 1 term drop karo:
v esc 2 = GM ( r 2 − 0 ) = r 2 GM = 7000 2 × 398600 = 113.89.
v esc ≈ 10.67 km/s .
Verify: Yahan circular speed hai 398600/7000 = 56.94 = 7.546 km/s, aur 2 × 7.546 = 10.67 km/s ✓ — Escape Velocity se exactly match karta hai.
Worked example Ex 5 — Hyperbolic orbit,
a < 0 trap (Cell C6)
Ek interstellar object Sun ke paas se hyperbolic orbit par guzarta hai jisme a = − 50 × 1 0 6 km (negative — yeh unbound hai). r = 150 × 1 0 6 km (1 AU) par yeh kitni fast move kar raha hai? Use karo G M ⊙ = 1.327 × 1 0 11 km 3 / s 2 .
Forecast: Hyperbola ke liye v 2 > 2 GM / r (escape speed aur usse bhi zyada). Toh local escape speed se upar ek speed expect karo.
Step 1. a ko negative rakhho — uska sign mat badlo .
Yeh step kyun? Convention "a < 0 for hyperbolas" exactly wohi hai jo − a 1 ko + ∣ a ∣ 1 bana deti hai, energy add karti hai taaki speed escape se zyada ho. Agar aapne sign "fix" kiya toh physics destroy ho jaati hai.
Step 2. a 1 = − 50 × 1 0 6 1 = − 2.0 × 1 0 − 8 ke saath substitute karo:
v 2 = 1.327 × 1 0 11 ( 150 × 1 0 6 2 − ( − 2.0 × 1 0 − 8 ) ) .
= 1.327 × 1 0 11 ( 1.3333 × 1 0 − 8 + 2.0 × 1 0 − 8 ) = 1.327 × 1 0 11 × 3.3333 × 1 0 − 8 .
v 2 = 4423.3 ⇒ v ≈ 66.51 km/s .
Verify: Local escape speed hai 2 GM / r = 2 × 1.327 × 1 0 11 / ( 150 × 1 0 6 ) = 1769.3 = 42.06 km/s. Hamara 66.5 > 42.06 ✓ — genuinely hyperbolic, jaise forecast tha.
a ke liye ulta solve karo (Cell C7)
Ek spacecraft r = 8000 km par hai aur v = 9.0 km/s move karta hua measure kiya gaya. Semi-major axis a kya hai, aur kya orbit bound hai ya unbound?
Forecast: 9.0 km/s — kya yeh 8000 km par escape se upar hai ya niche? Escape hai 2 × 398600/8000 = 9.98 km/s. Hum usse niche hain, toh orbit bound hona chahiye (positive a ).
Step 1. Vis-viva ko a 1 ke liye rearrange karo:
a 1 = r 2 − GM v 2 .
Yeh step kyun? Humein v aur r pata hai; sirf unknown a hai, toh ise isolate karo. Yeh "orbit diagnose karo" wala move hai.
Step 2. Plug in karo:
a 1 = 8000 2 − 398600 9. 0 2 = 2.5000 × 1 0 − 4 − 2.0321 × 1 0 − 4 = 4.679 × 1 0 − 5 .
a = 4.679 × 1 0 − 5 1 ≈ 21372 km .
Step 3. Classify karo: a > 0 aur finite ⇒ bound ellipse .
Verify: Is a se v recompute karo: v 2 = 398600 ( 2/8000 − 1/21372 ) = 398600 ( 2.5 × 1 0 − 4 − 4.679 × 1 0 − 5 ) = 398600 × 2.032 × 1 0 − 4 = 81.00 , toh v = 9.00 km/s ✓ round-trips. Aur a > 0 bound forecast se match karta hai. Dekho Orbital Elements .
Worked example Ex 7 — Real-world word problem: Hohmann
Δ v (Cell C8)
Ek satellite r 1 = 7000 km par circular orbit mein hai. Hum ise Hohmann Transfer Orbit use karke r 2 = 42000 km par circular orbit mein raise karna chahte hain. Pehla burn Δ v 1 dhundho — low circle chhodne aur transfer ellipse mein enter karne ke liye needed speed boost.
Forecast: Transfer ellipse ka perigee low circle par baithta hai, aur uski perigee speed (Ex 2 se, 9.88 km/s) low circular speed se fast hai. Toh Δ v 1 ek modest positive number hai.
Step 1. r 1 = 7000 par low circular speed (Cell C1 machinery):
v c 1 = 7000 398600 = 56.943 = 7.546 km/s .
Yeh step kyun? Burn se pehle satellite circular speed par move kar rahi hai; yeh hamara starting point hai.
Step 2. Transfer-ellipse perigee speed. Transfer ellipse ka r p = 7000 , r a = 42000 ⇒ a = 24500 km hai — Ex 2 ke same. Toh uski perigee speed v p = 9.883 km/s hai (Ex 2 se).
Yeh step kyun? Burn transfer ellipse ke perigee par hota hai, jo low circle ke saath coincide karta hai, toh hume us ellipse ki perigee speed chahiye.
Step 3. Burn sirf difference poori karta hai:
Δ v 1 = v p − v c 1 = 9.883 − 7.546 = 2.337 km/s .
Verify: Δ v 1 > 0 (upar climb karne ke liye speed-up ) ✓, aur itne transfer ke liye yeh realistic few-km/s figure hai ✓. Units km/s ✓.
Worked example Ex 8 — Exam twist: degenerate / sign trap (Cell C9)
"Ek student r = 30000 km par ek satellite ki speed compute karta hai jisme orbit ka a = 24500 km hai aur paata hai v 2 = 398600 ( 2/30000 − 1/24500 ) . Woh panic karta hai kyunki dono fractions close hain. Phir ek doosra student, jise ek hyperbolic orbit diya gaya hai a = 24500 km ke saath (woh sign bhool gaya), same tarah compute karta hai. Kya galat hai, aur pehle (legitimate) case mein correct speed kya hai?"
Forecast: Pehli calculation theek hai (bound ellipse, r < r a = ? ). Doosra classic sign mistake hai. Guess karo: legit answer couple of km/s hai.
Step 1. Pehla student — check karo ki orbit real hai. Yahan a = 24500 Ex 2 orbit hai, jiska apogee r a = 2 a − r p = 2 ( 24500 ) − 7000 = 42000 km hai. Kyunki 7000 ≤ 30000 ≤ 42000 , point r = 30000 actually orbit par exist karta hai. Achha.
Yeh step kyun? Vis-viva silently ek value return karta hai even un r ke liye jo orbit kabhi reach nahin karta; aapko pehle confirm karna chahiye r p ≤ r ≤ r a , warna "speed" fictitious hai.
Step 2. Compute karo:
v 2 = 398600 ( 30000 2 − 24500 1 ) = 398600 ( 6.6667 × 1 0 − 5 − 4.0816 × 1 0 − 5 ) = 398600 × 2.585 × 1 0 − 5 .
v 2 = 10.305 ⇒ v ≈ 3.210 km/s .
Step 3. Doosre student ki galti: hyperbola ke liye a negative hona chahiye. a = + 24500 use karna ek bound ellipse describe karta hai, woh unbound path nahin jो unka matlab tha — sign cosmetic nahin hai, yeh − a 1 ko energy add karne se energy remove karne mein flip kar deta hai. Fix hai a = − 24500 , jo bahut bada v 2 = 398600 ( 2/30000 + 1/24500 ) deta.
Verify: Pehla case: 3.21 km/s perigee 9.88 aur apogee 1.65 ke beech hai? Nahi — ruko, 3.21 < 1.65 false hai; 1.65 < 3.21 < 9.88 ✓ (yeh beech mein hai, apogee ke kareeb kyunki 30000 far side ke paas hai). Consistent ✓.
Recall Main kaun se cell mein hoon? (quick reflex)
Ek orbit diya gaya, pehle uska type batao ::: circle→a = r ; ellipse→0 < a finite; parabola→a → ∞ ; hyperbola→a < 0 .
v aur r diya gaya, main type kaise dhundhunga ::: 1/ a = 2/ r − v 2 / GM compute karo; a ka sign batata hai (+ bound, ∞ parabolic, − hyperbolic).
Ellipse par computed v par trust karne se pehle kya check karna zaroori hai ::: ki r p ≤ r ≤ r a ho, warna point orbit par nahin hai.
a ka sign orbit ki shape hai
Plus = penned in (ellipse, bound). Infinity = the edge (parabola, escape). Minus = missile (hyperbola, flies off). Kabhi negative a "fix" mat karo — woh aapko sach bata raha hai.