MÓDULO 13

REACTOR BATCH

Consideremos un reactor discontinuo mostrado en la figura 13.1. El reactor se carga con el reactante
y luego se sella. El contenido del reactor se calienta hasta la temperatura de reacción con vapor
saturado suministrado a la chaqueta a la temperatura TS . El vapor condensado se recoge en una
trampa de vapor. La mezcla de reacción se agita continuamente a una rpm fija. El reactor está
equipado con un serpentín de enfriamiento a través de la cual pasa agua enfriada y se elimina el calor
exotérmico para monitorear la temperatura del sistema a lo largo de una trayectoria predeterminada
temperatura-tiempo deseada. Para seguir la trayectoria prescrita, son necesarios tanto el
calentamiento como el enfriamiento de la unidad de proceso.
A guade enf r iam i ento( ent ra da)

contr ol d
e fl u
jo

cont rol de Ag ua de en fri amie nto(sal ida ),T c

presión

s umini stro
d ev apo r
serpe ntin de re fri gerac ión

Ag it ador
cha queta

v apor con den sad o

Figura 13.1 Representación esquemática para el ejemplo de reactor batch.

Considere que ocurren en el reactor las siguientes reacciones consecutivas:

A  k1
 B k2
C
Se asume que la reacción A  B tiene una cinética de segundo orden, mientras que la reacción
B  C tiene una cinética de primer orden. En esta operación batch, B es el producto deseado. Si la
reacción procede por un largo tiempo, la producción del producto indeseable C será alta. Si
detenemos la reacción demasiado pronto, la conversión de A puede ser muy baja. En ambas
situaciones, los rendimientos del producto deseado B son bajos. Por lo tanto, hay un tiempo óptimo
en el cual podemos parar la operación y obtener la máxima cantidad de producto deseado. Este
tiempo es conocido como batch time.

13.1 Modelo matemático

Suposiciones
Hacemos las siguientes suposiciones para desarrollar el modelo del proceso:
-Ambas reacciones químicas son exotérmicas.

XIII.1
-El reactor está bien aislado de modo que es despreciable el intercambiador de calor entre el reactor
y el ambiente externo. Algunas partes de las líneas de suministro de H 2O de enfriamiento y de
vapor, que se exponen al medio ambiente, también están aisladas.
-Las temperaturas de entrada y salida del refrigerante no varían mucho debido a la velocidad de
flujo suficientemente rápida de esta corriente. Por lo tanto en el modelo matemático se usa la
temperatura promedio de H O de enfriamiento TC .
2
-El vapor condensado es descargado a su temperatura de saturación (no hay subenfriamiento).
-Asumimos volumen constante del reactor (V), con capacidad calorífica media constante Cp. Tenga
en cuenta que el volumen se refiere al volumen de la mezcla de reacción, no al volumen del reactor.
-El contenido del reactor está perfectamente mezclado y no exhibe gradientes significativos de
concentraciones o temperatura, en ninguna parte del reactor.
U ,
-El coeficiente global de transferencia de calor para la chaqueta de calentamiento j se asume
constante; pero el coeficiente global de transferencia de calor para el refrigerante U C , se asume que
es función de la velocidad de flujo FC del refrigerante, como:
1 1 1
 0,8

U C 4550 FC 10,8

Desarrollo del modelo

El modelamiento del reactor se basa en un conjunto de ecuaciones de balance fundamentales. Estas

ecuaciones son declaraciones matemáticas de leyes físicas que requieren la conservación masa,
energía, momento y las cantidades de las especies químicas. La ecuación básica de balance es:

 Velocidad de   Velocidad de   Velocidad de flujo   Velocidad de generación 

 acumulación dentro  = flujo queingresa  - desalida desde el  +  por reacción química 
       
 delsistema   alsistema  sistema  dentro delsistema 

Ecuación de continuidad total

En la operación batch:
Velocidad de flujo neto de masa en el reactor =0
Velocidad de generación de masa dentro del reactor =0
d  V 

