Module 28: The Kepler Problem: Planetary Mechanics

28.1 Introduction Keplers Laws:1

1. Each planet moves in an ellipse with the sun at one focus. 2. The radius vector from the sun to a planet sweeps out equal areas in equal time. 3. The period of revolution T of a planet about the sun is related to the major axis A of the ellipse by
T 2 = k A3

(28.1.1)

where k is the same for all planets.

28.2 Planetary Orbits: The Kepler Problem

Since Johannes Kepler first formulated the laws that describe planetary motion, scientists endeavored to solve for the equation of motion of the planets. In his honor, this problem has been named The Kepler Problem. When there are more than two bodies, the problem becomes impossible to solve exactly. The most important three-body problem at the time involved finding the motion of the moon, since the moon interacts gravitationally with both the sun and the earth. Newton realized that if the exact position of the moon were known, the longitude of any observer on the earth could be determined by measuring the moons position with respect to the stars. In the eighteenth century, Leonhard Euler and other mathematicians spent many years trying to solve the three-body problem, and they raised a deeper question. Do the small contributions from the gravitational interactions of all the planets make the planetary system unstable over long periods of time? At the end of 18th century, Pierre Simon Laplace and others found a series solution to this stability question, but it was unknown whether or not the series solution converged after a long period of time. Henri Poincar proved that the series actually diverged. Poincar went on to invent new mathematical methods that produced the modern fields of differential geometry and topology in order to answer the stability question using geometric arguments, rather than analytic methods. Poincar and others did manage to show that the three-body problem was indeed stable, due to the existence of periodic Solution. Just as in the time of Newton and Leibniz and the invention of calculus,

As stated in An Introduction to Mechanics, Daniel Kleppner and Robert Kolenkow, McGraw-Hill, 1973, p 401.

unsolved problems in celestial mechanics became the experimental laboratory for the discovery of new mathematics.

28.3 Reducing the Two-Body Problem into a One-Body Problem

We shall begin our solution of the two-body problem by showing how the motion of two bodies interacting via a gravitational force (two-body problem) is mathematically equivalent to the motion of a single body with a reduced mass given by

m1 m2 m1 + m2

(28.3.1)

that is acted on by an external central gravitational force. Once we solve for the motion of the reduced body in this equivalent one-body problem, we can then return to the real twobody problem and solve for the actual motion of the two original bodies. The reduced mass was introduced in Section 10.7 of these notes. That section used similar but different notation from that used in this chapter. Consider the gravitational force between two bodies with masses m1 and m2 as shown in Figure 28.1.

Figure 28.1 Gravitational force between two bodies.

! Choose a coordinate system with a choice of origin such that body 1 has position r1 and ! ! body 2 has position r2 (Figure 28.2). The relative position vector r pointing from body 2 ! ! ! ! ! to body 1 is r = r1 ! r2 . We denote the magnitude of r by r = r , where r is the distance

is the unit vector pointing from body 2 to body 1, so that between the bodies, and r ! r=r r
(28.3.2)

Figure 28.2 Coordinate system for the two-body problem. The force on body 1 (due to the interaction of the two bodies) can be described as
! mm = !G 1 2 2 r . F1, 2 = ! F1, 2 r r

(28.3.3)

Recall that Newtons Third Law requires that the force on body 2 is equal in magnitude and opposite in direction to the force on body 1,

! ! F1, 2 = !F2,1 .
Newtons Second Law can be applied individually to the two bodies:

(28.3.4)

! ! d 2r1 F1, 2 = m1 2 , dt ! ! d 2r2 F2,1 = m2 2 . dt

Dividing through by the mass in each of Equations (28.3.5) and (28.3.6) yields
! F1, 2 ! d 2r1 = 2 , m1 dt

(28.3.5)

(28.3.6)

(28.3.7)

! F2,1

! d 2r2 = 2 . m2 dt

(28.3.8)

Subtracting the expression in Equation (28.3.8) from that in Equation (28.3.7) gives

! F1, 2

! ! ! d 2r1 d 2r2 d 2r ! = 2 ! 2 = 2 . m1 m2 dt dt dt

! F2,1

(28.3.9)

Using Newtons Third Law as given in Equation (28.3.4), Equation (28.3.9) becomes
! ! ! 1 1 " d 2r F1, 2 # + $= 2 . % m1 m2 & dt

(28.3.10)

Using the reduced mass , as defined in Equation (28.3.1) ,

1 1 1 = + , m1 m2
Equation (28.3.10) becomes

(28.3.11)

! F1, 2

! d 2r = 2 dt ! ! d 2r F1, 2 = 2 dt

(28.3.12)

! where F1, 2 is given by Equation (28.3.3).

Our result has a special interpretation using Newtons Second Law. Let be the ! with respect to an origin O , reduced mass of a reduced body with position vector r = r r is the unit vector pointing from the origin O to the reduced body. Then the where r equation of motion, Equation (28.3.12), implies that the body of reduced mass is under the influence of an attractive gravitational force pointing toward the origin. So, the original two-body gravitational problem has now been reduced to an equivalent one-body problem, involving a reduced body with reduced mass under the influence of a central . Note that in this reformulation, there is no body located at the central point force ! F1, 2 r (the origin O ). The parameter r in the two-body problem is the relative distance between the original two bodies, while the same parameter r in the one-body problem is the distance between the reduced body and the central point.

28.4 Energy and Angular Momentum, Constants of the Motion

Consider the reduced body with reduced mass given by Equation (28.3.1), orbiting about a central point under the influence of a radially attractive force given by Equation (28.3.3). The equivalent one-body problem has two constants of the motion, energy E and the angular momentum L about the origin O . Energy is a constant because there are 4

no external forces acting on the reduced body, and angular momentum is constant about the origin because the only force is directed towards the origin, and hence the torque about the origin due to that force is zero (the vector from the origin to the reduced body is anti-parallel to the force vector and sin ! = 0 ). Since angular momentum is constant, the orbit of the reduced body lies in a plane with the angular momentum vector pointing perpendicular to this plane. Choose polar coordinates ( r , ! ) for the reduced body (see Figure 28.3), where r is the distance of the reduced body from the central point that is now taken as the origin, ! is the angle that the reduced body makes with respect to a! chosen direction, and which . increases positively in the counterclockwise direction, and Fgrav = ! F1, 2 r

Figure 28.3 Coordinate system for the orbit of the reduced body. Since the force is conservative, the potential energy with choice of zero reference point U (!) = 0 is given by
U (r ) = ! G m1 m2 . r

(28.4.1)

The total energy E is constant, and the sum of the kinetic energy and the potential energy is
E= 1 2 G m1 m2 v ! . 2 r

(28.4.2)

The kinetic energy term, v 2 / 2 , has the reduced mass and the relative speed v of the two bodies. The velocity is cylindrical coordinates is given by (add link)
! , + vtan! v = vrad r v = vrad + vtan
2 2 2

" dr % " d! % = $ ' + $r ' # dt & # dt &

(28.4.3)

where vrad = dr / dt and vtan = r (d! / dt ) . Equation (28.A.2) then becomes

