3.2.10Orbital Mechanics & Astrodynamics

Vis-viva equation v² = GM(2 - r − 1 - a) — derivation

1,797 words8 min readdifficulty · medium6 backlinks

WHAT is it?

  • Circle: r=av2=GM/rr = a \Rightarrow v^2 = GM/r (constant speed).
  • Ellipse: aa finite, vv changes with rr (fast at perigee, slow at apogee).
  • Parabola (escape): av2=2GM/ra \to \infty \Rightarrow v^2 = 2GM/r.
  • Hyperbola: a<0v2>2GM/ra < 0 \Rightarrow v^2 > 2GM/r (excess speed).

WHY does total energy depend ONLY on aa?


HOW to derive it — from first principles

Step 1 — Write the total energy

The specific energy (per unit mass) of an orbiting body: ε=12v2kinetic    GMrpotential\varepsilon = \underbrace{\frac{1}{2}v^2}_{\text{kinetic}} \;-\; \underbrace{\frac{GM}{r}}_{\text{potential}}

Why this step? Gravitational PE per unit mass is GM/r-GM/r (zero at infinity, negative when bound). KE per unit mass is 12v2\tfrac12 v^2. Their sum is conserved because gravity does no non-conservative work.

Step 2 — Evaluate ε\varepsilon at perigee and apogee

At the apsides (perigee and apogee) the velocity is purely tangential (perpendicular to rr), so v=vpv = v_p or v=vav = v_a cleanly. Distances: rp=a(1e),ra=a(1+e)r_p = a(1-e), \qquad r_a = a(1+e)

Why this step? Energy is the same everywhere, but the two apsides give two easy equations we can solve together. Choosing them removes the radial velocity component.

\varepsilon = \frac12 v_p^2 - \frac{GM}{r_p} = \frac12 v_a^2 - \frac{GM}{r_a} \tag{1}

Step 3 — Use conservation of angular momentum

At the apsides vr\vec v \perp \vec r, so the specific angular momentum is simply h = r_p v_p = r_a v_a \;\Rightarrow\; v_a = v_p\frac{r_p}{r_a} \tag{2}

Why this step? Angular momentum is also conserved. With two conserved quantities (energy + LL) and two unknowns (vp,vav_p, v_a) we can solve everything.

Step 4 — Solve for vp2v_p^2

Substitute (2) into (1): 12vp2GMrp=12vp2rp2ra2GMra\frac12 v_p^2 - \frac{GM}{r_p} = \frac12 v_p^2\frac{r_p^2}{r_a^2} - \frac{GM}{r_a}

Group the vp2v_p^2 terms and the GMGM terms: 12vp2(1rp2ra2)=GM(1rp1ra)\frac12 v_p^2\left(1 - \frac{r_p^2}{r_a^2}\right) = GM\left(\frac{1}{r_p} - \frac{1}{r_a}\right)

Factor the left bracket as a difference of squares and the right side over a common denominator: 12vp2(rarp)(ra+rp)ra2=GMrarprpra\frac12 v_p^2 \cdot \frac{(r_a - r_p)(r_a + r_p)}{r_a^2} = GM\cdot\frac{r_a - r_p}{r_p r_a}

Cancel (rarp)(r_a - r_p): 12vp2ra+rpra2=GMrpra\frac12 v_p^2\,\frac{r_a + r_p}{r_a^2} = \frac{GM}{r_p r_a}

So vp2=2GMrarp(ra+rp)v_p^2 = \frac{2GM\, r_a}{r_p (r_a + r_p)}

Why this step? Pure algebra, but the difference-of-squares cancellation is what makes it collapse to something clean.

Step 5 — Plug in ra+rp=2ar_a + r_p = 2a

Since rp=a(1e)r_p = a(1-e) and ra=a(1+e)r_a = a(1+e), their sum is 2a2a. Then vp2=2GMrarp2a=GMraarpv_p^2 = \frac{2GM\,r_a}{r_p \cdot 2a} = \frac{GM\,r_a}{a\,r_p}

Step 6 — Get the energy, then generalize

Now compute ε\varepsilon at perigee: ε=12vp2GMrp=GMra2arpGMrp=GMrp(ra2a1)\varepsilon = \frac12 v_p^2 - \frac{GM}{r_p} = \frac{GM\,r_a}{2a\,r_p} - \frac{GM}{r_p} = \frac{GM}{r_p}\left(\frac{r_a}{2a} - 1\right)

Use ra=2arpr_a = 2a - r_p: ε=GMrp(2arp2a1)=GMrp(rp2a)=GM2a\varepsilon = \frac{GM}{r_p}\left(\frac{2a - r_p}{2a} - 1\right) = \frac{GM}{r_p}\left(-\frac{r_p}{2a}\right) = -\frac{GM}{2a}

