Vis-viva equation · Lablypedia

Related on Lably

Through Orbital mechanics, Gravitational field and Work (physics)

Lab · Class 11Orbital Mechanics Lab · Class 8–11Gravitation Lab · Class 11Black Holes Lab · Class 8–12Satellite Systems

Vis viva (Latin for "living force") is a term from the history of mechanics and the name given to the orbital equation originally derived by Isaac Newton. It represents the principle that the difference between the total 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.

Formulation

For any Keplerian orbit (elliptic, parabolic, hyperbolic, or radial), the vis-viva equation is as follows: $v^{2} = G M (\frac{2}{r} - \frac{1}{a})$ where:

v is the relative speed of the two bodies
r is the distance between the two bodies' centers of mass
a is the length of the semi-major axis (a > 0 for ellipses, a = ∞ or 1/a = 0 for parabolas, and a < 0 for hyperbolas)
G is the gravitational constant
M is the mass of the central body

The product of GM can also be expressed as the standard gravitational parameter using the Greek letter μ.

2 more chapters are free with a Lably account.

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; and
the specific orbital energy $ε$ , 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).

The formula for escape velocity can be obtained from the Vis-viva equation by taking the limit as $a$ approaches $\infty$ : $v_{e}^{2} = G M (\frac{2}{r} - 0) \to v_{e} = \sqrt{\frac{2 G M}{r}}$ For a given orbital radius, the escape velocity will be $\sqrt{2}$ times the orbital velocity.

Derivation for elliptic orbits (0 ≤ eccentricity < 1)

Specific total energy, $ε = \frac{v^{2}}{2} - \frac{G M}{r}$ is constant throughout the orbit. Using the subscripts a for apoapsis (apogee) and p for periapsis (perigee), the constant energy at two points gives: $\frac{v_{a}^{2}}{2} - \frac{G M}{r_{a}} = \frac{v_{p}^{2}}{2} - \frac{G M}{r_{p}}$

Rearranging, $\frac{v_{a}^{2}}{2} - \frac{v_{p}^{2}}{2} = \frac{G M}{r_{a}} - \frac{G M}{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 specific angular momentum $h = r_{p} v_{p} = r_{a} v_{a} = constant$ , thus $v_{p} = \frac{r_{a}}{r_{p}} v_{a}$ : $\frac{1}{2} (1 - \frac{r_{a}^{2}}{r_{p}^{2}}) v_{a}^{2} = \frac{G M}{r_{a}} - \frac{G M}{r_{p}}$ $\frac{1}{2} (\frac{r_{p}^{2} - r_{a}^{2}}{r_{p}^{2}}) v_{a}^{2} = \frac{G M}{r_{a}} - \frac{G M}{r_{p}}$

Isolating the kinetic energy at apoapsis and simplifying, $\begin{matrix} \frac{1}{2} v_{a}^{2} & = (\frac{G M}{r_{a}} - \frac{G M}{r_{p}}) \cdot \frac{r_{p}^{2}}{r_{p}^{2} - r_{a}^{2}} \\ \frac{1}{2} v_{a}^{2} & = G M (\frac{r_{p} - r_{a}}{r_{a} r_{p}}) \frac{r_{p}^{2}}{r_{p}^{2} - r_{a}^{2}} \\ \frac{1}{2} v_{a}^{2} & = G M \frac{r_{p}}{r_{a} (r_{p} + r_{a})} \end{matrix}$

From the geometry of an ellipse, $2 a = r_{p} + r_{a}$ where a is the length of the semimajor axis. Thus, $\frac{1}{2} v_{a}^{2} = G M \frac{2 a - r_{a}}{r_{a} (2 a)} = G M (\frac{1}{r_{a}} - \frac{1}{2 a}) = \frac{G M}{r_{a}} - \frac{G M}{2 a}$

Substituting this into our original expression for specific orbital energy, $ε = \frac{v^{2}}{2} - \frac{G M}{r} = \frac{v_{p}^{2}}{2} - \frac{G M}{r_{p}} = \frac{v_{a}^{2}}{2} - \frac{G M}{r_{a}} = - \frac{G M}{2 a}$

Thus, $ε = - \frac{G M}{2 a}$ and the vis-viva equation may be written $\frac{v^{2}}{2} - \frac{G M}{r} = - \frac{G M}{2 a}$ or $v^{2} = G M (\frac{2}{r} - \frac{1}{a})$

Therefore, the conserved angular momentum L = mh can be derived using $r_{a} + r_{p} = 2 a$ and $r_{a} r_{p} = b^{2}$ , where a is semi-major axis and b is semi-minor axis of the elliptical orbit, as follows: $v_{a}^{2} = G M (\frac{2}{r_{a}} - \frac{1}{a}) = \frac{G M}{a} (\frac{2 a - r_{a}}{r_{a}}) = \frac{G M}{a} (\frac{r_{p}}{r_{a}}) = \frac{G M}{a} {(\frac{b}{r_{a}})}^{2}$ and alternately, $v_{p}^{2} = G M (\frac{2}{r_{p}} - \frac{1}{a}) = \frac{G M}{a} (\frac{2 a - r_{p}}{r_{p}}) = \frac{G M}{a} (\frac{r_{a}}{r_{p}}) = \frac{G M}{a} {(\frac{b}{r_{p}})}^{2}$

Therefore, specific angular momentum $h = r_{p} v_{p} = r_{a} v_{a} = b \sqrt{\frac{G M}{a}}$ , and

Total angular momentum $L = m h = m b \sqrt{\frac{G M}{a}}$

More in Astronomy

EarthThird planet from the Sun SunStar at the centre of the Solar System MoonNatural satellite orbiting Earth AstronomyScientific study of celestial objects Solar SystemPlanetary system consisting of the Sun and objects orbiting it Romanization of GreekRendering Greek in Latin characters