Derivation of Kepler's Third Law

1. Jan 2, 2009

ximath

1. The problem statement, all variables and given/known data

I am studying gravitation and I have been trying to derive Kepler's Third Law using Kepler's Second Law.

2. Relevant equations

The second law :

$$\frac{dA}{dt}$$ = $$\frac{L}{2m}$$

3. The attempt at a solution

To start with, I thought if we take $$\int$$ $$\frac{L}{2m}dt$$ from 0 to T; we would find Area of the ellipse.

We also know that $$A = \Pi * a * b$$. So, to find T, I could use these two equations.

Moreover, to eliminate b, eccentricity could be used.
$$e^{2} = 1 - \frac{b^{2}}{a^{2}}$$

substituting b would make

$$A = \Pi * a^{2} * \sqrt{1 - e^{2}}$$

$$\frac{L.T}{2m} = \Pi * a^{2} * \sqrt{1 - e^{2}}$$

and we know that

$$L = m.Vtan.R$$

where Vtan is tangential component of V. (to the radius R)

And we also know that since T is period, $$R = ea + a$$

Thus we can write

$$\frac{m.Vtan.(ea + a)}{2m} .T= \Pi * a^{2} * \sqrt{1 - e^{2}} (1)$$

Now, we need to express Vtan. If we write centripetal force, we can conclude that

$$Fnet = m.\frac{Vtan^{2}}{R} = \frac{G.m.Msun}{R^{2}}$$

where

$$R = ea + a$$

and Msun is the mass of the sun.

$$Vtan = \sqrt{\frac{G Msun}{ea + a}$$

and finally, we can rearrange equation (1) as

$$T = \frac{2.\Pi.a^{3/2}\sqrt{1-e^{2}}}{\sqrt{GMsun}.\sqrt{1+e}}$$

simplfying a little bit would yield to the result

$$T = \frac{2.\Pi.a^{3/2}\sqrt{1-e}}{\sqrt{GMsun}}$$

I seems to be close, but wrong!

Kepler's third law :

$$T = \frac{2.\Pi.a^{3/2}}{\sqrt{GMsun}}$$

So the term with squareroot and 1-e shouldn't have been there... I have been trying to find what I'm doing wrong, and I thought my assumption that

$$Vtan = \sqrt{\frac{G Msun}{ea + a}$$

seems to be wrong. Although it sounds logical to me... I am stuck here and need your help.

Last edited: Jan 2, 2009
2. Jan 2, 2009

LowlyPion

http://astro.berkeley.edu/~converse/Lagrange/Kepler%27sThirdLaw.htm [Broken]

Last edited by a moderator: May 3, 2017
3. Jan 2, 2009

ximath

I have checked it out. But I would like to correct my attempt.

Expressing Vtan in terms of centripetal force didn't work... I wonder why it is incorrect. Please help me expressing Vtan somehow that would make my attempt correct.

4. Jan 2, 2009

LowlyPion

Well for one thing your L = mvr is for circular orbits isn't it?

5. Jan 2, 2009

ximath

L = m V r is incorrect, what I assume is;

L = m Vtan r

where Vtan is the component of the velocity vector V which is momentarily perpendicular to the radius; since each planet can actually make circular motion, momentarily.

Moreover, L = mVtanr is used to derive kepler's second law, so I don't think that one is problematic.

6. Jan 3, 2009

7. Jan 4, 2009

ximath

I am desperately waiting for some answers...

8. Jan 29, 2010

fgb

Hi!

Wow, I was solving the same problem and I ended up here in this forum, so i joined just to be able to answer your question. It was the same problem I had and it took me some research to figure it out.

Indeed, your expression for Vtan is not quite right, as the one you used for the centripetal force, and that happens for a very subtle reason.

First of all, as you're dealing with the general case of an elliptical orbit, you have to remember that the distance R from the planet to the Sun IS NOT a constant, and therefore you can't really say R = ea + a because that is not true through the entire orbit. (unless you have the special case of a circle, and you can easily check that your derivation indeed holds for kepler's for e = 0).

Now, the actual Vtan for the generalized elliptical case is a bit trickier. As the distance planet-Sun (which I will call r) changes, there's an extra r'' term for the centripetal acceleration besides V²/r. [in fact, since for an ellipse V is not always perpendicular to R, the correct formula would be V²/(radius of curvature) instead of simply r.].

(r" stands for the second time derivative of r). I will also write V²/r as r*$$\omega^2$$, where $$\omega$$ stands for the angular velocity (which is not a constant).

$$a = - \frac{GM}{r^2}$$

$$\frac{d^2r}{dt^2} - r*\omega^2$$ $$= - \frac{GM}{r^2}$$

This is a Non-linear second order ordinary differential equation. Its solution, in fact, leads to Kepler's first law (if you solve it in terms of r as a function of $$\theta$$ you will find the equation of an ellipse).

I will not solve it here, but you can find its solution on wikipedia if I'm not mistaken, just search Kepler's Laws. From it's solution, you'd find that

$$L^2 = m^2*a*(1 - \epsilon^2)*GM$$ which will lead to the correct result.

I hope I was able to help :) English is not really my first language, so i'm not sure i was clear. Please let me know if there's anything you'd like me to rephrase.

9. Jan 29, 2010

jasper10

GMm/r^2 = mrw^2

GM/r^3 = 4pi^2 / T^2

T^2 / r^3 = 4pi^2/GM