Introduction to Orbital Mechanics

May 31, 2001

Introduction

In introductory physics classes, we learn Newtons laws of motion and his law of gravity, as well as Kepler,s laws for planetary motions. However, we usually make some simplications which do not necessarily reect on the true picture. Prominent among these simplications are: 1. All bodies are represented as point masses. 2. The gravitational eld of a body is spherically symmetric. 3. The central body (Sun or Earth) is xed in space. 4. There are only two bodies in any given problem, one of which doesnt move. 5. Masses dont change. The object of this class is to investigate these restrictions in the setting of satellites going around earth, however, we could easily extend this to the solar system. One other limitation of the way we study orbital motion in elementary physics is, that we look at the planet circling the the sun from an outside vantage point as if it where a clock work model. However, in reality we are sitting on one of the bodies and are a moving observer. The connections between the simple motions of earth and its satellite and the observations will also be studied.

Newtons Laws of motion and Gravity Keplers Laws, simple examples

2.1

Newtons Laws

All of you have previously encountered Newtons laws of motion. These laws form the basis of the part of mechanics which is known as dynamics. They are as follows: 1. A body does not change its state of rest or uniform straight-line motion unless it is compelled by some force to change that state. 2. The change of motion is proportional to the force and takes place in the same direction as the force. 3. Action is always contrary and equal to reaction. The rst and second laws of motion are usually combined in a single vector equation: F = ma = m dv , dt (1)

where F is the force, the scalar m is the mass of the body, and a its acceleration. This simplied version assumes tacitly that the mass is a constant. Often, in our daily life experiences, this is a valid assumption or at least a good approximation. But even, when you press the accelerator in a car the mass of the car changes due to the burning of fuel and rubber. However, this change is insignicant compared with the mass of the vehicle and can be neglected. Accelerating rockets is always accompanied by huge changes of the mass of the rocket itself. Furthermore, since space travel is often done at very high velocities, relativity eects can play a role, albeit a very minor one in the context of this class. To account for this we use the following form of Newtons law: dmv F= , (2) dt where v is the velocity of the body. It is easy to see that (1.2) implies (1.1) for a constant mass. Newton also published his law of gravity F= GM m r, |r|3 (3)

where F is the attracting force between the two bodies of mass M and mass m and r is a vector pointing from one body to the other. The constant G is the Gravitational constant which is measured at G = 6.672 1011 m3 /kgs2 . 2

If the bigger of the two bodies were xed at the origin of a coordinate system, and the mass m of the smaller body were xed, we could describe the motion of the smaller body by m GM m dv = r dt |r|3 dr = v dt (4) (5)

a coupled system of two rst order ordinary dierential equations. Since r and v are vectors, it is actually a coupled system of six ordinary dierential equations. Given initial conditions r0 and v0 we can describe all future positions of the body of mass m. This tells us that we usually need seven variables to describe the motion completely, time, three spatial components and three velocity components. How does the above model comply with Newtons third law? The above model is really a mode for the motion of a single body in a eld which has spherical symmetry. The second body of mass M is completely ignored. However, if we have two bodies, which are connected by a symmetric force (which Newtons third law implies) any motion of the smaller body will be accompanied by a motion of the larger body. We will, however, prove that in many cases this motion of the larger body can be ignored. But we will also give several signicant examples, where it cant be.

2.2

Keplers laws

Newton discovered his law of gravity around 1665. This gave nally an explanation of the motion of planets around the sun. Before Newton several quantitative descriptive models were given. First there was the Greek geocentric model using epicycles, a very complicated model, but highly successful in predicting the positions of planets for over 2000 years. Copernicus gave a heliocentric model, still insisting that planets move on circles. The most concise of these descriptive models was given around 1600 by Johannes Kepler. Precise observations by Tycho Brahe in the second half of the 16th century dispelled myth of circular motion. Based on these observations Kepler spelled out the following laws: 1. Planets move around the sun in elliptic orbits with the sun located at one focus of the ellipse. 2. As the planet moves in its orbit around the sun, equal areas as measured from the focus are swept out in equal times.