Velocidad de acumulación de masa dentro del reactor dt
Desde que  es la densidad de la mezcla y V es el volumen de la mezcla, V es la masa total
dentro del reactor. Sustituyendo equivalencias en la ecuación básica de balance, resulta
d  V 
0
dt (13.1)
Hemos asumido que V es constante y por lo tanto de acuerdo a la ecuación de continuidad,  ,
permanece invariable (constante).

Ecuaciones de continuidad de componentes

Para el componente A (reactante)
Velocidad de flujo neto del componente A en el reactor = 0

XIII.2
  rAV     rA  V
Velocidad de generación del componente A por reacción química
dN A

Velocidad de acumulación del componente A dentro del reactor dt

Comentario
Por convención  rA es la velocidad de desaparición de la especie A; consecuentemente,  A  es la
r

velocidad de formación de A. Aquí  A  representa el número de moles de la especie A en el

N
sistema en el tiempo t. Ahora la ecuación básica de balance toma la forma:
dN A
 0    rA  V
dt (13.2)
es decir:
1 dN A d  N A / V  dC A
      rA 
V dt dt dt (13.3)
C
donde A es la concentración de la especie A. Desde que la reacción A 
k1
 B tiene una cinética
de segundo orden, de modo que
rA  k1C A 2 (13.4)
k
donde 1 representa la constante de velocidad de reacción.
Sustituyendo la ecuación (13.4) en la ecuación (13.3) se obtiene:
dC A
 k1C A2
dt (13.5)
Esta es la ecuación de continuidad o ecuación de balance de masa del componente A.

Para el componente B (producto deseado)

Velocidad de flujo neto del componente B en el reactor  0

    rB  V
Velocidad de generación del componente B por reacción química
dN B

Velocidad de acumulación del componente B dentro del reactor dt
Comentario
N
Aquí, ( rB ) , es la velocidad de desaparición de la especie B; B , denota el número de moles de la
especie B en el sistema en el tiempo t. Reemplazando equivalencias en la ecuación de balance, se
tiene:
dN B
   rB  V
dt (13.6)
rB  k1C A 2  k2C B   rB  k2CB  k1C A2
1 dN B
    rB 
V dt
dC B
 k1C A2  k2 CB
dt (13.7)
Esta es la ecuación de continuidad o de balance de masa para la especie B.

XIII.3
Constante de velocidad de reacción
De acuerdo a la ley de velocidad de Arrhenius, la constante de velocidad de reacción, k (también
conocida como constante de velocidad específica de reacción) se expresa por
 E 
k   exp  
 RT  (13.8)
donde:  : Factor de frecuencia o factor pre-exponencial.
E : Energía de activación.
R : Constante universal de los gases.
La función exponencial dependiente de la temperatura en la ecuación cinética 13.8, representa una
de las severas no linealidades en los sistemas de ingeniería química. Para el ejemplo del reactor
batch, la ecuación de Arrhenius da las siguientes formas.
 E 
k1  10 exp  1 
 RT  (13.9)
 E 
k2   20 exp  2 
 RT  (13.10)
Las dos ecuaciones anteriores se utilizan en la ecuación de continuidad del componente y de
energía.

Ecuación de balance de energía

Para el reactor batch, la ecuación general de balance, puede expresarse como:

 Velocidad de   Velocidad neta   Velocidad de adición 

acumulación   de entrada de    de energía por 
     
de energía  energía   reacción exotérmica 
(13.11)
Entonces,
 Velocidad entrada   Velocidad de salida 
 Velocidad de     
 acumulación    de calor de chaqueta   de calor del 
   de calentamiento   reactor al agua de 
 de energía     
 al reactor   refrigeración 
 U j Aj  TS  T   U C AC  T  TC 

