# 6. General planar motion#

Although Newton’s laws of motion, the various force laws, and the three conservation laws we have derived, are all valid in three dimensions, we have so far restricted our study of motion almost exclusively to two special cases: linear motion in one dimension, and rotational motion in a plane, where the radius of the rotation is constant. Although for the second case we do need two directions to describe it, the motion itself is constricted to a circle, and thus essentially one-dimensional. In this section, we’ll look at general motion in a plane - which turns out to capture a large number of important nontrivial cases.

## 6.1. Projectile motion#

The simplest case of two-dimensional motion occurs when a particle experiences a force only in one direction. The prime example of this case is the motion of a projectile in Earth’s (or any other planet’s) gravitational field as locally described by Galilean gravity (equation (2.8)): \(\bm{F} = m \bm{g}\). Once a projectile has been fired with a certain initial velocity \(\bm{v}_0\), we can find its trajectory by solving the equation of motion that follows from Newton’s second law: \(m \bm{g} = m \ddot{\bm{r}}\). We can decompose \(\bm{r}\) and \(\bm{v}_0\) in horizontal (\(x\)) and vertical (\(z\)) components; each of them has its own one-dimensional equation of motion, which we already solved in Section 2.3. The horizontal component experiences no force and thus executes a simple linear motion with uniform velocity \(v_0 \cos(\theta_0)\), where \(\theta_0 = \arccos(\bm{v}_0 \cdot \bm{\hat{x}}) / v_0\) is the angle with the horizontal under which the projectile was fired and \(v_0 = |\bm{v}_0|\) the initial speed. Likewise, because the acceleration due to gravitation is constant, our projectile will execute a uniformly accelerated motion in the vertical direction with initial velocity \(v_0 \sin(\theta_0)\). If the projectile’s initial position is \((x_0, z_0)\), its motion is thus described by:

We can find the trajectory of the projectile through space by eliminating the time from equation (6.1). We have \(t = (x - x_0) / v_0 \cos(\theta_0)\), which gives for the \(z\) coordinate as a function of \(x\):

Equation (6.2) describes the well-known parabolic trajectory of a projectile under the force of gravity.

## 6.2. General planar motion in polar coordinates#

Although in principle all planar motion can be described in Cartesian coordinates, they are not always the easiest choice. Take, for example, a central force field (a force field whose magnitude only depends on the distance to the origin, and points in the radial direction), as we’ll study in the next section. For such a force field polar coordinates are a more natural choice than Cartesians. However, polar coordinates do carry a few subtleties not present in the Cartesian system, because the direction of the axes depends on position. We will therefore first derive the relevant expressions for the position, velocity and acceleration vector, as well as the components of the force vector, in polar coordinates for the general case.

As we already know (see Section 15.1.2), the position vector \(\vec{r} = x\bm{\hat{x}} + y\bm{\hat{y}}\) has a particularly simple expression in polar coordinates: \(\bm{r} = r \bm{\hat{r}}\), where \(r = \sqrt{x^2 + y^2}\). To find the velocity and acceleration vectors in polar coordinates, we take time derivatives of \(\bm{r}\). Note that because the orientation of the polar basis vectors depends on the position in space, the time derivative acts on both the distance to the origin \(r\) and the basis vector \(\bm{\hat{r}}\). Because the two polar basis vectors are each other’s derivatives with respect to \(\theta\) (see equation (15.4)), we find for their time derivatives:

For the velocity and acceleration vectors we then find:

Note that equations (5.3) and (5.6) are the special cases of equations (6.4) and (6.5) for which both the radius \(r\) and the angular velocity \(\omega = \dot{\theta}\) are constant.

Using equation (6.5) for \(\ddot{\bm{r}}\) in Newton’s second law, we get an expression decomposing the net force \(\bm{F}\) into a radial and an angular part, each of which consists of two terms:

The two terms in \(F_r\) are readily identified as the radial acceleration \(\ddot{r}\) (acting along the line through the origin) and the centripetal force (which causes objects to rotate around the origin, see equation (5.10)). The first term \(r \ddot{\theta}\) in \(F_\theta\) is the tangential acceleration \(\alpha\) of a rotating object whose angular velocity is changing (equation (5.8)). The last term in \(F_\theta\) we have not encountered before; it is known as the *Coriolis force*

