Linearizing a Nonlinear Mixing Tank System: Finding Steady-State Concentrations

[danago]

I have a problem to solve where i am required to linearise a nonlinear system. Basically, the system involves a mixing tank of height H(t) and cross section area A with two input streams with flow rates Q₁(t) and Q₂(t) with concentrations C₁ and C₂ of some species (concentrations are constant). There is one outlet stream with flow rate Q(t) and concentration C(t). We are assuming that the tank is always perfectly mixed and that the outlet flow rate is proportional to \sqrt{H(t)}.

So what i need to do is express the concentration of the outlet stream as some nominal steady-state value plus an additional term for deviations from this nominal value

i.e. C(t) = C₀ + c(t)

Where c(t) is a linear approximation. My first thought is to write a material balance on the species in the fluid:

i.e. \frac{d}{{dt}}\left( {CAH(t)} \right) = {C_1}{Q_1} + {C_2}{Q_2} - CQ

It wouldn't be too hard to solve this DE for C(t), but I am not sure where to go from there. I guess i could then differentiate C(t) with respect to each variable to linearise it, but this seems a little tedious so I am not sure if there is something else i should be doing?

Any input is greatly appreciated.

Thanks,
Dan.

 

[CEL]

I have a problem to solve where i am required to linearise a nonlinear system. Basically, the system involves a mixing tank of height H(t) and cross section area A with two input streams with flow rates Q₁(t) and Q₂(t) with concentrations C₁ and C₂ of some species (concentrations are constant). There is one outlet stream with flow rate Q(t) and concentration C(t). We are assuming that the tank is always perfectly mixed and that the outlet flow rate is proportional to \sqrt{H(t)}.

So what i need to do is express the concentration of the outlet stream as some nominal steady-state value plus an additional term for deviations from this nominal value

i.e. C(t) = C₀ + c(t)

Where c(t) is a linear approximation.


My first thought is to write a material balance on the species in the fluid:

i.e. \frac{d}{{dt}}\left( {CAH(t)} \right) = {C_1}{Q_1} + {C_2}{Q_2} - CQ

It wouldn't be too hard to solve this DE for C(t), but I am not sure where to go from there. I guess i could then differentiate C(t) with respect to each variable to linearise it, but this seems a little tedious so I am not sure if there is something else i should be doing?

Any input is greatly appreciated.

Thanks,
Dan.

Expand your nonlinear equation in a Taylor series and take the first two terms.

 

[danago]

Are you referring to the analytical solution of the differential equation? If i solve the DE (i used the integrating factor method), i get:

C(t) = \frac{{\int {({C_1}{Q_1} + {C_2}{Q_2}){e^{\int {\frac{k}{{A\sqrt H }}dt} }}dt + \alpha } }}{{AH{e^{\int {\frac{k}{{A\sqrt H }}dt} }}}}

where alpha is a constant of integration.

Im not sure if this is correct, and even if it is, it seems a little difficult to linearise, especially since i don't know the exact form of the flow rate functions Q₁ and Q₂

Thanks for the reply by the way :)

 

[CEL]

Are you referring to the analytical solution of the differential equation? If i solve the DE (i used the integrating factor method), i get:

C(t) = \frac{{\int {({C_1}{Q_1} + {C_2}{Q_2}){e^{\int {\frac{k}{{A\sqrt H }}dt} }}dt + \alpha } }}{{AH{e^{\int {\frac{k}{{A\sqrt H }}dt} }}}}

where alpha is a constant of integration.

Im not sure if this is correct, and even if it is, it seems a little difficult to linearise, especially since i don't know the exact form of the flow rate functions Q₁ and Q₂

Thanks for the reply by the way :)

Do you need an analytical solution, or a numerical solution will do?

If you can use a numerical solution, you can start from

\frac{d}{dt}C(t) = f(C, t)

and

C(t_0) = C_0

and calculate

C(T) = C_0 + \frac{d}{dt}C(t_0) T

then calculate C(2T) from C(T) and so on.

If you need an analytical solution, you can make

f(C, t) = f(C_0, t_0) + \frac{df}{dC}(C_0, t_0) \right( C - C_0 \left )
and integrate.

 

[danago]

Do you need an analytical solution, or a numerical solution will do?

If you can use a numerical solution, you can start from

\frac{d}{dt}C(t) = f(C, t)

and

C(t_0) = C_0

and calculate

C(T) = C_0 + \frac{d}{dt}C(t_0) T

then calculate C(2T) from C(T) and so on.

If you need an analytical solution, you can make

f(C, t) = f(C_0, t_0) + \frac{df}{dC}(C_0, t_0) \right( C - C_0 \left )
and integrate.

I think i understand what you are saying; using the product rule of differentiation with the species material balance gives me:

C\frac{{d(AH)}}{{dt}} + AH\frac{{dC}}{{dt}} = {C_1}{Q_1} + {C_2}{Q_2} - CQ

I can then express the derivative of AH as a sum of flow rates into and out of the tank:

C({Q_1} + {Q_2} - Q) + AH\frac{{dC}}{{dt}} = {C_1}{Q_1} + {C_2}{Q_2} - CQ

Rearranging this and substituting the relevant values for t=0 (denoted by a zero in the subscript) gives:

{\left. {\frac{{dC}}{{dt}}} \right|_{t = 0}} = \frac{{{C_{1}}{Q_{10}} + {C_{2}}{Q_{20}} - {C_0}({Q_{10}} + {Q_{20}})}}{{A{H_0}}}

I can then use this to give me the linearised form:

C(t) = {C_0} + {\left. {\frac{{dC}}{{dt}}} \right|_{t = 0}}t

Have i taken the right approach with this?

 

[CEL]

I think i understand what you are saying; using the product rule of differentiation with the species material balance gives me:

C\frac{{d(AH)}}{{dt}} + AH\frac{{dC}}{{dt}} = {C_1}{Q_1} + {C_2}{Q_2} - CQ

I can then express the derivative of AH as a sum of flow rates into and out of the tank:

C({Q_1} + {Q_2} - Q) + AH\frac{{dC}}{{dt}} = {C_1}{Q_1} + {C_2}{Q_2} - CQ

Rearranging this and substituting the relevant values for t=0 (denoted by a zero in the subscript) gives:

{\left. {\frac{{dC}}{{dt}}} \right|_{t = 0}} = \frac{{{C_{1}}{Q_{10}} + {C_{2}}{Q_{20}} - {C_0}({Q_{10}} + {Q_{20}})}}{{A{H_0}}}

I can then use this to give me the linearised form:

C(t) = {C_0} + {\left. {\frac{{dC}}{{dt}}} \right|_{t = 0}}t

Have i taken the right approach with this?

Yes, it is correct.

 

[danago]

Alright, thanks very much Your help is much appreciated!