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 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.
Agua de enfriamiento(entrada)

control de flujo

control de Agua de enfriamiento(salida),Tc

presión

suministro
de vapor serpentin de refrigeración

Agitador
chaqueta

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

Considere que ocurren en el reactor las siguientes reacciones consecutivas:

Se asume que la reacción tiene una cinética de segundo orden, mientras que la reacción
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 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 H2O de enfriamiento
-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.
-El coeficiente global de transferencia de calor para la chaqueta de calentamiento se asume
constante; pero el coeficiente global de transferencia de calor para el refrigerante se asume que
es función de la velocidad de flujo del refrigerante, como:

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:

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

Velocidad de acumulación de masa dentro del reactor

Desde que es la densidad de la mezcla y es el volumen de la mezcla, es la masa total
dentro del reactor. Sustituyendo equivalencias en la ecuación básica de balance, resulta

(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
Velocidad de generación del componente A por reacción química

Velocidad de acumulación del componente A dentro del reactor

Comentario
Por convención es la velocidad de desaparición de la especie A; consecuentemente, es la
velocidad de formación de A. Aquí representa el número de moles de la especie A en el
sistema en el tiempo t. Ahora la ecuación básica de balance toma la forma:

(13.2)
es decir:

(13.3)
donde es la concentración de la especie A. Desde que la reacción tiene una cinética
de segundo orden, de modo que
(13.4)
donde representa la constante de velocidad de reacción.
Sustituyendo la ecuación (13.4) en la ecuación (13.3) se obtiene:

(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

Velocidad de generación del componente B por reacción química

Velocidad de acumulación del componente B dentro del reactor

Comentario
Aquí, , es la velocidad de desaparición de la especie 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:

(13.6)

(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

(13.8)
donde: : Factor de frecuencia o factor pre-exponencial.
: Energía de activación.
: 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.

(13.9)

(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:

(13.11)
Entonces,

Aquí, y 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:

XIII.4
Comentario
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, Los subíndices 1 y 2 se usan para las reacciones y
respectivamente. Aquí , 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:

Velocidad de acumulación de energía

Aquí, son constantes; es la entalpía . Sustituyendo términos en la ecuación 13.11,

obtenemos la siguiente forma:

(13.12)

Reordenando,

(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 ( ), 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, . Realmente el controlador manipula ambas variables, y
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:

XIII.5
 (13.14)
En la práctica general, una variable de salida (T) es controlada utilizando solamente la variable .
Si usamos otra variable o (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 definida de la siguiente manera:
(13.15)
(13.16)
donde los valores máximo y mínimo de y son elegidos en función del proceso o tomando en
cuenta límites de seguridad. Evidentemente representa el máximo enfriamiento del sistema y
el máximo calentamiento. Esta variable paramétrica es utilizada como una variable única a
ser 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

(13.17)
donde

(13.18)

(13.19)

(13.20)

(13.21)

(13.22)

(13.23)
Ahora 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:

(13.24)
donde:
: Señal de polarización del control.
: Ganancia proporcional (ajuste de parámetro).
: Constante de la integral de tiempo.
: Error del controlador.

Los valores de los parámetros de ajuste son elegidos arbitrariamente como y

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

XIII.6
 concentración de la especie A en el arranque, kmol/m 3 1,0
concentración de la especie B en el arranque, kmol/m 3 0,00

factor de frecuencia de la reacción ,m3/kmol s 1,1

factor de frecuencia de la reacción ,s-1 172,2

energía de activación para la reacción ,kJ/kmol 2,09×104

energía de activación para la reacción ,kJ/kmol 4,18×104

calor de reacción para la reacción , kJ/kmol 4,18×104

calor de reacción para la reacción , kJ/kmol 8,36×104

capacidad calorífica, kJ/kg ºC 1,0
densidad, kg/m3 1000,0
temperatura del refrigerante, ºC 25,0
coeficiente global de transferencia de calor de la chaqueta, kJ/m 3ºC s 1,16
4,42
valor máximo de ,kJ/m2 ºC s
1,39
valor mínimo de ,kJ/ m2 ºC s
valor máximo de T ,ºC 125,0
valor mínimo de T,ºC 25,0
150,0
valor máximo de , ºC
70,0
valor mínimo de , ºC
constante universal de los gases, kJ/kmol K 8,314
m2/m3 17,0
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.

XIII.7
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 .La figura debe mostrar los
valores de las dos valores de las variables manipuladas: Velocidad de flujo de refrigerante y
temperatura del vapor.
Obviamente:

De modo que podemos escribir:

XIII.8
Solución

Codificación del programa

XIII.9
XIII.10
Ejecución del programa

XIII.11