EduKG
← Back

Universal gravitation and orbital motion

Learning objectives

  • Outcome

    Calculate the gravitational force between two point masses using F=Gm1m2r2F = G\dfrac{m_1 m_2}{r^2} and state its direction.

  • Outcome

    Derive v=GM/rv = \sqrt{GM/r} and T=2πr3/(GM)T = 2\pi\sqrt{r^3/(GM)} for a circular orbit from Newton's second law and apply them to numerical data.

  • Outcome

    Determine orbital radius, altitude above the surface, or the mass of a central body from a known orbital period.

Understand the idea

Universal Gravitation and Orbital Motion

Text

Any two point masses attract each other along the line joining them with a force F=Gm1m2r2F = G\dfrac{m_1 m_2}{r^2}, where G=6.674×10−11 N⋅m2/kg2G = 6.674\times10^{-11}\ \text{N·m}^2/\text{kg}^2 and rr is the distance between the centres of the bodies.

If there are several bodies, the net force on each is found by vector addition of the individual forces, following the superposition principle, just as in applying Newton's second law.

For a satellite in a circular orbit, gravity is the only force and it also plays the role of the centripetal force: GMmr2=mv2r.G\dfrac{Mm}{r^2} = m\dfrac{v^2}{r}.

The satellite mass mm cancels, so the speed and the period depend only on the mass of the central body MM and the orbital radius rr measured from the centre, not from the surface: v=GMr,T=2πr3GM.v = \sqrt{\dfrac{GM}{r}}, \qquad T = 2\pi\sqrt{\dfrac{r^3}{GM}}.

This gives Kepler's third law: for all satellites of one central body, T2∝r3T^2 \propto r^3.

Worked examples

Direct

Problem

The Earth has mass M=5.97×1024M = 5.97\times10^{24} kg. Find the gravitational force between the Earth and a satellite of mass m=1000m = 1000 kg located at a distance r=7.0×106r = 7.0\times10^{6} m from the Earth's centre.

Solution

Substitute the data into the law of universal gravitation: F=GMmr2=6.674×10−11⋅5.97×1024⋅1000(7.0×106)2.F = G\dfrac{Mm}{r^2} = 6.674\times10^{-11}\cdot\dfrac{5.97\times10^{24}\cdot 1000}{(7.0\times10^{6})^2}.

Numerator: 6.674×10−11⋅5.97×1027≈3.98×1017 N⋅m26.674\times10^{-11}\cdot 5.97\times10^{27} \approx 3.98\times10^{17}\ \text{N·m}^2. Denominator: (7.0×106)2=4.9×1013 m2(7.0\times10^{6})^2 = 4.9\times10^{13}\ \text{m}^2.

Therefore F=3.98×10174.9×1013≈8.1×103 N.F = \dfrac{3.98\times10^{17}}{4.9\times10^{13}} \approx 8.1\times10^{3}\ \text{N}. The force is about 8.1 kN, directed from the satellite toward the Earth's centre, i.e. it is an attractive force.

Application

Problem

A satellite moves in a circular orbit of radius r=7.0×106r = 7.0\times10^{6} m around the Earth (M=5.97×1024M = 5.97\times10^{24} kg). Find the satellite's speed and its orbital period.

Solution

Set the gravitational force equal to the centripetal force: GMmr2=mv2rG\dfrac{Mm}{r^2} = m\dfrac{v^2}{r}.

The satellite mass cancels, so v=GMr=6.674×10−11⋅5.97×10247.0×106=5.69×107≈7.54×103 m/s.v = \sqrt{\dfrac{GM}{r}} = \sqrt{\dfrac{6.674\times10^{-11}\cdot 5.97\times10^{24}}{7.0\times10^{6}}} = \sqrt{5.69\times10^{7}} \approx 7.54\times10^{3}\ \text{m/s}.

The period equals the circumference divided by the speed: T=2πrv=2π⋅7.0×1067.54×103≈5.8×103 s≈97 min.T = \dfrac{2\pi r}{v} = \dfrac{2\pi\cdot 7.0\times10^{6}}{7.54\times10^{3}} \approx 5.8\times10^{3}\ \text{s} \approx 97\ \text{min}.

Check: a period of about 97 minutes agrees with known low Earth orbits.

Transfer

Problem

A geostationary satellite must have a period T=86164T = 86164 s (one sidereal day) and move in the equatorial plane. Determine its orbital radius and its altitude above the Earth's surface. Given: M=5.97×1024M = 5.97\times10^{24} kg, Earth radius RE=6.37×106R_E = 6.37\times10^{6} m.

Solution

From the period formula T=2πr3GMT = 2\pi\sqrt{\dfrac{r^3}{GM}}, solve for the radius: r=(GMT24π2)1/3.r = \left(\dfrac{GMT^2}{4\pi^2}\right)^{1/3}.

Compute the numerator: GMT2=6.674×10−11⋅5.97×1024⋅(86164)2≈2.96×1024GMT^2 = 6.674\times10^{-11}\cdot 5.97\times10^{24}\cdot (86164)^2 \approx 2.96\times10^{24}. Divide by 4π2≈39.484\pi^2 \approx 39.48: r3≈7.5×1022 m3,r≈4.22×107 m.r^3 \approx 7.5\times10^{22}\ \text{m}^3, \qquad r \approx 4.22\times10^{7}\ \text{m}.