A A
Aquí, j y C son las áreas de transferencia de calentamiento de la chaqueta, y del sistema de
refrigeración respectivamente. Puesto que hemos asumido que el reactor está bien aislado, al igual
que las líneas de suministro de vapor y refrigerante; concluimos que no hay intercambio de calor
entre el reactor y el medio ambiente.

De nuevo, la velocidad de adición de energía por las reacciones exotérmicas, es:

 Velocidad deadición  C
 deenergía por 
   H i   rV
i    H1  k1C A V   H 2  k 2C BV
2

 reacción exotérmica  i 1

Comentario

XIII.4
No se debe confundir, k conocida como la velocidad de reacción o velocidad específica de
reacción, con r que es la velocidad de reacción.
Claramente, C  2. Los subíndices 1 y 2 se usan para las reacciones A   B y B  C
k1 k2

respectivamente. Aquí 
H1  y  H 2 
, son los calores exotérmicos de reacción para las dos
reacciones. Se conoce bien que el calor de reacción es negativo para la reacción exotérmica y
positivo para la reacción endotérmica. Ahora la:
d  Vh  dT
  VCP
Velocidad de acumulación de energía dt dt

Aquí,
 ,V , C
P son constantes; h es la entalpía C T P . Sustituyendo términos en la ecuación 13.11,
obtenemos la siguiente forma:
dT
VCP  U j Aj  TS  T   U C AC  T  TC    H 1  k1C A 2V   H 2  k 2CBV
dt (13.12)
Reordenando,
dT U j Aj U A  H1  k C 2 V   H 2  k C V
  TS  T   C C  T  TC   1 A 2 B
dt VC P VCP  V CP  V CP
dT  H1   H 2  k C  U j Aj T  T  U C AC T  T
 k1C A2  2 B  S   C
dt  CP CP VC P VCP (13.13)
Esta es la ecuación de balance de energía para el ejemplo del reactor batch.

Aplicación de un algoritmo de control.

El control de un sistema de reactor batch, es necesario porque el sistema no opera en estado
estacionario. De la figura 13.1 es claro que la velocidad de enfriamiento se controla mediante la
manipulación del flujo de agua ( FC ), a través del aceite de enfriamiento.
El control de la velocidad de calentamiento se ejecuta regulando la presión de vapor y por lo tanto
la temperatura del vapor, TS . Realmente el controlador manipula ambas variables, FC y TS
simultáneamente para regular la temperatura del reactor T.
Con el fin de alcanzar un máximo rendimiento del producto siempre es necesario crear una reacción
amigable con el medio ambiente en el reactor. Para el presente caso la temperatura del set point en
el reactor debe ser cambiada precisamente y la dinámica del set point está dada como una
trayectoria deseada temperatura- tiempo:
Td  t   54  71exp  2,5  103 t 
(13.14)
En la práctica general, una variable de salida (T) es controlada utilizando solamente la variable TS .
Si usamos otra variable FC o U C (coeficiente global del refrigerante) simultáneamente para manejar
T, el sistema de control estará determinado. Una forma de resolver esta dificultad es introducir una
variable paramétrica u definida de la siguiente manera:
TS   TS ,max  TS ,min  u  TS ,min
(13.15)
U C   U C ,min  U C ,max  u  U C ,max
(13.16)
donde los valores máximo y mínimo de TS y U C son elegidos en función del proceso o tomando en
cuenta límites de seguridad. Evidentemente u  0 representa el máximo enfriamiento del sistema y
u  1, el máximo calentamiento. Esta variable paramétrica es utilizada como una variable única a ser

XIII.5
manipulada en el problema de control batch indicado. Después de sustituir en la ecuación de
balance de energía las ecuaciones (13.15) y (13.16) obtenemos
dT
  1k1C A2   2 k2C B   a1  a2T    b1  b2T  u
dt (13.17)
donde

1 
 H1 
 CP (13.18)

2 
 H 2 
 CP (13.19)
 
a1  U j A jTS ,min  U C ,max AC TC  C PV
(13.20)

