Mechanical Engineer's Handbook

Gravitational Energy

In near-Earth applications, potential energy induced by gravity can be calculated using U=mghU = mgh. This, however, gets less accurate the further you get from Earth as gravity isn't constant.

Deriving ΔU\Delta U and U(r)U(r)

Recall that work is defined as the integral of force with respect to displacement:

WAB=ABFdrW_{AB} = \int_{A}^{B} \vec{F} \cdot d \vec{r}

Also recall that change in potential energy is equal to the negative of work done by a conservative force:

ΔU=WAB=ABFdr\Delta U = -W_{AB} = -\int_{A}^{B} \vec{F} \cdot d \vec{r}

Evaluating this integral from initial distance r1r_1 to final distance r2r_2, we get:

ΔU=GMEm(1r11r2)\Delta U = GM_{E}m \left( \frac{1}{r_1} - \frac{1}{r_2} \right)
(1)

This formula for ΔU\Delta U provides the change in potential energy of an object with mass mm moving from distance r1r_1 to distance r2r_2 relative to Earth's center.

To define absolute potential energy UU at a single position rr, rather than a difference between two points, we must establish a reference point where potential energy is zero. We can define U=0U = 0 at an infinite distance away (rr \to \infty).

Expanding ΔU=UfUi\Delta U = U_f - U_i where the initial position is r1=rr_1 = r and the final position is sent out to infinity (r2=r_2 = \infty, so Uf=0U_f = 0):

UfUi=GMEm(1r11r2)U_f - U_i = GM_{E}m \left( \frac{1}{r_1} - \frac{1}{r_2} \right)
0U=GMEm(1r1)0 - U = GM_{E}m \left( \frac{1}{r} - \frac{1}{\infty} \right)
U=GMEm(1r0)-U = GM_{E}m \left( \frac{1}{r} - 0 \right)
U=GMEmrU = -\frac{GM_{E}m}{r}
(2)

Notice that in Equation 2, potential energy is negative. This is a direct consequence of choosing U=0U = 0 at rr \to \infty.

Because gravity is an attractive force, it does positive work as an object falls inward toward Earth. Since a loss in potential energy corresponds to positive work (ΔU=W\Delta U = -W), moving closer to Earth decreases potential energy further below zero.

An object at a finite distance rr is bound by Earth's gravitational field. It would require an input of positive energy equal to GMEmr\frac{GM_E m}{r} to pull it completely free to infinity.

Example 1

A 3000 kg\text{3000 kg} mass is boosted from the surface of Earth to a height of 500 km\text{500 km}. How much energy was expended?

We must first recognize this problem is asking for ΔU\Delta U, not an absolute UU. This is because the mass starts on the surface, therefore Ui0U_i \neq 0. We say UiU_i is the potential energy of the mass on the surface of Earth, and UfU_f is the potential energy at a height of 500 km\text{500 km}.

ΔU=GMEm(1rE1rE+r)\Delta U = GM_{E}m \left( \frac{1}{r_E} - \frac{1}{r_E + r} \right)

Plugging in our known values for GG, MEM_E, mm, rEr_E, and rr:

ΔU=(6.671011Nm2kg2)(5.971024kg)(3000kg)(16.37106m16.37106m+5105m)\Delta U = \left( 6.67 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2} \right) (5.97 \cdot 10^{24} kg) (3000kg) \left( \frac{1}{6.37 \cdot 10^6 m} - \frac{1}{6.37 \cdot 10^6 m + 5 \cdot 10^5m} \right)
ΔU1.361010J\Delta U \approx 1.36 \cdot 10^{10} J

Conservation of Energy

Because gravity is a conservative force, we can apply the conservation of energy laws to gravity at extended distances.

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

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

Also recall that for conservation of energy, the initial total mechanical energy must equal the final total mechanical energy:

Ki+Ui=Kf+UfK_i + U_i = K_f + U_f

Substituting in the kinetic energy for a mass in motion, as well as our definition of UU from Equation 2:

12mvi2GMmri=12mvf2GMmrf\frac{1}{2} m v_i^2 - \frac{GMm}{r_i} = \frac{1}{2}mv_f^2 - \frac{GMm}{r_f}
(3)

Note that MEM_E from Equation 2 was replaced with a generic MM, as the equation doesn't strictly apply to Earth. It is still assumed that m<Mm<M.

Bound vs. Unbound Systems

An object's trajectory is determined by its total mechanical energy E=K+UE = K + U:

  • If E<0E < 0, the object does not have enough kinetic energy to break away from the gravitational potential energy. It will remain in a closed orbit.
  • If E=0E = 0, the object has just enough kinetic energy to reach rr \to \infty where it comes to a complete stop (vf=0v_f = 0). This represents the threshold for escape.
  • If E>0E > 0, the object escapes to infinity and still retains leftover kinetic energy (vf>0v_f > 0), following a hyperbolic trajectory.