and is associated with a velocity in both the radial and the angular direction. It is fairly weak on everyday length scales, but gets large on global length scales. In particular, if you move over the surface of the Earth (necessarily with a nonzero angular component of your velocity), it tends to deflect you from a straight path. On the Northern hemisphere, if you move horizontally, it tends to push you to the right; it also pushes you west when going up, and east when going down. Coriolis forces are responsible for the rotational movement of air around high and low pressure zones, causing respectively clockwise and counterclockwise currents around them on the Northern hemisphere (see Fig. 6.1). We’ll encounter the Coriolis force again in the more general three-dimensional setting in Section 7.2.

## 6.3. Motion under the action of a central force#

A *central force* is a force that points along the (positive or negative) radial direction \(\bm{\hat{r}}\), and whose magnitude depends only on the distance \(r\) to the origin - so \(\bm{F}(\bm{r}) = F(r) \bm{\hat{r}}\). Central forces can be defined in both two and three dimensions, with the three-dimensional concept of the radial distance (to the origin) and direction (direction of increasing \(r\)) completely analogous to the two-dimensional case. Two important examples of central forces are (general) Newtonian gravity (2.9) and the Coulomb force (2.10) between two charged objects. Although these forces are three-dimensional examples, discussing them here is appropriate, as the following theorem shows.

The motion of a particle under the action of a central force takes place in a plane.

Proof. We first note that a central force can exert no torque on an object: \(\bm{\tau} = \bm{r} \times \bm{F} = F(r) (\bm{r} \times \bm{\hat{r}}) = 0\). Consequently, under the action of a central force, angular momentum is conserved. Moreover, we have \(\bm{r} \cdot \bm{L} = \bm{r} \cdot (\bm{r} \times \bm{p}) = 0\) and \(\bm{v} \cdot \bm{L} = \bm{v} \cdot (\bm{r} \times m\bm{v}) = 0\). Both the position vector \(\bm{r}\) and the velocity vector \(\bm{v}\) thus lie in the plane perpendicular to \(\bm{L}\). As \(\bm{L}\) is conserved \(\bm{r}\) and \(\bm{v}\) must be confined to the plane perpendicular to \(\bm{L}\) and through the origin.

Applying the results of the previous section to the motion of a single particle under the action of a central force, we find (for the plane in which the particle moves):

where we used that for a single particle, the magnitude of the angular momentum is given by \(L = m r^2 \dot{\theta}\). Rewriting equation (6.10) gives

where \(F_\mathrm{cf}\) is known as the *centrifugal force*, as it tends to move our particle away from the origin. We can write the centrifugal force as the derivative of a potential:

Writing the original force as the derivative of a potential \(U(r)\) as well, we can write down an equation for the total energy of the system:

For both Newtonian gravity and the Coulomb force, the potential can be written as \(U(r) = - \alpha/r\), where \(\alpha = G m_1 m_2\) for gravity and \(\alpha = -k_\mathrm{e} q_1 q_2\) for Coulomb’s law. We can then rewrite the energy equation as a differential equation for \(r(t)\):

To describe the motion of the particle, rather than specifying \(r(t)\) and \(\theta(t)\), we would like to express \(r\) as a function of \(\theta\). We can rewrite equation (6.14) to a differential equation for \(r(\theta)\) by invoking the chain rule:

where we again used that \(L = m r^2 \dot{\theta}\). Equation (6.14) now becomes:

Equation (6.16) isn’t pretty, but can be simplified a little by absorbing the \(1/r^2\) on the left hand side into the derivative, and completing the square on the right hand side:

We can simplify equation (6.17) further by introducing a new variable, \(z=\frac{1}{r} - m\alpha / L^2\). We also introduce a dimensionless constant \(\varepsilon = \sqrt{1 + 2 E L^2 / m\alpha^2}\) and an inverse length \(q = m\alpha \varepsilon / L^2\). With these substitutions, our equation becomes:

We can solve equation (6.18) by separation of variables:

Taking the reference angle \(\theta_0\) (our integration constant) to be zero, we find \(z(\theta) = q \cos(\theta)\). Translating back to \(r(\theta)\), we obtain a fairly simple solution:

What the solution (6.20) (the *orbit* of our particle under the action of the central force) actually looks like, depends on the value of our dimensionless variable \(\varepsilon\), known as the *eccentricity* of the orbit. To find out which orbits we can get, we translate equation (6.20) back to Cartesian coordinates, using \(x = r \cos\theta\). Defining \(k = L^2/m\alpha\), we get \(k = r + \varepsilon r \cos \theta = r + \varepsilon x\), or \(r = k - \varepsilon x\). Now using \(r^2 = x^2 + y^2\), we get

