3.3 Rate Problems
A surprising number of engineering models start from one sentence: the rate of change of a quantity equals the rate at which it comes in minus the rate at which it goes out. Write that sentence in symbols and you have a first-order differential equation. In this section we apply it to radioactive decay and population growth, to a hot part cooling in air, to salt in a stirred tank, to a population that runs out of room, and to a simple chemical reaction.
The mathematics is always a separable or a linear equation from Chapter 2. What changes from problem to problem is the modelling: deciding what the rates are, and using measured data to find the constants.
- Name the unknown and its units, and the time unit.
- Write the balance , with units checked term by term.
- Solve the equation (separable or linear) with the initial condition.
- Use extra data (a second measurement) to find any unknown constant.
- Answer the question and check that the result makes physical sense, especially as .
Exponential growth and decay
Many quantities change at a rate proportional to the amount present. In a radioactive sample, each nucleus has the same small chance of decaying per second, so the number of decays per second is proportional to the number of nuclei. In a bacterial culture with plenty of food, each cell divides at about the same rate, so the growth rate is proportional to the number of cells. Both lead to
with for growth and for decay. Separating variables () and using the initial condition gives
The constant is a rate per unit time: means that about of the current amount is lost per hour, at that instant. A more intuitive way to describe it is the time needed to halve or double.
For decay, with , the half-life is the time after which half of any amount remains:
For growth, with , the doubling time is .
The half-life does not depend on or on when you start counting. After one half-life remains, after two , after half-lives . So a convenient form of the solution is
Technetium-99m, widely used in medical imaging, has a half-life of about h. A hospital receives mg. Find the decay constant, the amount left after one day, and the time when only mg remains.
Solution
- Decay constant. .
- Solution. mg, in hours.
- After 24 h. h is four half-lives: mg.
- When is ? . Take logarithms:
The short half-life is the point of this isotope: it lasts long enough for imaging and then decays away within a couple of days.
Often is not given. Instead we have two measurements, and we must find both and from them.
A bacterial culture grows at a rate proportional to its size. Two hours after the start of an experiment it contains cells, and six hours after the start cells. How many cells were there at the start? What is the doubling time? When will the culture reach cells?
Solution
- Model. , in hours.
- Divide the two data. . Dividing eliminates . So
- Initial size. cells.
- Doubling time. h.
- 20 000 cells. Start from the measurement at : , so and h.
In the calculator below, enter the amount at and one later measurement. It finds , the half-life or doubling time, and answers a query (an amount at a given time, or the time for a given amount). The preset is Example 3.3.2 with the clock started at the first measurement.
Radioactive decay follows the exponential law extremely well. Population growth follows it only while nothing limits it. A single bacterium that doubled every 3 hours, without limit, would produce more mass than the whole Earth in less than three weeks. Real populations slow down; the logistic model later in this section captures that.
Newton's law of cooling
Put a hot metal part on a bench. It loses heat to the surrounding air, quickly at first and then more and more slowly. Newton's law of cooling states that the rate of change of the body's temperature is proportional to the difference between its temperature and the ambient (surrounding) temperature :
The minus sign makes physical sense: if , the temperature decreases; if , it increases (the same law describes warming).
The law holds well when heat leaves mainly by convection, the temperature inside the body is nearly uniform, and the temperature difference is moderate. For very hot bodies, radiation (which grows like the fourth power of absolute temperature) becomes important, and the law is only a first approximation.
Solving the equation
- New unknown. Let (the "excess temperature"). Since is constant, , and the equation becomes .
- Exponential law. , where .
- Back to .
Every body approaches the ambient temperature, and the excess temperature decays exponentially. The constant depends on the body (its surface area, mass, heat capacity) and on the air flow. It is almost always found from a measurement.
An aluminium casting leaves the mould at C and is left to cool in a workshop at C. After min its temperature is C. Find its temperature after half an hour, and the time when it is cool enough to handle (C).
Solution
- Model. , in minutes.
- Find . gives , so and
- After 30 min. No need to compute numerically: , so
- Safe to handle. gives , so
A temperature sensor stored at C is brought into a laboratory at C. After min it reads C. Assuming Newton's law, how long until it reads within C of the room temperature?
Solution
- Model. .
- Find . gives : the gap halves every minutes.
- Within 1 °C. We need , i.e. the gap must shrink by a factor :
Engineers call the time constant of the sensor (here min). A sensor needs several time constants to "catch up" with a sudden change. This lag must be taken into account in any temperature control loop.
Try your own data in the cooling calculator. Give the ambient temperature, the initial temperature and one later reading; you can ask for the temperature at a given time or the time at which a given temperature is reached.
Mixing problems
A tank holds a solution of salt (or a pollutant, a dye, a chemical) in water. A solution of known concentration flows in, the tank is kept well stirred, and the mixture flows out (Figure 3.3.3). How much salt is in the tank at time ?
Let be the amount of salt (kg) and the volume of liquid (L). Salt enters with the inflow and leaves with the outflow:
- rate in kg/min;
- rate out , because "well stirred" means the outflow has the same concentration as the whole tank.
The volume changes at the net flow rate, so . The balance is
This is a linear first-order equation in . When , the volume is constant and the equation is also separable.
The outflow carries the concentration of the tank, , not the inflow concentration . And when the flow rates differ, is a function of ; writing the initial volume in the denominator is wrong.
A L tank contains brine with kg of dissolved salt. Brine containing kg of salt per litre flows in at L/min, and the well-stirred mixture flows out at the same rate. Find the amount of salt after one hour, the long-term amount, and the time when there are kg of salt in the tank.
Solution
- Rates. In: kg/min. The volume stays L, so out: kg/min.
- Equation. , .
- Solve. Write it as . This has the same form as Newton's law of cooling with "ambient value" , so
- After 60 min. kg.
- Long term. kg. This is : in the end the tank has the inflow concentration.
- 40 kg. , so min.
A tank with a capacity of L initially holds L of water in which kg of salt is dissolved. Brine with kg/L flows in at L/min, and the stirred mixture flows out at L/min. How much salt is in the tank at the moment it starts to overflow, and what is the concentration then?
Solution
- Volume. The tank gains L/min, so . It is full when , i.e. at min.
- Rates. In: kg/min. Out: kg/min.
- Equation. , . This is linear.
- Integrating factor. .
- Integrate.
- Initial condition. At : , so .
- Solution.
- At overflow (). kg, and the concentration is kg/L.
The concentration is close to, but still below, the inflow value kg/L, as expected. After min the model no longer applies: the outflow rate changes when the tank overflows.
Use the mixing simulator to explore. With the volume is constant; otherwise it shows and the linear equation with its integrating factor.
Logistic population growth
Exponential growth cannot last. Food, space or nutrients run short, and the growth rate per individual falls as the population grows. The simplest model of this lets the per-capita growth rate decrease linearly, from a value when the population is tiny to zero when it reaches a level :
This is the logistic equation. The constant is the intrinsic growth rate and is the carrying capacity, the largest population the environment can sustain. The same equation describes the spread of a new technology through a market, or of an infection through a closed group.
Solving the logistic equation
Before solving, look at the right side. It is zero at and at , so these constants are equilibrium solutions. For the right side is positive (growth); for it is negative (decline). Now let us solve with , .
- Rewrite and separate. , so
- Partial fractions.
- Integrate. where from the initial condition (this keeps the sign right in both cases and ).
- Solve for . , so and
For , and , the solution of , , is
It satisfies as . If , the graph is S-shaped, with an inflection point where ; there the growth rate takes its largest value, .
Where does the inflection come from? Differentiate the equation once more with the chain rule: . While we have , so changes sign exactly where . Growth speeds up below half the capacity and slows down above it.
A lake is stocked with fish. Ecologists estimate that the lake can support fish and that the intrinsic growth rate is per year. Using the logistic model, find , the population after years, the time when the population reaches and when it reaches , and the largest yearly growth rate.
Solution
- Solution. , so
- After 10 years. fish.
- 4000 fish. , so and years.
- 7000 fish. , so and years.
- Largest growth rate. At : fish per year.
Fisheries use exactly this fact: harvesting is most productive when the stock is kept near , where the population replaces itself fastest.
The logistic simulator draws for your own , and and marks the inflection point.
Chemical reaction rates
In chemical kinetics, the rate of a reaction is the rate at which reactants are used up (or products form). Two simple cases lead to equations we can solve.
First-order reactions. If a substance A converts to products at a rate proportional to its own concentration, . This is the decay law again: , with a half-life that does not depend on the starting concentration. The decomposition of many compounds, and the decay of nuclei, behave this way.
Second-order reactions. In a reaction , one molecule of A meets one molecule of B. The rate is then proportional to the product of the two concentrations. Let the initial concentrations be and (mol/L), and let be the concentration of C formed. Each mole of C uses one mole of A and one of B, so the concentrations left are and , and
The equation is separable. For , the partial-fraction identity
gives, after integration,
As , approaches the smaller of and : the reaction stops when the limiting reactant is used up.
In the reaction , the initial concentrations are mol/L and mol/L. After min, mol/L of C has formed. Find the rate constant, , the concentration of C after min, and its final value.
Solution
- Equation. , . Here , so .
- Integrate. . At : . Subtract and divide by 5:
- Rate constant. At , : , so and
- Solve for . Write . Then , which rearranges to
- After 30 min. , so mol/L.
- Final value. As , and mol/L: all of B, the limiting reactant, is used up.
Exercises
- Carbon-14 has a half-life of about years. A wooden beam from an archaeological site contains of the carbon-14 found in living wood. Estimate its age.
Answer
years. - The concentration of a drug in the blood decays exponentially. It is mg/L two hours after an injection and mg/L six hours after it. Find the initial concentration and the half-life.
Answer
mg/L; h. - A sample of transformer oil cools from C to C in min in a room at C. Find its temperature after min and the time when it reaches C.
Answer
; C; min. - A thermometer at C is placed in a bath at C and reads C after s. When does it read C?
Answer
; s. - A body cooling in a room of unknown temperature reads , and C at , and min. Assuming Newton's law, find the room temperature and . (Hint: with , show .)
Answer
C; , . - A L tank is full of pure water. Brine with kg/L enters at L/min and the mixture leaves at the same rate. Find the salt after min and in the long run.
Answer
; kg; kg. - A L tank holds water with kg of dissolved salt. Pure water flows in at L/min and the mixture flows out at L/min. When will only kg of salt remain?
Answer
; min. - A L tank holds water with kg of salt. Pure water flows in at L/min and the mixture is drained at L/min. Find , and the amount and concentration of salt after min.
Answer
; for ; kg, concentration kg/L. - A population obeys the logistic law with per year, and . Find and the time when .
Answer
; years. - An invasive water plant covers m of a m lake. Five weeks later it covers m. Assuming logistic growth with m, find and the time when half the lake is covered.
Answer
with per week; half covered at weeks. - In with mol/L and L/(mol·min), the product concentration satisfies . Find and the time when half of A has reacted.
Answer
, so ; half reacted () at min. - In a first-order reaction, of the reactant is converted in min. When will be converted?
Answer
; min.