Escape Velocity

The escape velocity (vescv_{esc}) is the minimum speed an object must be given at initial radius rr to escape a central mass MM completely without further propulsion.

This occurs precisely at the boundary condition where total mechanical energy is zero (E=0E = 0).

Deriving Escape Velocity

Setting total energy E=0E = 0 at distance rr:

K+U=0K + U = 0
12mvesc2GMmr=0\frac{1}{2}m v_{esc}^2 - \frac{GMm}{r} = 0

Solving for vescv_{esc}:

12mvesc2=GMmr\frac{1}{2}m v_{esc}^2 = \frac{GMm}{r}
vesc2=2GMrv_{esc}^2 = \frac{2GM}{r}
vesc=2GMrv_{esc} = \sqrt{\frac{2GM}{r}}
(3)

Example 2

Calculate the escape velocity from the surface of the Moon given MM=7.351022 kgM_M = 7.35 \cdot 10^{22} \text{ kg} and rM=1.74106 mr_M = 1.74 \cdot 10^{6} \text{ m}.

Using Equation 3:

vesc=2GMMrM=2(6.671011Nm2kg2)(7.351022 kg)1.74106 m2.38103 m/sv_{esc} = \sqrt{\frac{2GM_M}{r_M}} = \sqrt{\frac{2 \left(6.67 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2}\right) \left(7.35 \cdot 10^{22} \text{ kg}\right)}{1.74 \cdot 10^{6}\text{ m}}} \approx 2.38 \cdot 10^3 \text{ m/s}

Example 3

A rocket is launched vertically from Earth's surface with an initial speed vi=8000 m/sv_i = 8000 \text{ m/s}. Ignoring air resistance, what maximum height above Earth's surface does it reach?

Use conservation of energy. Remember that when a rocket reaches its maximum height (apex), vf=0v_f = 0.

Start with energy conservation:

12mvi2GMEmri=12mvf2GMEmrf\frac{1}{2} m v_i^2 - \frac{GM_E m}{r_i} = \frac{1}{2}mv_f^2 - \frac{GM_E m}{r_f}

Cancel out mass mm:

12vi2GMEri=12vf2GMErf\frac{1}{2} v_i^2 - \frac{GM_E}{r_i} = \frac{1}{2}v_f^2 - \frac{GM_E}{r_f}

Since vf=0v_f = 0 at the apex, omit the final kinetic energy term:

12vi2GMEri=GMErf\frac{1}{2} v_i^2 - \frac{GM_E}{r_i} = - \frac{GM_E}{r_f}

Rearrange to isolate 1rf\frac{1}{r_f}:

GMEri12vi2=GMErf\frac{GM_E}{r_i} - \frac{1}{2} v_i^2 = \frac{GM_E}{r_f}
1rivi22GME=1rf\frac{1}{r_i} - \frac{v_i^2}{2GM_E} = \frac{1}{r_f}

Combine terms on the left:

2GMErivi22riGME=1rf\frac{2GM_E - r_i v_i^2}{2 r_i GM_E} = \frac{1}{r_f}

Invert to solve for rfr_f:

rf=2riGME2GMErivi2r_f = \frac{2 r_i GM_E}{2 GM_E - r_i v_i^2}

Plugging in known values (ri=6.371106 mr_i = 6.371 \cdot 10^6 \text{ m}, G=6.671011Nm2kg2G = 6.67 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2}, ME=5.9721024 kgM_E = 5.972 \cdot 10^{24} \text{ kg}, vi=8000 m/sv_i = 8000 \text{ m/s}):

rf=2(6.371106 m)(6.671011Nm2kg2)(5.9721024 kg)2(6.671011Nm2kg2)(5.9721024 kg)(6.371106 m)(8000 m/s)2r_f = \frac{2 (6.371 \cdot 10^6 \text{ m})(6.67 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2})(5.972 \cdot 10^{24} \text{ kg})}{2 (6.67 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2})(5.972 \cdot 10^{24} \text{ kg}) - (6.371 \cdot 10^6 \text{ m})(8000 \text{ m/s})^2}
rf1.304107 mr_f \approx 1.304 \cdot 10^7 \text{ m}

Finally, calculate the height hh above Earth's surface:

h=rfri=1.304107 m6.371106 m6.67106 mh = r_f - r_i = 1.304 \cdot 10^7 \text{ m} - 6.371 \cdot 10^6 \text{ m} \approx 6.67 \cdot 10^6 \text{ m}

Example 4