We can now distinguish four possibilities (see Fig. 6.2a):

\(\varepsilon = 0\): In this case, equation (6.21) becomes \(x^2 + y^2 = k^2\), so our orbit is a

*circle*with the origin at its center.\(0 < \varepsilon < 1\): For this case, with some algebra (see Exercise 6.2), we can rewrite equation (6.21) as \(((x-x_0)/a)^2 + (y/b)^2 = 1\), where \(a = k / (1-\varepsilon^2)\), \(x_0 = - \varepsilon a\), and \(b = k / \sqrt{1-\varepsilon^2}\). These orbits are

*ellipses*, with the center of the ellipse at \((x_0, 0)\), semi-major axis \(a\), semi-minor axis \(b\), and focal length \(f = \sqrt{a^2 - b^2} = k \varepsilon / (1-\varepsilon^2) = - x_0\). One of the foci thus lies at the origin.\(\varepsilon = 1\): Equation (6.21) now becomes \(y^2 = k^2 - 2kx\), which is the equation for a

*parabola*(extending along the negative \(x\)-axis) with its ‘top’ (in this case, rightmost point) at \((k/2, 0)\) and focal length \(k/2\), so the (single) focus is again located at the origin.\(\varepsilon > 1\): This case again requires some algebra to rewrite equation (6.21) in a recognizable standard form: \(((x-x_0)/a)^2 - (y/b)^2 = 1\), where \(a = k / (\varepsilon^2 - 1)\), \(x_0 = \varepsilon a\) and \(b = k / \sqrt{\varepsilon^2 - 1}\). These orbits are

*hyperbola*, crossing the \(x\)-axis at \((x_0, 0)\), and approaching asymptotes \(y = \pm b ((x/a) - \varepsilon)\), which meet at \((x_0+a,0)\). The focal length is now \(f = \sqrt{a^2 + b^2} = k \varepsilon / (\varepsilon^2 - 1) = \varepsilon a = x_0 + a\), so the focus of the hyperbola is also located at the origin.

In mathematics, these four possible types of orbits are known as *conic sections*: the curves you can get by intersecting a cone with a plane. Specifically, when the central force is gravity, such as in the solar system (where the sun is so much heavier than everything else combined that to good approximation we can describe orbits as being determined by the sun’s gravitational field alone), the four cases listed above are the only possible orbits bodies can have. The planets, dwarf planets, asteroids, and many minor objects in our solar system all follow elliptical orbits around the sun, albeit with different eccentricity^{[2]} - Earth’s is almost zero (0.017), but that of Mars is significantly higher (0.09), and of Pluto much higher still (0.25). Comets, on the other hand, typically parabolic or hyperbolic orbits. Spacecraft such as the Voyager and New Horizons probes are often put on trajectories past planets that are not their final destination, to pick up (or lose) speed through a gravitational assist (in which they take a little bit of momentum from the planet’s orbit); those paths past planets are typically hyperbola. Getting a spacecraft to orbit another planet (i.e., in a bound, so elliptical) orbit is actually much harder, but again, the resulting orbit is described by the maths presented above.

From the dimensionless eccentricity of an orbit, we can work out the total energy of that orbit easily. As \(\varepsilon = \sqrt{1 + 2 E L^2 / m \alpha^2}\), we find

For elliptical orbits (where \((1-\varepsilon^2) = k/a = L^2 / (\alpha m a)\), with \(a\) the semi-major axis), this expression gives a negative total energy that depends on the length of the semi-major axis only

which is exactly half the potential energy of such an orbit. For a parabolic orbit (the boundary case between a closed and an open orbit), \(\varepsilon = 1\), and so the total energy of the orbit vanishes; for a hyperbolic orbit we have \((\varepsilon^2 - 1) = k/a\) and we get a positive energy total energy.

## 6.4. Kepler’s laws#

The fact that the planets move in elliptical orbits was first discovered by Kepler, based on observational data alone (he didn’t have the benefit, as we do, of living after Newton and thus knowing about Newton’s law of gravity). Kepler summarized his observational facts in three laws, which we can, with the benefit of hindsight, prove to be corollaries of Newton’s laws.

Johannes Kepler (1571-1630)

