Reactor Quı́mico.

1st Luis Carlos Bonilla Góngora 2nd Cristian Mateo Castro Perfetti
Universidad del Cauca Universidad del Cauca
Facultad de Ingenierı́a en Facultad de Ingenierı́a en
Electrónica y Telecomunicaciones Electrónica y Telecomunicaciones
Popayán, Cauca Popayán, Cauca
[email protected] [email protected]

Abstract—This document presents the development of a math- D. Obtención de matrices Jacobianas

ematical model for a simple chemical reactor, its linearization and
subsequent simulation using Matlab/Simulink tools. The system
Para continuar con el proceso de linealización es necesario
studied comprises an ideal reactor. The differential equations tener las anteriores ecuaciones (1) y (2) despejadas con solo
that describe the dynamic behavior of the system are detailed, el término con la derivada a un lado de la ecuación, obte-
considering mass and energy balances that may be presented in niendo las nuevas ecuaciones de balance de masa y energı́a
the provided model. Subsequently, the model is linearized around respectivamente;
four stable operating points through the Jacobian linearization
method, obtaining a representation in matrices. Finally, the
linearized model is implemented in Simulink. d qCAi (t) −E qCA (t)
Index Terms—Reactor Quı́mico, modelo, MATLAB, lineal- CA (t) = − Ko e RT (t)) CA (t) − (4)
ización, espacio de estados. dt V V

I. I NTRODUCTION −E
d qρCp Ti (t) Ko e RT (t)) CA (t)∆Hr
Jacobianos Balance de energı́a y masa derivada de expo- T (t) = − −
dt V ρCv ρCv
nencial qCp T (t)
− U AT (t) + U ATs − (5)
II. D ESARROLLO V Cv
A. Balance de Masa Después realizamos una factorización que nos permita agru-
par términos iguales, de la siguiente manera;
Para el desarrollo del ejercicio en general es de vital
importancia el desarrollo de las ecuaciones de masa, donde
se obtiene la siguiente ecuación; d qCAi (t) h −E qi
CA (t) = + CA (t) −Ko e RT (t)) − (6)
dt V V
d
V CA (t) = qCAi (t) − V Ra (t) − qCA (t) (1)
dt −E
B. Balance de Energı́a d qρCp Ti (t) Ko e RT (t)) CA (t)∆Hr
T (t) = − −
Igualmente es de importancia el desarrollo de las ecuaciones dt V ρCv ρCv
 
de balance de enrgı́a en el sistema, donde se obtiene la qCp
+ U ATs + T (t) −U A − (7)
siguiente ecuación; V Cv
Ahora debemos definir los puntos de linealización, que son
d aquellos valores de las variables de entrada y de salida en los
V ρCv T (t) = qρCp Ti (t)− cuales la planta va a estar trabajando normalmente.
dt
− V Ra (t)∆Hr − U A [T (t) − Ts ] − qρCp T( t) (2)
C¯A C¯
   
, Ai
 
x̄ ū → (8)
C. Fórmula de la Reacción T̄ T̄i
Resulta de en este tipo de plantas o procesos, el reactor Ya que tenemos las funciones y los puntos de linealización
contiene mezclas de componentes o quı́micos, los cuales definidos, es necesario armar cada matriz para después de-
pueden generar calor al entrar en contacto entre sı́; por lo tanto sarrollarla individualmente. Como sabemos, la derivada de un
se provee una ecuación que describe dicha energı́a y masa que estado se compone de un componente A multiplicado por el
se produce en la mezcla, la cual se ilustra a continuación; estado sin derivar más la entrada U por la componente B; en
−E
este caso A y B serı́an matrices de 2x2 y resultarı́an de la
Ra (t) = Ko e RT (t)) CA (t) (3) siguiente forma;
 " # Reemplazando dichas variables las ecuaciones resultarı́an
