Mechanical Engineer's Handbook

Satellite Orbits and Energy

For simplicity, this chapter only covers orbits along a circular path.

Deriving Orbital Velocity

Consider a satellite of mass mm orbiting Earth (mass MEM_E) in a circular path at a distance rr from Earth's center. The only force acting on the satellite is gravity, which acts as the required centripetal force holding the satellite in orbit:

GmMEr2=mac\frac{GmM_E}{r^2} = ma_c

Recall that centripetal acceleration is defined as ac=v2ra_c = \frac{v^2}{r}. Substituting this expression yields:

GmMEr2=mv2r\frac{GmM_E}{r^2} = m \frac{v^2}{r}

Canceling the satellite's mass mm and multiplying both sides by rr isolates the squared velocity:

v2=GMErv^2 = \frac{GM_E}{r}

Taking the square root gives the orbital speed:

vorbit=GMErv_{\text{orbit}} = \sqrt{\frac{GM_E}{r}}
(1)

Applying Orbital Velocity

Using the expression for orbital velocity, we can determine the kinetic energy (KK) of a satellite in a circular orbit:

K=12mv2=12m(GMEr)=GMEm2rK = \frac{1}{2}m v^2 = \frac{1}{2} m \left(\frac{GM_E}{r}\right) = \frac{GM_E m}{2r}
(2)

We can also apply orbital velocity to find the period of circular orbit. If we say the length of the circular orbit is 2πr2 \pi r and lasts TT seconds, we can derive TT with the following:

vorbit=2πrTv_{\text{orbit}} = \frac{2 \pi r}{T}
GMEr=2πrT\sqrt{\frac{GM_E}{r}} = \frac{2 \pi r}{T}

Inverting both sides yields:

rGME=T2πr\sqrt{\frac{r}{GM_E}} = \frac{T}{2 \pi r}
T=2πrrGMET = 2 \pi r \sqrt{\frac{r}{GM_E}}

We can further simplify by bringing rr inside of the radical:

T=2πr2rGMET = 2 \pi \sqrt{r^2} \sqrt{\frac{r}{GM_E}}
T=2πr3GMET = 2 \pi \sqrt{\frac{r^3}{GM_E}}
(3)

Total Mechanical Energy in Circular Orbits

Now that we have formulas for kinetic and potential energy for objects in orbit we can find the total mechanical energy:

Recall that the total mechanical energy for an object is the sum of its kinetic and potential energies:

Emech=K+UE_{\text{mech}} = K + U

We can apply known KK from Equation 2 with our known UU from the prior chapter to get EmechE_{\text{mech}}:

Emech=GMEm2rGMEmrE_{\text{mech}} = \frac{GM_Em}{2r} - \frac{GM_Em}{r}

Simplifying, we get the following:

Emech=GMEm2r2GMEm2rE_{\text{mech}} = \frac{GM_Em}{2r} - \frac{2GM_Em}{2r}
Emech=GMEm2GMEm2rE_{\text{mech}} = \frac{GM_Em - 2GM_Em}{2r}
Emech=GMEm2rE_{\text{mech}} = -\frac{GM_Em}{2r}
(4)

Notice that for a circular orbit, Emech=K=12UE_{\text{mech}} = -K = \frac{1}{2}U. Because total energy is negative (Emech<0E_{\text{mech}} < 0), the satellite is in a bound orbit.

Energy Required to Reach Orbit

Launching an object into orbit requires supplying energy for two separate things:

  1. Gravitational Potential Energy: Raising the object from Earth's surface (r=REr = R_E) to orbital altitude (r=RE+hr = R_E + h).
  2. Kinetic Energy: Accelerating the object up to orbital speed (vorbitv_{\text{orbit}}).

Because the object is at rest on the surface before launch, its initial energy is purely potential:

Einitial=GMEmREE_{\text{initial}} = -\frac{GM_E m}{R_E}

Once in orbit, its total energy is:

Efinal=GMEm2rE_{\text{final}} = -\frac{GM_E m}{2r}

The minimum work required to launch the satellite into orbit is:

W=EfinalEinitial=GMEm2r+GMEmRE=GMEm(1RE12r)W = E_{\text{final}} - E_{\text{initial}} = -\frac{GM_E m}{2r} + \frac{GM_E m}{R_E} = GM_E m \left( \frac{1}{R_E} - \frac{1}{2r} \right)
(5)

