iLabs Forces

Forces / Gravitation

Gravitation

Gravitation lab

Planets on ellipses that sweep out equal areas, Newton's cannon and the escape speed, how gravity changes above and inside a planet, and two stars circling their centre of mass. Change a number and watch the orbit change; the worked explanation below does the same sums with your values. (Free fall and projectiles are in the Mechanics lab.)

t = 0.00 s

Graph

Graph

The physics behind it
Formulas used in this lab

Law of gravitation Class 9 & 11

Newton's law
F = Gm1m2/r2, with G = 6.674 × 10−11 N·m2/kg2
Gravitational field
g = GM/r2 (outside); at a surface g0 = GM/R2
Weight
W = mg (mass does not change from place to place, weight does)
Potential energy
U = −GMm/r; gravitational potential V = −GM/r (zero at infinity)
Near the surface
ΔU = mgh for small h

Variation of g Class 11

At height h
gh = g0R2/(R + h)2 ≈ g0(1 − 2h/R) for h << R
At depth d (uniform density)
gd = g0(1 − d/R) = g0r/R; zero at the centre
Shell theorem
a spherical shell has no gravitational effect inside it; outside it acts as if its mass were at the centre
Latitude and rotation
geff = g − ω2R cos2λ: least at the equator, equal to g at the poles
Potential inside
V = −GM(3R2 − r2)/2R3; at the centre 1.5 times its surface value

Satellites and escape Class 11

Orbital speed
v = √(GM/r) = √(g0R2/r)
Period
T = 2πr/v = 2π√(r3/GM)
Energies of a circular orbit
KE = GMm/2r, PE = −GMm/r, total E = −GMm/2r = −KE
Escape speed
ve = √(2GM/r) = √2 vc = √(2gR) = 11.2 km/s from the Earth's surface
Geostationary orbit
T = 24 h, r = 42 164 km (h = 35 786 km), v = 3.07 km/s, above the equator
Orbit shape from energy
E < 0: ellipse or circle; E = 0: parabola; E > 0: hyperbola

Kepler's laws Class 11

First law
orbits are ellipses with the Sun at one focus; rmin = a(1 − e), rmax = a(1 + e)
Second law
equal areas in equal times: dA/dt = L/2m = constant (angular momentum is conserved); vprp = vara
Third law
T2 = (4π2/GM)a3; T2 = a3 for T in years and a in AU
Energy of an ellipse
E = −GMm/2a, with the speed from v2 = GM(2/r − 1/a)
Eccentricity
e = (rmax − rmin)/(rmax + rmin); from the energy and angular momentum per unit mass, e = √(1 + 2EL2/(GM)2)

Two-body motion Class 11

Centre of mass
m1r1 = m2r2, r1 + r2 = d
Period
T2 = 4π2d3/G(m1 + m2)
Reduced mass
μ = m1m2/(m1 + m2)
Speeds
v1/v2 = m2/m1 = r1/r2