Physics 320
Astrophysics I:  Lecture #7
Prof. Dale E. Gary

Orbital Mechanics II

B: Orbits and Energy

1. Conservation of Energy

Recall the concepts of potential energy U, and kinetic energy K, the sum of which gives the total energy
E = K + U
where K = 1/2 mv2.  But what is the potential energy appropriate to our planetary system?

Remember that only differences in potential energy are important.  If we raise a mass at the surface of the Earth by a distance h, we do work against gravity

W F . dr  = -mgDr = -mgh

and the negative of this work is the change of potential energy

DU = -W = - F . dr  = mgh
Raising a mass from the floor to the desk raises the potential energy by mgh, and raising it by the same height from the desk to a shelf also raises U by mgh -- only the difference DU matters. Thus, we are free to choose our zero of potential energy anywhere we wish.

We will choose DU = 0 at r = infinity, which means the potential energy anywhere in the system is negative.  Now what is the potential energy at position r?  We first place a test mass of mass m at infinity, and then move it radially inward to distance r from the center (from mass M) with the force of gravity acting all the way along this path.  We have

r GMm
U = - -
r .dr

which evaluates to

r dr GMm | r
U = GMm
 = -
| inf

or finally,

U = -GMm/r         (Gravitational potential energy)
As advertised, the potential energy is negative for any r, and approaches zero as r approached infinity. The total energy for an elliptical orbit, then, is
E = K + U = 1/2 mv2- GMm/r
    = 1/2 mG(M+m)[2/r- 1/a] - GMm/r

E = - GMm/2a

where we have used the vis viva equation (which is why this is valid only for elliptical orbits), and we have made the approximation that M >> m.

Note that the total energy is negative, and is just a constant.  Thus, energy is conserved along the orbit, as of course it must be.

2. Total Energy for any Orbit
Now we will examine the total energy for any orbit, not limited to elliptical orbits.  To do this, we need to use conservation of angular momentum.  It is possible to show, although we will not do so, that the angular momentum and eccentricity are related by
L2 / GMm2r = 1 + e cos q,
so solving for r, we get
r = (L2 / GMm2) / (1 + e cos q).
Note that this has the same form as the general expression for the polar equation for a conic section.  Let us now repeat the calculation of total energy in the same way as before:
vr = dr/dt = (L2 / GMm2) d/dt [(1 + e cos q)]-1
    = (L2 / GMm2) e sin q dq/dt / (1 + e cos q)2
Recall that |L| = |r x p| = mrvq, so that rvq = L/m. But since vq = r dq/dt, we have r2dq/dt = L / m,  so. Putting this into the above equation,
vr = (GMm / L) (e sin q).
The corresponding perpendicular component of the velocity is even simpler to derive
vq = r dq/dt = r2dq/dt / r = L / rm = (GMm / L) (1 + e cos q).
So the total velocity is
v2 = vr2 + vq2 = (GMm / L)2 (1 + 2e cos q + e2)
and finally the kinetic energy is
= 1/2 mv2 = 1/2 m(GMm / L)2 (1 + 2e cos q + e2)
while the potential energy is
U = -GMm/r = -(GMm)2m/L2(1 + e cos q)
    = -1/2 m(GMm / L)2 (2 + 2e cos q)
We finally come to the rather simple expression for total energy, for an orbit of the form of any conic section:
E = K + U = 1/2 m(GMm / L)2 (e2 - 1)            (Total energy for any orbit)
It is instructive to solve this equation for the eccentricity, to get
e = [1 + 2L2E / (GMm)2m]1/2
In particular, note that
3. Planetary Motion and Effective Potential
We can consider the kinetic energy as having a radial part and an angular part:
= 1/2 mv2 = 1/2 mvr2 + 1/2 mvq2
=                Kr    +     Kq
where vq2 = L2 / m2r2    ==>    Kq = 1/2 mvq2 =1/2 L2 / mr2.  Consider a body moving inward through the solar system.  As r decreases,  the angular kinetic energy increases for smaller r due to conservation of angular momentum.  Since the total energy is conserved, this increasing angular kinetic energy comes in part from the radial kinetic energy.  The system acts as though there were a radial force opposing the inward motion.  Let us write the total energy as
Kr+ Kq + U(r) = 1/2 mvr2 + 1/2 L2 / mr2 - GMm/r = Kr+ Ueff
so we call this combination an effective potential.  Graphically, we have:
There is a nice way to interpret orbits using this diagram.  Since E = constant on any orbit, if E < 0 we have a bound orbit, with two turning points where Kr= 0,
rmin = rperi
rmax = rap
Note that any orbit with E < 0 has two turning points.  For E = 0, there is only one turning point, and the object reaches infinity with zero energy.  Finally, for E > 0, there is again one turning point, but the object reaches infinity with energy left over.

Effective potential in general relativity (strongly curved space-time).