Example 1

The International Space Station (ISS) orbits at an altitude of approximately h=400 kmh = 400 \text{ km} above Earth's surface. Given ME=5.97×1024 kgM_E = 5.97 \times 10^{24} \text{ kg} and RE=6.37×106 mR_E = 6.37 \times 10^6 \text{ m}, calculate the orbital speed (vorbitv_{\text{orbit}}) and the orbital period (TT) of the ISS.

The radius of the orbit rr is measured from the center of Earth: r=RE+hr = R_E + h. Use vorbit=GMErv_{\text{orbit}} = \sqrt{\frac{GM_E}{r}} and T=2πr3GMET = 2\pi \sqrt{\frac{r^3}{GM_E}}.

First, convert altitude to meters and calculate the total orbital radius from Earth's center:

r=RE+h=6.37×106 m+0.40×106 m=6.77×106 mr = R_E + h = 6.37 \times 10^6\text{ m} + 0.40 \times 10^6\text{ m} = 6.77 \times 10^6\text{ m}

To find orbital velocity, substitute rr and Earth's mass into Equation 1:

vorbit=GMEr=(6.674×1011Nm2/kg2)(5.97×1024 kg)6.77×106 m7.67×103 m/sv_{\text{orbit}} = \sqrt{\frac{GM_E}{r}} = \sqrt{\frac{(6.674 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2)(5.97 \times 10^{24}\text{ kg})}{6.77 \times 10^6\text{ m}}} \approx 7.67 \times 10^3\text{ m/s}

To find the period, substitute rr into Equation 3:

T=2πr3GME=2π(6.77×106 m)3(6.674×1011Nm2/kg2)(5.97×1024 kg)5550 sT = 2\pi \sqrt{\frac{r^3}{GM_E}} = 2\pi \sqrt{\frac{(6.77 \times 10^6\text{ m})^3}{(6.674 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2)(5.97 \times 10^{24}\text{ kg})}} \approx 5550\text{ s}

Converting the period into minutes:

T=5550 s60 s/min92.5 minutesT = \frac{5550\text{ s}}{60\text{ s/min}} \approx 92.5\text{ minutes}

Example 2

A 1000 kg1000\text{ kg} satellite is launched from rest on Earth's surface into a circular orbit at an altitude of h=400 kmh = 400\text{ km} (r=6.77×106 mr = 6.77 \times 10^6\text{ m}). How much work must be done to place this payload into orbit?

The required work equals the total change in mechanical energy: W=EfinalEinitial=GMEm2r(GMEmRE)W = E_{\text{final}} - E_{\text{initial}} = -\frac{GM_E m}{2r} - \left(-\frac{GM_E m}{R_E}\right).

Using the work-energy relation for orbital insertion (Equation 5):

W=GMEm(1RE12r)W = GM_E m \left( \frac{1}{R_E} - \frac{1}{2r} \right)

First, calculate the common product GMEmGM_E m:

GMEm=(6.674×1011Nm2/kg2)(5.97×1024 kg)(1000 kg)3.984×1017 Nm2GM_E m = (6.674 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2)(5.97 \times 10^{24}\text{ kg})(1000\text{ kg}) \approx 3.984 \times 10^{17}\text{ N}\cdot\text{m}^2

Now substitute RE=6.37×106 mR_E = 6.37 \times 10^6\text{ m} and r=6.77×106 mr = 6.77 \times 10^6\text{ m}:

W=(3.984×1017 Nm2)(16.37×106 m12(6.77×106 m))W = (3.984 \times 10^{17}\text{ N}\cdot\text{m}^2) \left( \frac{1}{6.37 \times 10^6\text{ m}} - \frac{1}{2(6.77 \times 10^6\text{ m})} \right)
W=(3.984×1017 Nm2)(1.570×107 m10.739×107 m1)W = (3.984 \times 10^{17}\text{ N}\cdot\text{m}^2) \left( 1.570 \times 10^{-7}\text{ m}^{-1} - 0.739 \times 10^{-7}\text{ m}^{-1} \right)
W3.31×1010 JW \approx 3.31 \times 10^{10}\text{ J}