3. The square of the period of the orbit is proportional to the cube of the semimajor axis of the ellipse. We will prove these laws from Newtons laws later in the course. Keplers laws are purely descriptive and make no attempt to explain the motion. For this reason they are a more practical tool to study the motions of bodies in a gravitational eld. Keplers laws still allow for circular motions, however, they are not restricted to these motions. Nevertheless, most of the planets, the moon and many articial satellites actually travel along nearly circular orbits.

2.3

Circular Orbits

In this section we assume that the central body of mass ME is the earth and is xed (or moving uniformly). And the second body of mass m is a small articial satellite. The Earths mass ME = 5.9742 1024 kg Sputnik I on the other hand weighed 83.6 kg. So its gravitational pull on the Earth can probably be neglected. (This is certainly not true for its eect on world history or American politics). We will also assume that the satellite moves on a circular orbit. And that earth is at the center of this orbit. Let r be the radius of the circular orbit. The distance traveled after a time t is given by s(t) = r((t) (0)). Hence, v = r = r The area swept in a time t is given by A(t) = r2 (t). Keplers second law implies that this area is proportional to t. It follows that (t) = t for some constant angular velocity . Thus we have the rule: Satellites on circular orbits must move at a constant angular velocity. In polar coordinates we have then r(t) = and therefore r (t) = 2 r cos t r sin t 4 = 2 r(t). r cos t r sin t ,

(4) can thus be written as m 2 r(t) = or simply as 2 r3 = GME . If T is the period of this motion then T = 2. substituting this into (6) we get r3 GME = , (7) 2 T 4 2 this is of course Keplers third law in the special case of circular orbits. The expression GME is called the geocentric gravitational constant. Its value is GME = 3986004 108 m3 /s2 . Example: A satellite on circular orbit, has a period of 90 min. What is radius of the orbit and the altitude of the satellite. From (7) we get: r=
3

GME m r(t), r3 (6)

GME T 2 = 6652555.5m = 6652.6km. 4 2

The radius of the Earth is about 6360 km, hence the altitude is about 292 km.

2.4

Homework:

Communications satellites should stay in or close to the same place of the sky all the time. Otherwise one would have to constantly move the satellite dish to get good reception. Such orbits are called geo-stationary. Assume a circular geo-stationary orbit in the equatorial plane of the earth. 1. Find the period and the radius of this orbit. 2. Assume that the satellite is at a xed longitude of 118 west (same as LA) in what direction do you have to point your dish if you live in LA (32 north). 3. For the same satellite what direction do you have to point your dish if you live in Las Vegas (33 north and 112 west.)?

3
3.1

The two body problem

The object of this section is to correct some of the aspects of the computations in the preceding section. We rst drop the assumption that the central body of mass M is at rest or moves uniformly. We replace with the following more sensible assumption. Assumption: The center of mass of the two body system consisting is at rest or moves uniformly. One can always arrive at this assumption by decomposing the motion of the two body system into two pieces, a motion of the center of mass and a motion of the two components with respect to the center of mass. If we assume that there is no external forces (except the gravitational forces between the two bodies), Newtons rst law tells us that the center of mass is either at rest or moves uniformly. In the case of a satellite moving around Earth, this assumption is not quite correct, since the center of mass moves in the gravitational eld of the sun, hence it is not a uniform motion with constant speed. However, on the surface of the Earth the acceleration due to the gravitational pull of the sun is approximately .0059m/s2 compared to Earth gravity of 9.81m/s2 , and we will ignore it for the time being. We center our coordinate system at the center of gravity, which sits on the line connecting the two bodies. In circular motion the distance between the two bodies is constant r. We have r = r1 + r2 , where r1 is the distance of the body of mass M from the center of mass, and r2 is the distance of the body of mass m from the center of mass. We have from the denition of the center of mass: M r1 = mr2 . Furthermore, we have the following balance of centrifugal forces: Mm mr2 2 = G 2 r M m M r1 2 = G 2 r Adding these equations and using the property of the center of mass we get: r3 2 = (r1 + r2 )3 2 = G(M + m), m (8) M This is the corrected form of Keplers third law for circular orbits. One sees that this is now dependent on the mass of the satellite. In the case Sputnik I, the correction factor m 83.6 = = 1.44 1023 M 5.9742 1024 r3 2 = GM 1 + 6 or

