5. Applications & Modelling0%

Differential Equations · Topic 5 of 5

5. Applications & Modelling

Video coming soon3 worked examples

Theory

Differential equations earn their keep in modelling. Exam questions rarely hand you the equation — they describe a situation in words and expect you to form it first. The key is that the phrase “rate of change” means a derivative with respect to time, and “proportional to” introduces a constant:

  • “the rate of increase of PP is proportional to PP    dPdt=kP\implies \dfrac{dP}{dt} = kP
  • “the rate of decrease of NN is proportional to NN    dNdt=kN\implies \dfrac{dN}{dt} = -kN
  • “the rate of cooling is proportional to the excess temperature above the surroundings θ0\theta_0    dθdt=k(θθ0)\implies \dfrac{d\theta}{dt} = -k(\theta - \theta_0)

That last one is Newton's law of cooling. The substitution u=θθ0u = \theta - \theta_0 turns it into the same separable equation as the others.

The workflow is always the same: form the equation, separate the variables and integrate, then use the given data to pin down the constants. There are usually two unknowns — the constant of integration and the constant of proportionality kk — so expect two pieces of information.

The Golden Rule: translate the words into an equation before doing any calculus, and count your unknowns against the data you have been given. If you have two unknowns and only one condition, you have mis-formed the equation.

⚠️ Common Examiner Traps

  • Sign of kk: for decay or cooling, write the minus sign into the equation and keep kk positive. Doing both — a minus sign and a negative kk — reverses the model.
  • Using conditions in the wrong order: the initial condition (at t=0t=0) gives the constant of integration most easily, because the exponential term becomes 11. Find that first, then use the second condition for kk.
  • Leaving the answer in log form: convert lnP=kt+C\ln P = kt + C to P=AektP = Ae^{kt} before substituting — it is far less error-prone.
  • Forgetting the surrounding temperature: in a cooling problem the variable that decays exponentially is the excess θθ0\theta - \theta_0, not θ\theta itself. Add θ0\theta_0 back at the end.

Worked examples

Example 1

A population grows at a rate proportional to its current size. Initially it is 500500, and after 3 hours it has reached 20002000. Find the population after 5 hours.

Step 1: Translate the description into a differential equation:

dPdt=kP\frac{dP}{dt} = kP

Step 2: Separate the variables and integrate:

1PdP=kdt    lnP=kt+C\int \frac{1}{P}\,dP = \int k\,dt \implies \ln P = kt + C

Step 3: Convert to exponential form immediately, writing A=eCA = e^{C}:

P=AektP = Ae^{kt}

Step 4: Use the initial condition P=500P = 500 when t=0t = 0. Since e0=1e^0 = 1, this gives AA at once:

500=A    P=500ekt500 = A \implies P = 500e^{kt}

Step 5: Use the second condition, P=2000P = 2000 when t=3t = 3, to find kk:

2000=500e3k    e3k=4    k=ln432000 = 500e^{3k} \implies e^{3k} = 4 \implies k = \frac{\ln 4}{3}

Step 6: Evaluate at t=5t = 5:

P=500e5k=500×45/3=500×10.0795040P = 500e^{5k} = 500 \times 4^{5/3} = 500 \times 10.079 \approx 5040

Example 2

A body at 90C90^\circ\text{C} is left to cool in a room kept at 20C20^\circ\text{C}. After 10 minutes its temperature is 60C60^\circ\text{C}. Using Newton's law of cooling, find its temperature after 25 minutes.

Step 1: Newton's law of cooling gives, with room temperature 2020:

dθdt=k(θ20)\frac{d\theta}{dt} = -k(\theta - 20)

Step 2: Substitute u=θ20u = \theta - 20. Since 2020 is constant, dudt=dθdt\frac{du}{dt} = \frac{d\theta}{dt}, so the equation becomes separable in uu:

dudt=ku    u=Aekt\frac{du}{dt} = -ku \implies u = Ae^{-kt}

Step 3: At t=0t = 0 the temperature is 9090, so the excess is u=70u = 70:

A=70    u=70ektA = 70 \implies u = 70e^{-kt}

Step 4: At t=10t = 10 the temperature is 6060, so the excess is 4040:

40=70e10k    e10k=47    k=110ln ⁣(74)=0.0559640 = 70e^{-10k} \implies e^{-10k} = \frac{4}{7} \implies k = \frac{1}{10}\ln\!\left(\frac{7}{4}\right) = 0.05596

Step 5: Find the excess at t=25t = 25:

u=70e25(0.05596)=70e1.399=70×0.2468=17.28u = 70e^{-25(0.05596)} = 70e^{-1.399} = 70 \times 0.2468 = 17.28

Step 6: Add the room temperature back to get the actual temperature:

θ=20+17.28=37.3C (1 d.p.)\theta = 20 + 17.28 = 37.3^\circ\text{C} \ \text{(1 d.p.)}

Example 3

A tank leaks so that the rate of decrease of its volume is proportional to the square root of the volume remaining. Initially it holds 100100 litres, and after 5 minutes it holds 6464 litres. Find the volume after 10 minutes.

Step 1: The rate is decreasing and proportional to V\sqrt{V}, so:

dVdt=kV\frac{dV}{dt} = -k\sqrt{V}

Step 2: Separate the variables, writing the left side with a negative index:

V1/2dV=kdt    2V1/2=kt+C\int V^{-1/2}\,dV = \int -k\,dt \implies 2V^{1/2} = -kt + C

Step 3: Apply the initial condition V=100V = 100 at t=0t = 0:

2100=C    C=202\sqrt{100} = C \implies C = 20

Step 4: Apply the second condition, V=64V = 64 at t=5t = 5, to find kk:

264=5k+20    16=5k+20    k=0.82\sqrt{64} = -5k + 20 \implies 16 = -5k + 20 \implies k = 0.8

Step 5: The model is therefore 2V=200.8t2\sqrt{V} = 20 - 0.8t. Substitute t=10t = 10:

2V=208=12    V=6    V=362\sqrt{V} = 20 - 8 = 12 \implies \sqrt{V} = 6 \implies V = 36

The tank holds 36 litres after 10 minutes.