a2   U j Aj  U C ,max AC  CPV  (13.21)
b1  U j Aj  TS ,max  TS ,min    U C ,max  U C ,min  AC TC   CPV
(13.22)
b2   U C ,max  U C ,min  AC  CPV
(13.23)
Ahora u , puede obtenerse empleando un algoritmo de control. Un controlador clásico PI será
implementado para regular la temperatura en el reactor batch. La ecuación del controlador PI es:
 1 
u  t   uS  kC  e  t    e  t  dt 
 i  (13.24)
donde:
uS : Señal de polarización del control.
K C : Ganancia proporcional (ajuste de parámetro).
 i : Constante de la integral de tiempo.
e  t   Td  t   T  t 
: Error del controlador.
kC  0,1 ºC 1 
Los valores de los parámetros de ajuste son elegidos arbitrariamente como y
 i  360s

Condiciones de arranque o puesta en marcha y características del sistema

C A0 concentración de la especie A en el arranque, kmol/m 3 1,08

CB 0 concentración de la especie B en el arranque, kmol/m 3 0,00
10 factor de frecuencia de la reacción A  B ,m3/kmol s 1,1
 20 factor de frecuencia de la reacción B  C ,s-1 172,2
E1 energía de activación para la reacción A  B ,kJ/kmol 2,09×104
E2 energía de activación para la reacción B  C ,kJ/kmol 4,18×104

 H1  calor de reacción para la reacción A  B , kJ/kmol 4,18×104

 H 2  calor de reacción para la reacción B  C , kJ/kmol 8,36×104

CP capacidad calorífica, kJ/kmolºC 1,0

XIII.6
 densidad, kg/m3 1000,0
TC temperatura del refrigerante, ºC 25,0
Uj coeficiente global de transferencia de calor de la chaqueta,kJ/m 3ºC s 1,16

U C ,max U C ,kJ/m3 ºC s 4,42

valor máximo de
U C ,min U 1,39
valor mínimo de C ,kJ/ m3 ºC s
Tmax valor máximo de T ,ºC 125,0
Tmin valor mínimo de T,ºC 25,0
TS ,max TS , ºC 150,0
valor máximo de
TS ,min T 70,0
valor mínimo de S , ºC
R constante universal de los gases, kJ/kmol K 8,314
AC / V m2/m3 17,0
Aj / V m2/m3 30,0
tiempo batch, h 1
intervalo de tiempo de integración, s 0,1

Simulación dinámica: Se necesita desarrollar el simulador dinámico basado en el modelo derivado

del proceso. El modelo del reactor batch comprende las ecuaciones diferenciales ordinarias
formuladas (dos ecuaciones de continuidad): Una para el componente A y otro para el componente
B y una ecuación de balance de energía; soportadas por la forma algebraica de las ecuaciones de
velocidad de reacción (dos ecuaciones una para cada reacción). Las ecuaciones diferenciales
ordinarias pueden resolverse utilizando el método de Runge Kutta de cuarto orden. Otra forma es
resolver las ecuaciones diferenciales ordinarias utilizando el enfoque de Euler con el fin de observar
la respuesta del lazo cerrado del controlador proporcional también el esquema de control sugiere
que este sea implementado en el simulador.

Se sugiere que el simulador de este reactor batch sea programado utilizando el lenguaje MATLAB.

La ejecución del programador debe obtener la representación del programa debe obtener una figura,
el seguimiento del rendimiento proporcionado del controlador PI diseñado.
 1 
u  t   uS  kC  e  t    e  t  dt 
 i 
El seguimiento de la temperatura es obtenida mediante el ajuste de la velocidad de calentamiento y
enfriamiento a través de la manipulación de las variables paramétricas u .La figura debe mostrar los
valores de las dos valores de las variables manipuladas: Velocidad de flujo de refrigerante y
temperatura del vapor.

XIII.7
Solución

Codificación del programa

XIII.8
XIII.9
Ejecución del programa

XIII.10
XIII.11