Systems of Differential Equations — Question 3

PDF ↗

Question 3

Two perfectly mixed tanks have constant volumes 100100 L and 200200 L. Fresh water enters tank 11 at 33 L/min. Tank 11 sends 55 L/min to tank 22; tank 22 returns 22 L/min to tank 11 and drains 33 L/min outside. Let x,yx,y be the salt amounts in grams, with x(0)=100x(0)=100, y(0)=0y(0)=0. You may use the inward-boundary criterion: for a locally Lipschitz vector field, the nonnegative quadrant is forward invariant if its boundary derivatives point inward or are tangent.

Tasks

  1. Check both volume balances and derive the first-order system for salt amounts. Identify the concentration used in each outflow term.

  2. Write the amount system in matrix form. Then derive the system for concentrations c1=x/100c_1=x/100, c2=y/200c_2=y/200, including initial data.

  3. Find the derivative of total salt and explain why the internal transfer rates cancel even though the volumes differ. Find every equilibrium.

  4. Prove that nonnegative initial amounts remain nonnegative and that total salt cannot increase. Give explicit bounds for this initial state without solving the coupled system.

Original worksheet page 1: question and worked solution for 5-4-003
Show solutionHide solution

Question 3 – Solution

Strategy. Multiply each flow rate by the concentration in its source tank, then check the total balance.

Step 1: Balance volumes and salt separately. Tank 11 receives 3+2=53+2=5 L/min and sends 55; tank 22 receives 55 and sends 2+3=52+3=5. Their volumes are constant. The return concentration is y/200y/200 and the forward concentration is x/100x/100, giving x′=−x/20+y/100,y′=x/20−y/40.\boxed{x'=-x/20+y/100,\qquad y'=x/20-y/40.} Fresh water supplies no salt. Each coefficient has units min−1^{-1}.

Step 2: Change from amounts to concentrations. For X=(x,y)TX=(x,y)^T, X′=(−1/201/1001/20−1/40)X,X(0)=(100,0)T.X'=\begin{pmatrix}-1/20&1/100\\1/20&-1/40\end{pmatrix}X, \qquad X(0)=(100,0)^T. Dividing each balance by its own tank volume gives c′=(−1/201/501/40−1/40)c,c(0)=(1,0)T.\boxed{c'=\begin{pmatrix}-1/20&1/50\\1/40&-1/40\end{pmatrix}c, \qquad c(0)=(1,0)^T.} This is a constant invertible change of state; the two matrices need not match.

Step 3: Check the external balance and equilibria. Adding the amount equations gives (x+y)′=−3y/200,\boxed{(x+y)'=-3y/200}, exactly the outside drain rate times its source concentration. At an equilibrium this forces y=0y=0, and then y′=x/20=0y'=x/20=0 forces x=0x=0. The only equilibrium is the salt-free state.

Step 4: Prove positivity and bounds. The field is linear and therefore locally Lipschitz. On x=0x=0, y≥0y\ge 0, we have x′=y/100≥0x'=y/100\ge 0; on y=0y=0, x≥0x\ge 0, we have y′=x/20≥0y'=x/20\ge 0. The given inward-boundary criterion therefore proves that the nonnegative quadrant is forward invariant. The constant linear system exists for all forward times, so the conclusion holds for every t≥0t\ge 0. Now the total balance implies 0≤x+y≤1000\le x+y\le 100. Hence 0≤x,y≤1000\le x,y\le 100 g, 0≤c1≤10\le c_1\le 1 and 0≤c2≤1/20\le c_2\le 1/2 g/L. These are guaranteed bounds, not claims that every bound is attained.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 5-4-003

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