13.5 Kepler’s Laws of Planetary Motion

Kepler’s First Law

The prevailing view during the time of Kepler was that all planetary orbits were circular. The data for Mars presented the greatest challenge to this view and that eventually encouraged Kepler to give up the popular idea. Kepler’s first law states that every planet moves along an ellipse, with the Sun located at a focus of the ellipse. An ellipse is defined as the set of all points such that the sum of the distance from each point to two foci is a constant. Figure shows an ellipse and describes a simple way to create it.

For elliptical orbits, the point of closest approach of a planet to the Sun is called the perihelion. It is labeled point A in Figure. The farthest point is the aphelion and is labeled point B in the figure. For the Moon’s orbit about Earth, those points are called the perigee and apogee, respectively.

An ellipse has several mathematical forms, but all are a specific case of the more general equation for conic sections. There are four different conic sections, all given by the equation

The variables r and [latex]\theta[/latex] are shown in Figure in the case of an ellipse. The constants [latex]\alpha[/latex] and e are determined by the total energy and angular momentum of the satellite at a given point. The constant e is called the eccentricity. The values of [latex]\alpha[/latex] and e determine which of the four conic sections represents the path of the satellite.

One of the real triumphs of Newton’s law of universal gravitation, with the force proportional to the inverse of the distance squared, is that when it is combined with his second law, the solution for the path of any satellite is a conic section. Every path taken by m is one of the four conic sections: a circle or an ellipse for bound or closed orbits, or a parabola or hyperbola for unbounded or open orbits. These conic sections are shown in Figure.

If the total energy is negative, then [latex]0\le e \lt 1[/latex], and Figure represents a bound or closed orbit of either an ellipse or a circle, where [latex]e=0[/latex]. [You can see from Figure that for [latex]e=0[/latex], [latex]r=\alpha[/latex], and hence the radius is constant.] For ellipses, the eccentricity is related to how oblong the ellipse appears. A circle has zero eccentricity, whereas a very long, drawn-out ellipse has an eccentricity near one.

If the total energy is exactly zero, then [latex]e=1[/latex] and the path is a parabola. Recall that a satellite with zero total energy has exactly the escape velocity. (The parabola is formed only by slicing the cone parallel to the tangent line along the surface.) Finally, if the total energy is positive, then [latex]e \gt 1[/latex] and the path is a hyperbola. These last two paths represent unbounded orbits, where m passes by M once and only once. This situation has been observed for several comets that approach the Sun and then travel away, never to return.

We have confined ourselves to the case in which the smaller mass (planet) orbits a much larger, and hence stationary, mass (Sun), but Figure also applies to any two gravitationally interacting masses. Each mass traces out the exact same-shaped conic section as the other. That shape is determined by the total energy and angular momentum of the system, with the center of mass of the system located at the focus. The ratio of the dimensions of the two paths is the inverse of the ratio of their masses.

Orbital Transfers

People have imagined traveling to the other planets of our solar system since they were discovered. But how can we best do this? The most efficient method was discovered in 1925 by Walter Hohmann, inspired by a popular science fiction novel of that time. The method is now called a Hohmann transfer. For the case of traveling between two circular orbits, the transfer is along a “transfer” ellipse that perfectly intercepts those orbits at the aphelion and perihelion of the ellipse. Figure shows the case for a trip from Earth’s orbit to that of Mars. As before, the Sun is at the focus of the ellipse.

For any ellipse, the semi-major axis is defined as one-half the sum of the perihelion and the aphelion. In Figure, the semi-major axis is the distance from the origin to either side of the ellipse along the x-axis, or just one-half the longest axis (called the major axis). Hence, to travel from one circular orbit of radius [latex]{r}_{1}[/latex] to another circular orbit of radius [latex]{r}_{2}[/latex], the aphelion of the transfer ellipse will be equal to the value of the larger orbit, while the perihelion will be the smaller orbit. The semi-major axis, denoted a, is therefore given by [latex]a=\frac{1}{2}({r}_{1}+{r}_{2})[/latex].

Let’s take the case of traveling from Earth to Mars. For the moment, we ignore the planets and assume we are alone in Earth’s orbit and wish to move to Mars’ orbit. From Figure, the expression for total energy, we can see that the total energy for a spacecraft in the larger orbit (Mars) is greater (less negative) than that for the smaller orbit (Earth). To move onto the transfer ellipse from Earth’s orbit, we will need to increase our kinetic energy, that is, we need a velocity boost. The most efficient method is a very quick acceleration along the circular orbital path, which is also along the path of the ellipse at that point. (In fact, the acceleration should be instantaneous, such that the circular and elliptical orbits are congruent during the acceleration. In practice, the finite acceleration is short enough that the difference is not a significant consideration.) Once you have arrived at Mars orbit, you will need another velocity boost to move into that orbit, or you will stay on the elliptical orbit and simply fall back to perihelion where you started. For the return trip, you simply reverse the process with a retro-boost at each transfer point.

To make the move onto the transfer ellipse and then off again, we need to know each circular orbit velocity and the transfer orbit velocities at perihelion and aphelion. The velocity boost required is simply the difference between the circular orbit velocity and the elliptical orbit velocity at each point. We can find the circular orbital velocities from Figure. To determine the velocities for the ellipse, we state without proof (as it is beyond the scope of this course) that total energy for an elliptical orbit is

where [latex]{M}_{\text{S}}[/latex] is the mass of the Sun and a is the semi-major axis. Remarkably, this is the same as Figure for circular orbits, but with the value of the semi-major axis replacing the orbital radius. Since we know the potential energy from Figure, we can find the kinetic energy and hence the velocity needed for each point on the ellipse. We leave it as a challenge problem to find those transfer velocities for an Earth-to-Mars trip.

