I2, I3 - Initial conditions and models

I2, I3 - Initial conditions and models

A differential equation gives a rate law; an initial condition selects the solution that describes a particular situation. I2 develops this choice of constants, and I3 turns statements about motion and temperature into first-order models. You will interpret the resulting functions, including their units, valid times and limitations.

Selecting a particular solution

A drink starts at 80 ∘C80\,^{\circ}\mathrm{C} in a room at 20 ∘C20\,^{\circ}\mathrm{C}. If its cooling rate is proportional to its temperature difference from the room, will it ever reach room temperature exactly? We need both a rate law and a starting temperature to answer that question.

For example, a proposed rate law y′=−0.1(y−20)y'=-0.1(y-20) has the general solution y=20+Ce−0.1ty=20+Ce^{-0.1t}. Different values of CC describe different starting temperatures; the rate law alone does not choose one. The condition y(0)=80y(0)=80 gives C=60C=60. We will interpret that solution in the last part.

Particular solution

A particular solution is a solution with its arbitrary constants fixed by stated conditions. An initial condition specifies a value of the unknown function, or a relevant derivative, at a stated value of the independent variable. Substitute into the solution, not just into one isolated term.

For a first-order linear equation with continuous coefficients on an interval, one value y(x0)=y0y(x_0)=y_0 inside that interval selects exactly one solution. Indeed, two solutions with the same value have difference C/μC/\mu; at x0x_0, this forces C=0C=0, so the solutions coincide throughout that interval.

Worked example 1: a condition at zero

Question. Solve y′+2y=8y'+2y=8, given y(0)=7y(0)=7.

Plan. Find the general solution using an integrating factor, then substitute the complete initial condition.

(e2xy)′=8e2xmultiply by e2xe2xy=4e2x+Cintegratey=4+Ce−2x7=4+Cuse x=0C=3\begin{aligned} (e^{2x}y)'&=8e^{2x} && \text{multiply by }e^{2x}\\ e^{2x}y&=4e^{2x}+C && \text{integrate}\\ y&=4+Ce^{-2x}\\ 7&=4+C && \text{use }x=0\\ C&=3 \end{aligned}

Answer. y=4+3e−2xy=4+3e^{-2x}.

Check. At zero, y=7y=7. Also y′+2y=−6e−2x+8+6e−2x=8y'+2y=-6e^{-2x}+8+6e^{-2x}=8.

Keep the differential equation and the value 77, but move the condition to x=ln⁡2x=\ln2. The exponential no longer equals one at the specified point.

Worked example 2: a condition away from zero

Question. Solve y′+2y=8y'+2y=8, given y(ln⁡2)=7y(\ln2)=7.

Plan. Use the same general solution, evaluating the exponential at the actual stated point.

(e2xy)′=8e2xy=4+Ce−2xintegrate and divide7=4+Ce−2ln⁡2=4+C4C=12\begin{aligned} (e^{2x}y)'&=8e^{2x}\\ y&=4+Ce^{-2x} && \text{integrate and divide}\\ 7&=4+Ce^{-2\ln2}=4+\frac C4\\ C&=12 \end{aligned}

Answer. y=4+12e−2xy=4+12e^{-2x}.

Check. e−2ln⁡2=1/4e^{-2\ln2}=1/4 gives y(ln⁡2)=7y(\ln2)=7. Direct substitution gives y′+2y=−24e−2x+8+24e−2x=8y'+2y=-24e^{-2x}+8+24e^{-2x}=8.

The second curve has y(0)=16y(0)=16, which is consistent with reaching 77 later. The same value at different times need not select the same curve.

Watch out

Wrong: y(ln⁡2)=7⇒C=7−4y(\ln2)=7\Rightarrow C=7-4. That would use x=0x=0. Substitute the stated independent-variable value into every term before solving for CC.

Select a new curve with a condition away from zero.

Formulating rate and motion models

A model starts with a signed statement about change. For an amount, write “rate in minus rate out”. For straight-line motion, choose a positive direction and use dv/dt=acceleration\mathrm{d}v/\mathrm{d}t=\text{acceleration}; a force balance then gives m dv/dt=resultant forcem\,\mathrm{d}v/\mathrm{d}t=\text{resultant force}. The mass-times-acceleration rule is enough here; no additional motion formula is needed.

