Step 1: Always check volume balance first. For the concentration formulas \( Q_{i}/V_{i} \) to stay valid, each tank's total inflow rate must equal its total outflow rate.
- Tank 1: in \( =6\,\text{L/min} \) (external) \( +\,3\,\text{L/min} \) (from T2) \( =9\,\text{L/min} \) ; out \( =9\,\text{L/min} \) . Balanced. ✓
- Tank 2: in \( =4\,\text{L/min} \) (external) \( +\,6\,\text{L/min} \) (from T1) \( =10\,\text{L/min} \) ; out \( =10\,\text{L/min} \) (of which \( 3\,\text{L/min} \) returns to Tank 1 and \( 7\,\text{L/min} \) leaves the system). Balanced. ✓
Both volumes are constant, so we may use the simple concentration formulas \( Q_{1}/25 \) and \( Q_{2}/30 \) throughout.
Step 2: Write each rate equation as (rate in) \( - \) (rate out). For each salt stream, the mass rate is (concentration) \( \,\times \, \) (volumetric flow rate).
Tank 1 derivative. Salt enters from the external pipe and from Tank 2; it leaves through the single \( 9\,\text{L/min} \) outflow.
\[ Q_{1}^{'}=\underbrace{\underbrace{(3\,\text{g/L})(6\,\text{L/min})}}+\underbrace{\underbrace{\left( \dfrac{Q_{2}}{30} \right)(3\,\text{L/min})}}-\underbrace{\underbrace{\left( \dfrac{Q_{1}}{25} \right)(9\,\text{L/min})}}=18+\dfrac{Q_{2}}{10}-\dfrac{9Q_{1}}{25} \]
Tank 2 derivative. Salt enters from the external pipe and from Tank 1; it leaves through the total \( 10\,\text{L/min} \) outflow (we do not need to split it into the " \( 3 \) to T1" and " \( 7 \) out" parts — the concentration leaving is the same regardless of destination).
\[ Q_{2}^{'}=\underbrace{\underbrace{(5\,\text{g/L})(4\,\text{L/min})}}+\underbrace{\underbrace{\left( \dfrac{Q_{1}}{25} \right)(6\,\text{L/min})}}-\underbrace{\underbrace{\left( \dfrac{Q_{2}}{30} \right)(10\,\text{L/min})}}=20+\dfrac{6Q_{1}}{25}-\dfrac{Q_{2}}{3} \]
Step 3: Collect in matrix-vector form. Read the coefficients of \( Q_{1} \) and \( Q_{2} \) off each equation to fill in the \( 2\times 2 \) matrix \( Q \) , and the constants into the forcing vector \( \mathbf{g} \) .
\[ \begin{bmatrix} Q_{1}^{'} \\ Q_{2}^{'} \end{bmatrix}=\begin{bmatrix} -\dfrac{9}{25} & \dfrac{1}{10} \\ \dfrac{6}{25} & -\dfrac{1}{3} \end{bmatrix}\begin{bmatrix} Q_{1} \\ Q_{2} \end{bmatrix}+\begin{bmatrix} 18 \\ 20 \end{bmatrix} \]
Append the initial conditions \( Q_{1}(0)=80\,\text{g} \) and \( Q_{2}(0)=50\,\text{g} \) to complete the IVP.
\[ \boxed{{\mathbf{Q}}^{'}(t)=\begin{bmatrix} -9/25 & 1/10 \\ 6/25 & -1/3 \end{bmatrix}\mathbf{Q}(t)+\begin{bmatrix} 18 \\ 20 \end{bmatrix},\,\mathbf{Q}(0)=\begin{bmatrix} 80 \\ 50 \end{bmatrix}} \]
↑ Back to question index