Derivation of Kepler's Third Law

[ximath]

Homework Statement

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

Homework Equations

The second law :

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

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.

Thanks in advance.

 

[LowlyPion]

Perhaps this link will help?
http://astro.berkeley.edu/~converse/Lagrange/Kepler%27sThirdLaw.htm

 

[ximath]

Hello, thanks for your help.

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.

 

[LowlyPion]

Hello, thanks for your help.

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.

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

 

[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.

 

[ximath]

Anyone, please ?

 

[ximath]

I am desperately waiting for some answers...

 

[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.

 

[jasper10]

GMm/r^2 = mrw^2

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

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