**Johannes Kepler** (1571-1630) was a German astronomer and mathematician who made major contributions to understanding the motion of the planets. Copernicus had published his heliocentric (rather than geocentric) view of the universe posthumously in 1543, but the two systems were still heavily debated in Kepler’s time. Having been convinced that Copernicus was right, Kepler set out to construct a geometric description of the solar system. He initially tried to do so using polyhedra and Platonic solids, but found that these could not accurately describe the data. In 1600, Kepler met with astronomer Tycho Brahe, who had made meticulous observations of the known planets, and, having been convinced of Kepler’s skills in mathematics, shared his data with him. After Tycho’s death in 1601, Kepler succeeded him as imperial mathematician in Prague, where he developed his laws over the next decade. Unfortunately, Kepler’s Calvinist views got him in trouble frequently with both the Catholic and the Lutheran church, which led to his excummunication, but he managed to avoid further persecution by moving frequently, and he always could continue his scientific work. The Kepler spacecraft and mission, launched in 2009 to hunt for extrasolar terrestial planets, is named in his honor.

(Kepler’s first law)

The planets move in elliptical orbits, with the sun at one of the foci.

(Kepler’s second law)

A line segment joining a planet and the sun sweeps out equal areas during equal intervals of time.

Proof. This law is nothing but a special case of conservation of angular momentum. Consider a small piece of the orbit, in which the planet moves a distance \(\mathrm{d}x\). The lines connecting the initial and final points of this piece of orbit with the sun make an angle \(\mathrm{d}\theta\). If the initial distance from the planet to the sun was \(r\), and the final distance \(r + \mathrm{d}r\), we have, to first order, \(\mathrm{d}x = r \mathrm{d}\theta\). The infinitesimal area the planet has swiped out is then given by (area of a triangle): \(\mathrm{d}A = \frac12 r \mathrm{d}x = \frac12 r^2 \mathrm{d}\theta\). If we want to know how much area was swept out over an amount of time, we need to know the time derivative of \(A\), which is thus given by \(\mathrm{d}A / \mathrm{d}t = \frac12 r^2 \mathrm{d}\theta / \mathrm{d}t\). Now using that the angular momentum of the planet is given by \(L = m r^2 \dot{\theta}\), we find

which is constant if \(L\) is conserved.

(Kepler’s third law)

The square of the period \(T\) of an orbit is proportional to the cube of its semi-major-axis length \(a\):

where \(M_\odot\) is the mass of the sun.

Proof. We integrate equation (6.24) over the period of a whole orbit, which gives \(A = LT / 2m\). By Kepler’s first law, the orbit is an ellipse, so its area equals \(A = \pi a b\), with \(a\) and \(b\) the ellipse’s semi-major and semi-minor axes. The two axes are related by \(b = a \sqrt{1-\varepsilon^2}\), with \(\varepsilon\) again the eccentricity of the ellipse. Making these substitutions and squaring the resulting relation, we get:

Using \(k = L^2 / m\alpha\), like in equation (6.21), and the observation that for an elliptical orbit \(k/(1-\varepsilon^2) = a\), we get \(L^2/m(1-\varepsilon^2) = \alpha a\). Now for orbits in the solar system, \(\alpha = G M_\odot m\), so we arrive at equation (6.25).

Kepler’s third law allows us to calculate the mass of any celestial body that is orbited by a (natural or artificial) satellite, simply by measuring the period and semi-major axis of the satellite’s orbit.

## 6.5. Problems#

(Geosynchronous orbit)

A particularly useful orbit for satellites is the *geosynchronous* one: the orbit in which the satellite rotates around the Earth in exactly one day, so with respect to the ground, it is always in the same position. Find the altitude (i.e., distance above the Earth’s surface) for a circular geosynchronous orbit.

(Properties of elliptical orbits)

In Section 6.3, we found four possibilities for the orbits under a central force, all described by equation (6.20): \(r(\theta) = k / (1+\varepsilon \cos\theta)\). For \(0 < \varepsilon < 1\) these orbits are ellipses. In this problem, we’ll find expressions for the quantities that we usually use to parametrize an ellipse (see Fig. 6.2b) in terms of \(\varepsilon\) and \(k\).

The easiest parameter to find is the length \(a\) of the semi-major axis, as twice the semi-major axis equals the distance between \(r(0)\) and \(r(\pi)\). From this equality, show that \(a = k / (1-\varepsilon^2)\).

The focal distance \(f\) is the distance from the center of the ellipse to the focal point. As we put the focal point in the origin, we have \(f = a - r(0)\). Show that \(f = \varepsilon a= k\varepsilon / (1-\varepsilon^2)\).

Find \(x_0\) (no extra calculation required).