Altitude above the surface: h=r−RE≈4.22×107−6.37×106≈3.58×107 m,h = r - R_E \approx 4.22\times10^{7} - 6.37\times10^{6} \approx 3.58\times10^{7}\ \text{m}, that is about 35,800 km. Note that the satellite mass does not enter the calculation: for a given period, the orbital radius is determined only by the Earth's mass.

Practice

Text

Self-check with a simulation. Real orbital launches are unavailable and hazardous for independent study, so use a model: the PhET simulation "Gravity and Orbits" in "Model" mode.

  1. Turn on the gravity force and velocity vectors. Confirm that the force points toward the central body and the velocity is tangent to the orbit.
  2. With the central body's mass fixed, record the circular orbital radius rr and the measured period TT for 4–5 different radii.
  3. Plot T2T^2 against r3r^3. The points should fall on a straight line through the origin.
  4. Estimate the slope of the line and compare it with the theoretical value 4π2/(GM)4\pi^2/(GM).
  5. Double the mass of the central body and verify that at the same radius the period decreases by a factor of 2\sqrt{2}.

Visual reference

Purpose

Show that gravity acts as the centripetal force in a circular orbit, and how the speed and period depend on the orbital radius.

Representation

Diagram and graph.

  • Diagram: a massive body MM at the centre; a satellite mm at distance rr from the centre. From the satellite toward the centre, a force vector F⃗\vec{F} of magnitude GMmr2G\dfrac{Mm}{r^2}; along the orbit, a velocity vector v⃗\vec{v} perpendicular to F⃗\vec{F}. Dashed lines show two orbits of different radii labelled r1<r2r_1 < r_2.
  • Graph: r3r^3 on the horizontal axis, T2T^2 on the vertical axis; plotted data points and a straight line through the origin with slope 4π2/(GM)4\pi^2/(GM). Caption: "T2∝r3T^2 \propto r^3 — Kepler's third law".

Real-world examples

The ISS: orbital speed and period

Scenario

The International Space Station flies about 400400 km above Earth's surface. Earth's radius is R=6,371×106R = 6{,}371\times10^{6} m, so the distance from Earth's centre to the station is r=6,77×106r = 6{,}77\times10^{6} m. Earth's mass is M=5,97×1024M = 5{,}97\times10^{24} kg and G=6,67×10−11 N⋅m2/kg2G = 6{,}67\times10^{-11}\ \text{N}\cdot\text{m}^2/\text{kg}^2. The calculation gives a speed v=GM/r≈7,7v=\sqrt{GM/r}\approx 7{,}7 km/s and a period T=2πr/v≈5,5×103T=2\pi r/v\approx 5{,}5\times10^{3} s, that is about 9292 minutes.

Connection

The central part of the theory is at work here: the only force on the station is Earth's pull, F=GMm/r2F=GMm/r^2, and it plays the role of the centripetal force F=mv2/rF=mv^2/r. Setting them equal gives v=GM/rv=\sqrt{GM/r}: the station's mass mm cancels, so at a given altitude every body moves with the same speed.

A GPS satellite in your phone

Scenario

Your phone fixes its position from signals of GPS satellites flying at an altitude of 20 20020\,200 km, that is at a distance r=2,66×107r = 2{,}66\times10^{7} m from Earth's centre. With M=5,97×1024M = 5{,}97\times10^{24} kg we get v=GM/r≈3,9v=\sqrt{GM/r}\approx 3{,}9 km/s, and the period T=2πr/v≈4,3×104T=2\pi r/v\approx 4{,}3\times10^{4} s ≈12\approx 12 h. The satellite makes two revolutions per day, so it is above the same point of Earth at the same time each day.

Connection

The example tests how orbital quantities depend on radius: comparing the GPS satellite with the ISS, we see that when rr grows by a factor of 3,93{,}9, the speed falls by about a factor of 22 (7,7/3,9≈27{,}7/3{,}9\approx 2) while the period grows. This is Kepler's third law T2∝r3T^2\propto r^3, which follows directly from F=GMm/r2F=GMm/r^2 and F=mv2/rF=mv^2/r.

The weight of a high jumper

Scenario

A high jumper of mass 7070 kg takes off at 5,05{,}0 m/s and rises 1,31{,}3 m. At Earth's surface g=GM/R2=6,67×10−11⋅5,97×1024/(6,371×106)2≈9,8g=GM/R^2=6{,}67\times10^{-11}\cdot5{,}97\times10^{24}/(6{,}371\times10^{6})^2\approx 9{,}8 m/s$^2$, so their weight is mg≈687mg\approx 687 N. During the jump the distance to Earth's centre changes by 1,31{,}3 m, that is by the fraction 1,3/6,371×106≈2×10−71{,}3/6{,}371\times10^{6}\approx 2\times10^{-7}; by the inverse-square law the weight changes by only 70⋅9,8⋅4×10−7≈3×10−470\cdot9{,}8\cdot4\times10^{-7}\approx 3\times10^{-4} N.

Connection

The example shows that "weight" is the same universal-gravitation force F=GMm/R2F=GMm/R^2, only applied at a planet's surface, and that the free-fall acceleration g=GM/R2g=GM/R^2 does not depend on the jumper's mass. The 1/R21/R^2 proportionality explains why the change in gg can be neglected at jump height but not at ISS altitude.