Skip to main content

CLP-4 Vector Calculus

Section 1.9 Optional — Central Forces

One of the great triumphs of Newtonian mechanics was the explanation of Kepler's laws 1 , which said
  1. The planets trace out ellipses about the sun as focus.
  2. The radius vector \(\vr\) sweeps out equal areas in equal times.
  3. The square of the period of each planet is proportional to the cube of the major axis of the planet's orbit.
Newton showed that all of these behaviours follow from the assumption that the acceleration \(\va(t)\) of each planet obeys the law of motion \(m\va =\vF\) where \(m\) is the mass of the planet and
\begin{equation*} \vF = -\frac{GMm}{r^3}\vr \end{equation*}
is the “gravitational force” applied on the planet by the sun. Here \(G\) is a constant 2 , called the “gravitational constant” or the “universal gravitational constant”, \(M\) is the mass of the sun, \(\vr\) is the vector from the sun to the planet and \(r=|\vr|\text{.}\)
In this section, we'll show that some of these properties follow from the weaker assumption that the acceleration \(\va(t)\) of each planet obeys the law of motion \(m\va =\vF\) with \(\vF\) being a central force. That is, the assumption that \(\vF\) is parallel to \(\vr\text{.}\) The verification that the other properties follow from the specific form of the gravitational force, proportional to \(r^{-2}\text{,}\) will be delayed until the optional §1.10.
So, in this section, we assume that we have a parametrized curve \(\vr(t)\) and that this curve obeys
\begin{equation*} m\difftwo{\vr}{t}(t) = \vF\big(\vr(t)\big) \end{equation*}
where, for all \(\vr\in\bbbr^3\text{,}\) \(\vF(\vr)\) is parallel to \(\vr\text{.}\) We shall show that
  1. The path \(\vr(t)\) lies in a plane through the origin and that
  2. the radius vector \(\vr\) sweeps out equal areas in equal times.
We'll start by trying to guess what the plane is. Pretend that we know that \(\vr(t)\) lies in a fixed plane through the origin. Then \(\vv(t)=\diff{\vr}{t}(t)\) lies in the same plane and \(\vr(t)\times\vv(t)\) is perpendicular to the plane. If our path really does lie in a fixed plane, \(\vr(t)\times\vv(t)\) cannot change direction — it must always be parallel to the normal vector to the plane. So let's define
\begin{equation*} \vOm(t) = \vr(t)\times\vv(t) \end{equation*}
and check how it depends on time. By the product rule,
\begin{align*} \diff{\vOm}{t}(t) &=\diff{\ }{t}\big(\vr(t)\times\vv(t)\big) =\vv(t)\times\vv(t) + \vr(t)\times\va(t)\\ &=\frac{1}{m}\vr(t)\times\vF\big(\vr(t)\big)\\ &=\vZero \tag{A} \end{align*}
because \(\vr(t)\) and \(\vF\big(\vr(t)\big)\) are parallel. So \(\vOm(t)\) is 3  in fact independent of \(t\text{.}\) It is a constant vector that we'll just denote \(\vOm\text{.}\)
As \(\vr(t)\times\vv(t)=\vOm\text{,}\) we have that \(\vr(t)\) is always perpendicular to \(\vOm\) and
\begin{equation*} \vr(t)\cdot\vOm =0 \end{equation*}
  • If \(\vOm\ne \vZero\text{,}\) this is exactly the statement that \(\vr(t)\) always lies in the plane through the origin with normal vector \(\vOm\text{.}\)
  • If \(\vOm=\vZero\text{,}\) then \(\vr(t)\) is always parallel to \(\vv(t)\) and there is some function \(\alpha(t)\) such that
    \begin{equation*} \diff{\vr}{t}(t) = \vv(t) = \alpha(t)\,\vr(t) \end{equation*}
    This is a first order, linear, ordinary differential equation that we can solve by using an integrating factor. Set
    \begin{equation*} \beta(t) = \int_0^t\alpha(t)\ \dee{t} \end{equation*}
    Then
    \begin{align*} \diff{\vr}{t}(t) = \alpha(t)\,\vr(t) &\iff e^{-\beta(t)} \diff{\vr}{t}(t) -\alpha(t)e^{-\beta(t)}\,\vr(t)=0\\ &\iff \diff{\ }{t}\big[e^{-\beta(t)}\vr(t)\big] = 0\\ &\iff e^{-\beta(t)}\vr(t) = \vr(0)\\ &\iff \vr(t) = e^{\beta(t)}\vr(0) \end{align*}
    so that \(\vr(t)\) lies on a line through the origin. This makes sense — the particle is always moving parallel to its radius vector.
