Unit 3, lesson 15

Exponential Growth, Decay, and Newton Cooling

How a rate proportional to the amount gives exponential growth and decay, what the half-life and the time constant measure, and how Newton's law of cooling moves a temperature toward a constant ambient temperature.

  • Video 15
  • Growth, Decay and Population Models

A radioactive sample loses the same fraction of its remaining atoms in every hour, and a cup of coffee in a quiet room loses roughly the same fraction of its excess over room temperature in every minute. How long until half the sample is gone, and when will the coffee be cool enough to drink? Both questions lead to a rate proportional to the quantity that is changing. This lesson solves that equation, reads off the half-life and the time constant, and applies the same idea to a temperature difference in Newton’s law of cooling.

Key terms

Proportional rate of change

Definition With tt as the independent variable and Q′=dQ/dtQ' = dQ/dt, a quantity QQ has a rate of change proportional to itself if

Q′=kQQ' = kQ

for a constant kk. The constant kk has units of inverse time, and wherever Q≠0Q \ne 0 it equals the relative rate Q′/QQ'/Q. Every solution is defined for all tt and has the form

Q(t)=Q0ekt,Q(t) = Q_0 e^{kt},

where Q0=Q(0)Q_0 = Q(0). No solution is missed, because (e−ktQ)′=e−kt(Q′−kQ)=0(e^{-kt}Q)' = e^{-kt}(Q' - kQ) = 0, so e−ktQ(t)e^{-kt}Q(t) is constant on any interval where QQ solves the equation.

In plain words The rate is a fixed multiple of the current amount. Twice as much material changes twice as fast, so the same fraction of the amount changes in every unit of time.

Example Separation gives the same formula. On an interval where Q≠0Q \ne 0, dividing by QQ gives (ln⁡∣Q∣)′=k(\ln\lvert Q \rvert)' = k, so ln⁡∣Q∣=kt+C\ln\lvert Q \rvert = kt + C and Q=±eCektQ = \pm e^{C}e^{kt}. The zero function, lost in the division, is the case Q0=0Q_0 = 0. For instance, the problem Q′=0.03 QQ' = 0.03\,Q with tt in years and Q(0)=500Q(0) = 500 has the solution Q(t)=500e0.03tQ(t) = 500e^{0.03t}, and Q(10)=500e0.3≈674.93Q(10) = 500e^{0.3} \approx 674.93.

Exponential growth and decay

Definition Let Q′=kQQ' = kQ with Q0>0Q_0 > 0. If k>0k > 0, the solution Q0ektQ_0e^{kt} increases without bound, which is exponential growth. If k<0k < 0, it decreases toward 00 as t→∞t \to \infty but stays positive at every finite time, which is exponential decay. If k=0k = 0, it is constant. In every case, for all tt and every time step Δt\Delta t,

Q(t+Δt)=ek Δt Q(t),Q(t + \Delta t) = e^{k\,\Delta t}\,Q(t),

so equal time steps multiply the amount by equal factors. On the phase line of Q′=kQQ' = kQ, the equilibrium Q=0Q = 0 is unstable for k>0k > 0 and asymptotically stable for k<0k < 0.

In plain words Exponential change multiplies by the same factor over each equal stretch of time, while linear change adds the same amount. Because the amount changes, the absolute rate kQkQ changes with it.

Example Take Q0=10Q_0 = 10 and k=ln⁡3k = \ln 3 per hour, so Q(t)=10e(ln⁡3)t=10⋅3tQ(t) = 10e^{(\ln 3)t} = 10 \cdot 3^{t} and the amount triples every hour. A linear model with the same increase in the first hour adds 2020 every hour.

Hour00112233
Exponential101030309090270270
Linear1010303050507070

The exponential rate Q′=(ln⁡3)QQ' = (\ln 3)Q is about 10.9910.99 per hour when Q=10Q = 10 and about 98.8898.88 per hour when Q=90Q = 90.

Half-life

Definition For decay write k=−λk = -\lambda with λ>0\lambda > 0, so Q′=−λQQ' = -\lambda Q and Q(t)=Q0e−λtQ(t) = Q_0e^{-\lambda t} with Q0>0Q_0 > 0. The half-life t1/2t_{1/2} is the time in which the amount falls to half its value. It solves e−λt1/2=12e^{-\lambda t_{1/2}} = \tfrac{1}{2}, so

t1/2=ln⁡2λ.t_{1/2} = \frac{\ln 2}{\lambda}.

The same duration halves the amount from any starting time, since Q(t+t1/2)=12Q(t)Q(t + t_{1/2}) = \tfrac{1}{2}Q(t) for every tt. After nn half-lives the amount is Q0/2nQ_0/2^n, and in general Q(t)=Q0 2−t/t1/2Q(t) = Q_0\,2^{-t/t_{1/2}}. For growth with k>0k > 0 the doubling time (ln⁡2)/k(\ln 2)/k plays the same role.

In plain words Each half-life removes half of what is left, not half of the original amount. So the amount keeps halving and never reaches zero at a finite time.

