Skip to content

Problem Irodov 1.209

A planet AA travels along an elliptical orbit about the Sun. At some instant its distance from the Sun is r0r_0, its speed is v0v_0, and the angle between the radius vector (Sun to planet) and the velocity vector at that instant is α\alpha. Find the maximum and minimum distances the planet ever attains from the Sun over the course of its orbit.

The planet A at the given instant, a distance r_0 from the Sun m_S, its velocity v_0 tangent to the orbit at angle α to the radius vector; r_min and r_max mark perihelion and aphelion.

The planet AA at the given instant, a distance r0r_0 from the Sun mSm_S, its velocity v0\vect v_0 tangent to the orbit at angle α\alpha to the radius vector; rminr_{\mathrm{min}} and rmaxr_{\mathrm{max}} mark perihelion and aphelion.

Given
  • r0r_0
  • v0v_0
  • α\alpha
Find
rmax⁡r_{\max} and rmin⁡r_{\min}

Solution

Step 1Angular momentum from the given point

Angular momentum about the Sun is conserved, so its value at the given point equals its value at either turning point, a distance rr away with speed vv there.

mr0v0sin⁡α=mrvm r_0 v_0\sin\alpha = m r v
v=r0v0sin⁡αrv = \frac{r_0 v_0\sin\alpha}{r}
(1)
Step 2Energy conservation between the two points

Mechanical energy is conserved along the whole orbit, so it too is the same at the given point and at the turning point; dividing through by the planet’s mass mm turns this into a statement about the specific energy alone.

12v02−γmSr0=12v2−γmSr(2)\frac12 v_0^2 - \frac{\gamma m_S}{r_0} = \frac12 v^2 - \frac{\gamma m_S}{r} \tag{2}
Step 3A quadratic for the turning radius

Two shorthands carry the rest of the algebra: ρ\rho for the turning radius measured in units of the given distance, and η\eta for the ratio of the given kinetic energy scale to the gravitational one at that same distance.

ρ=rr0,η=r0v02γmS\rho = \dfrac{r}{r_0},\qquad \eta = \dfrac{r_0v_0^2}{\gamma m_S}
(3)

Substituting (1) for vv into the energy balance (2) and rewriting it in these two shorthands turns it into a quadratic in ρ\rho.

12v02−γmSr0=12v02sin⁡2αρ2−γmSρr01−2η=sin⁡2αρ2−2ηρρ2(1−2η)=sin⁡2α−2ρη\begin{align*}\frac12 v_0^2 - \frac{\gamma m_S}{r_0} &= \frac12\frac{v_0^2\sin^2\alpha}{\rho^2} - \frac{\gamma m_S}{\rho r_0} \\ 1 - \frac{2}{\eta} &= \frac{\sin^2\alpha}{\rho^2} - \frac{2}{\eta\rho} \\ \rho^2\left(1-\frac{2}{\eta}\right) &= \sin^2\alpha - \frac{2\rho}{\eta}\end{align*}

Multiplying through by η\eta and collecting every term on one side,

(2−η)ρ2−2ρ+ηsin⁡2α=0(2-\eta)\rho^2 - 2\rho + \eta\sin^2\alpha = 0
(4)
Step 4Solve the quadratic

The quadratic formula applied to (4) gives both roots, one for each turning point, at once.

ρ=2±4−4(2−η)ηsin⁡2α2(2−η)=1±1−(2−η)ηsin⁡2α2−η\rho = \frac{2\pm\sqrt{4-4(2-\eta)\eta\sin^2\alpha}}{2(2-\eta)} = \frac{1\pm\sqrt{1-(2-\eta)\eta\sin^2\alpha}}{2-\eta}

Multiplying back by r0r_0, with the plus sign giving the farther turning point (rmax⁡r_{\max}) and the minus sign the nearer one (rmin⁡r_{\min}); the book itself reports both distances as the two signs of this one expression, so a single boxed formula carries both asked quantities.

rmax⁡/min⁡=r02−η[1±1−(2−η)ηsin⁡2α],η=r0v02γmS\begin{gathered}r_{\max/\min} = \frac{r_0}{2-\eta}\left[1\pm\sqrt{1-(2-\eta)\eta\sin^2\alpha}\right],\\ \eta=\frac{r_0v_0^2}{\gamma m_S}\end{gathered}
(5)
Check

If the given point is itself a turning point, α=90∘\alpha=90^\circ and sin⁡2α=1\sin^2\alpha=1; the radicand becomes 1−(2−η)η=(1−η)21-(2-\eta)\eta=(1-\eta)^2, so the square root is ∣1−η∣|1-\eta|, and for η<1\eta<1 the plus sign in the boxed formula gives r=r0(2−η)/(2−η)=r0r=r_0(2-\eta)/(2-\eta)=r_0, exactly recovering the starting point as one of its own turning points, as it must.

Remark

The denominator 2−η2-\eta vanishes at η=2\eta=2, i.e. at v02=2γmS/r0v_0^2=2\gamma m_S/r_0, which is precisely the escape speed at r0r_0: rmax⁡r_{\max} diverges there because the trajectory is no longer a closed ellipse but a parabola, consistent with η<2\eta<2 being the actual condition for a bound orbit to exist at all.