Modeling with First Order DE’s — Question 10

PDF ↗

Question 10

A reservoir initially contains 1010 L of usable liquid. During a four-minute operation, a pump supplies a constant rate u≥0u\ge 0 L/min, scheduled demand removes liquid at rate d(t)=2+td(t)=2+t L/min, and leakage removes liquid at a rate equal to 0.500.50 times the current volume per minute. Here tt is in minutes; the slope in d(t)d(t) is 11 L/min2^2.

The reservoir must retain at least 44 L at every time 0≤t≤40\le t\le 4. Assume no capacity restriction, no other flows, and that the stated demand and leakage laws apply while liquid remains.

Tasks

  1. Derive and solve the volume IVP in terms of uu.

  2. Find the smallest constant pump rate that satisfies the reserve requirement throughout the operation. Justify feasibility over the whole interval, not only at its endpoint.

  3. For this rate, calculate the cumulative leaked volume by t=4t=4 using an overall volume balance.

  4. Compare the answer with the minimum pump rate obtained by ignoring leakage. Sketch the optimal volume and reserve level, and explain why the no-leakage design cannot meet the original requirement.

Original worksheet page 1: question and worked solution for 2-7-010
Show solutionHide solution

Question 10 – Solution

Strategy. Solve the storage balance, obtain a necessary pump rate from the final reserve, and then verify that this candidate satisfies the full time constraint.

Step 1: Model and solve the storage balance. Input minus demand minus leakage gives S′=u−(2+t)−0.50S,S(0)=10,S'=u-(2+t)-0.50S,\qquad S(0)=10, with all rates in L/min. A particular solution is 2u−2t2u-2t; adding Ae−t/2Ae^{-t/2} and applying S(0)=10S(0)=10 gives S(t)=2u−2t+(10−2u)e−t/2.\boxed{S(t)=2u-2t+(10-2u)e^{-t/2}}. It has initial value 1010, and differentiation gives the required balance on any time interval on which the physical operation remains valid.

Step 2: Find and verify the minimum rate. Since S(4)≥4S(4)\ge 4 is necessary, u≥u*=12−10e−22(1−e−2)=6+1e2−1≈6.157 L/min.\boxed{u\ge u_*= \frac{12-10e^{-2}}{2(1-e^{-2})} =6+\frac 1{e^2-1}\approx 6.157\text{ L/min}}. For this candidate, 6<u*<76<u_*<7 and S′(t)=−2+(u*−5)e−t/2<0(0≤t≤4).S'(t)=-2+(u_*-5)e^{-t/2}<0\qquad(0\le t\le 4). Thus SS decreases from 1010 to exactly 44 L and never violates the reserve. The candidate is feasible, and every smaller rate fails the necessary final condition, proving optimality.

See the diagram in the original worksheet below.

Step 3: Recover cumulative leakage. The scheduled demand totals ∫04(2+t)dt=16\int_0^4(2+t)\,dt=16 L. Initial storage plus pumped volume equals final storage plus demand plus leakage, so Lleak=10+4u*−16−4=4u*−10≈14.63 L.\boxed{L_{\mathrm{leak}}=10+4u_*-16-4=4u_*-10 \approx 14.63\text{ L}}.

Step 4: Compare the incomplete model. Ignoring leakage gives S0=10+(u−2)t−t2/2S_0=10+(u-2)t-t^2/2. This concave function has its minimum on [0,4][0,4] at an endpoint; S0(4)≥4S_0(4)\ge 4 gives the no-leakage minimum u=2.5u=2.5 L/min. That rate is below u*u_* and fails the actual reserve constraint. It may eventually exhaust the actual reservoir, so negative values of its unconstrained storage formula would not describe physical stored liquid. The optimal design above stays positive and avoids that limitation.

Original worksheet page 2: question and worked solution for 2-7-010

Original worksheet layout. Use Enlarge or open the PDF for a closer view.