Respuesta de una reactor químico a un cambio de concentración de entrada

Por: publicado el Oct 6, 2022

Problema 5.5 (Process Systems Analysis and Control - Coughanowr, LeBlanc)

gráfica problema 5.5
Considere el tanque agitado mostrado en la figura. La reacción que ocurrre es:

ABA \to B\\

Con una velocidad de reacción igual a:

r=kC0r=kC_0

Donde

r = (mol A)/(volumen)/(tiempo)k = constante de velocidad de reaccioˊnC0(t) = concentracioˊn de A en el reactor en el tiempo t (mol A)/(volumen)V = volumen de la mezcla en el reactorF = caudal de alimentacioˊn constante (volumen)/(tiempo)Ci(t)concentracioˊn de A en la entrada (mol A)/(volumen)\begin{array}{l} r \text{ = (mol A)/(volumen)/(tiempo)}\\ k \text{ = constante de velocidad de reacción}\\ C_0(t)\text{ = concentración de A en el reactor en el tiempo t } (mol\space A)/(volumen)\\ V\text{ = volumen de la mezcla en el reactor}\\ F\text{ = caudal de alimentación constante (volumen)/(tiempo)}\\ C_i(t)\text{concentración de A en la entrada (mol A)/(volumen)} \end{array}

Asumiendo densidad y volumen constante V, derive la función de tranferencia, relacionando la concentración en el reactor y la concentración de entrada. Dibuje la respuesta del reactor para un cambio tipo paso unitario en la concentración de entrada.

Resolviendo

Escribiendo nuestro balance de materia, sabiendo que no=C0Vn_o = C_0V

CiFC0FkC0V=dnodt=d(C0V)dtC_iF-C_0F-kC_0V=\frac{dn_o}{dt} = \frac{d(C_0V)}{dt}

CiFC0FkC0V=Vd(C0)dt     (1)C_iF-C_0F-kC_0V=V \frac{d(C_0)}{dt} \space\space\space\space\space\textbf{(1)}

Realizando el balance en estado estacionario

CisFC0sFkC0sV=0     (2)C_{is}F-C_{0s}F-kC_{0s}V=0 \space\space\space\space\space\textbf{(2)}

Restado (2) de (1) y tranformando a variables desviación

CiFCisF(C0FC0sF)k(C0VC0sV)=Vd(C0C0s)dtC_iF-C_{is}F-(C_0F-C_{0s}F)-k(C_0V-C_{0s}V)=V \frac{d(C_0-C_{0s})}{dt}

CiFC0FkC0V=Vd(C0)dtC'_iF-C'_0F-kC'_0V=V \frac{d(C'_0)}{dt}

Aplicando la transformada de Laplace y despejando y sabiendo que C0(t=0)=0C'_0(t=0) = 0

Ci(s)FC0(s)FkC0(s)V=V(sC0(s)C0(t=0))C'_i(s)F-C'_0(s)F-kC'_0(s)V=V (sC'_0(s)-C'_0(t=0))

Ci(s)FC0(s)FkC0(s)V=VsC0(s)C'_i(s)F-C'_0(s)F-kC'_0(s)V=V sC'_0(s)

Obteniendo nuestra función transferencia

C0(s)Ci(s)=FVs+F+kV\mathbf{\frac{C'_0(s)}{C'_i(s)}=\frac{F}{Vs+F+kV}}

Para poder hacer la gráfica con la variación de la concentración de entrada, reordenemos nuestra función.

C0(s)Ci(s)=F/(F+kV)Vs/(F+kV)+1\frac{C'_0(s)}{C'_i(s)}=\frac{F/(F+kV)}{Vs/(F+kV)+1}

Haciendo un cambio de variable

Kp=F/(F+kV)τ=V/(F+kV)K_p = F/(F+kV)\\ \tau = V/(F+kV)\\

C0(s)Ci(s)=Kpτs+1\frac{C'_0(s)}{C'_i(s)}=\frac{K_p}{\tau s+1}

Para un cambio en la concentración de entrada tipo paso unitario Con A como una constante cualquiera Ci=A/sC'_i=A/s.

C0(s)=AsKpτs+1C'_0(s)=\frac{A}{s}\frac{K_p}{\tau s+1}

Reordenando para realizar la antitransformada

C0(s)=AKp+AKpτsAKpτss(τs+1)C'_0(s)=\frac{A\cdot K_p+A\cdot K_p \cdot\tau s-A\cdot K_p \cdot\tau s}{s(\tau s+1)}

C0(s)=AKpsAKpττs+1C'_0(s)=\frac{A\cdot K_p}{s}-\frac{A\cdot K_p \cdot\tau}{\tau s+1}

C0(s)=AKpsAKps+1/τC'_0(s)=\frac{A\cdot K_p}{s}-\frac{A\cdot K_p}{s+1/\tau}

Antitransformando

C0(t)=AKp(1et/τ)C'_0(t) = A\cdot K_p(1-e^{-t/\tau})

Graficando esta respuesta

respuesta del sistema a un cambio en la concentración de entrada

Referencias

  • Coughanowr, D. R.; LeBlanc, S. E. (2009). Process Systems Analysis and Control (3rd edition). McGraw-Hill. ISBN 978-0-07-339789-4.