Physics+ Kepler Third Law: UNIZOR.COM - Physics+ 4 All - Laws of Newton

Notes to a video lecture on UNIZOR.COM

Laws of Newton -
Kepler's Third Law

Kepler's Third Law states that for all objects moving around a fixed source of gravitational field along elliptical orbits the ratio of a square of their period of rotation to a cube of a semi-major axis is the same.

As in the case of the Kepler's First Law, this Third Law has been based on numerous experiments and years of observation.

Based on all the knowledge conveyed in previous lectures on Kepler's Laws, we will derive this Third Law theoretically.

Let's make a simple derivation of Kepler's Third Law in case of a circular orbit.

In this case the velocity vector of an object circulating around a central point is always perpendicular to a position vector from a center to an object.
Since the gravitational force is collinear with a position vector, it is also perpendicular to velocity, which is tangential to a circular orbit. Therefore, gravitational force makes no action along a velocity vector which makes the magnitude of the velocity vector constant.

Let's introduce the following characteristics of motion:
t - absolute time,
r - radius of a circular orbit of a moving object,
F - vector of gravity,
M - mass of the source of gravitational field,
m - mass of object moving in the gravitational field,
r - position vector from the source of gravitational field to a moving object,
r'=v - velocity vector of a moving object,
r"=v'=a - acceleration vector of a moving object,
T - period of circulation,
ω=2π/T - scalar value of angular velocity,
Here bold letters signify vectors, regular letters signify scalars and magnitudes of corresponding vectors, single and double apostrophes signify first and second derivative by time.

Constant magnitude v of velocity vector means constant angular velocity ω and obvious equality v=r·ω.

Magnitude a of an acceleration vector can be simply found by representing a position vector as a pair of Cartesian coordinates (x,y):
x = r·cos(ωt)
y = r·sin(ωt)
x' = −r·ω·sin(ωt)
y' = r·ω·cos(ωt)
x" = −r·ω²·cos(ωt)=−ω²·x
y" = −r·ω²·sin(ωt)=−ω²·y
and, therefore,
a = r" = −ω²·r
(collinear with r and F)
from which follows
a = |a| = |−ω²·r| = ω²·r

According to the Newton's Second Law,
F = m·a
According to the Universal Law of Gravitation,
F = G·M·m/r²
Therefore,
a = ω²·r = G·M/r²
from which follows
ω² = G·M/r³

Since ω=2π/T,
4π²/T² = G·M/r³
T²/r³ = 4π²/(G·M) - constant
End of proof for circular orbit.

Let's prove it in a more complicated general case of any elliptical orbit.
We will use the First and the Second Kepler's Laws as well as the results presented in the previous lecture Planet Orbit Geometry to derive this Third Law.

Recall the Kepler's Second Law (see the lecture Kepler's Second Law in this course).
We have introduced a function A(t) that represents an area of a sector bounded by r(0), r(t) and a trajectory from a planet's position P(0) at time t=0 to its position at any moment of time P(t).
Then the area swept by position vector r(t) during the object's motion from time t₁ to t₂ equals
ΔA_[t1,t2] = A(t₂) − A(t₁)

Using the above symbols, the Kepler's Second Law can be formulated as
If t₂−t₁ = t₄−t₃ then
A(t₂)−A(t₁) = A(t₄)−A(t₃)
The above condition is equivalent to a statement that
dA(t)/dt is constant or, equivalently, that A(t) is a linear function of time t with A(0)=0.

In the same lecture we have proven that
dA(t)/dt = ½|L|/m
where L is an angular momentum of a moving object (constant in a central force field) and m is object's mass.

Therefore,

A(t) = t·½|L|/m

Assume, we want to know how much area is swept by position vector r(t) during a period T of a complete round movement of a planet around the Sun.
Obviously, it's
A(T) = T·½|L|/m
At the same time, A(T) is an area of an elliptical orbit of a planet, and we know that the area of an ellipse along which a planet moves equals to
A(T) = π·a·b
where a is semi-major and b is semi-minor axes (see the lecture More on Ellipse in this course).

Therefore,
π·a·b = T·½|L|/m

Let's assume that at time t=0 a planet is at the furthest from the Sun point (aphelion), it's initial position vector is r₀ and its velocity vector is v₀.
At this initial point on an orbit the position vector r₀, lying along the major axis of an ellipse, and a tangential to an ellipse vector of velocity v₀ are perpendicular to each other.

Therefore, the magnitude of the angular momentum |L| equals to a product of planet's mass, a magnitude of its position vector (that is, a distance of the Sun) and a magnitude of its velocity vector:
|L| = m·|r₀|·|v₀| or L/m = r₀·v₀

Now the above formula for a period T of a planet's rotation around the Sun is
π·a·b = ½T·r₀·v₀

In the previous lecture Planet Orbit Geometry we have derived the expressions for major and minor axes of an elliptical orbit of a planet in terms of its initial position and velocity at aphelion:

a =

r0 2−β

b =	r0√β·(2−β) 2−β

where β=r₀·v₀²/(G·M)

Let's substitute these expression into a formula connecting a period T with an area of an ellipse:
πr₀²√β·(2−β)/(2−β)²=½T·r₀·v₀

Let's square both sides to get rid of a radical:
π²r₀⁴β/(2−β)³=¼T²·r₀²·v₀²

Next is just technicality.
Cancel one r₀ from both sides
π²r₀³β/(2−β)³=¼T²·r₀·v₀²

Replace r₀³/(2−β)³ with a³ (see formula above)
π²a³β=¼T²·r₀·v₀²

Replace β with r₀·v₀²/(G·M) (see formula above)
π²a³r₀·v₀²/(G·M)=¼T²·r₀·v₀²

Cancel r₀·v₀² on both sides
π²a³/(G·M)=¼T²

Final result:

T² a3

=

4π² G·M

The right side is a constant that contains no planet-specific parameters (like initial position and velocity), which means that any planet has the following property (Kepler's Thirt Law).
The ratio of a square of the period of a planet's rotation around the Sun to a cube of a semi-major axis of its elliptical orbit is constant that depends only on a mass of the Sun.