EduKG
← Back

gravitation Kepler laws and cosmic velocities

Learning objectives

  • Outcome

    Compute the attraction force between two point masses with F=Gm1m2r2F = G\frac{m_1 m_2}{r^{2}} and state its direction.

  • Outcome

    Find the circular-orbit speed and the first cosmic velocity from v=GM/r=gRv = \sqrt{GM/r} = \sqrt{gR}.

  • Outcome

    Apply Kepler's third law T2/a3=constT^{2}/a^{3} = \mathrm{const} to calculate an orbital period or radius.

  • Outcome

    Derive the second cosmic velocity v2=2GM/Rv_2 = \sqrt{2GM/R} from the zero-total-energy condition.

Understand the idea

Gravitation, Kepler's Laws and Cosmic Velocities

Text

Newton showed that two point masses attract each other along the line joining them with the force F=G m1m2r2,G≈6.67×10−11 m3 kg−1 s−2.F = G\,\frac{m_1 m_2}{r^{2}},\qquad G \approx 6.67\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}}.

For a body on a circular orbit of radius rr this force plays the role of the centripetal force, so GMm/r2=mv2/rG M m / r^{2} = m v^{2}/r, hence

v=GMr.v = \sqrt{\frac{GM}{r}}.

Newton's law of gravitation applied to general (not necessarily circular) two-body motion together with Newton's second law yields Kepler's three laws: a planet moves along an ellipse with the Sun at one focus; the radius vector sweeps equal areas in equal times; the squares of the periods are proportional to the cubes of the semi-major axes, T2/a3=constT^{2}/a^{3} = \mathrm{const}. The circular-orbit relation above is a special case: it reproduces the third law only for circular motion (T2∝r3T^{2} \propto r^{3}), while the first and second laws require non-circular motion.

Kepler's second law is conservation of angular momentum: at perihelion and aphelion v r=constv\,r = \mathrm{const}, so a comet moves faster when it is close to the Sun.

Equating gravitation to the "weight" at a planet's surface gives the first cosmic velocity: v1=gR=GM/R≈7.9 km/sv_1 = \sqrt{gR} = \sqrt{GM/R} \approx 7.9\ \mathrm{km/s} for the Earth.

The second cosmic velocity follows from the zero-total-energy condition 12v2−GM/R=0\frac{1}{2}v^{2} - GM/R = 0, hence v2=2GM/R=2 v1≈11.2 km/sv_2 = \sqrt{2GM/R} = \sqrt{2}\,v_1 \approx 11.2\ \mathrm{km/s}.

Worked examples

Direct

Problem

A satellite of mass m=500m = 500 kg moves in a circular orbit of radius r=1.28×107r = 1.28\times10^{7} m around the Earth. Taking GM=3.986×1014 m3/s2GM = 3.986\times10^{14}\ \mathrm{m^3/s^2}, find the gravitational force acting on the satellite and the gravitational field strength (free-fall acceleration) at that height.

Solution

Given: m=500m = 500 kg, r=1.28×107r = 1.28\times10^{7} m, GM=3.986×1014 m3/s2GM = 3.986\times10^{14}\ \mathrm{m^3/s^2}.

  1. Law of universal gravitation: F=GMmr2.F = \frac{GMm}{r^{2}}.
  2. Square of the radius: (1.28×107)2=1.638×1014 m2(1.28\times10^{7})^{2} = 1.638\times10^{14}\ \mathrm{m^2}.
  3. Numerator: GMm=3.986×1014⋅500=1.993×1017GMm = 3.986\times10^{14}\cdot 500 = 1.993\times10^{17}.
  4. Force: F=1.993×10171.638×1014≈1.2×103F = \dfrac{1.993\times10^{17}}{1.638\times10^{14}} \approx 1.2\times10^{3} N ≈1.2\approx 1.2 kN.
  5. Field strength: g=GMr2=3.986×10141.638×1014≈2.4 m/s2g = \dfrac{GM}{r^{2}} = \dfrac{3.986\times10^{14}}{1.638\times10^{14}} \approx 2.4\ \mathrm{m/s^2} — about four times smaller than at the surface, because r=2Rr = 2R.

Application

Problem

Find the radius of the geostationary orbit, on which a satellite circles the Earth with period T=8.64×104T = 8.64\times10^{4} s (24 h), and its speed. Use Kepler's third law and GM=3.986×1014 m3/s2GM = 3.986\times10^{14}\ \mathrm{m^3/s^2}.

Solution

