About Opentopia
Opentopia Directory Encyclopedia Tools

Kepler's laws of planetary motion

Encyclopedia : K : KE : KEP : Kepler's laws of planetary motion

 

Johannes Kepler's primary contributions to astronomy/astrophysics were his three laws of planetary motion. Kepler, a nearly blind though brilliant German mathematician, derived these laws, in part, by studying the observations of the keen-sighted Danish astronomer Tycho Brahe. The article on Johannes Kepler gives a less mathematical description of the laws, as well as a treatment of their historical and intellectual context.

Sir Isaac Newton's later discovery of the laws of motion and universal gravitation depended strongly on Kepler's work. Although from the modern point of view, Kepler's laws can be seen as a consequence of Newton's laws, historically, it was the other way around: Kepler provided a kinematic mathematical model of the empirical observations, which Newton then interpreted using calculus and his new physics.

Contents

Kepler's first law

Kepler's first law
The orbit of a planet about a star is an ellipse with the star at one focus.

There is no object at the other focus of a planet's orbit. The semimajor axis, a, is half the major axis of the ellipse. In some sense it can be regarded as the average distance between the planet and its star, but it is not the time average in a strict sense, as more time is spent near apocentre than near pericentre.

Kepler's second law

Kepler's second law

''A line joining a planet and its star sweeps out equal areas during equal intervals of time This is also known as the law of equal areas.

