Guía 1 - Ejercicio 3: Tanques Mezcladores

Ejercicio 3

Planteo: Tanques A y B con volumen VV, llenos con agua en t0=0t_0=0 / t>0t>0, solución de volumen vv y masa mm fluye a A por seg / mezcla fluye A →\to B a la misma tasa / mezcla fluye de B a la misma tasa.

dσAdt+vVσA=mV,dσBdt+vVσB=vVσA  ⟹  mB(t)=mVv(1−e−vt/V)−mte−vt/V\frac{d\sigma_A}{dt} + \frac{v}{V}\sigma_A = \frac{m}{V} \quad , \quad \frac{d\sigma_B}{dt} + \frac{v}{V}\sigma_B = \frac{v}{V}\sigma_A \implies m_B(t) = \frac{mV}{v} \left(1 - e^{-vt/V}\right) - mte^{-vt/V}

Resolución para el tanque A:

σA′+vVσA=mV\sigma_A' + \frac{v}{V}\sigma_A = \frac{m}{V} σA′=mV−vVσA\sigma_A' = \frac{m}{V} - \frac{v}{V}\sigma_A dσAdt=mV−vVσA\frac{d\sigma_A}{dt} = \frac{m}{V} - \frac{v}{V}\sigma_A ∫dσAmV−vVσA=∫dt\int \frac{d\sigma_A}{\frac{m}{V} - \frac{v}{V}\sigma_A} = \int dt −Vv∫−vVdσAmV−vVσA=t+c-\frac{V}{v} \int \frac{-\frac{v}{V} d\sigma_A}{\frac{m}{V} - \frac{v}{V}\sigma_A} = t + c −Vvln⁡(σA−vV)=t+c  ⟹  ln⁡(...)=−vVt+c~-\frac{V}{v} \ln\left(\sigma_A - \frac{v}{V}\right) = t + c \implies \ln(...) = -\frac{v}{V}t + \tilde{c} σA(t)−mv=Ke−vt/V\sigma_A(t) - \frac{m}{v} = K e^{-vt/V} σA(t)=Ke−vt/V+m/v\sigma_A(t) = K e^{-vt/V} + m/v

Condición inicial: σA(0)=0=K+m/v  ⟹  K=−m/v\sigma_A(0) = 0 = K + m/v \implies K = -m/v

  ⟹  σA(t)=mv(1−e−vt/V)\implies \sigma_A(t) = \frac{m}{v} \left(1 - e^{-vt/V}\right)

Resolución para el tanque B:

  ⟹  σB′+vVσB=vVmv(1−e−vt/V)=mV(1−e−vt/V)\implies \sigma_B' + \frac{v}{V}\sigma_B = \frac{v}{V} \frac{m}{v} \left(1 - e^{-vt/V}\right) = \frac{m}{V}\left(1 - e^{-vt/V}\right)

Recordatorio de factor integrante usado: x′+P(t)x=Q(t)x' + P(t)x = Q(t), μ(t)=e∫P(t)dt\mu(t) = e^{\int P(t) dt}   ⟹  μx′+μPx=μQ  ⟹  (μx)′=μQ  ⟹  x(t)=1μ(t)(∫μ(t)Q(t)dt+C)\implies \mu x' + \mu P x = \mu Q \implies (\mu x)' = \mu Q \implies x(t) = \frac{1}{\mu(t)} \left( \int \mu(t)Q(t) dt + C \right)

  ⟹  μ(t)=e∫vVdt=evt/V\implies \mu(t) = e^{\int \frac{v}{V} dt} = e^{vt/V} σB(t)=1evt/V(∫evt/VmV(1−e−vt/V)dt+c)\sigma_B(t) = \frac{1}{e^{vt/V}} \left( \int e^{vt/V} \frac{m}{V} \left(1 - e^{-vt/V}\right) dt + c \right) =e−vt/V(∫mVevt/Vdt−∫mVdt)= e^{-vt/V} \left( \int \frac{m}{V} e^{vt/V} dt - \int \frac{m}{V} dt \right) =e−vt/Vc+e−vt/V(mVVvevt/V−mVt)= e^{-vt/V} c + e^{-vt/V} \left( \frac{m}{V} \frac{V}{v} e^{vt/V} - \frac{m}{V} t \right) =c~e−vt/V+mv−mVte−vt/V= \tilde{c} e^{-vt/V} + \frac{m}{v} - \frac{m}{V} t e^{-vt/V}

Condición inicial: σB(0)=0=c~+mv−0  ⟹  c~=−m/v\sigma_B(0) = 0 = \tilde{c} + \frac{m}{v} - 0 \implies \tilde{c} = -m/v

  ⟹  σB(t)=−mve−vt/V+mv−mVte−vt/V\implies \sigma_B(t) = -\frac{m}{v} e^{-vt/V} + \frac{m}{v} - \frac{m}{V} t e^{-vt/V}