Example Start with 8080 units and a half-life of 44 hours. Then, with tt in hours,

λ=ln⁡24 h−1≈0.1733 h−1,Q(t)=80e−(ln⁡2)t/4=80⋅2−t/4.\begin{gathered} \lambda = \frac{\ln 2}{4}\ \mathrm{h}^{-1} \approx 0.1733\ \mathrm{h}^{-1}, \\[0.4em] Q(t) = 80e^{-(\ln 2)t/4} = 80 \cdot 2^{-t/4}. \end{gathered}

So Q(4)=40Q(4) = 40, Q(8)=20Q(8) = 20 and Q(12)=10Q(12) = 10. Differentiating gives Q′=−λQQ' = -\lambda Q, and Q(0)=80Q(0) = 80, so the formula satisfies both the rate law and the starting amount.

Time constant

Definition For Q′=−λQQ' = -\lambda Q with λ>0\lambda > 0, the time constant is τ=1/λ\tau = 1/\lambda, measured in units of time, so Q(t)=Q0e−t/τQ(t) = Q_0e^{-t/\tau}. Over each time constant the amount is divided by ee,

Q(t+τ)=Q(t)e≈0.368 Q(t),Q(t + \tau) = \frac{Q(t)}{e} \approx 0.368\,Q(t),

so about 36.836.8 percent remains and about 63.263.2 percent is lost. The two time scales are related by t1/2=τln⁡2≈0.693 τt_{1/2} = \tau\ln 2 \approx 0.693\,\tau, so the time constant is about 1.4431.443 half-lives. The tangent line to the decay curve at any time reaches zero exactly one time constant later. For growth with k>0k > 0, the same name is used for 1/k1/k, the time in which the amount grows by the factor ee.

In plain words The time constant is the natural clock of the decay. If the amount kept falling at its current rate, it would be gone after one time constant, but the rate shrinks along with the amount, so about 3737 percent is still there.

Example For the 44-hour half-life above, τ=4/ln⁡2≈5.771\tau = 4/\ln 2 \approx 5.771 hours, which is longer than the half-life, and Q(τ)=80/e≈29.43Q(\tau) = 80/e \approx 29.43 units. Measured in minutes, the same decay has λ=(ln⁡2)/240≈0.002888\lambda = (\ln 2)/240 \approx 0.002888 per minute and τ≈346.2\tau \approx 346.2 minutes, the same duration. The product λt\lambda t must have no units, so tt and λ\lambda have to use the same time unit.

Newton’s law of cooling

Definition Let T(t)T(t) be the temperature of an object at time tt, with T′=dT/dtT' = dT/dt, and let AA be the constant ambient temperature of its surroundings. Newton’s law of cooling is the model

T′=−k(T−A),T' = -k(T - A),

where k>0k > 0 is a constant with units of inverse time. If T>AT > A the rate is negative and the object cools, and if T<AT < A the rate is positive and it warms. The difference D=T−AD = T - A satisfies D′=−kDD' = -kD, so with T0=T(0)T_0 = T(0)

T(t)=A+(T0−A)e−kt.T(t) = A + (T_0 - A)e^{-kt}.

The model assumes that the object has nearly the same temperature throughout, that heat flows at a rate proportional to the temperature difference, and that the surroundings stay at AA.

In plain words The excess of the object’s temperature over the room’s decays exponentially. A hot object cools and a cold one warms, both toward AA, and the time constant 1/k1/k belongs to that excess, not to the temperature itself.

Example Put an object at 90∘C90^{\circ}\mathrm{C} in a room at 20∘C20^{\circ}\mathrm{C}, with k=15k = \tfrac{1}{5} per minute. The initial difference is 70∘C70^{\circ}\mathrm{C}, and with tt in minutes

T(t)=20+70e−t/5,T(t) = 20 + 70e^{-t/5},

so T(0)=90T(0) = 90. Differentiating gives T′=−14e−t/5T' = -14e^{-t/5}, which equals −15(T−20)=−15⋅70e−t/5-\tfrac{1}{5}(T - 20) = -\tfrac{1}{5} \cdot 70e^{-t/5}, so the formula solves the equation. After one time constant, 55 minutes, the excess is 70/e≈25.75∘C70/e \approx 25.75^{\circ}\mathrm{C} and T(5)≈45.75∘CT(5) \approx 45.75^{\circ}\mathrm{C}. The object reaches 40∘C40^{\circ}\mathrm{C} when 70e−t/5=2070e^{-t/5} = 20, at t=5ln⁡3.5≈6.26t = 5\ln 3.5 \approx 6.26 minutes.

Constant ambient temperature

Definition The ambient temperature is the temperature of the surroundings. Newton’s formula assumes it is a constant AA, which is reasonable when the surroundings are large enough that the object does not noticeably heat or cool them. Then the equilibrium T=AT = A is asymptotically stable, and every solution tends to it as t→∞t \to \infty. If the ambient temperature changes with time, A=A(t)A = A(t) with AA differentiable, the difference D=T−AD = T - A satisfies