This completes the verification that \(\vr(t)\) lies in a plane through the origin.
Now we show that the radius vector \(\vr(t)\) sweeps out equal areas in equal times. In other words, we now verify that the rate at which \(\vr(t)\) sweeps out area is independent of time. To do so we rewrite the statement that \(|\vr(t)\times\vv(t)\big|\) is constant in polar coordinates. Writing \(\vr(t) = r(t)\hat\vr\big(\theta(t)\big)\) and then applying Lemma 1.8.2.b gives that
\begin{align*} \text{constant} = \big|\vr\times\vv\big| &= \Big|r\hat\vr\times\Big(\diff{r}{t}\ \hat\vr + r\ \diff{\theta}{t}\ \hat\vth\Big)\Big| =r^2\diff{\theta}{t}\\ &\text{since}\quad |\hat\vr\times\hat\vr|=0,\ |\hat\vr\times\hat\vth|=1 \end{align*}
is constant. It now suffices to observe that \(r(t)^2\diff{\theta}{t}(t)\) is exactly twice the rate at which \(\vr(t)\) sweeps out area. To see this, just look at the figure below. The shaded area is essentially a wedge of a circular disk of radius \(r\text{.}\) (If \(r(t)\) were independent of \(t\text{,}\) it would be exactly a wedge of a circular disk.) Its area is the fraction \(\frac{\dee{\theta}}{2\pi}\) of the area of the full disk, which is
\begin{equation*} \frac{\dee{\theta}}{2\pi}\ \pi r^2 = \frac{1}{2}r^2\,\dee{\theta} \end{equation*}

Exercises Exercises

Exercise Group.

Exercises — Stage 3
1. (✳).
Let \(\vr(t) = x(t)\,\hi + y (t)\,\hj + z(t)\,\hk\) be the position of a particle at time \(t\) . Suppose the motion of the particle satisfies the differential equation \(\difftwo{\vr}{t} = f (r) \vr\) where \(r = |\vr|\) .
  1. Suppose \(f(r)\) is an arbitrary function of \(r\) . Prove or disprove each of the following statements.
    1. The motion of the particle is planar.
    2. The path of the particle sweeps out equal areas in equal times.
  2. Find all forms of \(f(r)\) for which the motion of the particle always lies on a straight line.
  3. Give a specific form of \(f(r)\) for which the motion of the particle could lie on an ellipse.
2. (✳).
An object moves along a curve in the \(xy\)-plane having polar equation \(r=\frac{1}{\theta+\al}\) (where \(\al\) is a constant) under the influence of a central force so that the object has no transverse acceleration.
  1. Verify that \(r^2\dot\theta=h\) remains constant as the object moves.
  2. Express the magnitude of the acceleration of the object as a function of \(r\) and \(h\text{.}\)
The German astronomer Johannes Kepler (1571–1630) developed these laws during the course of an attempt to relate the five extraterrestrial planets then known to the five Platonic solids. He based the laws on a great number of careful measurements made by the Danish Astronomer Tycho Brahe (1546–1601). Then Isaac Newton (English, 1642–1727) provided the explanation in 1687. Kepler also wrote a paper entitled “On the Six-Cornered Snowflake”. Tycho Brahe lost his nose in a sword duel and wore a prosthetic nose from then on. The story is that Brahe died from a burst bladder that resulted from his refusing to leave the dinner table before his host.
Its value is about \(6.67408\times10^{-11} \text{m}^3\,\text{kg}^{-1}\,\text{sec}^{-2}\text{.}\)
Physicists call \(m\,\vOm(t)\) the angular momentum at time \(t\) and refer to (A) as (an example of) conservation of angular momentum. Conservation of angular momentum is exploited in gyro-compasses and by ice skaters (to spin faster/slower).