2 2 1 (! dr $ ! d' $ + G m1 m2 E = *# & + # r & - . . 2 *" dt % " dt % r ) ,

(28.4.4)

The magnitude of the angular momentum with respect to the center of mass is
L = r vtan = r 2 d! . dt

(28.4.5)

28.5 The Orbit Equation for the Reduced Body

There are two approaches to describing the motion of the reduced body. We can try to find both the distance from the origin, r (t ) and the angle, ! (t ) , as functions of the parameter time, but in most cases explicit functions cant be found analytically. We can also find the distance from the origin, r (! ) , as a function of the angle ! . This second approach offers a spatial description of the motion of the reduced body. We begin with Newtons Second Law

! ! F = a !G m1m2 r2
2 # d 2r # d" & & . = % 2 ! r% ( ( r r ( % dt $ ' dt $ '

(28.5.1)

Note that the motion is no longer circular, the radius is not constant so there is a term in the radial acceleration d 2 r / dt 2 along with the usual centripetal acceleration term ! r ( d" / dt )2 . Setting the components equal, using the constant of motion
L = r 2 ( d! / dt ) and rearranging,

d 2r L2 Gm1m2 . 2 = ! dt r3 r2

(28.5.2)

What we will do is use the same substitution u = 1/ r and change the independent variable from t to r , using the chain rule twice, since Equation (28.5.2) is a second-order equation. That is, the first time derivative is
dr dr du dr du d! = = . dt du dt du d! dt

(28.5.3)

From r = 1/ u we have dr / du = !1/ u 2 . Combining with the angular velocity d! / dt = L / r 2 in terms of L and u , d! / dt = L u 2 / , Equation (28.5.3) becomes

dr 1 du Lu 2 du L =" 2 =" , dt u d! d!
Taking the second derivative with respect to t ,
d 2 r d " dr # d " dr # d! = $ %= $ % dt 2 dt & dt ' d! & dt ' dt d 2u L " L # = ( 2 $ u2 % d! & ' =( d 2u 2 L2 u . d! 2 2

(28.5.4)

(28.5.5)

Substituting into Equation (28.5.2), and using r = 1/ u ,

d 2u 2 L2 L2 3 " 2u = u " Gm1m2 u 2 . d!

Canceling the common factor of u 2 and rearranging,

(28.5.6)

Gm1m2 d 2u + u = . d! 2 L2

(28.5.7)

Equation (28.5.7) is mathematically equivalent to the Harmonic Oscillator Equation with a constant term. The solution consists of two parts: the angle-independent solution
u0 =

Gm1m2 L2

(28.5.8)

and a sinusoidally varying term of the form

uH (! ) = A cos(! " ! 0 ) ,

(28.5.9)

where A and ! 0 are constants determined by the form of the orbit. The expression in Equation (28.5.8) is the inhomogeneous solution and represents a circular orbit. The expression in Equation (28.5.9) is the homogeneous solution (as hinted by the subscript) and must have two independent constants. Define
r0 ! L2 / G m1 m2 .

(28.5.10)

This is called the semilatus rectum. Then u0 = 1 / r0 and we can write the solution to Equation (28.5.7) as the sum of the inhomogeneous (Eq. (28.5.8)) and homogenous (Eq. (28.5.9)) pieces 7

u = u0 + uH =

1 (1 + r0 A(cos(! " ! 0 ))) r0

r0 r= . 1 + r0 A cos(! " ! 0 )

(28.5.11)

We have two constants to choose. Define A ! " / r0 where

! = 1 + 2 Er0 2 / L2 = 1 + 2 EL2 / (G m1 m2 )2

(28.5.12)

a quantity called eccentricity. Choose ! 0 = " then cos(! " # ) = " cos(! ) . Thus Equation (28.5.11) can be written in the form r0 r= . (28.5.13) 1 ! " cos# The two constants of the motion in terms of r0 and ! are

L = ( G m1 m2 r0 ) E= 2 r0

1 2

G m1 m2 (! 2 " 1)

(28.5.14)

The orbit equation as given in Equation (28.5.13) is a general conic section and is perhaps somewhat more familiar in Cartesian coordinates. Let x = r cos! and y = r sin ! , with r 2 = x 2 + y 2 . The orbit equation can be rewritten as
r = r0 + ! r cos" .

(28.5.15)

Using the Cartesian substitutions for x and y , rewrite Equation (28.5.15) as

( x 2 + y 2 )1/ 2 = r0 + ! x .

(28.5.16)

Squaring both sides of Equation (28.5.16),

x 2 + y 2 = r0 2 + 2 ! x r0 + ! 2 x 2 .

(28.5.17)

After rearranging terms, Equation (28.5.17) is the general expression of a conic section with axis on the x -axis,
x 2 (1 ! " 2 ) ! 2 " x r0 + y 2 = r0 2

(28.5.18)

(we now see that the dotted axis in Figure 28.3 can be taken to be the x -axis). For a given r0 > 0 , corresponding to a given nonzero angular momentum as in Equation (28.5.10), there are four cases determined by the value of the eccentricity. Case 1: When ! = 0 , E = Emin < 0 and r = r0 . Equation (28.5.18) is the equation for a circle,
x 2 + y 2 = r0 2

(28.5.19)

Case 2: When 0 < ! < 1 , Emin < E < 0 and Equation (28.5.18) describes an ellipse,
y 2 + A x2 ! B x = k

(28.5.20)

where A > 0 and k is a positive constant. (Appendix 17.C shows how this expression may be expressed in the more traditional form involving the coordinates of the center of the ellipse and the semimajor and semiminor axes.) Case 3: When ! = 1 , E = 0 and Equation (28.5.18) describes a parabola,
y 2 r0 ! . x= 2r0 2

(28.5.21)

Case 4: When ! > 1 , E > 0 and Equation (28.5.18) describes a hyperbola,

y 2 ! A x2 ! B x = k

(28.5.22)

where A > 0 and k is a positive constant.

28.6 Energy Diagram, Effective Potential Energy, and Orbits

The energy (Equation (28.A.7)) of the reduced body moving in two dimensions can be reinterpreted as the energy of a reduced body moving in one dimension, the radial direction r , in an effective potential energy given by two terms,
U eff L2 G m1 m2 . = ! 2 2 r r

(28.6.1)

The total energy is still the same, but our interpretation has changed;
E = K eff + U eff 1 ! dr " L2 G m1 m2 , = $ % + # 2 2 & dt ' 2 r r
2

(28.6.2)

where the effective kinetic energy K eff associated with the one-dimensional motion is
K eff 1 ! dr " = # $ . 2 % dt &
2

(28.6.3)

The graph of U eff as a function of r = r / r0 , where r0 as given in Equation (28.5.10), is shown in Figure 28.4. The upper curve (red, if you can see this in color) is proportional to L2 / (2 r 2 ) ! 1 / 2r 2 . The lower blue curve is proportional to !Gm1m2 / r ! !1/ r . The sum U eff is represented by the green curve. The minimum value of U eff is at r = r0 , as will be shown analytically below. The vertical scale is in units of !U eff (r0 ).

Figure 28.4 Graph of effective potential energy. Whenever the one-dimensional kinetic energy is zero, K eff = 0 , the energy is equal to the effective potential energy,
E = U eff = L2 G m1 m2 . ! 2 2 r r

(28.6.4)

Recall that the potential energy is defined to be the negative integral of the work done by the force. For our reduction to a one-body problem, using the effective potential, we will introduce an effective force such that
B B ! ! U eff , B ! U eff , A = ! # F eff " dr = ! # Freff dr A A

(28.6.5)

10

The fundamental theorem of calculus (for one variable) then states that the integral of the derivative of the effective potential energy function between two points is the effective potential energy difference between those two points,
B

U eff , B ! U eff , A = "

dU eff dr dr A

(28.6.6)

Comparing Equation (28.6.6) to Equation (28.6.5) shows that the radial component of the effective force is the negative of the derivative of the effective potential energy,
Freff = ! dU eff dr

(28.6.7)

The effective potential energy describes the potential energy for a reduced body moving in one dimension. (Note that the effective potential energy is only a function of the variable r and is independent of the variable ! ). There are two contributions to the effective potential energy, and the total radial component of the force is
F
eff r

d d ! L2 G m1 m2 " = # U eff = # $ # % 2 dr dr & 2 r r '

(28.6.8)

Thus there are two forces acting on the reduced body,

Freff = Fr , centifugal + Fr , gravity ,

(28.6.9)

with an effective centrifugal force given by

Fr ,centrifugal = # d ! L2 " L2 $ %= dr & 2 r 2 ' r 3

(28.6.10)

and the conventional gravitational force

Fr , gravity = ! G m1 m2 . r2

(28.6.11)

With this nomenclature, lets review the four cases presented in Section 17.3. Case 1: Circular Orbit E = Emin

11

The lowest energy state, Emin , corresponds to the minimum of the effective potential energy, Emin = (U eff )min . When this condition is satisfied the effective kinetic energy is zero since E = K eff + U eff . The condition
K eff 1 ! dr " = # $ =0 2 % dt &
2

(28.6.12)

implies that the radial velocity is zero, so the distance r from the central point is a constant. This is the condition for a circular orbit. The condition for the minimum of the effective potential energy is dU eff L2 G m1m2 . (28.6.13) 0= =! 3 + dr r r2 We can solve Equation (28.6.13) for r ,
L2 , r ! r0 = G m1m2

(28.6.14)

reproducing Equation (28.5.10). Case 2: Elliptic Orbit Emin < E < 0 When K eff = 0 , the mechanical energy is equal to the effective potential energy, E = U eff , as in Equation (28.6.4). Having dr / dt = 0 corresponds to a point of closest or furthest approach as seen in Figure 28.4. This condition corresponds to the minimum and maximum values of r for an elliptic orbit,
E= L2 G m1 m2 ! 2 r 2 r

(28.6.15)

Equation (28.6.15) is a quadratic equation for the distance r ,

G m1 m2 L2 r + r! =0 E 2 E
2

(28.6.16)

with two roots

2 G m1 m2 ! ! G m1 m2 " L2 " r=# + $$ % . % % $ & 2E ' 2 E 2E & ' 1/ 2

(28.6.17)

Equation (28.6.17) may be simplified somewhat as 12

1/ 2 2 ! " " 2 L E G m1 m2 ! #1 #1 + r=% $ $ 2 $ $ 2E # # ( G m m ) 1 2 ' ' & &

(28.6.18)

Recall from Equation (28.5.12), the square root is the eccentricity ! , given by

! = 1 + 2 Er0 2 / L2 = 1 + 2 EL2 / (G m1 m2 )2
Thus Equation (28.6.18) becomes
r=" G m1 m2 (1 ! ). 2E

(28.6.19)

(28.6.20)

A little algebra shows that

r0 = 1" ! 2 L2 / G m1 m2 # $ 2 L2 E 1 " %1 + % (G m1 m2 ) 2 & & ' ( L2 / Gm1m2 = "2 L2 E / (Gm1m2 ) 2 G m1 m2 . 2E

(28.6.21)

="

Substituting the last expression in (28.6.21) into Equation (28.6.20) gives an expression for the points of closest and furthest approach,
r= r0 (1 ! ) . 1" ! 2

(28.6.22)

The minus sign corresponds to the distance of closest approach,

r " rmin = r0 1+ !

(28.6.23)

and the plus sign corresponds to the distance of furthest approach,

r " rmax = r0 . 1# !

(28.6.24)

Case 3: Parabolic Orbit E = 0

13

The effective potential energy, as given in Equation (28.6.1), approaches zero ( U eff ! 0 ) when the distance r approaches infinity ( r ! " ). Since the total energy is zero, when r ! " the kinetic energy also approaches zero, K eff ! 0 . This corresponds to a parabolic orbit (see Equation (28.5.21)). Recall that in order for a body to escape from a planet, the body must have a total energy E = 0 (we set U eff = 0 at infinity). This escape velocity condition corresponds to a parabolic orbit. For a parabolic orbit, the body also has a distance of closest approach. This distance rpar can be found from the condition
E = U eff = G m1 m2 L2 ! =0. 2 2 r r

(28.6.25)

Solving Equation (28.6.25) for r yields

rpar = L2 1 = r0 ; 2 G m1 m2 2

(28.6.26)

the fact that the minimum distance to the origin (the focus of a parabola) is half the semilatus rectum is a well-known property of a parabola. Case 4: Hyperbolic Orbit E > 0 When E > 0 , in the limit as r ! " the kinetic energy is positive, K eff > 0 . This corresponds to a hyperbolic orbit (see Equation (28.5.22)). The condition for closest approach is similar to Equation (28.6.15) except that the energy is now positive. This implies that there is only one positive solution to the quadratic Equation (28.6.16), the distance of closest approach for the hyperbolic orbit
rhyp = r0 . 1+ !

(28.6.27)

The constant r0 is independent of the energy and from Equation (28.5.12) as the energy of the reduced body increases, the eccentricity increases, and hence from Equation (28.6.27), the distance of closest approach gets smaller.

28.7 Orbits of the Two Bodies

The orbit of the reduced body can be circular, elliptical, parabolic or hyperbolic, depending on the values of the two constants of the motion, the angular momentum and the energy. Once we have the explicit solution (in this discussion, r (! ) ) for the reduced body, we can find the actual orbits of the two bodies. 14

Choose a coordinate system as we did for the reduction of the two-body problem (Figure 28.5).

Figure 28.5 Center of mass coordinate system. The center of mass of the system is given by

! ! ! m1 r1 + m2 r2 . R cm = m1 + m2

(28.7.1)

! ! Let r1! be the vector from the center of mass to body 1 and r2! be the vector from the center of mass to body 2. Then, by the geometry in Figure 28.5,
! ! ! ! ! ! r = r1 " r2 = r1! " r2

(28.7.2)

and hence
! ! ! ! ! ! ! ! m1 r1 + m2 r2 m2 (r1 ! r2 ) ! r1" = r1 ! R cm = r1 ! = = r. m1 + m2 m1 + m2 m1

(28.7.3)

A similar calculation shows that

! ! ! = " r. r2 m2

(28.7.4)

Thus each body undergoes a motion about the center of mass in the same manner that the reduced body moves about the central point given by Equation (28.5.13). The only difference is that the distance from either body to the center of mass is shortened by a factor / mi . When the orbit of the reduced body is an ellipse, then the orbits of the two bodies are also ellipses, as shown in Figure 28.6.

15

Figure 28.6 The elliptical motion of bodies under mutual gravitation. When one mass is much smaller than the other, for example m1 ! m2 , then the reduced mass is approximately the smaller mass,

m1 m2 mm ! 1 2 = m1 m1 + m2 m2

(28.7.5)

The center of mass is located approximately at the position of the larger mass, body 2 of mass m2 . Thus body 1 moves according to

! ! ! r1! = r " r m1
and body 2 is approximately stationary,

(28.7.6)

! ! m ! ! ! =" r" 1 r #0 r2 m2 m2

(28.7.7)

16

28.8 Keplers Laws

Elliptic Orbit Law Each planet moves in an ellipse with the sun at one focus. When the energy is negative, E < 0 , and according to Equation (28.5.12),
1

" #2 2 E L2 ! = $1 + % $ (G m m )2 % 1 2 ' &

(28.8.1)

and the eccentricity must fall within the range 0 " ! < 1 . These orbits are either circles or ellipses. Note the elliptic orbit law is only valid if we assume that there is only one central force acting. We are ignoring the gravitational interactions due to all the other bodies in the universe, a necessary approximation for our analytic solution. Equal Area Law The radius vector from the sun to a planet sweeps out equal areas in equal time. Using analytic geometry, the sum of the areas of the triangles in Figure 28.7 is given by

"A =
in the limit of small "! .

1 (r "! ) "r "A = 1 r "! r + (r "! ) "r (r "! )r + ( ) 2 2 2 2

(28.8.2)

Figure 28.7 Keplers equal area law. The average rate of the change of area, !A , in time, !t , is given by

"A =

1 (r "! )r (r "! ) "r + . 2 "t 2 "t

(28.8.3)

17

In the limit as !t " 0 , "! # 0 , this becomes

dA 1 2 d! = r dt 2 dt

(28.8.4)

Recall that according to Equation (28.A.6) (reproduced below as Equation (28.8.5)), the angular momentum is related to the angular velocity d! / dt by

d! L = dt r 2
and Equation (28.8.4) is then

(28.8.5)

dA L = . dt 2

(28.8.6)

Since L and are constants, the rate of change of area with respect to time is a constant. This is often familiarly referred to by the expression: equal areas are swept out in equal times (see Keplers Laws at the beginning of this chapter). Period Law The period of revolution T of a planet about the sun is related to the major axis A of the ellipse by T 2 = k A3 where k is the same for all planets. When Kepler stated his period law for planetary orbits based on observation, he only noted the dependence on the larger mass of the sun. Since the mass of the sun is much greater than the mass of the planets, his observation is an excellent approximation. Equation (28.8.6) can be rewritten in the form
2 dA =L. dt

(28.8.7)

Equation (28.8.7) can be integrated as

orbit

2 dA = ! L dt
0

(28.8.8)

where T is the period of the orbit. For an ellipse, 18

Area =

orbit

"

dA = ! ab

(28.8.9)

where a is the semimajor axis and b is the semiminor axis. Thus we have
T= 2 ! ab . L

(28.8.10)

Squaring Equation (28.8.10) then yields

T2 =

4! 2 2 a 2b 2 . L2

(28.8.11)

In Appendix 17.B, the angular momentum is given in terms of the semimajor axis and the eccentricity by Equation (28.B.10). Substitution for the angular momentum into Equation (28.8.11) yields
T2 = 4! 2 2 a 2 b2 . Gm1 m2 a(1 " # 2 )

(28.8.12)

In Appendix 17.B, the semi-minor axis is given by Equation (28.B.22) which upon substitution into Equation (28.8.12) yields
T2 = 4! 2 2 a 3 . Gm1 m2

(28.8.13)

Using Equation (28.3.1) for reduced mass, the square of the period of the orbit is proportional to the semi-major axis cubed,
T2 = 4! 2 a 3 . G (m1 + m2 )

(28.8.14)

19

28.9 Worked Examples

28.9.1 Example Satellite Orbits A satellite of mass ms is initially in a circular orbit of radius r0 around the earth. The earth has mass me >> ms and radius Re . Let G denote the universal gravitational constant. Express all your answers in terms of Re , me , ms , G and r0 as needed.

a) Find an expression for the speed v0 of the satellite when it is in the circular orbit. b) Find an expression for the mechanical energy E0 of the satellite when it is in the circular orbit. Take U (r ) ! 0 as r ! " . As a result of an orbital maneuver the satellite trajectory is changed to an elliptical orbit. This is accomplished by firing a rocket for a short time interval thus increasing the tangential speed of the satellite. The apogee (farthest distance from earth) of the elliptical orbit is three times the closest approach (perigee),
ra = 3rp = 3r0 .

(In the figure, the small red circle represents the earth.) c) Use conservation of energy and angular momentum for the elliptic orbit to find an expression for the speed of the satellite, vp , immediately after the rocket has finished firing.

Solution:

20

a) This preliminary part should be found directly from Newtons Second Law and the Universal Law of Gravitation. The magnitude of the acceleration for the circular orbit is 2 v0 / r0 , and so

ms

2 v0 mm = G s2 e r0 r0

v0 = Gme / r0 .
b) The total mechanical energy is the sum of the kinetic energy and the gravitational potential energy,

E0 = !G

ms me 1 2 + msv0 r0 2

1 mm =! G s e. 2 r0
c) Since rp = r0 , and ra = 3r0 , the condition that angular momentum is constant
rp vp = ra va becomes r0 vp = (3r0 )va , so va = vp / 3 . The condition that the mechanical

energy is constant then becomes,

1 1 mm mm 2 2 ms vp ! G s e = ms va !G s e r0 2 ra 2 1 " vp % mm = ms $ ' ! G s e 3r0 2 # 3& 4 2 2 me vp = G 3 r0 9 vp = (3 / 2)Gme / r0 .

As a simple check, note that vp > v0 . As a further check, some minor algebra shows that after the rocket burn, the final mechanical energy is Ef = !Gms me / (4 r0 ) = !Gms me / A , where A = r0 + 3r0 is the major axis of the ellipse.
2

21

28.9.2 Example Halleys Comet The equation for any orbit in an inverse square gravitational field is given by
r= r0 1 # ! cos "

(28.9.1)

where

r0 =

L2 . Gm1m2

(28.9.2)

In Equation (28.9.2), L is the angular momentum, = m1m2 (m1 + m2 ) is the reduced mass, and ! is the eccentricity of the orbit. When 0 < ! < 1 , the orbit is an ellipse (in the above figure, ! = 3/ 5 ). The period, T , depends only on the length of the major axis, A , of the ellipse, which is given by
# 2T 2G (m1 + m2 ) $ A=% & !2 ( )
1/ 3

2r0 1' " 2

(28.9.3)

. Halleys Comet is in an elliptic orbit about the sun. The eccentricity of the orbit is ! = 0.967 and the period is T = 76 y . The mass of the sun is msun = 1.99 !1030 kg . The mass of Halleys Comet is negligible compared to the sun. a) Using this data, determine the distance of Halleys Comet at closest approach rp (perihelion) to the sun, and furthest distance ra (aphelion) from the sun. b) What is the speed vp of Halleys Comet when it is closest to the sun? Solution:

22

Before diving into the numerical calculations, lets do some preliminary math. First, note that when mcomet ! msun (note that mcomet is not given in the problem), m1 + m2 ! msun and

msun mcomet 1 = mcomet ! mcomet . msun + mcomet 1 + mcomet / msun

(28.9.4)

Next, from Equation (28.9.1), we have

rp = r0 , 1+ ! ra = r0 1" !

(28.9.5)

and combining with Equation (28.9.3) gives

rp = A r0 (1 " ! ), 2 ra = A r0 (1 + ! ). 2

(28.9.6)

As a quick check, note that rp + ra = A , the major axis. Next, anticipating part (b), we expect to use angular momentum considerations. At perihelion, the comet must be moving perpendicular to the vector from the sun to the comet, and so the magnitude of the angular momentum in terms of vp , rp and mcomet is
L = mcomet vp rp .

(28.9.7)

Combining with Equation (28.9.2) (with m1 = msun , m2 = mcomet , or vice versa) and the simplification for as given in Equation (28.9.4) and solving for vp gives
2 Gr0 msun mcomet Gr0 msun G (1 + ! )msun vp = . = = = mcomet rp mcomet rp rp rp

(28.9.8)

The mass of the comet does indeed drop out of this problem if mcomet ! msun . Its time to do the numbers. The calculations presented here were done by computer, keeping almost arbitrary precision in the intermediate calculations, and rounded to three figures (even though the period is given to only two figures). For the Newtonian gravitational constant, G = 6.67 "10!11 N # m 2 # kg !2 was used.

23

a) Since the major axis A is used for determining both rp and ra , find that quantity first. From Equation (28.9.3), with m1 + m2 = msun ,

# 2T 2Gmsun $ A=% & !2 ' (

1/ 3

# 2 (76 y ) 3.16 )107 s * y "1 )(6.67 )10"11 N * m 2 * kg "2 )( 1.99 )1030 kg )$ & =% % & !2 ' ( 12 = 5.37 )10 m
from which Equation (28.9.6) gives
rp = 8.86 !1010 m, ra = 5.28 !1012 m .

1/ 3

(28.9.9)

(28.9.10)

These results are roughly half and thirty times the earth-sun distance, respectively; rp is roughly the distance from the sun to the ex-planet Pluto. If fact, the period of Halleys Comet is roughly 1/ 8 the period of Plutos orbit, consistent with Equation (28.9.3). A graph of the orbit is shown here:

The tiny dot (red, if viewed in color) represents the sun and is not to scale; a circle representing the sun to scale on this scale is too small to be seen ( rp > 100 Rsun ). b) Equation (28.9.8) then gives, with the result in Equation (28.9.10),
vp = G (1 + ! )msun = 5.43 #104 m $ s "1 . rp

(28.9.11)

This is essentially (but of course smaller than) the escape velocity from the sun. In fact 1 2 2 / vescape = (1 + ! ) , which is 0.9835 for the eccentricity of its not too hard to show that vp 2 this orbit.

24

28.9.3 Example Lunar Orbit Collision A lunar mapping satellite of mass ms is in a circular orbit around the moon, and the orbit has radius R0 = 1.5 Rm where Rm is the radius of the moon. A repair robot of mass mr < ms is injected into that orbit, but due to a NASA sign error it orbits in the opposite direction. The two collide and stick together in a useless metal mass. The point of this problem is to find whether they create more junk orbiting the moon or crash into the lunar surface. The mass of the moon is denoted by mm . The universal gravitational constant is denoted by G.

a. What is the initial orbital velocity of the mapping satellite, v0 ? Express your answer in terms of R0 , mm , and G. b. What is the speed of the space junk (satellite and robot) immediately after the collision? Write it as f v0 , where you must determine the number f . Express your answer in terms of ms and mr . c. After the collision, the orbit of the space junk has changed. Use conservation of energy and angular momentum to find an equation for the closest approach rp of the space junk to the moon. d. Solve your equation in part c for the number f assuming that the closest approach rp = Rm , the space junk hits the moon. Solution: a) The speed of the mapping satellite undergoing uniform circular motion can be found from the force equation,

Gms mm R0 2

=!

ms v0 2 R0

(28.9.12)

25

So the speed is
v0 = Gmp / R0

(28.9.13)

b) The momentum flow diagram for the collision is shown below.

Because there are no external forces, the momentum is constant and so

ms v0 ! mr v0 = ( ms + mr )v f .
Thus the speed after the collision is
vf = ms ! mr ms + mr v0 = fv0 ,

(28.9.14)

(28.9.15)

where the ratio of speed after the collision to the speed before the collision is given by the number

f =

vf v0

ms ! mr . ms + mr

(28.9.16)

c) After the collision, the energy equation is given by

G (ms + mr )m p 1 G (ms + mr )m p 1 .(28.9.17) (ms + mr )v f 2 ! = (ms + mr )v p 2 ! 2 R0 2 rp

Setting v f = fv0 and simplifying yields

1 2 2 Gm p 1 2 Gm p . f v0 ! = vp ! 2 R0 2 rp

(28.9.18)

The angular momentum equation is

26

(ms + mr ) R0 v f = (ms + mr )rp v p .

(28.9.19)

Again setting v f = fv0 , Eq. (28.9.19) becomes

R0 fv0 = rp v p

(28.9.20)

Eq. (28.9.20) implies that v p = R0 fv0 / rp which we can substitute into Eq. (28.9.18) yielding
1 2 2 Gm p 1 ! R0 fv0 " Gm p . f v0 # = $ # % % R0 2$ r r 2 p p & '
2

(28.9.21)

Collecting terms yields

! 1 1" 1 2 2 ! R0 2 " f v0 $1 # 2 % = Gm p $ # % . $ % $ % 2 rp ' & R0 rp ' &

(28.9.22)

If we assume that at the closest approach rp = Rm , the space junk hits the moon. Then using the values rp = Rm , v0 = Gm p / R0 , and R0 = (3 / 2) Rm , Eq. (28.9.22) becomes
Gm p ! 9 RM 2 " ! 2 1 " f = Gm # $1 # % $ %. p 3RM & 4 RM 2 ' & 3RM RM '
2

(28.9.23)

We can solve Eq. (28.9.23) for f :

f = 4/5 .

(28.9.24)

27

28.9.4 Example Lowest Energy Solution

! is given by The effective potential energy for a linear restoring central force F = !kr r

U effective =

L2 1 + kr 2 . 2 2 r 2

Find the radius and the energy for the motion with the lowest energy. What type of motion does this correspond to? If the energy is slightly greater than the lowest energy, what type of motion would that correspond to?

Solution: Taking the derivative of the effective potential with respect to the radius r and setting the derivative equal to zero at r0 ,

L2 + k r0 = 0 r03
1/ 4

" L2 # r0 = $ % . & k '

The minimum energy is
L2 ! k " U eff (r0 ) = # $ 2 % L2 & ! L2 k " =# $ . % &
1/ 2 1/ 2 1/ 2

k ! L2 " + # $ 2 % k &

The orbit at this radius is a circle.

28

For a slightly larger energy, the orbit will oscillate about the minimum radius. Although it takes a bit more work to show, and is not part of this problem, the orbit will be an ellipse with the origin at the center (not at a focus, as in a Keplerian orbit).

29

28.9.5 Example Effective Potential Energy A system of two particles with a reduced mass interacts via an attractive central force ! , where r is the relative distance between the two particles. The magnitude of F = !ar 3 r the angular momentum for the equivalent one particle problem about the origin is L . a) Find an expression for the effective potential energy and make a graph of the effective potential energy as a function of r . Find an expression for r that minimizes the effective potential energy. b) Indicate by a horizontal line on the graph sketch the total energy that would correspond to a constant relative distance r between the two particles. Find an expression for this energy. What type of motion does this correspond to for the equivalent one-particle problem? Solution: a) The potential energy associated with the given force, denoted U with no subscript, is given by

