Skip to content

Problem Irodov 1.208

A planet of mass mm moves along an elliptical orbit about the Sun, of mass mSm_S, under gravity alone. Using the conservation laws, show that the planet’s total mechanical energy EE depends only on the orbit’s semi-major axis aa, and not on its eccentricity. Find EE as a function of aa.

The planet's elliptical orbit about the Sun at a focus. At perihelion (distance r_p) the velocity v_p is purely tangential, and the same is true of v_a at aphelion (distance r_a); these are the two points the derivation compares.

The planet’s elliptical orbit about the Sun at a focus. At perihelion (distance rpr_p) the velocity vp\vect v_p is purely tangential, and the same is true of va\vect v_a at aphelion (distance rar_a); these are the two points the derivation compares.

Given
  • mm
  • mSm_S
  • aa
Find
E(a)E(a)

Solution

Step 1The turning points of the orbit

At perihelion and aphelion the radial velocity vanishes, so the velocity there is purely tangential; both points lie on the major axis, on opposite sides of the occupied focus, at distances rpr_p and rar_a from it. By the definition of the semi-major axis aa, the major axis has length 2a2a, so

rp+ra=2ar_p+r_a = 2a
(1)
Step 2Angular momentum ties the two speeds together

Angular momentum about the Sun is conserved along the orbit, and at the two apsides the velocity is purely tangential, so evaluating it there needs no projection; equating the two apsidal values relates the two speeds.

mvprp=mvara(2)m v_p r_p = m v_a r_a \tag{2}

so that

va=vp rprav_a = v_p\,\frac{r_p}{r_a}
(3)
Step 3Energy conservation between the two points

Only gravity does work on the planet, so its mechanical energy is the same at perihelion and at aphelion.

12mvp2−γmmSrp=12mva2−γmmSra(4)\frac12 m v_p^2 - \frac{\gamma m m_S}{r_p} = \frac12 m v_a^2 - \frac{\gamma m m_S}{r_a} \tag{4}
Step 4Eliminate the speeds

Substituting (3) for vav_a in (4) and dividing through by mm leaves an equation in vpv_p, rpr_p and rar_a alone.

12vp2−γmSrp=12vp2rp2ra2−γmSra12vp2(1−rp2ra2)=γmS(1rp−1ra)12vp2 (ra−rp)(ra+rp)ra2=γmS ra−rprpra\begin{align*} \frac12 v_p^2 - \frac{\gamma m_S}{r_p} &= \frac12 v_p^2\frac{r_p^2}{r_a^2} - \frac{\gamma m_S}{r_a} \\ \frac12 v_p^2\left(1-\frac{r_p^2}{r_a^2}\right) &= \gamma m_S\left(\frac{1}{r_p}-\frac{1}{r_a}\right) \\ \frac12 v_p^2\,\frac{(r_a-r_p)(r_a+r_p)}{r_a^2} &= \gamma m_S\,\frac{r_a-r_p}{r_p r_a} \tag{5}\end{align*}

Since ra≠rpr_a\neq r_p on a genuine ellipse, the common factor (ra−rp)(r_a-r_p) on both sides of (5) cancels; using (1) for ra+rpr_a+r_p then isolates vp2v_p^2.

vp2=2γmSrp⋅rara+rp=γmS raa rp(6)v_p^2 = \frac{2\gamma m_S}{r_p}\cdot\frac{r_a}{r_a+r_p} = \frac{\gamma m_S\,r_a}{a\,r_p} \tag{6}
Step 5The energy collapses to a function of aa alone

Substituting (6) into the perihelion energy per unit mass, then rearranging (1) to trade ra−2ar_a-2a for rpr_p, eliminates rpr_p as well.

Em=12vp2−γmSrp=γmS ra2a rp−γmSrp=γmSrp(ra2a−1)=γmSrp⋅ra−2a2ara−2a=ra−(rp+ra)=−rpEm=γmSrp⋅−rp2a\begin{align*}\frac{E}{m} &= \frac12 v_p^2 - \frac{\gamma m_S}{r_p} \\ &= \frac{\gamma m_S\,r_a}{2a\,r_p} - \frac{\gamma m_S}{r_p} \\ &= \frac{\gamma m_S}{r_p}\left(\frac{r_a}{2a}-1\right) \\ &= \frac{\gamma m_S}{r_p}\cdot\frac{r_a-2a}{2a} \\ r_a - 2a &= r_a - (r_p + r_a) = -r_p \\ \frac{E}{m} &= \frac{\gamma m_S}{r_p}\cdot\frac{-r_p}{2a}\end{align*}
E=−γm mS2aE = -\dfrac{\gamma m\, m_S}{2a}
(7)
Check

The distance rpr_p canceled completely: the same algebra started from the aphelion side would have left rar_a instead, and both must and do give the identical result, which is exactly the demonstration asked for. As an independent check, a circular orbit (e→0e\to0, rp→ra→ar_p\to r_a\to a) satisfies Newton’s second law directly, mv2/a=γmmS/a2mv^2/a=\gamma m m_S/a^2, so v2→γmS/av^2\to\gamma m_S/a without appeal to the elliptical algebra above, and E/m=v2/2−γmS/a→−γmS/(2a)E/m=v^2/2-\gamma m_S/a\to-\gamma m_S/(2a), matching (7).

Remark

Two orbits sharing the same aa but very different eccentricities, one nearly circular and one a thin sliver of an ellipse, carry exactly the same total energy; only their angular momentum differs, since the quantity conserved in (2) depends on rpr_p and rar_a individually, not merely on their sum. So aa alone fixes the energy, and by Kepler’s third law the period too, while the eccentricity is free to redistribute the same energy between kinetic and potential along the orbit however it likes.