# Vis-viva equation

In astrodynamics, the vis-viva equation, also referred to as orbital-energy-conservation equation, is one of the equations that model the motion of orbiting bodies. It is the direct result of the law of conservation of energy, which requires that the sum of kinetic and potential energy is constant at all points along the orbit.

Vis viva (Latin for "live force") is a term from the history of mechanics, and it survives in this sole context. It represents the principle that the difference between the aggregate work of the accelerating forces of a system and that of the retarding forces is equal to one half the vis viva accumulated or lost in the system while the work is being done.

## Equation

For any Kepler orbit (elliptic, parabolic, hyperbolic, or radial), the vis-viva equation[1] is as follows:

${\displaystyle v^{2}=GM\left({{2 \over {r}}-{1 \over {a}}}\right)}$

where:

The product of GM can also be expressed using the Greek letter μ.

## Derivation

In the vis-viva equation the mass m of the orbiting body (e.g., a spacecraft) is taken to be negligible in comparison to the mass M of the central body (e.g., the Earth). In the specific cases of an elliptical or circular orbit, the vis-viva equation may be readily derived from conservation of energy and momentum.

Specific total energy is constant throughout the orbit. Thus, using the subscripts a and p to denote apoapsis (apogee) and periapsis (perigee), respectively,

${\displaystyle \epsilon ={\frac {v_{a}^{2}}{2}}-{\frac {GM}{r_{a}}}={\frac {v_{p}^{2}}{2}}-{\frac {GM}{r_{p}}}}$

Rearranging,

${\displaystyle {\frac {v_{a}^{2}}{2}}-{\frac {v_{p}^{2}}{2}}={\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}}$

Recalling that for an elliptical orbit (and hence also a circular orbit) the velocity and radius vectors are perpendicular at apoapsis and periapsis, conservation of angular momentum requires ${\displaystyle h=r_{p}v_{p}=r_{a}v_{a}={\text{constant}}}$, thus ${\displaystyle v_{p}={\frac {r_{a}}{r_{p}}}v_{a}}$:

${\displaystyle {\frac {1}{2}}\left(1-{\frac {r_{a}^{2}}{r_{p}^{2}}}\right)v_{a}^{2}={\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}}$
${\displaystyle {\frac {1}{2}}\left({\frac {r_{p}^{2}-r_{a}^{2}}{r_{p}^{2}}}\right)v_{a}^{2}={\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}}$

Isolating the kinetic energy at apoapsis and simplifying,

${\displaystyle {\frac {1}{2}}v_{a}^{2}=\left({\frac {GM}{r_{a}}}-{\frac {GM}{r_{p}}}\right)\left({\frac {r_{p}^{2}}{r_{p}^{2}-r_{a}^{2}}}\right)}$
${\displaystyle {\frac {1}{2}}v_{a}^{2}=GM\left({\frac {r_{p}-r_{a}}{r_{a}r_{p}}}\right)\left({\frac {r_{p}^{2}}{r_{p}^{2}-r_{a}^{2}}}\right)}$
${\displaystyle {\frac {1}{2}}v_{a}^{2}=GM\left({\frac {r_{p}}{r_{a}(r_{p}+r_{a})}}\right)}$

From the geometry of an ellipse, ${\displaystyle 2a=r_{p}+r_{a}}$ where a is the length of the semimajor axis. Thus,

${\displaystyle {\frac {1}{2}}v_{a}^{2}=GM\left({\frac {2a-r_{a}}{r_{a}(2a)}}\right)}$

Substituting this into our original expression for specific orbital energy,

${\displaystyle \epsilon ={\frac {v_{a}^{2}}{2}}-{\frac {GM}{r_{a}}}=GM\left({\frac {2a-r_{a}}{2ar_{a}}}\right)-{\frac {GM}{r_{a}}}}$
${\displaystyle \epsilon =GM\left({\frac {2a-r_{a}}{2ar_{a}}}-{\frac {1}{r_{a}}}\right)=-{\frac {GM}{2a}}}$

Thus, ${\displaystyle \epsilon =-{\frac {GM}{2a}}}$ and the vis-viva equation may be written

${\displaystyle {\frac {v^{2}}{2}}-{\frac {GM}{r}}=-{\frac {GM}{2a}}}$

or

${\displaystyle v^{2}=GM\left({\frac {2}{r}}-{\frac {1}{a}}\right)}$.

Therefore, the conserved angular momentum L = mh is derived below using -

${\displaystyle r_{a}=a\left(1+e\right)}$
${\displaystyle r_{p}=a\left(1-e\right)}$
${\displaystyle r_{a}r_{p}=b^{2}}$

where a, b, e and m are semi-major axis, semi-minor axis, eccentricity and mass

${\displaystyle v_{a}^{2}=GM\left({\frac {2}{r_{a}}}-{\frac {1}{a}}\right)={\frac {GM}{a}}\left({\frac {2}{1+e}}-{\frac {1}{1}}\right)={\frac {GM}{a}}\left({\frac {1-e}{1+e}}\right)={\frac {GM}{a}}\left({\frac {r_{p}}{r_{a}}}\right)={\frac {GM}{a}}\left({\frac {b^{2}}{r_{a}^{2}}}\right)}$
${\displaystyle L=mh=v_{a}r_{a}m=mb\,{\sqrt {\frac {GM}{a}}}}$

## Practical applications

Given the total mass and the scalars r and v at a single point of the orbit, one can compute r and v at any other point in the orbit.[2]

Given the total mass and the scalars r and v at a single point of the orbit, one can compute the specific orbital energy ${\displaystyle \epsilon \,\!}$, allowing an object orbiting a larger object to be classified as having not enough energy to remain in orbit, hence being "suborbital" (a ballistic missile, for example), having enough energy to be "orbital", but without the possibility to complete a full orbit anyway because it eventually collides with the other body, or having enough energy to come from and/or go to infinity (as a meteor, for example).

## References

1. Tom Logsdon, Orbital Mechanics: theory and applications, John Wiley & Sons, 1998
2. For the three-body problem there is hardly a comparable vis-viva equation: conservation of energy reduces the larger number of degrees of freedom by only one.