dF1 dF1
dCA dT de la siguiente forma;
A= dF2 dF2 (9)
dCA dT h qi hqi
" dF1 dF1
# C˙A = C¯A −C2 − + T̄ [−C1 ] + C¯Ai (23)
v v
dCAi dTi
B= dF2 dF2 (10)
dCAi dTi      
¯ ∆Hr ∆Hr ρCp
Ahora se procederá a detallar cada desarrollo de las Ṫ = CA −C2 + T̄ −C1 + T̄i (24)
ρCv ρCv V Cv
derivadas parciales mostradas anteriormente;
Por ultimo obtenemos las funciones de transferencia, despe-
dF1 q jando la variable de interés sin derivar y volviendo el coefi-
= −Ko e−E/RT̄ − (11)
dCA v ciente que la acompaña en 1, de tal forma que el coeficiente
dF1 E que acompaña a la variable de interés derivada es el tao de
= −Ko C¯A e−E/RT̄ ¯2 (12) dicha función y los coeficientes de las demás variables serı́an
dTi RT
las ganancias, de la siguiente manera para el balance de masa;
dF2 −Ko e−E/RT̄ ∆HR h qi hqi ¯
= (13) C˙A + C¯A C2 − = CAi − C1 T̄ (25)
dCA ρCv v v
dF2 −Ko C¯A ∆HR −E/RT̄ E
= e (14) 1 q C1
dT ρCv RT¯2 C˙A + C¯A = C¯Ai − C̄ (26)
q q
C2 + v vC2 + q C2 + v
dF1 q
= (15) Para esta función;
dCAi V
dF1 1
=0 (16) τ1 = q (27)
dTi C2 + v

dF2 q q
= (17) K1 = (28)
dCAi V vC2 + q

dF2 qCp C1
= (18) K2 = − q (29)
dTi V Cv C2 + v
Ahora mediante multiplicación de matrices por el escalar Obteniendo la siguiente función de transferencia;
que en este caso es el estado sin derivar y la entrada,
obtenemos dos nuevas funciones como a continuación; k1 k2
G1 (s) = + (30)
τ1 s + 1 τ1 s + 1
h qi Y para el balance de energı́a;
C˙A = C¯A −Ko e−E/RT̄ − +
 v   
E qCp      
+ T̄ −Ko C¯A e−E/RT̄ ¯2 + C¯Ai (19) Ṫ + T̄ C1
∆Hr
= C¯A −C2
∆Hr
+ T̄i
qCp
(31)
RT v ρCv ρCv V Cv
" −E # ρCv −C2 ρ2 Cp
Ko e RT (t)) ∆H
r Ṫ + T̄ = C̄A + T̄i (32)
Ṫ = C¯A − + C1 ∆Hr C1 vC1 ∆Hr
ρCv
" −E #   Para esta función;
Ko e RT (t)) CA (t)∆Hr E qCp
+ T̄ − + T̄i (20) ρCv
ρCv RT¯2 V Cv τ2 = (33)
−C1 ∆Hr
Las ecuaciones indicadas anteriormente son las ecuaciones
linealizadas, a las cuales se les va a realizar un reemplazo por −C2
K3 = (34)
otras variables para que al momento de realizar las operaciones C1
siguientes, no se tengan mayores problemas, donde dichas
variables quedan ası́; ρ2 C p
K4 = (35)
vC1 ∆Hr
Ko E C¯A −E/RT̄
C1 = e (21) Obteniendo la siguiente función de transferencia;
RT¯2
k3 k4
C2 = Ko e−E/RT̄ (22) G2 (s) = +
τ2 s + 1 τ2 s + 1
(36)
E. Simulación en Simulink/MATLAB del sistema sin linealizar como estabilidad, tiempo de respuesta, sobreimpulso y otros
En la siguiente imagen se puede ver como el sistema parámetros relevantes para el diseño y control del reactor.
presenta un corto tiempo de estabilización pero al terminarlo, IV. R ESULTADOS
la señal tiende a un valor contante
-

F. Simulación en Simulink/MATLAB del sistema linealizado

Teniendo en cuenta el sistema ya linealizado, se procede a
simularlo y se obtiene la siguiente gráfica, que nos permite
ver que desde el tiempo = 0 el valor es constante y no tiene
tiempos de estabilización, ya que los puntos de linealización
fueron aplicados correctamente.

III. C ONCLUSIONES
- El informe demuestra que el modelo matemático desarrol-
lado representa adecuadamente el comportamiento dinámico
del reactor sencillo bajo condiciones especı́ficas. La lineal-
ización en un punto de operación permite simplificar el análisis
y control del sistema. - Es importante destacar que la lineal-
ización es una aproximación válida solo en un rango limitado
alrededor del punto de operación. Para grandes perturbaciones,
el modelo linealizado puede no ser preciso y se requerirı́a un
modelo no lineal completo. - La simulación en MATLAB fa-
cilita el análisis de los resultados, permitiendo evaluar métricas