Modeling with First Order DE’s — Question 7

PDF ↗

Question 7

A well-mixed pond maintains volume 10001000 m3^3. Water enters and leaves at 100100 m3^3/day. The entering water has a constant pollutant concentration u≥0u\ge 0 kg/m3^3, which is to be selected. Pollutant also disappears inside the pond at a rate equal to 0.100.10 times its current mass per day. Initially the pond concentration is 0.080.08 kg/m3^3.

For this hypothetical design problem, require the pond concentration to stay at or below 0.100.10 kg/m3^3 throughout the first five days. Assume perfect mixing and that disappearance does not change the liquid volume.

Tasks

  1. Derive the pollutant mass balance and the concentration IVP, keeping flushing and internal disappearance distinct.

  2. Solve for the concentration in terms of uu and verify the initial condition.

  3. Determine the largest constant incoming concentration satisfying the five-day requirement, and justify checking the entire time interval.

  4. Find the largest constant incoming concentration that would satisfy the same bound for all future time. Explain why it differs from the finite-horizon answer.

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

Question 7 – Solution

Strategy. Separate input, outflow and internal loss, then use the solution’s monotonicity to turn a whole-interval constraint into a justified endpoint test.

Step 1: Write the balances. If MM is pollutant mass in kg and c=M/1000c=M/1000, then M′=100u−100M1000−0.10M,M(0)=80.M'=100u-100\frac{M}{1000}-0.10M,\qquad M(0)=80. Each term is in kg/day. Dividing by the constant volume gives c′=0.10u−0.20c,c(0)=0.08.c'=0.10u-0.20c,\qquad c(0)=0.08. The total loss coefficient combines flushing at 0.100.10 day−1^{-1} with internal disappearance at the same rate.

Step 2: Solve and identify monotonicity. An integrating factor e0.20te^{0.20t} gives c(t)=u2+(0.08−u/2)e−0.20t,t≥0.\boxed{c(t)=\frac u2+(0.08-u/2)e^{-0.20t}},\qquad t\ge 0. Its initial value is 0.080.08, and differentiating verifies c′=0.10u−0.20cc'=0.10u-0.20c. It is decreasing when u<0.16u<0.16, constant when u=0.16u=0.16, and increasing when u>0.16u>0.16.

Step 3: Enforce the five-day bound. For u≤0.16u\le 0.16, the maximum is the acceptable initial concentration. For u>0.16u>0.16, the maximum on [0,5][0,5] is c(5)c(5). Therefore the largest allowed input is u5=2(0.10−0.08e−1)1−e−1≈0.2233 kg/m3.\boxed{u_5=\frac{2(0.10-0.08e^{-1})}{1-e^{-1}} \approx 0.2233\text{ kg/m}^3}. This value exceeds 0.160.16, so its concentration rises monotonically from 0.080.08 to exactly 0.100.10. Any larger uu fails at day five, proving optimality rather than just feasibility.

Step 4: Compare a permanent operating requirement. The limiting pond concentration is u/2u/2. A bound valid for all t≥0t\ge 0 requires and is guaranteed by u≤0.20 kg/m3.\boxed{u\le 0.20\text{ kg/m}^3}. For these inputs, c(t)c(t) lies between its initial value and u/2u/2. If u>0.20u>0.20, the limit exceeds the bound, so it must eventually be violated. The larger finite-horizon input uses the initially cleaner pond as temporary dilution capacity; it cannot be continued indefinitely under the same requirement.

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

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