D′=−kD−A′(t),D' = -kD - A'(t),

so the constant-ambient formula no longer solves the equation unless A′A' is zero. The equation T′+kT=kA(t)T' + kT = kA(t) is still linear, and the integrating factor ekte^{kt} solves it.

In plain words A fixed room temperature is the target of the motion. If the target moves, the object chases it, and the simple exponential formula no longer describes the chase.

Example A can at 4∘C4^{\circ}\mathrm{C} sits in a room at 22∘C22^{\circ}\mathrm{C} with k=0.1k = 0.1 per minute. Here T0−A=−18T_0 - A = -18, so T(t)=22−18e−0.1tT(t) = 22 - 18e^{-0.1t}, which rises toward 22∘C22^{\circ}\mathrm{C}. The excess is negative and shrinks in size, and the can is halfway, at 13∘C13^{\circ}\mathrm{C}, after 10ln⁡2≈6.9310\ln 2 \approx 6.93 minutes. That halfway time does not depend on the starting temperature, because every 10ln⁡210\ln 2 minutes halves the excess, whatever its size.

Check yourself

  1. A 5050 gram sample decays with a half-life of 66 days. How much remains after 1515 days, and what is the time constant?
Show the answer

Fifteen days are 15/6=2.515/6 = 2.5 half-lives, so the amount left, in grams, is

Q(15)=50⋅2−2.5=5042≈8.84.\begin{aligned} Q(15) &= 50 \cdot 2^{-2.5} \\[0.4em] &= \frac{50}{4\sqrt{2}} \approx 8.84. \end{aligned}

The decay rate is λ=(ln⁡2)/6\lambda = (\ln 2)/6 per day, so the time constant is τ=6/ln⁡2≈8.66\tau = 6/\ln 2 \approx 8.66 days.

  1. Coffee at 90∘C90^{\circ}\mathrm{C} sits in a room at 20∘C20^{\circ}\mathrm{C} and has cooled to 60∘C60^{\circ}\mathrm{C} after 1010 minutes. Assuming Newton’s law of cooling, when does it reach 40∘C40^{\circ}\mathrm{C}? A classmate models the temperature itself as T=90e−ctT = 90e^{-ct}. What answer does that give, and what is wrong with it?
Show the answer

The excess over the room decays, so T−20=70e−ktT - 20 = 70e^{-kt}. At t=10t = 10 the excess is 4040, so e−10k=4/7e^{-10k} = 4/7 and k=ln⁡(7/4)/10≈0.0560k = \ln(7/4)/10 \approx 0.0560 per minute. The temperature is 40∘C40^{\circ}\mathrm{C} when the excess is 2020, so e−kt=2/7e^{-kt} = 2/7 and

t=10ln⁡3.5ln⁡1.75≈22.4t = \frac{10\ln 3.5}{\ln 1.75} \approx 22.4

minutes. The classmate’s model gives e−10c=2/3e^{-10c} = 2/3 and then e−ct=4/9=(2/3)2e^{-ct} = 4/9 = (2/3)^2, so t=20t = 20 minutes, which is too early. That model sends the temperature toward 0∘C0^{\circ}\mathrm{C}, below the room, while Newton’s law moves it toward 20∘C20^{\circ}\mathrm{C} and never past it. The quantity that decays is the difference T−AT - A, not TT.

  1. The room warms steadily, with A(t)=20+tA(t) = 20 + t in degrees Celsius and tt in minutes, and an object starts at 90∘C90^{\circ}\mathrm{C} with k=15k = \tfrac{1}{5} per minute. Show that T=A(t)+70e−t/5T = A(t) + 70e^{-t/5} is not a solution, then solve the problem and describe the long-run behavior.
Show the answer

For that guess, T′=1−14e−t/5T' = 1 - 14e^{-t/5} while −15(T−A)=−14e−t/5-\tfrac{1}{5}(T - A) = -14e^{-t/5}, so the two sides differ by A′(t)=1A'(t) = 1. The equation is linear,

T′+15T=15(20+t).T' + \tfrac{1}{5}T = \tfrac{1}{5}(20 + t).

Multiplying by et/5e^{t/5} gives (et/5T)′=15et/5(20+t)(e^{t/5}T)' = \tfrac{1}{5}e^{t/5}(20 + t). Integration by parts gives the antiderivative et/5(15+t)e^{t/5}(15 + t), so T=15+t+Ce−t/5T = 15 + t + Ce^{-t/5}, and T(0)=90T(0) = 90 gives C=75C = 75. The solution is

T(t)=15+t+75e−t/5.T(t) = 15 + t + 75e^{-t/5}.

Since T−A=−5+75e−t/5→−5T - A = -5 + 75e^{-t/5} \to -5, the object ends up trailing the room by 5∘C5^{\circ}\mathrm{C}. The lag is A′/k=5A'/k = 5, the warming rate times the time constant.

How this lesson connects