Given: T=8.64×104T = 8.64\times10^{4} s, GM=3.986×1014 m3/s2GM = 3.986\times10^{14}\ \mathrm{m^3/s^2}.

  1. Kepler's third law in the form T2a3=4π2GM\dfrac{T^{2}}{a^{3}} = \dfrac{4\pi^{2}}{GM} gives a=(GM T24π2)1/3.a = \left(\frac{GM\,T^{2}}{4\pi^{2}}\right)^{1/3}.
  2. T2=(8.64×104)2=7.47×109 s2T^{2} = (8.64\times10^{4})^{2} = 7.47\times10^{9}\ \mathrm{s^2}, so GM T2=3.986×1014⋅7.47×109=2.98×1024 m3GM\,T^{2} = 3.986\times10^{14}\cdot 7.47\times10^{9} = 2.98\times10^{24}\ \mathrm{m^3}.
  3. 4π2=39.54\pi^{2} = 39.5, hence a=(2.98×102439.5)1/3=(7.54×1022)1/3≈4.22×107a = \left(\dfrac{2.98\times10^{24}}{39.5}\right)^{1/3} = (7.54\times10^{22})^{1/3} \approx 4.22\times10^{7} m, that is about 42 200 km from the Earth's centre.
  4. Speed: v=GMa=3.986×10144.22×107≈3.07 km/sv = \sqrt{\dfrac{GM}{a}} = \sqrt{\dfrac{3.986\times10^{14}}{4.22\times10^{7}}} \approx 3.07\ \mathrm{km/s}; a check through the orbit length gives the same: v=2πaT≈3.07 km/sv = \dfrac{2\pi a}{T} \approx 3.07\ \mathrm{km/s}.
  5. Height above the surface: 42 200−6 400≈35 80042\,200 - 6\,400 \approx 35\,800 km.

Transfer

Problem

A comet moves around the Sun along an ellipse with perihelion distance rp=1.0×1011r_p = 1.0\times10^{11} m and aphelion distance ra=5.0×1011r_a = 5.0\times10^{11} m. Using Kepler's second law and conservation of energy, find the comet's speed at perihelion and at aphelion. Take GM⊙=1.33×1020 m3/s2GM_{\odot} = 1.33\times10^{20}\ \mathrm{m^3/s^2}.

Solution

Given: rp=1.0×1011r_p = 1.0\times10^{11} m, ra=5.0×1011r_a = 5.0\times10^{11} m, GM=1.33×1020 m3/s2GM = 1.33\times10^{20}\ \mathrm{m^3/s^2}.

  1. Kepler's second law (conservation of angular momentum): vprp=varav_p r_p = v_a r_a, so va=vprpra=vp5v_a = v_p\dfrac{r_p}{r_a} = \dfrac{v_p}{5}.
  2. Conservation of energy per unit mass: 12vp2−GMrp=12va2−GMra.\frac{1}{2}v_p^{2} - \frac{GM}{r_p} = \frac{1}{2}v_a^{2} - \frac{GM}{r_a}.
  3. Substitute va=vp/5v_a = v_p/5 and ra=5rpr_a = 5r_p: 12vp2(1−125)=GMrp(1−15).\frac{1}{2}v_p^{2}\left(1 - \frac{1}{25}\right) = \frac{GM}{r_p}\left(1 - \frac{1}{5}\right).
  4. Hence vp2=53 GMrpv_p^{2} = \dfrac{5}{3}\,\dfrac{GM}{r_p}. With GMrp=1.33×10201.0×1011=1.33×109 m2/s2\dfrac{GM}{r_p} = \dfrac{1.33\times10^{20}}{1.0\times10^{11}} = 1.33\times10^{9}\ \mathrm{m^2/s^2} we get vp2=2.22×109v_p^{2} = 2.22\times10^{9}, that is vp≈4.7×104 m/s≈47 km/sv_p \approx 4.7\times10^{4}\ \mathrm{m/s} \approx 47\ \mathrm{km/s}.
  5. va=vp/5≈9.4 km/sv_a = v_p/5 \approx 9.4\ \mathrm{km/s}.
  6. Check: the total energy 12vp2−GMrp≈−2.2×108 J/kg<0\frac{1}{2}v_p^{2} - \dfrac{GM}{r_p} \approx -2.2\times10^{8}\ \mathrm{J/kg} < 0, so the orbit is closed (the comet does not leave the Solar System).

Practice

Text