! ! ! U = ! $! F ( r " ) # d r " =
! r r0

r0

a r " 3 dr " =

a 4 ( r ! r04 ) 4

(28.9.25)

where r0 is the radius at which we choose to set the potential energy equal to zero. While we could choose any (finite) value for r0 , for the purposes of making the graph well choose r0 = 0 . The effective potential energy is then
U eff = U + L2 1 a 4 L2 1 . = r + 2 r2 4 2 r2

(28.9.26)

From Equation (28.9.26) we see that the effective potential has both an attractive part ( ! r 4 ) and the repulsive part ( ! !1/ r 2 ). A plot of U eff as a function of r , with, for graphing purposes, a = L = = 1 is shown below. Plots of U and the repulsive term are shown as well; U eff is the upper curve, blue if viewed in color.

30

Differentiating Equation (28.9.26) and setting the derivative equal to zero at r = rmin yields

ar

3 min

L2 1 ! =0 3 rmin " L2 # =$ % . &a'

1/ 6

(28.9.27)

rmin

In the above plot, with a = L = = 1 , rmin = 1 , consistent with the graph. b) Using the value of rmin from (28.9.27) in Equation (28.9.26) yields

a ! L2 " U eff (rmin ) = # $ 4% a& = =

2/3

L2 ! L2 " + $ # 2 % a &

