Tarea 2

Dinámica de Procesos Químicos

Septiembre del 2022

1
Problema
Se dispone de un reactor continuo perfectamente agitado en el que ocurre la reacción A → B ,
cuya cinética se describe por la siguiente ecuación:
k 1C A
−r A =
1+ k 2 C A
3
m
Con k 1=1.0 h−1 y k 2=1.0 .
mol
La reacción ocurre en condiciones isotérmicas y volumen constante ( 2 m3) y las propiedades físicas
3
m
se pueden supones constantes. En estado estacionario el flujo de entrada es 1 y la
h
mol
concentración de A en la entrada es 0.5 3 .
m
a) Formule el modelo dinámico del proceso.
b) Clasifique las variables del problema y determine los grados de libertad.
c) Determine la trayectoria dinámica de la concentración de A ante i) un cambio escalón en
mol m
3
la concentración de entrada a 1 3 y ii) un cambio escalón a 1.5 en el flujo de
m h
entrada.
d) Presente las funciones de transferencia para este proceso. Determine su orden dinámico y
los parámetros característicos. ¿Los parámetros son invariantes? Compárelos en el estado
estacionario inicial vs. el estado estacionario final. Justifique su respuesta.
e) Calcule analíticamente la respuesta ante las perturbaciones del punto c utilizando las
funciones de transferencia y la transformada inversa de Laplace. Compare con el modelo
no lineal (gráficamente) y discuta los resultados.

2
 Supuestos
q 0 ,C A 0 , ρ0 1. El reactor está
perfectamente agitado.
2. Rxn. ocurre en condiciones
isotérmicas.
3. El volumen es constante.
4. Propiedades físicas
constantes.

q ,C A , ρ

a) Formule el modelo dinámico del proceso.

Balance total de materia:
d (ρV ) ⇒
d (V ) d (V )
⇒
= ρ0 q0 ( t ) −ρq ( t ) por ( 4 ) ρ =ρ q 0 ( t )− ρq (t ) por ( 3 ) = ρq 0 ( t ) −ρq ( t )=0
dt dt dt
⇒
❑ q0 ( t ) =q ( t ) (1)

Balance de materia del componente A:

d ( V C A) ⇒ d (C A ) ⇒ d (C A )
=C A 0 q 0−C A q+r A V por ( 4 ) y ec . ( ¿ ) V =C A 0 q−C A q +r A V ❑ V =q (C ¿ ¿ A 0−C A
dt dt dt
⇒ d (C A ) q k1 C A
❑ = (C ¿ ¿ A 0−C A )− (2)¿
dt V 1+k 2 C A
- Se observa que la ec.(2) es una función NO LINEAL, esto debido a que está la
multiplicación de los términos q C A 0 y q C A .

b) Clasifique las variables del problema y determine los grados de libertad.

Parámetros: k 1 , k 2 , ρ0 , ρ

Variables: q 0 ,q ,C A , C A 0 ,V

Manipulables: q

Perturbaciones: q 0 ,C A 0 , V

3
 Respuestas: C A

Ecuaciones: 2

G.L = N V −N E =5−2=3 g .l

- Como se tienen 3 grados de libertad, se deben especificar 3 variables. Estas las da el

3
m mol
enunciado y son: V ( 2 m ) , q 0 (1
3
), C A 0 (0.5 3 ).
h m

c) Determine la trayectoria dinámica de la concentración de A ante i) un cambio escalón en

mol m
3
la concentración de entrada a 1 3 y ii) un cambio escalón a 1.5 en el flujo de
m h
entrada.

Para determinar la trayectoria dinámica de la concentración de A, se calcula la respuesta

estacionaria C AS. El régimen permanente inicial del reactor se obtiene para las cargas fijas
al proceso, las cuales son V , q s , C A 0 S .

qs qs k 1 C As
0=
V
( C A 0 s−C As )−r s= ( C A 0 s−C As )−
V 1+k 2 C As

Utilizando el programa Matlab con el archivo resolvedor_CSTR_estac.m se obtiene un

mol
valor de C As=0.1861 3 . Luego de obtener el valor de la concentración en estado
m
estacionario, se pude resolver el sistema dinámico de manera NO LINEAL para los distintos
casos, igualmente con el programa Matlab, pero con el archivo CSTR_c.m (los detalles
están dentro de este).