We end this discussion by pointing out a few important details. First, we have not accounted for the gravitational potential energy due to Earth and Mars, or the mechanics of landing on Mars. In practice, that must be part of the calculations. Second, timing is everything. You do not want to arrive at the orbit of Mars to find out it isn’t there. We must leave Earth at precisely the correct time such that Mars will be at the aphelion of our transfer ellipse just as we arrive. That opportunity comes about every 2 years. And returning requires correct timing as well. The total trip would take just under 3 years! There are other options that provide for a faster transit, including a gravity assist flyby of Venus. But these other options come with an additional cost in energy and danger to the astronauts.

Kepler’s Second Law

Kepler’s second law states that a planet sweeps out equal areas in equal times, that is, the area divided by time, called the areal velocity, is constant. Consider Figure. The time it takes a planet to move from position A to B, sweeping out area [latex]{A}_{1}[/latex], is exactly the time taken to move from position C to D, sweeping area [latex]{A}_{2}[/latex], and to move from E to F, sweeping out area [latex]{A}_{3}[/latex]. These areas are the same: [latex]{A}_{1}={A}_{2}={A}_{3}[/latex].

Comparing the areas in the figure and the distance traveled along the ellipse in each case, we can see that in order for the areas to be equal, the planet must speed up as it gets closer to the Sun and slow down as it moves away. This behavior is completely consistent with our conservation equation, Figure. But we will show that Kepler’s second law is actually a consequence of the conservation of angular momentum, which holds for any system with only radial forces.

Recall the definition of angular momentum from Angular Momentum, [latex]\mathbf{\overset{\to }{L}}=\mathbf{\overset{\to }{r}}\times \mathbf{\overset{\to }{p}}[/latex]. For the case of orbiting motion, [latex]\mathbf{\overset{\to }{L}}[/latex] is the angular momentum of the planet about the Sun, [latex]\mathbf{\overset{\to }{r}}[/latex] is the position vector of the planet measured from the Sun, and [latex]\mathbf{\overset{\to }{p}}=m\mathbf{\overset{\to }{v}}[/latex] is the instantaneous linear momentum at any point in the orbit. Since the planet moves along the ellipse, [latex]\mathbf{\overset{\to }{p}}[/latex] is always tangent to the ellipse.

We can resolve the linear momentum into two components: a radial component [latex]{\mathbf{\overset{\to }{p}}}_{\text{rad}}[/latex] along the line to the Sun, and a component [latex]{\mathbf{\overset{\to }{p}}}_{\text{perp}}[/latex] perpendicular to [latex]\mathbf{\overset{\to }{r}}[/latex]. The cross product for angular momentum can then be written as

[latex]\mathbf{\overset{\to }{L}}=\mathbf{\overset{\to }{r}}\times \mathbf{\overset{\to }{p}}=\mathbf{\overset{\to }{r}}\times ({\mathbf{\overset{\to }{p}}}_{\text{rad}}+{\mathbf{\overset{\to }{p}}}_{\text{perp}})=\mathbf{\overset{\to }{r}}\times {\mathbf{\overset{\to }{p}}}_{\text{rad}}+\mathbf{\overset{\to }{r}}\times {\mathbf{\overset{\to }{p}}}_{\text{perp}}[/latex].

The first term on the right is zero because [latex]\mathbf{\overset{\to }{r}}[/latex] is parallel to [latex]{\mathbf{\overset{\to }{p}}}_{\text{rad}}[/latex], and in the second term [latex]\mathbf{\overset{\to }{r}}[/latex] is perpendicular to [latex]{\mathbf{\overset{\to }{p}}}_{\text{perp}}[/latex], so the magnitude of the cross product reduces to [latex]L=r{p}_{\text{perp}}=rm{v}_{\text{perp}}[/latex]. Note that the angular momentum does not depend upon [latex]{p}_{\text{rad}}[/latex]. Since the gravitational force is only in the radial direction, it can change only [latex]{p}_{\text{rad}}[/latex] and not [latex]{p}_{\text{perp}}[/latex]; hence, the angular momentum must remain constant.

Now consider Figure. A small triangular area [latex]\Delta A[/latex] is swept out in time [latex]\Delta t[/latex]. The velocity is along the path and it makes an angle [latex]\theta[/latex] with the radial direction. Hence, the perpendicular velocity is given by [latex]{v}_{\text{perp}}=v\text{sin}\theta[/latex]. The planet moves a distance [latex]\Delta s=v\Delta t\text{sin}\theta[/latex] projected along the direction perpendicular to r. Since the area of a triangle is one-half the base (r) times the height [latex](\Delta s)[/latex], for a small displacement, the area is given by [latex]\Delta A=\frac{1}{2}r\Delta s[/latex]. Substituting for [latex]\Delta s[/latex], multiplying by m in the numerator and denominator, and rearranging, we obtain

The areal velocity is simply the rate of change of area with time, so we have

Since the angular momentum is constant, the areal velocity must also be constant. This is exactly Kepler’s second law. As with Kepler’s first law, Newton showed it was a natural consequence of his law of gravitation.

Kepler’s Third Law

Kepler’s third law states that the square of the period is proportional to the cube of the semi-major axis of the orbit. In Satellite Orbits and Energy, we derived Kepler’s third law for the special case of a circular orbit. Figure gives us the period of a circular orbit of radius r about Earth:

For an ellipse, recall that the semi-major axis is one-half the sum of the perihelion and the aphelion. For a circular orbit, the semi-major axis (a) is the same as the radius for the orbit. In fact, Figure gives us Kepler’s third law if we simply replace r with a and square both sides.

We have changed the mass of Earth to the more general M, since this equation applies to satellites orbiting any large mass.