1/ 3

a1/ 3 L4 / 3 ! 1 1 " # + $ 2/3 % 4 2 & 3 a1/ 3 L4 / 3 . 4 2/3

(28.9.28)

This value is shown as the black horizontal line in the plot below; with a = L = = 1 ,
rmin = 1 , U eff (rmin ) = 3/ 4 .

31

The motion is a circle of radius rmin = 1 .

32

28.9.6 Example Transfer Orbit A space vehicle is in a circular orbit about the earth. The mass of the vehicle is ms and the radius of the orbit is 2 Re . It is desired to transfer the vehicle to a circular orbit of radius 4 Re . The mass of the earth is M e .

a) What is the difference in energy between the outer and inner circular orbits? b) An efficient way to accomplish the transfer is to use an elliptical orbit from point A at the inner circular orbit at to point B at the outer circular orbit (known as a Hohmann transfer orbit). This is accomplished by firing a rocket for a short time interval during each change of orbit thus increasing the tangential speed of the satellite. What velocity changes are required at the points of intersection, A and B ? c) Assume that the rocket burns fuel at a steady rate and the exhaust speed relative to the rocket is u . Using our early results for the rocket equation, how much fuel is burned at each of the rocket firings at the points of intersection, A and B ? Assume that before the firing at A , the mass of the rocket is m0 .