Multiplicando por el volumen VV para hallar la masa:

  ⟹  mB(t)=VσB(t)=−Vmve−vt/V+Vmv−VmVte−vt/V\implies m_B(t) = V \sigma_B(t) = -V \frac{m}{v} e^{-vt/V} + V \frac{m}{v} - \frac{Vm}{V} t e^{-vt/V}   ⟹  mB(t)=mVv(1−e−vt/V)−mte−vt/V\implies m_B(t) = \frac{mV}{v} \left(1 - e^{-vt/V}\right) - mte^{-vt/V}

Mi resolución

Objetivo: Demostrar que la masa del soluto en el tanque B, dadas las ecuaciones de concentración y las condiciones iniciales, está dada por:

MB(t)=mVv(1−e−vt/V)−mte−vt/VM_B(t) = \frac{mV}{v}\left(1 - e^{-vt/V}\right) - mte^{-vt/V}

1. Modelado y Resolución del Tanque A

La dinámica de la concentración en el tanque A está gobernada por la siguiente Ecuación Diferencial Ordinaria (EDO) lineal de primer orden:

dσAdt+vVσA=mV\frac{d\sigma_A}{dt} + \frac{v}{V}\sigma_A = \frac{m}{V}

Dado que el tanque comienza lleno de agua pura en t=0t=0, establecemos nuestro Problema de Valor Inicial (PVI) con la condición:

σA(0)=0\sigma_A(0) = 0

Método de Resolución (Factor Integrante):

Buscamos un factor integrante μ(t)=e∫P(t)dt\mu(t) = e^{\int P(t)dt}, donde P(t)=vVP(t) = \frac{v}{V}.

μ(t)=evVt\mu(t) = e^{\frac{v}{V}t}

Multiplicamos toda la EDO por μ(t)\mu(t) para transformar el miembro izquierdo en la derivada de un producto:

evVtdσAdt+evVtvVσA=evVtmVe^{\frac{v}{V}t}\frac{d\sigma_A}{dt} + e^{\frac{v}{V}t}\frac{v}{V}\sigma_A = e^{\frac{v}{V}t}\frac{m}{V} ddt(σA⋅evVt)=mVevVt\frac{d}{dt}\left(\sigma_A \cdot e^{\frac{v}{V}t}\right) = \frac{m}{V}e^{\frac{v}{V}t}

Integramos ambos miembros con respecto al tiempo tt:

σA⋅evVt=∫mVevVtdt\sigma_A \cdot e^{\frac{v}{V}t} = \int \frac{m}{V}e^{\frac{v}{V}t} dt σA⋅evVt=mV(VvevVt)+C=mvevVt+C\sigma_A \cdot e^{\frac{v}{V}t} = \frac{m}{V} \left(\frac{V}{v}e^{\frac{v}{V}t}\right) + C = \frac{m}{v}e^{\frac{v}{V}t} + C

Despejamos la concentración σA(t)\sigma_A(t):

σA(t)=mv+Ce−vVt\sigma_A(t) = \frac{m}{v} + C e^{-\frac{v}{V}t}

Aplicamos la condición inicial σA(0)=0\sigma_A(0) = 0 para hallar la constante CC:

0=mv+Ce0  ⟹  C=−mv0 = \frac{m}{v} + C e^0 \implies C = -\frac{m}{v}

Por lo tanto, la función de concentración para el Tanque A es:

σA(t)=mv(1−e−vVt)\sigma_A(t) = \frac{m}{v}\left(1 - e^{-\frac{v}{V}t}\right)

2. Modelado y Resolución del Tanque B

La EDO para la concentración en el tanque B acopla la salida del tanque A como su término fuente:

dσBdt+vVσB=vVσA\frac{d\sigma_B}{dt} + \frac{v}{V}\sigma_B = \frac{v}{V}\sigma_A

Sustituimos la solución hallada σA(t)\sigma_A(t) en la EDO de B:

dσBdt+vVσB=vV[mv(1−e−vVt)]\frac{d\sigma_B}{dt} + \frac{v}{V}\sigma_B = \frac{v}{V} \left[ \frac{m}{v}\left(1 - e^{-\frac{v}{V}t}\right) \right] dσBdt+vVσB=mV−mVe−vVt\frac{d\sigma_B}{dt} + \frac{v}{V}\sigma_B = \frac{m}{V} - \frac{m}{V}e^{-\frac{v}{V}t}

Nuevamente, el tanque arranca con agua pura, por lo que la condición inicial es:

σB(0)=0\sigma_B(0) = 0

Resolución:

Como el lado izquierdo de la ecuación es idéntico al del tanque A, el factor integrante es el mismo: μ(t)=evVt\mu(t) = e^{\frac{v}{V}t}. Multiplicamos toda la ecuación:

evVtdσBdt+evVtvVσB=mVevVt−mVe−vVtevVte^{\frac{v}{V}t}\frac{d\sigma_B}{dt} + e^{\frac{v}{V}t}\frac{v}{V}\sigma_B = \frac{m}{V}e^{\frac{v}{V}t} - \frac{m}{V}e^{-\frac{v}{V}t}e^{\frac{v}{V}t} ddt(σB⋅evVt)=mVevVt−mV\frac{d}{dt}\left(\sigma_B \cdot e^{\frac{v}{V}t}\right) = \frac{m}{V}e^{\frac{v}{V}t} - \frac{m}{V}

Integramos respecto a tt:

σB⋅evVt=∫(mVevVt−mV)dt\sigma_B \cdot e^{\frac{v}{V}t} = \int \left( \frac{m}{V}e^{\frac{v}{V}t} - \frac{m}{V} \right) dt σB⋅evVt=mvevVt−mVt+K\sigma_B \cdot e^{\frac{v}{V}t} = \frac{m}{v}e^{\frac{v}{V}t} - \frac{m}{V}t + K

Despejamos σB(t)\sigma_B(t):

σB(t)=mv−mVte−vVt+Ke−vVt\sigma_B(t) = \frac{m}{v} - \frac{m}{V}t e^{-\frac{v}{V}t} + K e^{-\frac{v}{V}t}

Evaluamos con la condición inicial σB(0)=0\sigma_B(0) = 0 para hallar la constante KK:

0=mv−0+Ke0  ⟹  K=−mv0 = \frac{m}{v} - 0 + K e^0 \implies K = -\frac{m}{v}

Reemplazamos KK y agrupamos términos:

σB(t)=mv−mve−vVt−mVte−vVt\sigma_B(t) = \frac{m}{v} - \frac{m}{v}e^{-\frac{v}{V}t} - \frac{m}{V}t e^{-\frac{v}{V}t} σB(t)=mv(1−e−vVt)−mVte−vVt\sigma_B(t) = \frac{m}{v}\left(1 - e^{-\frac{v}{V}t}\right) - \frac{m}{V}t e^{-\frac{v}{V}t}

3. Demostración de la Masa en el Tanque B

La relación fundamental entre concentración, masa y volumen es:

σ(t)=M(t)V  ⟹  M(t)=σ(t)⋅V\sigma(t) = \frac{M(t)}{V} \implies M(t) = \sigma(t) \cdot V

Multiplicando nuestra expresión de σB(t)\sigma_B(t) por el volumen VV:

MB(t)=V[mv(1−e−vVt)−mVte−vVt]M_B(t) = V \left[ \frac{m}{v}\left(1 - e^{-\frac{v}{V}t}\right) - \frac{m}{V}t e^{-\frac{v}{V}t} \right]

Distribuyendo la constante VV, obtenemos la expresión final solicitada (Q.E.D.):

MB(t)=mVv(1−e−vt/V)−mte−vt/VM_B(t) = \frac{mV}{v}\left(1 - e^{-vt/V}\right) - mte^{-vt/V}

4. Análisis e Interpretación Física del Modelo

A. El Efecto de Retraso e Inercia

Podemos reescribir la concentración en B en función de la de A:

σB(t)=σA(t)−mVte−vVt\sigma_B(t) = \sigma_A(t) - \frac{m}{V}t e^{-\frac{v}{V}t}

El término restado es estrictamente positivo para todo t>0t > 0. Esto indica que σB(t)<σA(t)\sigma_B(t) < \sigma_A(t) durante todo el régimen transitorio. Físicamente, el tanque B exhibe una inercia frente al soluto: debe esperar a que la mezcla en el tanque A adquiera concentración antes de comenzar a recibir soluto significativamente.

B. Comportamiento Asintótico (Estado Estacionario)

Al evaluar el límite de ambas concentraciones cuando t→∞t \to \infty:

  • En el tanque A: lim⁡t→∞σA(t)=mv\lim_{t \to \infty} \sigma_A(t) = \frac{m}{v}

  • En el tanque B, debido a que la exponencial domina sobre la función lineal tt, el término inercial se anula: lim⁡t→∞σB(t)=mv\lim_{t \to \infty} \sigma_B(t) = \frac{m}{v}

Conclusión: A largo plazo (estado estacionario), el sistema "olvida" las condiciones iniciales (agua pura). El volumen total de agua es completamente lavado y reemplazado, alcanzando ambos tanques un equilibrio dinámico donde la concentración interna iguala a la concentración del flujo de entrada original (m/vm/v).