which can be denitely ignored. Examples: The ratio of the mass of the moon to the mass of the earth is ML = 0.0123. ME The radius of the moons nearly circular orbit is about 384400 km. This implies that r1 = 0.0123 384400 = 472 i.e the center of the rotation is 472 km above the center of the earth. The orbital period of the moon is 27.3 86400 = 2358720s i. e. = 2.664 106 s1 . This produces an acceleration pointed away from the center of the earth and opposite the center of the earth of 2 (RE + r1 ) = 4.86 105 m/s2 , which is responsible for the second tide. Another more exciting example is from the search for extra solar planets. Planets do not radiate any energy, and it is therefore impossible to directly observe planets circling around stars. However, if a large planet of mass mP circles a star of mass MS with period T we have the following: r = rS + rP MS rS = m P rP G r3 = (MS + mP ) 2 T 4 2 The star will radiate on a dierent frequency, depending on its velocity relative to the observer. If the star moves on a circular path (in the same plane with the observer) with speed v = rS = 2rS /T , the maximal dierence between velocities is 2v . This can be computed from the shift in frequencies. One can also directly observe the period T . From this one can compute the ratio of the two masses and rS . Often one can also estimate MS from the total radiation of the star. This then completely determines the orbit and size of the planet.

3.2

The General two Body Problem I Elementary Conservation Laws

In this section we will investigate the two body problem in more general terms. To do this consider the situation in the gure below. There two bodies of masses m1 and m2 are at 7

locations r1 and r2 . We also assume that the only acting force is the gravitational attraction between the two bodies: F21 = G m1 m2 m1 m2 (r1 r2 ) = G 3 (r1 r2 ), 3 |r1 r2 | r

is the force acting on the body m2 , and F12 = F21 is the force acting on m1 .

Newtons law implies m1 d2 r1 = F12 dt2 d2 r2 m2 2 = F21 dt (9) (10)

Adding these two equations we get: 0 = m1 or d2 r1 d2 r2 d + m = 2 2 2 dt dt dt P = m1 m1 dr1 dr2 + m2 dt dt ,

dr1 dr2 + m2 dt dt is constant. This is known as conservation of momentum. It can be further interpreted as d = 0, where dt m1 m2 = r1 + r2 , m1 + m2 m1 + m2 is the location of the center of mass of the two bodies. We have does established Conservation of Linear Momentum: The center of mass of a two body system with no exterior forces moves at a uniform speed. Since the motion of the center of mass is trivial we will ignore it for the remainder of this derivation. We will put the the origin of our reference frame at the center of mass. (See the Figure below.)

To continue we multiply the equations motion by m2 and m1 respectively and subtract the from each other to get d2 G(m1 + m2 ) (r1 r2 ) = (r1 r2 ) . 2 dt |r1 r2 |3 We can further simplify this by introducing r1 and r2 , where rj is the distance between mj and the center of mass and set r1 r2 = (r1 + r2 ) r Furthermore, we introduce the reduced masses: m2 2 m1 + m2 m2 1 MR (m2 ) = m1 + m2 MR (m1 ) = If we only consider the motion of the body m1 , we can write: m1 and if we consider the motion of m2 : d2 r2 Gm2 MR (m2 ) r m2 2 = 2 dt r2 If one considers the motion of a small satellite of mass m around earth,. The reduced mass MR is nearly equal to the mass of the earth and the problem reduces to a one body problem. To continue we consider the motion of the mass m1 : m1 d2 r1 Gm1 MR = r1 3 2 dt r1
dr1 dt

(11) (12)

d2 r1 Gm1 MR (m1 ) = r 2 2 dt r1

This equation can be skalar multiplied by m1 d 2 dt dr1 dr1 dt dt

to get Gm1 MR dr1 r1 3 r1 dt d 1 = Gm1 MR . dt r1 =

1 Gm1 MR 2 m1 v1 = K. (13) 2 r1 The quantity on the left of this equation is called the total energy, and we have established: Conservation of Total Energy: The total energy of the two body problem is constant. 9

or