A 1000 kg1000 \text{ kg} asteroid far outside the solar system approaches Earth from deep space with an initial speed of v=3000 m/sv_\infty = 3000 \text{ m/s}. How fast will it be moving right before it impacts Earth's surface (rREr \approx R_E)?

Because the asteroid originates at "deep space" (rr \to \infty), its initial potential energy Ui=0U_i = 0. Since total energy E>0E > 0, it follows an unbound hyperbolic trajectory.

Apply conservation of mechanical energy between deep space (rir_i \to \infty) and Earth's surface (rf=REr_f = R_E):

Ki+Ui=Kf+UfK_i + U_i = K_f + U_f
12mv20=12mvf2GMEmRE\frac{1}{2}m v_\infty^2 - 0 = \frac{1}{2}m v_f^2 - \frac{GM_E m}{R_E}

Cancel out mass mm from all terms:

12v2=12vf2GMERE\frac{1}{2} v_\infty^2 = \frac{1}{2} v_f^2 - \frac{GM_E}{R_E}

Multiply by 2 and solve for final velocity vfv_f:

vf2=v2+2GMEREv_f^2 = v_\infty^2 + \frac{2GM_E}{R_E}
vf=v2+2GMEREv_f = \sqrt{v_\infty^2 + \frac{2GM_E}{R_E}}

Note that 2GMERE\frac{2GM_E}{R_E} is equal to vesc2v_{esc}^2, showing that impact velocity for an unbound object is always vf=v2+vesc2v_f = \sqrt{v_\infty^2 + v_{esc}^2}.

Plugging in known values (v=3000 m/sv_\infty = 3000 \text{ m/s}, G=6.6741011Nm2kg2G = 6.674 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2}, ME=5.9721024 kgM_E = 5.972 \cdot 10^{24} \text{ kg}, RE=6.371106 mR_E = 6.371 \cdot 10^6 \text{ m}):

vf=(3000 m/s)2+2(6.6741011Nm2kg2)(5.9721024 kg)6.371106 mv_f = \sqrt{(3000 \text{ m/s})^2 + \frac{2 (6.674 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2})(5.972 \cdot 10^{24} \text{ kg})}{6.371 \cdot 10^6 \text{ m}}}
vf=9106 m2/s2+1.251108 m2/s21.16104 m/sv_f = \sqrt{9 \cdot 10^6 \text{ m}^2/\text{s}^2 + 1.251 \cdot 10^8 \text{ m}^2/\text{s}^2} \approx 1.16 \cdot 10^4 \text{ m/s}

The asteroid impacts at approximately 11.6 km/s11.6 \text{ km/s} (or 11,580 m/s11,580 \text{ m/s}).

Example 5

A 1500 kg1500 \text{ kg} space probe is launched from Earth's surface on a parabolic trajectory with total mechanical energy E=0E = 0.

a.) What launch velocity viv_i is required?

b.) What will its speed be when it is very far from Earth (rr \to \infty)?

A total mechanical energy of E=0E = 0 represents the exact boundary between bound (E<0E < 0) and unbound (E>0E > 0) states. This corresponds directly to the escape velocity equation.

Set the total initial energy to zero:

E=Ki+Ui=0E = K_i + U_i = 0
12mvi2GMEmRE=0\frac{1}{2}m v_i^2 - \frac{GM_E m}{R_E} = 0

Cancel mass mm and solve for viv_i:

12vi2=GMERE\frac{1}{2} v_i^2 = \frac{GM_E}{R_E}
vi=2GMEREv_i = \sqrt{\frac{2GM_E}{R_E}}

Substitute known values (G=6.6741011Nm2kg2G = 6.674 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2}, ME=5.9721024 kgM_E = 5.972 \cdot 10^{24} \text{ kg}, RE=6.371106 mR_E = 6.371 \cdot 10^6 \text{ m}):

vi=2(6.6741011Nm2kg2)(5.9721024 kg)6.371106 mv_i = \sqrt{\frac{2(6.674 \cdot 10^{-11} \frac{\text{N} \cdot \text{m}^2}{\text{kg}^2})(5.972 \cdot 10^{24} \text{ kg})}{6.371 \cdot 10^6 \text{ m}}}
vi1.12104 m/s(11.2 km/s)v_i \approx 1.12 \cdot 10^4 \text{ m/s} \quad (11.2 \text{ km/s})

As rr \to \infty, potential energy approaches zero (U=0U_\infty = 0).

By energy conservation:

E=K+U=0E = K_\infty + U_\infty = 0
12mv2+0=0    v=0 m/s\frac{1}{2}m v_\infty^2 + 0 = 0 \implies v_\infty = 0 \text{ m/s}

The probe comes to a rest relative to Earth as it approaches infinite distance.