Self-study work in a simulation (the safe route: real orbital measurements are out of reach).

  1. Open the PhET simulation "Gravity and Orbits" or an equivalent one.
  2. Keep the central mass fixed and create at least five circular orbits with different radii rr; record rr and the period TT for each.
  3. For each orbit compute r3r^{3} and T2T^{2} and plot T2T^{2} against r3r^{3}.
  4. Check that the points lie on a straight line through the origin and find the slope k=T2/r3k = T^{2}/r^{3}; compare it with 4π2/(GM)4\pi^{2}/(GM).
  5. Halve the central mass, repeat one measurement and verify that the slope doubles — this is Kepler's third law.

Visual reference

Purpose

Show the geometric meaning of Kepler's laws and that the third law T2∝a3T^{2} \propto a^{3} is verified by a straight line on a graph.

Representation

A two-panel diagram.

  • Left — an elliptical orbit with the Sun at one focus; two sectors swept by the radius vector in equal time intervals are shaded and have equal areas (the second law).
  • Right — a graph: a3a^{3} on the horizontal axis, T2T^{2} on the vertical axis; the points for several planets lie on a straight line through the origin (the third law).

Real-world examples

The ISS: one orbit in 92 minutes

Scenario

The International Space Station flies about 400400 km above the Earth's surface. With the Earth's radius R≈6400R \approx 6400 km the orbital radius is r≈6800r \approx 6800 km =6,8⋅106= 6{,}8\cdot10^{6} m, and GM=3,99⋅1014GM = 3{,}99\cdot10^{14} m$^3$/s$^2$. Then the station's speed is

v=GMr=3,99⋅10146,8⋅106≈7,7 km/s,v=\sqrt{\frac{GM}{r}}=\sqrt{\frac{3{,}99\cdot10^{14}}{6{,}8\cdot10^{6}}}\approx 7{,}7\ \text{km/s},

and the period is T=2πr/v≈5550T = 2\pi r/v \approx 5550 s ≈92\approx 92 min. In one day the crew sees about 16 sunrises.

Connection

This is the first cosmic (circular) velocity at work: the gravitational force F=GMm/r2F=GMm/r^2 plays the role of the centripetal force, giving v=GM/rv=\sqrt{GM/r}. It is a direct consequence of Newton's second law applied to circular motion.

A geostationary weather and TV satellite

Scenario

A satellite-TV dish points at a satellite that hangs in one fixed spot in the sky. Its period must equal the sidereal day: T=86 164T = 86\,164 s. From Kepler's third law

a=GMT24π23=3,99⋅1014⋅(86 164)24π23≈4,22⋅107 m,a=\sqrt[3]{\frac{GMT^{2}}{4\pi^{2}}}=\sqrt[3]{\frac{3{,}99\cdot10^{14}\cdot(86\,164)^{2}}{4\pi^{2}}}\approx 4{,}22\cdot10^{7}\ \text{m},

that is 42 20042\,200 km from the Earth's centre and an altitude of ≈35 800\approx 35\,800 km. The satellite's speed is v=2πa/T≈3,1v = 2\pi a/T \approx 3{,}1 km/s.

Connection

The example tests Kepler's third law T2=4π2a3/(GM)T^{2} = 4\pi^{2}a^{3}/(GM): only at one radius does the orbital period equal one day, so the geostationary orbit is unique. The factor GMGM is the same gravitational constant as in the law of universal gravitation.

A rocket to Mars and the cosmic velocities

Scenario

A launch vehicle accelerates its stage to 7,97{,}9 km/s to reach low orbit, and to 11,211{,}2 km/s to leave the Earth for good. For a flight to Mars the engine on the low orbit (h=200h=200 km) is fired again, adding ≈3,6\approx 3{,}6 km/s. The trip follows half of a Hohmann ellipse: the Earth at 1,001{,}00 AU, Mars at 1,521{,}52 AU, so the semi-major axis is a=(1,00+1,52)/2=1,26a = (1{,}00+1{,}52)/2 = 1{,}26 AU. By Kepler's third law T=a3/2=1,42T = a^{3/2} = 1{,}42 years, so the transfer takes T/2≈0,71T/2 \approx 0{,}71 years ≈259\approx 259 days.

Connection

The second cosmic velocity v2=2GM/R=11,2v_{2}=\sqrt{2GM/R}=11{,}2 km/s follows from energy conservation: the kinetic energy exactly cancels the gravitational potential energy −GMm/R-GMm/R. The shape of the trajectory and the flight time are set by Kepler's first and third laws; near perihelion (at the Earth) the craft moves faster than near aphelion — that is Kepler's second law.