33

Solution: a) Because ms << m p , the reduced mass ! ms . The total mechanical energy is the sum of the kinetic and potential energies,
E = K +U = m M 1 ms v 2 ! G s e . 2 r

(28.9.29)

For a circular orbit, the orbital speed and orbital radius must be related by the radial component of Newtons Second Law

!G
We can rewrite Eq. (28.9.30) as

ms M e r2

= ! ms

v2 . r

(28.9.30)

1 1 m M ms v 2 = G s e . 2 2 r
Substituting Eq. (28.9.31) into Equation (28.9.29) gives
m M 1 m M 1 m M 1 E = G s e ! G s e = ! G s e = U (r ) . 2 r r 2 r 2

(28.9.31)

(28.9.32)

In moving from a circular orbit of radius 2 Re to an orbit of radius 4 Re , the total energy increases, as the energy becomes less negative. The change in energy is

!E = E (r = 4 Re ) " E (r = 2 Re ) 1 m M # 1 m M $ = " G s e "%" G s e & 2 4 Re 2 Re ( ' 2 1 m M = G s e. 8 Re

(28.9.33)

b) The satellite must increase its speed at point A in order to move to the larger orbit radius and increase its speed again at point B to stay in a circular orbit. Denote the satellite speed at point A while in the circular orbit as v A, i and after the speed increase (a rocket burn) as v A, f . Similarly, denote the satellites speed when it first reaches point B as v B , i and the speed of the satellite in the circular orbit at point B as v B , f . The speeds v A, i and v B , f are given by Equation (28.9.31).