Caso 1: Sistema sin perturbaciones en la entrada.

d ( CA ) q0 k1 C A
Ecuación para resolver del sistema dinámico: = (C ¿ ¿ A 0−C A )− ¿
dt V 1+ k 2 C A

La figura 1 muestra el comportamiento del sistema sin perturbaciones a través del

tiempo.

4
Figura 1. Gráfico concentración de A sin perturbaciones en la entrada.

5
 mol
Caso 2: Sistema con cambio escalón en la concentración de entrada de 1 3 .
m
d ( CA ) q0 k1 C A
Ecuación para resolver del sistema dinámico: = (C ¿ ¿ A 1−C A )− ¿
dt V 1+k 2 C A
Con C A 1 la nueva concentración por perturbación en la entrada. La figura 2 muestra el
comportamiento del sistema con perturbación en la concentración de entrada en el
reactor a través del tiempo.

Figura 2. Gráfico concentración de A con perturbación en C de entrada.

3
m
Caso 3: Sistema con cambio escalón en flujo de entrada de 1.5 .
h
d ( CA ) q1 k1 C A
Ecuación para resolver del sistema dinámico: = (C ¿ ¿ A 0−C A )− ¿
dt V 1+ k 2 C A
Con q 1 el nuevo flujo por perturbación en la entrada. La figura 3 muestra el
comportamiento del sistema con perturbación en el flujo de entrada en el reactor a través
del tiempo.

6
 Figura 3. Gráfico concentración de A con perturbación en q de entrada.

d) Presente las funciones de transferencia para este proceso. Determine su orden dinámico y
los parámetros característicos. ¿Los parámetros son invariantes? Compárelos en el estado
estacionario inicial vs. el estado estacionario final. Justifique su respuesta.

El primer paso es linealizar la ecuación del sistema dinámico de la siguiente manera:

d (C A) ⇒ d ( C' A )
V =q0 (C ¿ ¿ A 0−C A )−r A V ❑ V =(q 0 C ¿ ¿ A 0−q 0 s C A 0 s)−(q 0 C ¿ ¿ A−q0 s C As)−(r ¿ ¿ A−
dt dt
¿

( ¿)
' ' ' dr
( q 0 C A ) =C As q ' 0+ q0 s C ' A ; ( q 0 C A 0 ) =C A 0 s q ' 0 +q 0 s C ' A 0 ; r =r A−r As= d C |
A s
( C A −C As) =
k1 C ' A
(1+k 2 C As)
2

⇒ d ( C' A ) k1 C ' A
❑V =C A 0 s q ' 0 + q0 s C ' A 0−C As q' 0−q 0 s C' A −
dt (1+k 2 C As )2

d ( C' A )
( )
⇒ k1 V
❑V =C A 0 s q' 0 +q 0 s C' A 0 −C As q' 0−C ' A q 0 s +
dt ( 1+k 2 C As )
2

d ( C' A )
( )
⇒ k 1V
❑V +C ' A q0 s + =q '0 (C A 0 s−C As)+q 0 s C ' A 0
dt ( 1+k 2 C As)
2

7
 ⇒
V d ( C' A ) (C A 0 s−C As) q0 s
❑ +C ' A=q' 0 + C 'A 0

( ) ( )( )
k1V dt k1V k1 V (3)
q0 s+ 2
q0 s+ 2
q0 s + 2
( 1+ k 2 C As ) ( 1+ k 2 C As ) ( 1+ k 2 C As )

Aplicamos la transformada de Laplace a la ecuación (3) y se obtiene lo siguiente:

V ' ' ' ( C A 0 s−C As ) q0s '

s C A (s)+C A (s )=Q 0 (s ) + C A 0 (s)

( ) ( )( )
k1 V k 1V k 1V
q0s+ 2
q0 s + 2
q0 s + 2
( 1+k 2 C As ) ( 1+k 2 C As) ( 1+k 2 C As )

Se definen los parámetros característicos como:

V ( C A 0 s −C As ) q0s
τ= ; K 1= ; K2=

( ) ( ) ( )
k1 V k1 V k1V
q0s+ 2
q0s+ 2
q0s+ 2
( 1+k 2 C As ) ( 1+ k 2 C As ) ( 1+ k 2 C As )

Con K 1 y K 2 las ganancias, la ecuación queda expresada como sigue:

' ' ' '

τ s C A ( s ) +C A ( s )=Q 0 ( s ) K 1 + K 2 C A 0 ( s )

' '
⇒
' Q 0 ( s) K 1 K 2 C A 0 ( s )
❑ C A ( s )= + (4)
τ s+1 τ s+ 1

K1 K2
Siendo B1= y B 2= las funciones de transferencia, entonces:
τ s +1 τ s +1
⇒
' ' '
❑ C A ( s )=Q 0 ( s ) B1 +C A 0 ( s ) B2 (5)

- En la ec.(4) se puede ver que el sistema es de primer orden, esto se debe a que la
ecuación presenta un polinomio de grado 1 en el denominador.

Variabilidad de los parámetros

Para determinar la variabilidad de los parámetros, se utiliza el archivo de Matlab
resolvedor_CSTR_estac.m para determinar los nuevos estados estacionarios cambiando
las variables en la entrada del reactor que presentan una perturbación y así comparar los
estados estacionarios iniciales con los finales. Los valores obtenidos se presentan en la
siguiente tabla:
Parámetros E.E in. E.E f. 1 E.E f. 2

8
 mol 0.5 1 0.5
C A0( 3
)
m
m
3
1 1 1.5
q0( )
h
mol 0.1861 0.4142 0.2410
C As( 3 )
m
−1
τ (h ) 0.8259 1 0.7146
K1 0.1296 0.2929 0.0925
K2 0.4129 0.5 0.5360

- Como se aprecia en la tabla, los estados estacionarios van variando según cambian los
parámetros. Esto nos dice que la ecuación NO ES LINEAL.

Diagrama de bloques para el sistema linealizado:

e) Calcule analíticamente la respuesta ante las perturbaciones del punto c utilizando las
funciones de transferencia y la transformada inversa de Laplace. Compare con el modelo
no lineal (gráficamente) y discuta los resultados.

Las graficas que se muestran fueron obtenidas a través de Matlab con el archivo
grafNLvsL.m.

Caso 2: Para el caso 2, en el que hay un cambio escalón en la concentración de entrada de

mol
1 3 , se tiene la siguiente expresión junto a la función escalón:
m
⇒
0.5
La anti transformada L−1 C ' A 0 ( s ) =
C ' A 0 (t)=0.5 H (t) s
¿
¿
'
Reemplazando en la ec.(4) y eliminando el término de Q 0 ( s ) debido a que no hay cambio
en este:
Q' 0 ( s ) K 1 K 2∗0.5
C ' A ( s )= +
τ s+1 s (τ s+1)

9
 0

( )
⇒ −t
' K 2∗0.5 −1 '
C A ( s )= La anti transformada L C A ( t )=C A −C As= K 2∗0.5 1−e τ
s (τ s +1)

( )+C
⇒ −t
τ
❑ C A =K 2∗0.5 1−e As (6)

En la figura 4, se presenta el gráfico de comparación entre el sistema NO LINEAL y el

sistema linealizado:

Figura 4. Gráfica comparativa entre ec. no lineal y lineal para caso 2.

3
m
Caso 3: Para el caso 2, en el que hay un cambio escalón en flujo de entrada de 1.5 , se
h
tiene la siguiente expresión junto a la función escalón:

⇒
0.5
' La anti transformada L−1 C ' A 0 ( s ) =
C A 0 (t)=0.5 H (t) s
¿
¿
'
Reemplazando en la ec.(4) y eliminando el término de C A 0 ( s ) debido a que no hay cambio
en este:
'
' 0.5∗K 1 K 2 C A 0 ( s )
C A ( s )= +
s (τ s +1) τ s+1
0

10
 =K ∗0.5 (1−e )+C
⇒ −t
τ
❑C A 1 As (7)

En la figura 5, se presenta el gráfico de comparación entre el sistema NO LINEAL y el

sistema linealizado:

Figura 5. Gráfica comparativa entre ec. no lineal y lineal para caso 3.

- Se puede concluir de los gráficos que el sistema LINEALIZADO se aproxima bien al

sistema NO LINEAL a puntos cercanos al estado estacionario inicial, pero a medida que
transcurre el tiempo se van separando las curvas una de otra. Esto se debe a que el
sistema NO LINEAL puede presentar más oscilaciones con distintos estados
estacionarios y la LINEALIZACIÓN no los toma todos.

11