MHB Alex's question at Yahoo Answers regarding a mixing problem

  • Thread starter Thread starter MarkFL
  • Start date Start date
  • Tags Tags
    Mixing
MarkFL
Gold Member
MHB
Messages
13,284
Reaction score
12
Here is the question:

How do I solve this linear equation problem?


A tank with a capacity of 1600 L is full of a mixture of water and chlorine with a concentration of 0.0125 g of chlorine per liter. In order to reduce the concentration of chlorine, fresh water is pumped into the tank at a rate of 16 L/s. The mixture is kept stirred and is pumped out at a rate of 40 L/s. Find the amount of chlorine in the tank as a function of time. (Let y be the amount of chlorine in grams and t be the time in seconds.)

I have posted a link there to this topic so the OP can see my work.
 
Mathematics news on Phys.org
Hello Alex,

We want to determine the amount $y(t)$ in grams of chlorine present in the tank at time $t$ in seconds. So, let's begin with:

$$\text{time rate of change of chlorine = time rate of chlorine coming in minus time rate of chlorine going out}$$

$$\text{time rate of chlorine coming in=concentration of solution coming in times the time rate of volume coming in}$$

$$\text{time rate of chlorine going out=concentration of solution going out times the time rate of volume going out}$$

Let:

$$C_I$$ = the constant concentration of solution coming in.

$$R_I$$ - the constant rate of incoming solution.

$$C_O$$ = the variable concentration of solution going out.

$$R_O$$ - the constant rate of outgoing solution.

Hence, stated mathematically, we may write:

$$\frac{dy}{dt}=C_IR_I-C_OR_O$$

Concentration equals amount per volume, and the volume of solution is a function of $t$, let's call it $V(t)$. So, we have:

$$\frac{dy}{dt}=C_IR_I-\frac{y(t)}{V(t)}R_O$$

To find the volume of solution present in the tank at time $t$, consider the IVP:

$$\frac{dV}{dt}=R_I-R_O$$ where $$V(0)=V_0$$

Separating variables, and integrating, we have:

$$\int\,dV=\left(R_I-R_O \right)\int\,dt$$

$$V(t)=\left(R_I-R_O \right)t+C$$

Using the initial values, we may determine the parameter $C$:

$$V(0)=C=V_0$$

Hence:

$$V(t)=\left(R_I-R_O \right)t+V_0$$

Thus, we may now write:

$$\frac{dy}{dt}=C_IR_I-\frac{y(t)}{\left(R_I-R_O \right)t+V_0}R_O$$

Arranging the ODE in standard linear form, we may model the amount of chlorine present in the tank at time $t$ with the following IVP:

$$\frac{dy}{dt}+\frac{R_O}{\left(R_I-R_O \right)t+V_0}y(t)=C_IR_I$$ where $y(0)=y_0$

Computing the integrating factor, we find:

$$\mu(t)=e^{\int \frac{R_O}{\left(R_I-R_O \right)t+V_0}\,dt}=\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}$$

Multiplying the ODE by this integrating factor, we obtain:

$$\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}\frac{dy}{dt}+R_O\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}-1}y(t)=C_IR_I\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}$$

Observing that the left side of the equation is now the differentiation of the product $\mu(t)\cdot y(t)$, we may write:

$$\frac{d}{dt}\left(\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}\cdot y(t) \right)=C_IR_I\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}$$

Integrating with respect to $t$, we have:

$$\int\,d\left(\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}\cdot y(t) \right)=C_IR_I\int\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}\,dt$$

$$\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}}\cdot y(t)=\frac{C_IR_I}{R_I-R_O}\left(\left(R_I-R_O \right)t+V_0 \right)^{\frac{R_O}{R_I-R_O}+1}+C$$

Solving for $y(t)$ we find:

$$y(t)=\frac{C_IR_I}{R_I-R_O}\left(\left(R_I-R_O \right)t+V_0 \right)+C\left(\left(R_I-R_O \right)t+V_0 \right)^{-\frac{R_O}{R_I-R_O}}$$

Using the initial conditions, we may determine the parameter $C$:

$$y(0)=\frac{C_IR_I}{R_I-R_O}\left(V_0 \right)+C\left(V_0 \right)^{-\frac{R_O}{R_I-R_O}}=y_0$$

$$C=V_0^{\frac{R_O}{R_I-R_O}}\left(y_0-\frac{C_IR_IV_0}{R_I-R_O} \right)$$

Thus, the solution satisfying the IVP is:

$$y(t)=\frac{C_IR_I}{R_I-R_O}\left(\left(R_I-R_O \right)t+V_0 \right)+V_0^{\frac{R_O}{R_I-R_O}}\left(y_0-\frac{C_IR_IV_0}{R_I-R_O} \right)\left(\left(R_I-R_O \right)t+V_0 \right)^{-\frac{R_O}{R_I-R_O}}$$

$$y(t)=\frac{C_IR_I}{R_I-R_O}\left(\left(R_I-R_O \right)t+V_0 \right)+\left(y_0-\frac{C_IR_IV_0}{R_I-R_O} \right)\left(\frac{V_0}{\left(R_I-R_O \right)t+V_0} \right)^{\frac{R_O}{R_I-R_O}}$$

Using the given data for this problem:

$$C_I=0\,\frac{\text{g}}{\text{L}},\,R_I=16\, \frac{\text{L}}{\text{s}},\,R_O=40\, \frac{\text{L}}{\text{s}},\,V_0=1600\text{ L},\,y_0=0.0125\, \frac{\text{g}}{\text{L}}\cdot1600\text{ L}=20\text{ g}$$

we have:

$$y(t)=20\left(1-\frac{3}{200}t \right)^{\frac{5}{3}}$$

Here is a plot of the solution on the relevant domain $$0\le t\le\frac{200}{3}$$:

View attachment 1601
 

Attachments

  • alexplot.jpg
    alexplot.jpg
    10.3 KB · Views: 92
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Fermat's Last Theorem has long been one of the most famous mathematical problems, and is now one of the most famous theorems. It simply states that the equation $$ a^n+b^n=c^n $$ has no solutions with positive integers if ##n>2.## It was named after Pierre de Fermat (1607-1665). The problem itself stems from the book Arithmetica by Diophantus of Alexandria. It gained popularity because Fermat noted in his copy "Cubum autem in duos cubos, aut quadratoquadratum in duos quadratoquadratos, et...
I'm interested to know whether the equation $$1 = 2 - \frac{1}{2 - \frac{1}{2 - \cdots}}$$ is true or not. It can be shown easily that if the continued fraction converges, it cannot converge to anything else than 1. It seems that if the continued fraction converges, the convergence is very slow. The apparent slowness of the convergence makes it difficult to estimate the presence of true convergence numerically. At the moment I don't know whether this converges or not.
Back
Top