# 3.2: Cooling and Mixing

- Page ID
- 98036

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)## Newton’s Law of Cooling

Newton’s law of cooling states that if an object with temperature \(T(t)\) at time \(t\) is in a medium with temperature \(T_m(t)\), the rate of change of \(T\) at time \(t\) is proportional to \(T(t)-T_m(t)\); thus, \(T\) satisfies a differential equation of the form

\[\label{eq:4.2.1} T'=-k(T-T_m).\]

Here \(k > 0\), since the temperature of the object must decrease if \(T > T_m\), or increase if \(T < T_m\). We’ll call \(k\) the *temperature decay constant of the medium*.

For simplicity, in this section we’ll assume that the medium is maintained at a constant temperature \(T_m\). This is another example of building a simple mathematical model for a physical phenomenon. Like most mathematical models it has its limitations. For example, it is reasonable to assume that the temperature of a room remains approximately constant if the cooling object is a cup of coffee, but perhaps not if it is a huge cauldron of molten metal. (For more on this see *Exercise 4.2.17*.)

To solve Equation \ref{eq:4.2.1}, we rewrite it as

\[T'+kT=kT_m. \nonumber\]

Since \(e^{-kt}\) is a solution of the complementary equation, the solutions of this equation are of the form \(T=ue^{-kt}\), where \(u'e^{-kt}=kT_m\), so \(u'=kT_me^{kt}\). Hence,

\[u=T_me^{kt}+c, \nonumber\]

so

\[T=ue^{-kt}=T_m+ce^{-kt}. \nonumber\]

If \(T(0)=T_0\), setting \(t=0\) here yields \(c=T_0-T_m\), so

\[\label{eq:4.2.2} T=T_m+(T_0-T_m)e^{-kt}. \]

Note that \(T-T_m\) decays exponentially, with decay constant \(k\).

Below is a video on Newton's Law of Cooling.

A ceramic insulator is baked at \(400^\circ\)C and cooled in a room in which the temperature is \(25^\circ\)C. After 4 minutes the temperature of the insulator is \(200^\circ\)C. What is its temperature after 8 minutes?

**Solution**

Here \(T_0=400\) and \(T_m=25\), so Equation \ref{eq:4.2.2} becomes

\[\label{eq:4.2.3} T=25+375e^{-kt}.\]

We determine \(k\) from the stated condition that \(T(4)=200\); that is,

\[200=25+375e^{-4k}; \nonumber\]

hence,

\[e^{-4k} = {175\over 375} = {7\over 15}. \nonumber\]

Taking logarithms and solving for \(k\) yields

\[k=-{1\over 4} \ln {7\over 15}={1\over 4}\ln {15\over 7}. \nonumber\]

Substituting this into Equation \ref{eq:4.2.3} yields

\[T=25+375 e^{-{t\over 4} \ln {15\over 7}} \nonumber\]

(Figure 3.2.1 ). Therefore the temperature of the insulator after 8 minutes is

\[\begin{array}{rl} T(8) & = 25+375 e^{-2 \ln {15\over 7}} \\[9pt] & = 25+375 \left({7\over 15}\right)^2 \approx 107^\circ \mbox{C}. \end{array}\nonumber \]

An object with temperature \(72^\circ\)F is placed outside, where the temperature is \(-20^\circ\)F. At 11:05 the temperature of the object is \(60^\circ\)F and at 11:07 its temperature is \(50^\circ\)F. At what time was the object placed outside?

**Solution**

Let \(T(t)\) be the temperature of the object at time \(t\). For convenience, we choose the origin \(t_0=0\) of the time scale to be 11:05 so that \(T_0=60\). We must determine the time \(\tau\) when \(T(\tau)=72\). Substituting \(T_0=60\) and \(T_m=-20\) into Equation \ref{eq:4.2.2} yields

\[T = -20+\bigl(60-(-20)\bigr)e^{-kt}\nonumber\]

or

\[\label{eq:4.2.4} T = -20+80e^{-kt}. \]

We obtain \(k\) from the stated condition that the temperature of the object is 50\(^\circ\)F at 11:07. Since 11:07 is \(t=2\) on our time scale, we can determine \(k\) by substituting \(T=50\) and \(t=2\) into Equation \ref{eq:4.2.4} to obtain

\[50 = -20+80e^{-2k}\nonumber\]

This is shown in Figure 3.2.2 .

Hence,

\[e^{-2k}={70\over 80}={7\over 8}.\nonumber\]

Taking logarithms and solving for \(k\) yields

\[k =-{1\over 2} \ln {7\over 8} = {1\over 2} \ln {8\over 7}.\nonumber\]

Substituting this into Equation \ref{eq:4.2.4} yields

\[T = -20+80 e^{-{t\over 2}\ln {8\over 7}},\nonumber\]

and the condition \(T(\tau)=72\) implies that

\[72 =-20+80 e^{-{\tau\over 2} \ln {8\over 7}};\nonumber\]

hence,

\[e^{-{\tau\over 2} \ln {8\over 7}} ={92\over 80} = {23\over 20}.\nonumber\]

Taking logarithms and solving for \(\tau\) yields

\[\tau = -{2 \ln {23\over 20}\over \ln {8\over 7}} \approx-2.09\ \mbox{min}.\nonumber\]

Therefore the object was placed outside about 2 minutes and 5 seconds before 11:05; that is, at 11:02:55.

## Mixing Problems

In the next two examples a saltwater solution with a given concentration (weight of salt per unit volume of solution) is added at a specified rate to a tank that initially contains saltwater with a different concentration. The problem is to determine the quantity of salt in the tank as a function of time. This is an example of a *mixing problem*. To construct a tractable mathematical model for mixing problems we assume in our examples (and most exercises) that the mixture is stirred instantly so that the salt is always uniformly distributed throughout the mixture. *Exercises 4.2.22* and *4.2.23* deal with situations where this isn’t so, but the distribution of salt becomes approximately uniform as \(t\to\infty\).

A tank initially contains 40 pounds of salt dissolved in 600 gallons of water. Starting at \(t_0 = 0\), water that contains 1/2 pound of salt per gallon is poured into the tank at the rate of 4 gal/min and the mixture is drained from the tank at the same rate (Figure 3.2.3 ).

- Find a differential equation for the quantity \(Q(t)\) of salt in the tank at time \(t > 0\), and solve the equation to determine \(Q(t)\).
- Find \(\lim_{t\to\infty}Q(t)\).

**Solution a**

To find a differential equation for \(Q\), we must use the given information to derive an expression for \(Q'\). But \(Q'\) is the rate of change of the quantity of salt in the tank changes with respect to time; thus, if *rate in* denotes the rate at which salt enters the tank and *rate out* denotes the rate by which it leaves, then

\[\label{eq:4.2.5} Q' = \mbox{rate in}-\mbox{rate out}.\]

The rate in is

\[\left({1\over 2}\ \mbox{lb/gal}\right) \times (4\ \mbox{gal/min}) = 2\ \mbox{lb/min}.\nonumber\]

Determining the rate out requires a little more thought. We’re removing 4 gallons of the mixture per minute, and there are always 600 gallons in the tank; that is, we are removing \(1/150\) of the mixture per minute. Since the salt is evenly distributed in the mixture, we are also removing \(1/150\) of the salt per minute. Therefore, if there are \(Q(t)\) pounds of salt in the tank at time \(t\), the rate out at any time \(t\) is \(Q(t)/150\). Alternatively, we can arrive at this conclusion by arguing that

\[\begin{align*} \mbox{rate out} & = (\mbox{concentration})\times(\mbox{rate of flow out}) \\[5pt] &= (\mbox{lb/gal})\times(\mbox{gal/min})\\[4pt] &= {Q(t)\over600} \times 4 \\[4pt] &={Q(t)\over150}. \end{align*}\nonumber \]

We can now write Equation \ref{eq:4.2.5} as

\[Q' = 2-{Q\over 150}. \nonumber\]

This first order equation can be rewritten as

\[Q'+{Q\over 150} = 2. \nonumber\]

Since \(e^{-t/150}\) is a solution of the complementary equation, the solutions of this equation are of the form \(Q=ue^{-t/150}\), where \(u'e^{-t/150}=2\), so \(u'=2e^{t/150}\). Hence,

\[u = 300e^{t/150}+c, \nonumber\]

so

\[\label{eq:4.2.6} Q=ue^{-t/150}=300+ce^{-t/150}\]

(Figure 3.2.4 ). Since \(Q(0)=40\), \(c=-260\); therefore,

\[Q=300-260e^{-t/150}.\nonumber\]

**Solution b**

From Equation \ref{eq:4.2.6}, we see that that \(\lim_{t \to \infty}Q(t)=300\) for any value of \(Q(0)\). This is intuitively reasonable, since the incoming solution contains 1/2 pound of salt per gallon and there are always 600 gallons of water in the tank.

Below is a video on mixing.

A 500-liter tank initially contains 10 g of salt dissolved in 200 liters of water. Starting at \(t_0=0\), water that contains 1/4 g of salt per liter is poured into the tank at the rate of 4 liters/min and the mixture is drained from the tank at the rate of 2 liters/min (Figure [figure:4.2.5}). Find a differential equation for the quantity \(Q(t)\) of salt in the tank at time \(t\) prior to the time when the tank overflows and find the concentration \(K(t)\) (g/liter) of salt in the tank at any such time.

**Solution**

We first determine the amount \(W(t)\) of solution in the tank at any time \(t\) prior to overflow. Since \(W(0) = 200\) and we are adding 4 liters/min while removing only 2 liters/min, there’s a net gain of 2 liters/min in the tank; therefore,

\[W(t) = 2t+200.\nonumber\]

Since \(W(150)=500\) liters (capacity of the tank), this formula is valid for \(0 \le t \le 150\).

Now let \(Q(t)\) be the number of grams of salt in the tank at time \(t\), where \(0 \le t \le 150\). As in Example 3.2.3

\[\label{eq:4.2.7} Q' = \mbox{rate in}-\mbox{rate out}.\]

The rate in is

\[\label{eq:4.2.8} \left({1\over 4}\ \mbox{g/liter}\,\right) \times (4\ \mbox{liters/min}\,) = 1\ \mbox{g/min}.\]

To determine the rate out, we observe that since the mixture is being removed from the tank at the constant rate of 2 liters/min and there are \(2t+200\) liters in the tank at time \(t\), the fraction of the mixture being removed per minute at time \(t\) is

\[{2\over 2t+200} = {1\over t+100}.\nonumber\]

We’re removing this same fraction of the salt per minute. Therefore, since there are \(Q(t)\) grams of salt in the tank at time \(t\),

\[\label{eq:4.2.9} \mbox{rate out} = {Q(t)\over t+100}.\]

Alternatively, we can arrive at this conclusion by arguing that

\[\begin{array}{lcl} \mbox{rate out} & = & (\mbox{concentration})\times(\mbox{rate of flow out}) =(\mbox{g/liter})\times(\mbox{liters/min})\\[10pt] &=&{Q(t)\over2t+200}\times 2={Q(t)\over t+100}. \end{array} \nonumber \]

Substituting Equation \ref{eq:4.2.8} and Equation \ref{eq:4.2.9} into Equation \ref{eq:4.2.7} yields

\[\label{eq:4.2.10} Q'=1-{Q\over t+100},\quad \text{so} \quad Q'+{1\over t+100} Q=1.\]

By separation of variables, \(1/(t+100)\) is a solution of the complementary equation, so the solutions of Equation \ref{eq:4.2.10} are of the form

\[Q={u\over t+100}, \quad \text{where} \quad {u'\over t+100 }=1, \quad \text{so} \quad u'=t+100.\nonumber\]

Hence,

\[\label{eq:4.2.11} u = {(t+100)^2\over 2}+c.\]

Since \(Q(0)=10\) and \(u=(t+100)Q\), Equation \ref{eq:4.2.11} implies that

\[(100)(10) = {(100)^2\over 2}+c,\nonumber\]

so

\[c=100(10)-{(100)^2\over 2} =-4000\nonumber\]

and therefore

\[u = {(t+100)^2\over 2} -4000.\nonumber\]

Hence,

\[Q = {u\over t+200}= {t+100\over 2}-{4000\over t+100}.\nonumber\]

Now let \(K(t)\) be the concentration of salt at time \(t\). Then

\[K(t) = {1\over 4}-{2000\over(t+100)^2}\nonumber\]

This is shown in Figure 3.2.6 .

Below is a video on salt in a tank and differential equations.