Physics 430 Lecture 19 Kepler Orbits Dale E

Physics 430: Lecture 19 Kepler Orbits Dale E. Gary NJIT Physics Department

8. 6 The Kepler Orbits q Last time we derived a general equation for the path of a body in the 2 -body central force problem: where, through a change of variables we substituted u = 1/r. q This is obeyed for any central force F(r), but let’s look specifically at the gravitational case (the Kepler problem), where, using g = Gm 1 m 2, we have q Inserting this into the path equation, we have the simpler, linear equation q The solution can be found by one last substitution, w(f) = u(f) - gm/l 2, which transforms the equation into our old friend again with solution w(f) = A cos(f - d). We will choose coordinates for which d = 0, so the final solution, then, is November 5, 2009

The Final Kepler Path q Finally, substituting for u = 1/r, we have Bounded Orbits q q q The dimensionless constant is going to play a big role in the shape of the orbit, depending on whether it is greater or less than 1. If e < 1, then the denominator is always positive for any value of f. If e > 1, there is a range of values of f for which the denominator vanishes, and r blows up (the object is unbound). So e = 1 is the demarcation between bound and unbound orbits. Because we want to talk about bound orbits, we will first take e < 1. In the above equation, as cos f oscillates between -1 and 1, the orbital distance r varies between rmin =perihelion (perigee) rmax = aphelion (apogee) November 5, 2009

Bounded Orbits, cont’d The shape of the orbit, then, looks like the figure at right. q We now want to prove that this shape is an ellipse, and to do that you will show in HW Prob. 8. 16 that planet q can be written in the form: star a c f d O b where The graphical meanings of a, b, c and d are shown in the figure. Here a is called the semi-major axis (half the longer axis) and b is the semi-minor axis. q The constant e is the eccentricity of the ellipse, and can be determined from q q Notice that as e 0, d goes to zero, a and b become equal, and the ellipse becomes a circle. As e 1 , d a , a and b/a 0, and the ellipse grows long and skinny (i. e. very eccentric). November 5, 2009

Example 8. 4: Halley’s Comet q Statement of the problem: n q Halley’s comet follows a very eccentric orbit, with e = 0. 967. Given that the closest approach to the Sun (perihelion) is 0. 59 AU (astronomical units), what is its greatest distance from the Sun? Solution: n Notice that rmax/rmin = (1 + e)/(1 – e). Therefore Orbital Period (Kepler’s third law) q Recall that Kepler’s second law (Chapter 3) states that the line between the Sun and a planet sweeps out equal areas in equal times, and is related to angular momentum l by q Since the total area of an ellipse is A = pab, the period is November 5, 2009

Orbital Period-2 q Squaring this and using the definitions of b and c given earlier: q Recall that, for the Sun, g = Gm 1 m 2 Gm. Msun, which then gives Kepler’s third law: q What is interesting is that this does not depend on the mass of the satellite, so the law is obeyed for all bodies (planets, comets, asteroids) so long as they do not get too massive relative to the Sun. Example 8. 5: Period of Low-Orbit Earth Satellite q Use Kepler’s third law to estimate the period of a satellite in a circular orbit close to Earth (a few 10’s of miles up). November 5, 2009

Relation between Energy and e q There is an important relation between the eccentricity and the energy of the orbit. To find this, think about the effective potential energy curve that we saw last time. At the closest approach (inner turning point), rmin, the total energy E = Ueff(rmin) (of course, we could alternatively use rmax): q Recalling that rmin = c/(1 + e), and substituting c = l 2/gm, we have: q Putting this into the above equation for energy, after some algebra: valid for any eccentricity November 5, 2009

8. 7 The Unbound Kepler Orbits q Going back to our original equation for the path , let’s now consider the case e 1 (which corresponds to E 0). In this case, the denominator blows up at some values of f, hence the orbits are unbound. q For the special case e = 1 (which corresponds to E = 0), we can convert the above to the cartesian form which is an equation for a parabola. For a parabola, the legs of the parabola eventually go parallel and at infinite distance they approach but never quite reach q For e > 1, the denominator blows up for some other value of f, such that q >1 =1 In this case, it can be shown that the cartesian form is a hyperbola: 0 where the legs go out at angles ±fmax. (the angles of the asymptotes). q The geometrical relationships are shown in the summary plot. November 5, 2009 <1 e

Summary of Kepler Orbits q Important relations of Kepler orbits are: e>1 path equation (hyperbola) e=1 (parabola) energy equation eccentricity energy orbit e=0 E<0 circle 0<e<1 E<0 ellipse e=1 E=0 parabola e>1 E>0 hyperbola e=0 (circle) e<1 (ellipse) scale factor for orbit November 5, 2009

8. 8 Changes of Orbit q Let’s first give the general approach to finding changes in the elliptical orbit of, say, a spacecraft orbiting Earth. The most general way of writing the path equation is where d is some inclination angle of the elliptical orbit, and the constants c, e and d are written with subscript 1 to indicate their initial values. q To change its orbit, a spacecraft can fire its engine in some particular direction for a brief time, thus causing an instantaneous change in velocity. From the change in velocity, we can calculate its new total energy and angular momentum, and thus calculate a new c 2 and e 2. The new orbit and the old orbit have to agree for some particular ro and fo, where the spacecraft was when its velocity changed, so we can calculate the remaining quantity d 2 by Hopefully you can see that this is straightforward, although tedious. q We can do a simple but interesting problem—a tangential thrust at perigee. Say the velocity changes from v 1 to v 2 = lv 1. q November 5, 2009

Tangential Thrust at Perigee q At perigee, cos(f – d) = 1, so the “continuity” between orbits gives q Because l is proportional to velocity, and the constant c is proportional to l, we have Therefore If l > 1 (increase in velocity), the new P orbit has a higher eccentricity, and a higher angular momentum and higher energy. q Likewise, if l < 1 (decrease in velocity), the new orbit is more circular, has a smaller angular momentum, and lower energy. q But what happens when the initial orbit is already circular ( e 1 = 0), and we decrease the velocity? q November 5, 2009 P

Example 8. 6: Changing between Circular Orbits If we want to go from, say, Earth to Mars, we have to boost out of an essentially circular Earth orbit onto an ellipse that takes us to Mars, and then go into another circular orbit to match that of Mars. The situation is shown in the figure at right. 2 q This is called a Hohman transfer ellipse, and is the transfer P orbit that takes the least energy. Let’s calculate the velocity P’ 1 changes needed to do this transfer. 3 q There are three orbits involved, two with zero eccentricity (e 1 = 0, e 3 = 0), with orbital radii R 1 and we’ll assume R 3 = 2 R 1. q The first match between orbits 1 and 2 requires q q The second match between orbits 2 and 3 requires q This is easily solved for l to give Need 15% of Earth orbital speed November 5, 2009