At the ‘top’ of the ellipse (where \(y=b\)), we have \(\cos(\theta - \pi) = f/r\) and \(r^2 = b^2 + f^2\). From these two equations, show that \(b = a \sqrt{1-\varepsilon^2} = k / \sqrt{1-\varepsilon^2}\).

*Hint: first show that \(r=a\) at this position; from that you also get \(f = \sqrt{a^2-b^2}\).*

(Velocity vector in an elliptical orbit)

Prove that for an elliptical orbit, the velocity of the orbiting object is purely tangential at both the pericenter and apocenter (the closest and furthest point from the origin in the orbit). *Hint: take the time derivative of equation (6.20) to find the radial speed.*

(Juno spacecraft)

At the time of writing, the Juno spacecraft is in a polar orbit around Jupiter. It’s perijove (closest distance to Jupiter) is at an altitude of \(4200\;\mathrm{km}\), or \(75600\;\mathrm{km}\) from Jupiter’s center; Juno’s apojove is \(8.1\cdot 10^6 \;\mathrm{km}\) from Jupiter’s center. The spacecraft weighs \(1593\;\mathrm{kg}\). Jupiter’s mass (known from the orbits of its moons) is \(1.90 \cdot 10^{27}\;\mathrm{kg}\).

Find the semi-major axis and eccentricity of Juno’s orbit around Jupiter.

How long does it take Juno to complete one orbit?

Find Juno’s maximum and minimum speed while traversing this orbit. You may find the result of Exercise 6.3 useful in answering this question.

Suppose NASA wants to make Juno’s orbit more circular. Should it have Juno fire its thruster rockets to accelerate or decelerate its angular velocity?

(Kepler’s laws and reduced mass)

Kepler’s laws apply to the case that an object with relatively small mass \(m\) orbits an object with large mass \(M\), which we assume stays fixed. Technically, this is incorrect: both objects rotate about their common center of mass. Fortunately, we can still use the expressions derived in this section, with a small modification. To see how this works, we write down the equations of motion for the two objects, due to the force they exert on each other:

where \(\bm{x}_1\) is the position of the object with mass \(m\), \(\bm{x}_2\) that of the object with mass \(M\), and \(\bm{r} = \bm{x}_2 - \bm{x}_1\) their separation. We denote the position of the center of mass of the system by \(\bm{R}\), so

As there is no external force acting on the system, the total momentum is conserved and therefore the center of mass cannot accelerate. Prove this statement by showing that

(6.28)#\[(m+M)\ddot{\bm{R}} = 0.\]From equation (6.26), derive an equation of motion for the separation \(\bm{r}\) between the two objects. The equations you found in (a) and (b) together are equivalent to the equations of motion in equation (6.26), but they are uncoupled: we don’t need to know the position of the center of mass to find the separation, and vice versa.

Show that you can re-write the equation of motion for the separation between the two objects as \(\bm{F}(\bm{r}) = \mu \ddot{\bm{r}}\), where \(\mu\) is the

*reduced mass*(which we also encountered when studying collisions in the center of mass frame), given by\[\mu = \frac{m M}{m+M}.\]Note that solving the final equation for the separation \(\bm{r}\) is entirely equivalent to solving the equation of motion of a single particle under the action of a central force, with the modification that the mass of the particle is replaced by the reduced mass. The central force \(\bm{F}(r)\) is still just the force that the two masses exert on each other. For the case that \(m \ll M\), the reduced mass is approximately equal to \(m\).

Calculate the reduced mass of the Earth-Moon two body problem. Can we state that the Moon revolves around the Earth?

Nowhere in the derivations in this problem did we assume that \(m \ll M\). The same rules apply to any two objects. Consider the opposite limit: two objects (these might for instance be binary stars) of equal mass \(M\) that rotate around their common center of mass. Show that for this case, for circular orbits the orbital period is given by

(6.29)#\[T^2 = \frac{2\pi^2 d^3}{GM},\]where \(d\) is the distance between the two objects.

(Walking on a merry-go-round)

A student with mass \(65.0\;\mathrm{kg}\) stands at the center of a simple merry-go-round that consists of a large disk of radius \(1.5\;\mathrm{m}\) and mass \(25\;\mathrm{kg}\) and is making a full rotation every \(2.0\;\mathrm{s}\). The student walks out to a distance of \(0.50\;\mathrm{m}\) from the center.

Find the rotational frequency of the merry-go-round with the student at this point.

What are the forces acting on the student at this point?