34

While the satellite is moving from point A to point B (that is, during the transfer, after the first burn and before the second), both mechanical energy and angular momentum are conserved. Conservation of energy relates the speeds and radii by

m M m M 1 1 ms (v A, f ) 2 ! G s e = ms (vB , i ) 2 ! G s e . 2 2Re 2 4Re

Conservation of angular momentum relates the speeds and radii by
ms v A, f (2 Re ) = ms vB , i (4 Re ) ! v A, f = 2 vB , i .

(28.9.34)

(28.9.35)

Substitution of Equation (28.9.35) into Equation (28.9.34) yields, after minor algebra, v A, f = Equation (28.9.31) gives v A, i = 1 G Me , vB, f = 2 Re 1 G Me . 4 Re (28.9.37) 2 G Me 1 GM e . , vB, i = 3 Re 6 Re (28.9.36)

Thus the change in speeds at the respective points is given by

# 2 1 & G Me !v A = v A, f " v A, i = % " ( Re 2' $ 3

(28.9.38)

and
# 1 1 & G Me !v B = v B , f " v B , i = % " . ( Re 6' $ 4

(28.9.39)

Note that at both points, the speed must increase. c) Recall from our rocket equation (add correct link W08D2 Worked Example Rocket Problem) that the change in the rocket speed is given by
$ ' m0 !vr " vr , f # vr ,0 = u ln & ) % m0 # !m f (

(28.9.40)

where the mass of the rocket after the burn is mr , f = m0 ! "m f with !m f > 0 the amount of fuel that is necessary to burn in order to increase the speed by !vr . We can write Eq. (28.9.40) as 35

! Exponentiating both sides yields

# m0 ! "m f & "vr = ln % (. m0 u $ '

(28.9.41)

e
We can now solve for !m f :

! "vr / u

m0 ! "m f m0
" !vr / u

(28.9.42)

!m f = m0 (1 " e
So for the first firing at point A we have that

).

(28.9.43)