Suppose a planet takes 1 day to travel from points A to B. During this time, an imaginary line, from the Sun '''to the planet, will sweep out a roughly triangular area. This same amount of area will be swept every day regardless of where in its orbit the planet is.

As a planet travels in its elliptical orbit, its distance from the Sun will vary. As an equal area is swept during any period of time and since the distance from a planet to its orbiting star varies, one can conclude that in order for the area being swept to remain constant, a planet must vary in speed. The physical meaning of the law is that the planet moves faster when it is closer to the sun. This is because the sun's gravity accelerates the planet as it falls toward the sun, and decelerates it on the way back out.

Kepler's third law

The squares of the orbital periods of planets are directly proportional to the cubes of the semi-major axis of the orbits.
[T^2 \propto a^3]
[T] = orbital period of planet
[a] = semimajor axis of orbit
So the expression T2a–3 has the same value for all planets in the solar system as it has for Earth.

That value is (with T in seconds, a in meters) [3.00\times 10^ \frac}} \pm \ 0.7%\, ].

Thus, not only does the length of the orbit increase with distance, the orbital speed decreases, so that the increase of the sidereal period is more than proportional.

See the actual figures: attributes of major planets.

This law is also known as the harmonic law.

Accuracy and limitations

As discussed below, Kepler's third law needs to be modified when the orbiting body's mass is not negligible compared to the mass of the body being orbited. However, the correction is fairly small in the case of the planets orbiting the Sun.

A more serious limitation of Kepler's laws is that they do not work well in a system consisting of more than two bodies. They were designed to give a description of the motion of the planets around the sun, which is satisfactory for most planets but a particularly bad approximation in the case of the Earth-Sun-Moon system. For calculations of the Moon's orbit, Kepler's laws are far less accurate than the empirical method invented by Ptolemy more than a thousand years before.

Newton generalized Kepler's first law with the realization that an object moving at greater than escape velocity (e.g., some comets) would have an open parabolic or hyperbolic orbit rather than a closed elliptical one. Thus, all of the conic sections are possible orbits. The second law is still valid for open orbits (since angular momentum is still conserved), but the third law is inapplicable because the motion is not periodic.

Because electrical forces also obey an inverse square law, Kepler's laws also apply to bodies interacting electrically. If the interaction is repulsive, as in Rutherford scattering, then the orbit is hyperbolic, with the center of repulsion outside the hyperbola. In the planetary model of the atom, the electrons were imagined to orbit the nucleus exactly like the planets orbiting the Sun, and Niels Bohr even went so far as to attempt to assign specific atomic energy levels to specific elliptical orbits. The Keplerian approximation is successful when the quantum-mechanical wavelengths are small compared to the size of the orbit (as in Rutherford scattering, and in some unusual artificially populated electron orbitals with very high values of the principal quantum number).

Kepler's laws do not consider the emission of radiation. The emission of gravitational radiation is negligible in our solar system, but important in some stellar systems containing black holes or neutron stars (any system involving large masses combined with large accelerations). In the planetary model of the atom, the emission of electromagnetic radiation should have led to the collapse of all atoms, and this was a hint that classical physics was in need of modification.

Kepler's laws don't incorporate relativity, so, for example, they don't correctly predict the precession of Mercury's orbit. The successful explanation of this non-Keplerian behavior was one of the most dramatic early successes of general relativity. Kepler's laws fail even more dramatically in the case of an object orbiting a black hole; here, general relativity predicts that many orbits —those that pass through the event horizon— are one-way streets to a dead-end region of spacetime.

Connection to Newton's laws and conservation laws

Kepler did not understand why his laws were correct; it was Isaac Newton who discovered the answer to this more than fifty years later. The second law can also be seen as a statement of conservation of angular momentum, which is a logical consequence of Newton's laws in the special case of a force that acts along the line connecting two objects.

Kepler's first law

Newton said that "every object in the universe attracts every other object along a line of the centres of the objects, proportional to each object's mass, and inversely proportional to the square of the distance between the objects."

The following assumes that acceleration a is of the form

[ \mathbf = \frac} = f(r)\hat}.]
Recall that in polar coordinates

[\frac} = \dot r\hat} + r\dot\theta\hat,]
[\frac} = (\ddot r - r\dot\theta^2)\hat} + (r\ddot\theta + 2\dot r \dot\theta)\hat.]
In component form, we have

[\ddot r - r\dot\theta^2 = f(r),]
[r\ddot\theta + 2\dot r\dot\theta = 0.]
Substituting for [\ddot \theta] and [\dot r] in the second equation, we have

[r + 2 \dot\theta = 0]
which simplifies to

[\frac = -2\frac.]
When integrated, this yields

[\log\dot\theta = -2\log r + \log\ell,]
[ \log\ell = \log r^2 + \log\dot\theta, ]
[\ell = r^2\dot\theta],
for some constant [\ell], which can be shown to be the specific angular momentum. Now we substitute. Let

[r = \frac,]
[\dot r = -\frac\dot u = -\frac\frac\frac= -\ell\frac,]
[\ddot r = -\ell\frac\frac = -\ell\dot\theta\frac= -\ell^2u^2\frac.]
The equation of motion in the [\hat}] direction becomes

[\frac + u = - \fracf\left(\frac\right).]
Newton's law of gravitation relates the force per unit mass to the radial distance as

[ f \left( \right) = f(r)= - \, = - GM u^2 ]
where G is the constant of universal gravitation and M is the mass of the star.

As a result,

[\frac + u = \frac]
This differential equation has the general solution:

[u = \frac \bigg[ 1 + ecos(theta-theta_0) bigg] .]
for arbitrary constants of integration e and θ0.

Replacing u with 1/r and letting θ0 = 0:

[r = = \frac ]
This is indeed the equation of a conic section with eccentricity e and the origin at one focus. Thus, Kepler's first law is a direct result of Newton's law of gravitation and Newton's second law of motion.

Kepler's second law

Assuming Newton's laws of motion, we can show that Kepler's second law is consistent. By definition, the angular momentum [\mathbf] of a point mass with mass [m] and velocity [\mathbf] is :

[\mathbf \equiv \mathbf \times \mathbf = \mathbf \times ( m \mathbf )].
where [\mathbf] is the position vector of the particle and [\mathbf = m \mathbf ] is the momentum of the particle.

By definition,

[\mathbf = \frac} ].
As a result, we have
[\mathbf = \mathbf \times m\frac}].
Taking the derivative of both sides with respect to time, we have

[\frac} = (\mathbf \times \mathbf) + \left( \frac} \times m\frac} \right)= ( \mathbf \times \mathbf ) + ( \mathbf \times \mathbf ) = 0 ]
since the cross product of parallel vectors is 0. Note that F is always parallel to r, since the force is entirely radial, and p is always parallel to v by definition. Therefore, we can now say that [|\mathbf|] is constant.

The area swept out by the line joining the planet and the sun is half the area of the parallelogram formed by [\mathbf] and [d\mathbf].

[dA = \begin\frac\end |\mathbf \times d\mathbf| = \begin\frac\end \left|\mathbf \times \frac}dt\right| = \fracdt]
Since [|\mathbf|] is constant, the area swept out by the planet in a time interval is also a constant. Q.E.D.

Kepler's third law

Newton used the third law as one of the pieces of evidence used to build the conceptual and mathematical framework supporting his law of gravity. If we take Newton's laws of motion as given, and consider a hypothetical planet that happens to be in a perfectly circular orbit of radius r, then we have [F=mv^2/r] for the sun's force on the planet. The velocity is proportional to r/T, which by Kepler's third law varies as one over the square root of r. Substituting this into the equation for the force, we find that the gravitational force is proportional to one over r squared. (Newton's actual historical chain of reasoning is not known with certainty, because in his writing he tended to erase any traces of how he had reached his conclusions.) Reversing the direction of reasoning, we can consider this as a proof of Kepler's third law based on Newton's law of gravity, and taking care of the proportionality factors that were neglected in the argument above, we have:

[T^2 = \frac \cdot r^3]
where: The same arguments can be applied to any object orbiting any other object. This discussion implicitly assumed that the planet orbits around the stationary sun, although in reality both the planet and the sun revolve around their common center of mass. Newton recognized this, and modified this third law, noting that the period is also affected by the orbiting body's mass. However typically the central body is so much more massive that the orbiting body's mass may be ignored. Newton also proved that in the case of an elliptical orbit, the semimajor axis could be substituted for the radius. The most general result is:

[T^2 = \frac \cdot a^3]
where: For objects orbiting the sun, it can be convenient to use units of years, AU, and solar masses, so that G, 4π2 and the various conversion factors cancel out. Also with m<<M we can set m+M = M, so we have simply [T^2=a^3]. Note that the values of G and planetary masses are not known with good accuracy; however, the products GM (the Keplerian attraction) are known to extremely high precision.

Proving Kepler's third law

Define point A to be the periapsis, and point B as the apoapsis of the planet when orbiting the sun.

Kepler's second law states that the orbiting body will sweep out equal areas in equal quantities of time. If we now look at a very small periods of time at the moments when the planet is at points A and B, then we can approximate the area swept out as a triangle with an altitude equal to the distance between the planet and the sun, and the base equal to the time times the speed of the planet.

[\begin\frac\end \cdot(1-\epsilon)a\cdot V_A\,dt= \begin\frac\end \cdot(1+\epsilon)a\cdot V_B\,dt]
[(1-\epsilon)\cdot V_A=(1+\epsilon)\cdot V_B]
[V_A=V_B\cdot\frac]
Using the law of conservation of energy for the total energy of the planet at points A and B,

[\frac-\frac =\frac-\frac]
[\frac-\frac =\frac-\frac]
[\frac=\frac\cdot \left ( \frac-\frac \right ) ]
[\frac\right ) ^2-V_B^2}=\frac\cdot \left ( \frac \right ) ]
[V_B^2 \cdot \left ( \frac\right ) ^2-V_B^2=\frac\cdot \left ( \frac \right ) ]
[V_B^2 \cdot \left ( \frac\right )=\frac ]
[V_B^2 \cdot \left ( \frac \right) =\frac ]
[V_B^2 \cdot 4\epsilon =\frac ]
[V_B =\sqrt}.]
Now that we have [V_B], we can find the rate at which the planet is sweeping out area in the ellipse. This rate remains constant, so we can derive it from any point we want, specifically from point B.

[\frac=\frac\cdot(1+\epsilon)a\cdot V_B \,dt}= \begin\frac\end \cdot(1+\epsilon)a\cdot V_B ]
:[= \begin\frac\end \cdot(1+\epsilon)a\cdot \sqrt} = \begin\frac\end \cdot\sqrt]
However, the total area of the ellipse is equal to [\pi a \sqrta]. (That's the same as [\pi a b], because [b=\sqrta]). The time the planet take out to sweep out the entire area of the ellipse equals the ellipse's area, so,

[T\cdot \frac=\pi a \sqrta]
[T\cdot \begin\frac\end \cdot\sqrt=\pi \sqrta^2]
[T=\fraca^2}} =\frac}=\frac}\sqrt]
[T^2=\fraca^3.]
However, if the mass m is not negligible in relation to M, then the planet will orbit the sun with the exact same velocity and position as a very small body orbiting an object of mass [M+m] (see reduced mass). To integrate that in the above formula, M must be replaced with [M+m], to give

[T^2=\fraca^3.]
Q.E.D.

Solution for the motion as a function of time

The Keplerian problem assumes an orbit with semimajor axis a, semiminor axis b, and eccentricity ε. To convert the laws into predictions, Kepler began by adding the orbit's auxiliary circle (that with the major axis as a diameter) and defined these points:

and three angles measured from perihelion:

Then

[\mboxcxz=\mboxcxs+\mboxsxz][=\mboxcxs+\mboxcyz]
[\frac2E=a\varepsilon\frac a2\sin E+\frac2M]
giving Kepler's equation

[M=E-\varepsilon\sin E].
To connect E and T, assume [r=\mboxsp] then

[a\varepsilon+r\cos T=a\cos E] and [r\sin T=b\sin E]
[r=\frac=\frac]
[\tan T=\frac=\frac ba\frac=\frac\sin E}]
which is ambiguous but useable. A better form follows by some trickery with trigonometric identities:

[\tan\frac T2=\sqrt\frac\tan\frac E2]
(So far only laws of geometry have been used.)

Note that [\mboxspz] is the area swept since perihelion; by the second law, that is proportional to time since perihelion. But we defined [\mboxspz=\frac ba\mboxcyz=\frac ba\frac2M] and so M is also proportional to time since perihelion—this is why it was introduced.

We now have a connection between time and position in the orbit. The catch is that Kepler's equation cannot be rearranged to isolate E; going in the time-to-position direction requires an iteration (such as Newton's method) or an approximate expression, such as

[E\approx M+\left(\varepsilon-\frac18\varepsilon^3\right)\sin M+\frac12\varepsilon^2\sin 2M+\frac38\varepsilon^3\sin 3M]
via the Lagrange reversion theorem. For the small ε typical of the planets (except Pluto) such series are quite accurate with only a few terms; one could even develop a series computing T directly from M.[link]

See also

External links

 

From Wikipedia, the Free Encyclopedia. Original article here. Support Wikipedia by contributing or donating.
All text is available under the terms of the GNU Free Documentation License See Wikipedia Copyrights for details.

Search Titles
0123456789
ABCDEFGHIJ
KLMNOPQRST
UVWXYZ?

E-mail this article to:

Personal Message: