Vis-Viva Equation

The vis-viva equation (named after Gottfried Wilhelm Leibniz's Latin term vis viva, or "living force," for twice the mechanical kinetic energy of a moving body) is one of the fundamental governing equations of celestial mechanics and spacecraft astrodynamics. In the unperturbed Keplerian two-body problem, it expresses the conservation of specific total mechanical energy (ε) by directly relating the instantaneous orbital speed v of a body at radial distance r from a central gravitating primary to the semi-major axis (a) of its conic-section trajectory:

v² = G M ( 2 / r − 1 / a ) = μ ( 2 / r − 1 / a )

where G = 6.67430 × 10⁻¹¹ m³ kg⁻¹ s⁻² is the Newtonian gravitational constant, M is the mass of the central body (or M + m for comparable masses), and μ = GM is the standard gravitational parameter (for the Sun, μ_☉ = 1.3271244 × 10²⁰ m³ s⁻²; for Earth, μ_⊕ = 3.986004418 × 10¹⁴ m³ s⁻²).

Because the vis-viva equation holds across all four non-degenerate conic sections—circles (e = 0), ellipses (0 < e < 1), parabolas (e = 1, a → ∞), and hyperbolas (e > 1, a < 0)—and depends only on the scalar radial distance r and the semi-major axis a (independent of orbital eccentricity e or flight-path angle γ at a given r and a), it serves as the workhorse relationship for calculating spacecraft orbital velocities, Hohmann transfer orbit impulses, planetary escape velocity thresholds, aerobraking energy dissipation, and hyperbolic gravity assist encounters across the Solar System.

🪐 Interactive Orbital Mechanics Simulation: Apply the Vis-Viva Equation in real time to compute periapsis speeds, transfer burns, and hyperbolic trajectories in Gravity Assist.

Conservation of Specific Orbital Energy (ε)

Consider a small body or spacecraft of mass m moving in the inverse-square gravitational field of a spherically symmetric primary body of mass M ≫ m. The total mechanical energy E of the orbiting body is the sum of its translational kinetic energy (E_t = ½*mv*²) and its Newtonian gravitational potential energy (U = −*GMm*/r, referenced to U(∞) = 0 at infinite separation). Dividing E by the spacecraft mass m defines the conserved specific mechanical energy (ε, measured in J kg⁻¹ or m² s⁻²):

ε = E / m = ½ v² − μ / r = constant

To evaluate the constant ε in terms of the geometric semi-major axis a of an elliptical orbit (0 ≤ e < 1), we examine the two apsides—periapsis (r_p = *a*[1 − e]) and apoapsis (r_a = *a*[1 + e])—where the radial velocity component vanishes (dr/dt = 0) and the velocity vector v is perpendicular to the position vector r. Because specific orbital angular momentum h = r × v is conserved along the orbit, its scalar magnitude at periapsis and apoapsis satisfies:

h = r_p v_p = r_a v_a ⇒ v_a = ( r_p / r_a ) v_p

Equating the specific mechanical energy ε at periapsis and apoapsis yields:

½ v_p² − μ / r_p = ½ v_a² − μ / r_a = ½ ( r_p / r_a )² v_p² − μ / r_a

Factoring out v_p² and using the difference of squares (r_a² − r_p²) = (r_a − r_p)(r_a + r_p) together with the geometric apsidal sum r_p + r_a = 2*a* gives:

½ v_p² [ ( r_a² − r_p² ) / r_a² ] = μ [ ( r_a − r_p ) / ( r_p r_a ) ] ⇒ ½ v_p² = μ r_a / [ r_p ( r_p + r_a ) ] = μ r_a / ( 2 a r_p )

Substituting ½*v*_p² back into the specific energy expression at periapsis reveals that specific orbital energy depends only on the central gravitational parameter μ and the semi-major axis a:

ε = μ r_a / ( 2 a r_p ) − μ / r_p = μ ( r_a − 2*a* ) / ( 2 a r_p ) = μ ( −*r*_p ) / ( 2 a r_p ) = − μ / ( 2 a )

Equating the instantaneous specific energy ε = ½*v*² − μ/r at any point r to −*μ*/(2*a*) and multiplying through by 2 recovers the vis-viva equation:

v² = μ ( 2 / r − 1 / a )

Circular, Elliptical, Parabolic, and Hyperbolic Solutions

By varying the semi-major axis a relative to the instantaneous radius r, the vis-viva equation partitions all two-body orbital trajectories into four dynamical classes:

1. Circular Orbit (r = a, e = 0, ε = −*μ*/[2*r*] < 0): Setting a = r everywhere along the trajectory gives a constant orbital speed equal to the local circular velocity (v_circ):

v_circ = √( μ / r ) For Low Earth Orbit at altitude h = 400 km (r = 6,778 km), v_circ = 7.67 km s⁻¹; for Earth around the Sun at r = 1.000 AU, v_circ = 29.78 km s⁻¹.

1. Elliptical Orbit (a > r/2, 0 < e < 1, ε = −*μ*/[2*a*] < 0): A bound closed orbit in which translational kinetic energy is traded back and forth with gravitational potential energy as the body oscillates between periapsis r_p = a(1 − e) and apoapsis r_a = a(1 + e). Evaluating the vis-viva equation at the apsides yields the maximum and minimum speeds along an ellipse:

v_p = √[ ( μ / a ) · ( 1 + e ) / ( 1 − e ) ] , v_a = √[ ( μ / a ) · ( 1 − e ) / ( 1 + e ) ]

1. Parabolic Escape Trajectory (a → ∞, e = 1, ε = 0): In the limit a → ∞, the term 1/a vanishes, yielding the local escape velocity (v_esc) required to reach r → ∞ with zero residual velocity:

v_esc = √( 2 μ / r ) = √2 · v_circ

1. Hyperbolic Flyby Trajectory (a < 0, e > 1, ε = +½*v*_∞² > 0): For an unbound hyperbolic orbit, the geometric semi-major axis is negative (a = −|a| = −*μ*/v_∞²). Substituting −1/a = +*v*_∞²/μ into the vis-viva equation gives:

v² = 2 μ / r + v_∞² = v_esc(r)² + v_∞² where v_∞ is the hyperbolic excess speed (the asymptotic speed remaining as r → ∞), and C₃ = v_∞² = 2*ε* is the characteristic launch energy.

Orbital Maneuvers, the Counterintuitive Braking Paradox, and the Oberth Effect

The vis-viva equation reveals a counterintuitive property of orbital mechanics known as the orbital speed–altitude paradox: in a stable circular orbit (v_circ = √[μ/r]), an spacecraft orbiting at a lower altitude (r₁ < r₂) travels faster than a spacecraft in a higher circular orbit, even though its total mechanical energy (ε = −*μ*/[2*r*]) is more negative (lower). Conversely, when a satellite in the upper thermosphere experiences atmospheric drag during aerobraking, aerodynamic drag removes total mechanical energy (dε/dt < 0), which decreases the semi-major axis a; as the orbit contracts to a smaller radius r, gravitational potential energy −*μ*/r drops twice as much as total energy ε = −*μ*/(2*a*) (by the virial theorem, ⟨*E*_t⟩ = −½⟨*U*⟩), causing the satellite's actual orbital speed v to increase as it spirals downward!

Furthermore, differentiating the vis-viva equation at a fixed radial distance r (dr = 0 during an impulsive rocket burn) yields the sensitivity of semi-major axis a to an instantaneous tangential velocity impulse Δ*v*:

v dv = ( μ / [ 2 a² ] ) da ⇒ da = ( 2 a² v / μ ) dv

Because da/dv is directly proportional to the instantaneous orbital speed v, firing a rocket thruster at periapsis (where r = r_p is smallest and v = v_p is largest) maximizes the change in semi-major axis a and orbital energy ε per unit of propellant expended—the mathematical basis of the Oberth effect.

Why does the vis-viva equation not include the orbital eccentricity (e) or the direction of the velocity vector?

Because Newtonian gravity (U = −*GMm*/r) is a spherically symmetric central potential, potential energy depends only on the scalar radial distance r = |r|, and translational kinetic energy (½*mv*²) depends only on the scalar magnitude of velocity v = |v|, not on the flight-path angle γ between r and v. Two spacecraft at the same radial distance r with the same speed v but launched at different elevation angles γ₁ ≠ γ₂ possess the exact same specific mechanical energy ε = ½*v*² − μ/r = −*μ*/(2*a*) and therefore trace out ellipses with the identical semi-major axis a (and identical orbital period T = 2*π*√[a³/μ]), even though their angular momenta (h = rv cos γ) and eccentricities (e = √[1 + 2*εh*²/μ²]) differ!

How do you use the vis-viva equation to calculate the Δ*v* required for a Hohmann transfer from Low Earth Orbit to Geostationary Orbit?

Consider a spacecraft in a circular Low Earth Orbit (LEO) at altitude h₁ = 400 km (r₁ = 6,778 km, where v_circ,1 = √[μ_⊕/r₁] = 7.669 km s⁻¹) transferring to Geostationary Earth Orbit (GEO) at r₂ = 42,164 km (v_circ,2 = √[μ_⊕/r₂] = 3.075 km s⁻¹). The elliptical Hohmann transfer orbit connecting r₁ and r₂ has semi-major axis a_H = ½(r₁ + r₂) = 24,471 km. Substituting r = r₁ and a = a_H into the vis-viva equation gives the required transfer speed at perigee: v_p = √[μ_⊕(2/r₁ − 1/a_H)] = 10.152 km s⁻¹ (requiring a first burn Δ*v*₁ = v_p − v_circ,1 = 2.483 km s⁻¹). At apogee (r = r₂), the vis-viva equation gives v_a = √[μ_⊕(2/r₂ − 1/a_H)] = 1.608 km s⁻¹, requiring a circularization burn Δ*v*₂ = v_circ,2 − v_a = 1.467 km s⁻¹ (total Δ*v* = 3.950 km s⁻¹).

When does the two-body vis-viva equation require corrections in real spaceflight?

The classical vis-viva equation assumes two point masses (or spherically symmetric masses) isolated from external forces. It requires perturbation corrections in four main regimes: (1) near oblate planets where equatorial bulges (J₂ harmonic, such as Earth's J₂ = 1.0826 × 10⁻³) add a quadrupole potential term; (2) in multi-body regimes near Lagrange points where the Jacobi integral in the rotating three-body frame replaces two-body energy conservation; (3) during low-thrust continuous propulsion (solar sail or Hall-effect ion thrusters) or atmospheric aerobraking where a(t) varies continuously; and (4) close to compact relativistic objects where general relativistic post-Newtonian corrections apply.

1. Bate, R. R., Mueller, D. D., White, J. E., & Saylor, W. W. (2020). Fundamentals of Astrodynamics (2nd ed.). Dover Publications. ISBN 978-0486497044
2. Vallado, D. A. (2013). Fundamentals of Astrodynamics and Applications (4th ed.). Microcosm Press. ISBN 978-1881883180
3. Curtis, H. D. (2019). Orbital Mechanics for Engineering Students (4th ed.). Butterworth-Heinemann / Elsevier. https://doi.org/10.1016/C2016-0-02107-1
4. Murray, C. D., & Dermott, S. F. (1999). Solar System Dynamics. Cambridge University Press. https://doi.org/10.1017/CBO9781139174817