!m f , A = m0 (1 " e

" !v A / u

(28.9.44)

where !v A is given by Eq. (28.9.38). The amount of fuel after the firing is then

m0, B = m0 e

! "v A / u

(28.9.45)

The amount of fuel burned during the second firing at point B is then

!m f , B = m0, B (1 " e
where !vB is given by Eq. (28.9.39).

" !v B / u

) = m0 e

" !v A / u

(1 " e

" !v B / u

(28.9.46)

36

Chapter 28 Appendices
28.A: Derivation of the Orbit Equation
Two presentations of the result of Equation (28.5.13)

28.B: Properties of Elliptical Orbits

The dynamical properties of objects in elliptical orbits

28.C: Analytic Geometric Properties of Ellipses

Demonstrating how the results of Appendices 17.A and 17.B are consistent with more familiar representations of ellipses

28.D: Even More on Kepler Orbits

Using geometry, vector algebra, but minimal calculus to find the orbit equation, introducing the Laplace-Runge-Lenz vector.

Appendix 28.A: Derivation of the Orbit Equation Using Energy Methods

Consider the reduced body with reduced mass given by Equation (28.3.1), orbiting about a central point under the influence of a radially attractive force given by Equation (28.3.3). Since the force is conservative, the potential energy with choice of zero reference point U (!) = 0 is given by
U (r ) = ! G m1 m2 . r

(28.A.1)

The total energy E is constant, and the sum of the kinetic energy and the potential energy is
E= 1 2 G m1 m2 v ! . 2 r

(28.A.2)

The kinetic energy term, v 2 / 2 , has the reduced mass and the relative speed v of the two bodies. The velocity is cylindrical coordinates is given by (add link)
! , + vtan! v = vrad r v = vrad + vtan
2 2 2

" dr % " d! % = $ ' + $r ' # dt & # dt &

(28.A.3)

where vrad = dr / dt and vtan = r (d! / dt ) . Equation (28.A.2) then becomes 37

2 2 1 (! dr $ ! d' $ + G m1 m2 E = *# & + # r & - . . 2 *" dt % " dt % r ) ,

(28.A.4)

The magnitude of the angular momentum with respect to the center of mass is
L = r vtan = r 2 d! . dt

(28.A.5)

We shall explicitly eliminate the ! dependence from Equation (28.A.4) by using our expression in Equation (28.A.5),

d! L . = dt r 2
The mechanical energy as expressed in Equation (28.A.4) then becomes
1 ! dr " 1 L2 G m1 m2 . E = $ % + # 2 & dt ' 2 r 2 r
2

(28.A.6)

(28.A.7)

Equation (28.A.7) is a separable differential equation involving the variable r as a function of time t and can be solved for the first derivative dr / dt ,
dr = dt 2 ! 1 L2 G m1 m2 " 2 + E # $ % . 2 r2 r & '
1

(28.A.8)

Instead of solving for the position of the reduced body as a function of time, we shall find a geometric description of the orbit by finding r (! ) . We first divide Equation (28.A.6) by Equation (28.A.8) to obtain
L d! d! r2 . = dt = 1 1 dr dr 2 2 2 "2# " 1 L G m1 m2 # + dt % & % E $ & 2 2 r r '( ' (

(28.A.9)

The variables r and ! are separable; d! = L (1 / r 2 )

1/ 2

2 # G m1 m2 & L2 E " + 2 % ( r ' 2 r $

dr .

(28.A.10)

38

What follows involves a good deal of hindsight, allowing selection of convenient substitutions in the math in order to get a clean result. First, note the many factors of the reciprocal of r . Make the substitution u = 1/ r , du = !(1/ r 2 ) dr , with the result
d! = " L du . 1/ 2 2 2 # $ L 2 u + G m1 m2 u & %E" 2 ' (

(28.A.11)

Experience in evaluating integrals suggests that we make the absolute value of the factor multiplying u 2 inside the square root equal to unity. That is, multiplying numerator and denominator by 2 / L ,
d! = " du . (2 E / L " u + 2( G m1 m2 / L2 )u )1/ 2
2 2

(28.A.12)

As both a check and a motivation for the next steps, note that the left side (d! ) of Equation (28.A.12) is dimensionless, and so the right side must be. This means that the factor of G m1 m2 / L2 in the square root must have the same dimensions as u , or
length !1 ; so, define
r0 ! L2 / G m1 m2 .

(28.A.13)

a quantity we previously encountered in Eq. (28.5.10) called the semilatus rectum. The differential equation then becomes d! = " =" du (2 E / L " u 2 + 2u / r0 )1/ 2
2

du (2 E / L + 1 / r0 " u 2 + 2u / r0 " 1 / r0 2 )1/ 2

2 2

du =" 2 (2 E / L + 1 / r0 2 " (u " 1 / r0 )2 )1/ 2 =" r0 du (2 Er0 / L + 1 " ( r0 u " 1)2 )1/ 2
2 2

(28.A.14)

Next, we note that the combination of terms 2 Er0 2 / L2 + 1 is dimensionless, so define

! = 1 + 2 Er0 2 / L2

(28.A.15)

39

a quantity we also previously encountered (Eq. (28.5.12)) called the eccentricity. The last expression in is then
d! = " r0 du (# " ( r0 u " 1)2 )1/ 2
2

(28.A.16)

From here, well combine a few calculus steps, going immediately to the substitution r0 u # 1 = ! cos " , r0 du = #! sin " d" , yielding d! = " We can now integrate and find that "# sin $ d$ = d$ . (# 2 " # 2 cos 2 $ )1/ 2 (28.A.17)

! = " d! = " d# = # + constant

(28.A.18)

We have a choice in selecting the constant, and if we pick ! = " $ # , ! = " + # , cos ! = # cos" , the result is
r= 1 r0 = u 1 # ! cos"

(28.A.19)

agreeing with our result in Equation (28.5.13). Note that if we chose the constant of integration to be zero, the result would be
r= 1 r0 = u 1 + ! cos"

(28.A.20)

which is the same trajectory reflected about the vertical axis in Figure 28.3, indeed the same as rotating by ! .

Appendix 28.B: Properties of an Elliptical Orbit

We now consider the special case of an elliptical orbit.

40

Figure 28.B.1 Ellipse. In Figure 28.B.1, let a denote the semimajor axis, b denote the semiminor axis and x0 denote the distance from the center of the ellipse to the origin of our coordinate system ( r ,! ) . We shall now express the parameters a , b and x0 in terms of the constants of the motion L , E , , m1 and m2 . The semimajor axis: See Equation (28.A.19) above. The major axis A = 2a is given by
A = 2 a = rmax + rmin

(28.B.1)

where the distance of furthest approach occurs when ! = 0 , hence

rmax = r (! = 0) =

r0 1" #

(28.B.2)

and the distance of nearest approach occurs when ! = " , hence

rmin = r (! = " ) =

r0 1+ #

(28.B.3)

Figure 28.B.2: nearest and furthest approach Thus 41

1" r r # r a= $ 0 + 0 %= 0 2 . 2 ' 1& ! 1+ ! ( 1& !

(28.B.4)

The semilatus rectum r0 can be re-expressed in terms of the semimajor axis and the eccentricity,
r0 = a(1 ! " 2 ) .

(28.B.5)

We can now express the distance of nearest approach, Equation (28.B.3), in terms of the semimajor axis and the eccentricity,

rmin

a(1 " ! 2 ) = = = a(1 " ! ) . 1+ ! 1+ ! r0

(28.B.6)

In a similar fashion the distance of furthest approach is

rmax =

r0 1! "

a(1 ! " 2 ) = a 1+ " . 1! "

(28.B.7)

Figure 17.B.2 shows the distances of nearest and furthest approach. Using our results for r0 and ! from Equations (17.3.12) and (17.3.13), we have for the semimajor axis

L2 1 a= Gm1 m2 (1 ! (1 + 2 E L2 / G m m 2 )) 1 2

(28.B.8)

=!

G m1 m2 2E

and so the energy is determined by the semimajor axis,

E=! G m1 m2 . 2a

(28.B.9)

The angular momentum is related to the semilatus rectum r0 by Equation (17.3.12). Using Equation (28.B.5) for r0 , we can express the angular momentum (28.B.3) in terms of the semimajor axis and the eccentricity, L = Gm1 m2 r0 = Gm1 m2 a(1 ! " 2 ) . Note that (28.B.10)

42

(1 ! " 2 ) =

Gm1 m2 a

(28.B.11)

Location x0 of the center of the ellipse: From Figure 28.B.1, the distance from a focal point to the center of the ellipse is
x0 = rmax ! a .

(28.B.12)

Using Equation (28.B.7) for rmax , we have that

x0 = a (1 + ! ) " ! = ! a.

(28.B.13)

Thus, from Equations (17.3.12), (28.B.13) and (28.B.8),

x0 = ! a = "
The semi-minor axis: From Figure 28.B.1,

G m1 m2 2E

(1 + 2 E L2 / (G m1 m2 )2 ) .

(28.B.14)

b = ( rb 2 ! x0 2 ) where

(28.B.15)

rb =
which can be rewritten as

r0 , 1 # ! cos" b

(28.B.16)

rb # rb ! cos" b = r0 .

(28.B.17)

Note that from Figure 28.B.1,

x0 = rb cos! b ,

(28.B.18)

so that
rb = r0 + ! x0 .

(28.B.19)

43

Substituting Equation (28.B.13) for Equation (28.B.19) yields

x0

and Equation (28.B.5) for

r0

into

rb = a(1 ! " 2 ) + a" 2 = a .

(28.B.20)

The fact that rb = a is a well-known property of an ellipse reflected in the geometric construction, that the sum of the distances from the two foci to any point on the ellipse is a constant. Thus the semi-minor axis b becomes b = rb 2 ! x0 2 = a 2 ! " 2 a 2 = a 1 ! " 2 . Using Equation (28.B.11) for 1 " ! 2 , we have for the semi-minor axis (28.B.21)

b = aL2 / Gm1 m2 .
We can now use Equation (28.B.8) for a in the above expression, yielding
b = aL2 / Gm1 m2 = L ! G m1 m2 1 / Gm1 m2 = L ! 2E 2 E

(28.B.22)

(28.B.23)

Speeds at nearest and furthest distances: At nearest approach the velocity vector is tangent to the orbit, so the angular momentum is
L = rmin v p

(28.B.24)

and the speed at nearest approach is

v p = L / rmin .

(28.B.25)

Figure 28.B.3 Speeds at nearest and furthest approach

44

Using Equation (28.B.10) for the angular momentum and Equation (28.B.6) for rmin , Equation (28.B.25) becomes vp = L = rmin

Gm1 m2 (1 ! " 2 ) a(1 ! " )

Gm1 m2 (1 ! " 2 )

a(1 ! " )2

Gm1 m2 (1 + " ) .(28.B.26) a(1 ! " )

A similar calculation show that the speed va at furthest approach, L = rmax

va =

Gm1 m2 (1 ! " 2 ) a(1 + " )

Gm1 m2 (1 ! " 2 )

a(1 + " )

Gm1 m2 (1 ! " ) .(28.B.27) a(1 + " )

Appendix 28.C: Analytic Geometric Properties of Ellipses

Consider Equation (28.5.18), and for now take ! < 1 , so that the equation is that of an ellipse. It takes some, but not a great deal, of algebra to put this into the more familiar form

(x ! x0 )
a2

y2 =1 b2

(28.C.1)

where the ellipse has axes parallel to the x and y coordinate axes, center at ( x0 , 0) , semimajor axis a and semiminor axis b . We rewrite Equation (28.5.18) as

x2 "
We next complete the square,

2! r0 y2 r02 x + = . 1" ! 2 1" ! 2 1" ! 2

(28.C.2)

2" r0 " 2 r02 y2 r02 " 2 r02 x ! x+ + = + 1! "2 (1 ! " 2 )2 1 ! " 2 1 ! " 2 (1 ! " 2 )2
2

" r0 & y2 r02 # = x ! + % ( $ 1 ! " 2 ' 1 ! " 2 (1 ! " 2 )2

We now multiply each side by the factor (1 ! " 2 )2 / r02 yielding

(28.C.3)

" r0 & (1 ! " 2 )2 y 2 (1 ! " 2 )2 # x ! + 2 =1 % ( $ r0 (1 ! " 2 ) 1! "2 ' r02

Eq. (28.C.4) simplifies to

(28.C.4)

45

" r0 & # x! % ( $ 1! "2 '

# r0 & % ( $ (1 ! " 2 ) '
2

y2 # % $ r0 & ( 1! "2 '

2

=1

(28.C.5)

The last expression in (28.C.3) is the equation of an ellipse

( x ! " a )2 y 2 + 2 =1 a2 b
with semimajor axis
a= r0 , 1" ! 2

(28.C.6)

(28.C.7)

semiminor axis

b=

r0 1" !
2

= a 1" ! 2

(28.C.8)

and center at (! r0 / (1 " ! 2 ), 0) = (! a, 0) , as found in Equation (28.B.13).

Appendix 28.D: Even More on Kepler Orbits

Weve seen so far that for a Kepler Orbit, we have two constants of the motion: the angular momentum and the total energy. Since the angular momentum is a vector with three components, these constitute a total of four scalar constants of the motion. The Kepler Problem has six degrees of freedom (three position, three velocity), and so we expect to be able to find two more scalar constants of the motion. We might expect to be able to further identify any orbit by a vector in the plane of the orbit, perpendicular to the angular momentum, and this is indeed the case. Symmetry suggests that this vector would be along the major axis, and well see that this is the case as well. What follows uses a good deal of vector algebra, but minimal calculus, and leads to the orbit equation in a surprisingly simple form. Well need two results from vector algebra that we havent had to use yet. Specifically, ! ! ! for vectors a , b and c , we have
! ! ! ! ! ! ! ! ! a ! (b ! c ) = b(a " c ) # c (a " b) ! ! ! ! ! ! ! ! ! a " (b ! c ) = b " ( c ! a ) = c " (a ! b).

(28.D.1)

These relations are not hard to derive in Cartesian coordinates; the derivations will not be reproduced here. As a check, however, note that the vector on the right side of the first 46

! ! ! ! ! ! relation is perpendicular to both a and b ! c . If a , b and c are non-coplanar, the common magnitude of the scalars in the second relation is the expression for the volume of a parallelepiped with the three vectors forming the sides.
Lets start with the known constant angular momentum,
! ! ! ! ! L = r !p = r ! v

(28.D.2)

and re-express this quantity in a way that will allow us to use Newtons Laws. Specifically, consider the velocity in terms of polar coordinates,
dr ! dr +r v= r dt dt

(28.D.3)

so that the angular momentum can be expressed as

! " " dr % % dr ) ! $ r ' = r2 $ r ! ' L = (r r dt & # dt & #

(28.D.4)

from which

! 1 L dr ! = 2 . r dt r

(28.D.5)

The advantage to this operation is that we now have an explicit scalar factor of 1/ r 2 , which can and should be related to the same factor that appears in Newtons Law of , Gravitation. However, in order to use that law, we need a vector relation involving r into both sides of Equation (28.D.5), yielding and so well cross r
! " r L dr ! #$r # % = 2 # r dt ' r & dr 1 ! d2 !" ! ( =( r % # L. $ dt Gm1m2 & dt 2 '

(28.D.6)

In the above, the first relation in Equation (28.D.1) was used to simplify the left side, and ! was used on the Newtons Law of Gravitation, in the form ( d 2r / dt 2 ) = !(Gm1m2 / r 2 )r right side. Note the cancellation of the factor of the reduced mass . Equation (28.D.6) may now be integrated to obtain

47

! ! v!L ! = r +A Gm1m2 ! ! ! v!L " A=r , Gm1m2

(28.D.7)

! ! is in the plane of the orbit, and L is where A is a constant vector. Since r ! perpendicular to the plane of the orbit, A must lie in the plane of the orbit, as indicated ! , and hence above. Further, by considering extreme points of the orbit, where v ! r ! ! ! , we see that at these points A is in the direction parallel to the major axis. v!L " r ! ! Since A is a constant vector, A must always be in this direction. By considering the ! vector A at perihelion (or at any point on the orbit of a circular orbit), we can see that the ! direction of A is that from the perihelion point to the focus; well need this result below, when we find the orbit equation. ! The magnitude of A is readily found by calculating
! ! !r " A2 = A ! A = r 2 1 ! ! ! ( v # L) + r Gm1m2 Gm1m2 ! ! ! ! ( v # L ) ! ( v # L ) .(28.D.8) 2

The first dot product is manifestly 1 . The middle term, the cross term, is found using the second relation in Equation (28.D.1),

! 1 ! ! 1 2 ! ! ! ! ! ( v " L ) = L ! (r " v) = L. r L ! (r " v ) = r r

(28.D.9)

! ! ! ! The third term is most easily evaluated by recalling that v ! L , so that v ! L = vL and ! ! ! ! ( v ! L ) " ( v ! L ) = v 2 L2 . Combining, we see that
A2 = 1 ! = 1+ = 1+ 2 1 2 1 L + v 2 L2 2 Gm1m2 r (Gm1m2 ) " 1 2 Gm1m2 % 2 L2 v ! ' r & (Gm1m2 )2 $ #2 2 L2 E = ( 2. (Gm1m2 )2 (28.D.10)

! Thus, the constant vector A is directed along the major axis and has magnitude equal to the eccentricity. ! ! The orbit equation is now found algebraically by taking the dot product of A and r ;
48

! ! ! !r # A ! r = Ar cos" = r

1 ! ! ! r ! ( v $ L) . Gm1m2

(28.D.11)

The first term is merely the magnitude r of the position vector. The second term, ! , is repeating the calculation of Equation (28.D.9) with r instead of r

1 L2 = r0 , Gm1m2
with the result
Ar cos! = " r cos! = r # r0 .

(28.D.12)

(28.D.13)

Solving for r gives the orbit equation in the form

r= r0 . 1 # ! cos"

(28.D.14)

! It should be noted that what we call the vector A is a negative scalar multiple of the Laplace-Runge-Lenz vector (yes, it took three people to come up with this). Specifically, the L-R -L vector is in many sources given as
! ! ! ! ! r A LRL = p ! L " Gm1m2 = " Gm1m2 A . r

(28.D.15)

! ! ! Our choice of the form for A allows A = ! and the direction of A to lead to
Equation (28.D.11) without introduction of extra minus signs.

49

MIT OpenCourseWare http://ocw.mit.edu

8.01SC Physics I: Classical Mechanics

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

50