This is the big result: ε=GM2a\boxed{\varepsilon = -\dfrac{GM}{2a}} — energy depends only on aa. ✅

Step 7 — Read off vis-viva

Since ε\varepsilon is constant, at any point rr with speed vv: 12v2GMr=GM2a\frac12 v^2 - \frac{GM}{r} = -\frac{GM}{2a}

Multiply by 2 and rearrange: v2=GM(2r1a)\boxed{v^2 = GM\left(\frac{2}{r} - \frac{1}{a}\right)}

Figure — Vis-viva equation v² = GM(2 - r − 1 - a) — derivation

Worked Examples


Common Mistakes


Recall Feynman: explain to a 12-year-old

Imagine a ball on a stretchy string going around you in an oval. When it swings in close it whips around fast; when it drifts far out it goes slow — but the total "oomph" (speed-energy + height-energy) never changes. The vis-viva equation is just a recipe: tell me how far the ball is right now and how big its whole loop is, and I'll tell you exactly how fast it's moving. The "size of the loop" is the only thing that sets the total oomph.


Flashcards

What does aa represent in the vis-viva equation?
The semi-major axis of the orbit (a fixed property of the whole orbit).
What does rr represent in the vis-viva equation?
The instantaneous distance from the central body's center to the orbiting body.
State the specific orbital energy in terms of aa.
ε=GM/(2a)\varepsilon = -GM/(2a) — depends only on orbit size.
Derive escape speed from vis-viva.
Set aa\to\infty so 1/a01/a\to0: v2=2GM/rv^2 = 2GM/r, hence vesc=2GM/rv_{esc}=\sqrt{2GM/r}.
Vis-viva for a circular orbit reduces to?
v=GM/rv=\sqrt{GM/r}, since r=ar=a.
Why is energy the same at perigee and apogee?
Gravity is conservative; total mechanical energy is conserved everywhere on the orbit.
What two conservation laws are used in the derivation?
Conservation of energy and conservation of angular momentum.
For a hyperbolic orbit, what is the sign of aa?
Negative, giving v2>2GM/rv^2 > 2GM/r (excess hyperbolic speed).
Relationship between escape and circular speed at the same radius?
vesc=2vcircv_{esc} = \sqrt{2}\,v_{circ}.
What is GMGM called and what symbol is often used?
The standard gravitational parameter, μ\mu.

Connections

Concept Map

implies

sums to

sums to

evaluated at

gives two eqs

relates v_a v_p

combine

yields

substitute back

special cases

Gravity is conservative

Specific energy const

KE half v squared

PE minus GM over r

Angular momentum const

Apsides perigee apogee

Energy at r_p and r_a

Solve for v_p squared

E equals minus GM over 2a

Vis-viva equation

Orbit types by a

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, vis-viva equation ka matlab simple hai: kisi bhi orbit me, agar tumhe pata hai ki abhi tum central body se kitni door ho (ye hai rr) aur tumhari poori orbit kitni badi hai (ye hai semi-major axis aa), to tum exact speed vv nikal sakte ho. Formula hai v2=GM(2/r1/a)v^2 = GM(2/r - 1/a). Yaad rakho — rr har second change hota hai, lekin aa poori orbit ka fixed size hai. Dono alag cheez hain, confuse mat hona.

Iske peeche ka asli jaadu hai energy conservation. Gravity ek conservative force hai, isliye total energy (kinetic 12v2\frac12 v^2 plus potential GM/r-GM/r) poori orbit me kabhi nahi badalti. Derivation me hum perigee aur apogee dono jagah energy aur angular momentum equal lagate hain, thoda algebra karte hain, aur magic se nikalta hai ε=GM/(2a)\varepsilon = -GM/(2a) — yaani total energy sirf orbit ke size par depend karti hai, shape ya position par nahi. Bas isi ko rearrange karo to vis-viva mil jaata hai.

Practical me ye bahut kaam ka hai. Satellite launch karna ho, Hohmann transfer se Mars bhejna ho, ya escape velocity nikalni ho — sab vis-viva se. Jaise escape ke liye orbit infinitely badi (aa \to \infty) kar do, 1/a1/a zero ho jaata hai, aur v=2GM/rv = \sqrt{2GM/r} aa jaata hai. Circle ke liye r=ar=a daalo, v=GM/rv=\sqrt{GM/r}. Ek hi formula, saari orbits cover. Isliye ise yaad rakhna sabse zyada "return on investment" wali cheez hai astrodynamics me.

Go deeper — visual, from zero

Test yourself — Orbital Mechanics & Astrodynamics

Connections