For a well-mixed tank of fixed volume VV, containing amount qq, an outflow of rr volume units per unit time removes material at rate rq/Vrq/V. If material enters at rate a(t)a(t), then q′=a(t)−(r/V)qq'=a(t)-(r/V)q. “Well mixed” is the assumption that makes the exiting concentration equal to q/Vq/V.

For cooling in an environment at temperature u(t)u(t), use y′=−k(y−u)y'=-k(y-u), where k>0k>0. Above the environment, y′<0y'<0; below it, y′>0y'>0. This sign check is more reliable than interpreting the word “rate” as automatically positive. The units of kk are reciprocal time.

Worked example 3: constant driving force

Question. A particle of mass 2 kg2\,\mathrm{kg} moves in a straight line. The forward force is 6 N6\,\mathrm{N}, and resistance is 2v N2v\,\mathrm{N}, opposing forward velocity v m s−1v\,\mathrm{m}\,\mathrm{s}^{-1}. It starts from rest at t=0t=0. Find v(t)v(t) and its displacement s(t)s(t) from the starting point.

Plan. Form the signed force equation, solve for velocity, then integrate velocity with s(0)=0s(0)=0.

2v′=6−2vforward minus resistancev′+v=3(etv)′=3etintegrating factor etv=3+Ce−tv(0)=0⇒C=−3s=3t+3e−t+Dintegrate velocitys(0)=0⇒D=−3\begin{aligned} 2v'&=6-2v && \text{forward minus resistance}\\ v'+v&=3\\ (e^tv)'&=3e^t && \text{integrating factor }e^t\\ v&=3+Ce^{-t}\\ v(0)=0&\Rightarrow C=-3\\ s&=3t+3e^{-t}+D && \text{integrate velocity}\\ s(0)=0&\Rightarrow D=-3 \end{aligned}

Answer. v=3−3e−tv=3-3e^{-t} and s=3t−3+3e−ts=3t-3+3e^{-t}, for t≥0t\ge0, in metres per second and metres respectively.

Check. 2v′=6e−t=6−2v2v'=6e^{-t}=6-2v, s′=vs'=v, and both starting values are zero. Velocity tends to 3 m s−13\,\mathrm{m}\,\mathrm{s}^{-1}, where the two forces balance.

Now change only the driving force to one that decreases with time. The same resistance law and starting conditions remain.

Worked example 4: a changing driving force

Question. In Worked example 3, replace the forward force by (6−2t) N(6-2t)\,\mathrm{N}, for 0≤t≤30\le t\le3. Find v(t)v(t) and s(t)s(t) on this interval.

Plan. Keep the time-dependent force inside the differential equation and integrate its product with the factor.

2v′=6−2t−2v(etv)′=(3−t)etnormalise and multiplyetv=(4−t)et+Cintegration by partsv=4−t+Ce−tv(0)=0⇒C=−4s=4t−12t2+4e−t+Ds(0)=0⇒D=−4\begin{aligned} 2v'&=6-2t-2v\\ (e^tv)'&=(3-t)e^t && \text{normalise and multiply}\\ e^tv&=(4-t)e^t+C && \text{integration by parts}\\ v&=4-t+Ce^{-t}\\ v(0)=0&\Rightarrow C=-4\\ s&=4t-\frac12t^2+4e^{-t}+D\\ s(0)=0&\Rightarrow D=-4 \end{aligned}

The integration step can be checked immediately: ((4−t)et)′=(3−t)et((4-t)e^t)'=(3-t)e^t.

Answer. v=4−t−4e−tv=4-t-4e^{-t}, s=4t−12t2+4e−t−4s=4t-\frac12t^2+4e^{-t}-4, for 0≤t≤30\le t\le3.

Check. v′+v=(−1+4e−t)+(4−t−4e−t)=3−tv'+v=(-1+4e^{-t})+(4-t-4e^{-t})=3-t; also s′=vs'=v and the two initial conditions hold. Here v≥0v\ge0: it rises then falls, with endpoint values 00 and 1−4e−3>01-4e^{-3}>0, so the assumed resistance direction is consistent throughout the interval.

The formula may exist for times beyond 33, but the stated force model has only been specified up to that time. A mathematical expression does not extend a model's licence to describe the world.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account