Enfriamiento de Reactor
Enfriamiento de Reactor
Enfriamiento de Reactor
Universidad de Carabobo
farteaga20@gmail.com
2000
Alfredo E. Varela V. / Daniel E. Machado L. / Richard J. Romero G.
DISEO DE UN SISTEMA DE ENFRIAMIENTO PARA UN REACTOR DE TANQUE
AGITADO
Ingeniera UC, junio, ao/vol. 7, nmero 001
Universidad de Carabobo
Valencia, Venezuela
En el presente trabajo se hace un diseo de un sistema de enfriamiento para un reactor de mezcla de una
empresa productora resinas.
Para enfriar las paredes del reactor desde temperaturas cercanas 150 C hasta 40 C se propone utilizar un aceite
trmico (terminol), cuya temperatura es controlada mediante un intercambiador de calor externo de tubos y
coraza. Para la simulacin del proceso se elabor un programa en Visual Basic bajo ambiente Excel 97. Se
emple el modelo de celdas (Gaddis y Schlnder, 1979) para el modelo del intercambiador de calor y dos
controladores del tipo proporcional para el control de la temperatura del aceite trmico y de la pared del reactor.
Los resultados de la simulacin indican que la temperatura de la pared del reactor se estabiliza en el valor
deseado en un tiempo de 40 minutos, que es mucho menor que el tiempo empleado si permite que el tanque se
enfre por si solo. Si este tiempo se empleara en produccin, incidira en un aumento del 6 % de su capacidad de
produccin mensual.
Abstracts
A cooling system is proposed for a mixed tank reactor in a resin plant. Using a thermal oil is suggested to cool
the reactor walls from about 150 C to 40 C. Temperature of the thermal oil is controlled using an external
shell-and-tube heat exchanger. Process simulation is performed with computer program in Visual Basic running
under Excel 97. A cell model (Gaddis y Schlnder, 1979) and two proportional controllers were used for the
heat exchanger and the control of thermal-oil temperature ad reactor wall temperature, respectively.
Simulation results suggest that the reactor wall temperature is reached in about 40 minutes which is much lower
than the time employed by allowing the tank to cool by itself. Using this time interval in production will
increase the monthly resin production in 6%.
INTRODUCCION
El equipo principal de una planta de produccin de resinas es un reactor de tanque agitado con una capacidad de
1500 galones. La planta fabrica cinco tipos de resinas: acrlicas, polister, alqudicas, epxica y nitrogenada. En
los diferentes procesos de elaboracin de resinas existen algunas que alcanzan temperaturas de reaccin de
aproximadamente 230C. Luego de retirada la carga de producto, es necesario enfriar las paredes del reactor
antes de alimentar los materiales de la siguiente carga. Esta etapa se viene efectuando por intercambio
espontneo con el ambiente o utilizando un procedimiento manual. La formulacin de tesinas acrlicas establece
que algunas materias primas slidas se deben cargar al reactor a una temperatura por debajo de los 40 C para
evitar su descomposicin, lo que puede producir gases inflamables e ignicin espontnea. Las resinas de
poliesteres y alqudicas, que en su formulacin utilizan materia prima slida, requieren temperaturas del reactor
por debajo de l0 80 C, a fin de evitar su descomposicin.
El tiempo empleado con el enfriamiento espontneo incide directamente en la rata de produccin, restando
versatilidad y seguridad al reactor. Este trabajo propone un sistema que permite enfriar el reactor en forma
uniforme y reducir el tiempo de procesamiento, contribuyendo al aumento de la produccin.
SISTEMA DE ENFRIAMIENTO Y ESPECIFICACIONES DE EQUIPOS
La Figura 1 muestra un esquema del sistema de enfriamiento sugerido. El cual consiste en un esquema de
controla en cascada. La seal de salida del primer controlador ser el punto de referencia de un segundo, y sta
se comparar con la seal de temperatura del aceite proveniente del transmisor colocado a la salida del
intercambiador. Una esquema en cascada asegurara una rpida correccin de la temperatura del aceite de
enfriamiento.
ms rpida que la del intercambiador, se puede considerar que los cambios en las posiciones de las vlvulas se
vern reflejados en cambios en escaln en la temperatura del agua de alimentacin al intercambiador.
El intercambiador de calor fue diseado siguiendo procedimientos conocidos (Kern, 1965) y las caractersticas
del mismo son mostradas mediante la Tabla 1. Las especificaciones de la chaqueta del reactor fueron tomadas
en la planta y son dadas mediante la Tabla 2. Los parmetros de diseos de estos equipos son usados para
simular el comportamiento del sistema.
MODELOS MATEMTICOS
A continuacin se plantean los modelos matemticos del intercambiador de calor, pared del reactor y chaqueta
del reactor. La dinmica correspondiente a otros equipos de proceso y conexiones se considera despreciable.
Intercambiador de calor
Se emple el modelo de celdas (Gaddis y Schlnder, 1979) para el intercambiador de calor. El nmero de celdas
es determinado por el nmero de deflectores en la coraza y el nmero de paso de tubo.
Tabla 1. Especificaciones del intercambiador de calor
Las Figura 2 muestra en el esquema un Intercambiador de calor 1-2. Las celdas pueden ser conectadas de
acuerdo con la estructura del intercambiador de calor y la disposicin de las conexiones fsicas de entrada y
salida.
Para el anlisis de la transferencia de calor en cada celda se hacen las siguientes suposiciones:
1.- La celda est compuesta de dos partes separadas por una superficie diatrmica comn A (i).
2.- Los fluidos se consideran perfectamente mezclados en cada celda, y por consiguiente tienen una
temperatura uniforme.
3.- Toda transferencia de calor entre los fluidos ocurren a travs del rea comn A (i), cualquier otra
superficie que limita la celda es adiabtica, o es un rea de entrada / salida para los fluidos.
4.- Todas las propiedades fsicas son independientes de temperatura.
5.- Las partes metlicas del intercambiador de calor son incorporadas en el modelo de celda bajo las
siguientes consideraciones:
La masa metlica de los tubos est a la temperatura media entre la temperatura del fluido que va
por los tubos y la del fluido del lado de la coraza.
La masa metlica de la coraza est a la temperatura del fluido del lado de la coraza.
6.- La conduccin de calor a travs de las paredes de tubos y la divisin entre celdas es despreciable.
La Figura 3 muestra una modelo de una celda con flujos en contracorriente, donde se identifican las siguientes
variables:
El subndice "s" significa que pertenece al lado de la coraza, y el "t" al lado de los tubos. Se asume tambin que
estas condiciones son vlidas para cualquier forma de flujo, es decir: concurrente, contracorriente o flujo
cruzado.
donde: t = densidad del fluido por los tubos, lb / pie3, Cpt = calor especifico del fluido por los tubos, Btu /
lb.F, Gt = flujo msico del fluido por los tubos, lb / h, U = coeficiente global de transferencia de calor, Btu /
h.ft2.F. En este balance de energa, el calor especfico de la celda i es definido como as
donde xt(i) es un parmetro que considera la masa metlica del tubo en la celda y cuya justificacin es propuesta
por Correa y Marchetti (1987). En forma adimensional, la Ecuacin (1) se expresa
Donde
Las ecuaciones (3) y (5) son validas para cualquiera de las NM celdas en las cuales sea dividido el
intercambiador de calor. Por lo tanto las 2NM ecuaciones diferenciales de primer orden describen el
comportamiento dinmico del intercambiador.
El procedimiento usado para integrar estas ecuaciones se basa en la aproximacin a la derivada mediante
Sustituyendo las Ecuaciones (7) y (8) en las ecuaciones (3) (5), y usando k y k + 1 para los instantes de tiempo 0
y 0 + 0 resulta:
donde
Para simplificar, los argumentos i y j no se incluyeron en los coeficientes a y b. Sin embargo, estos tienen
diferentes valores si el rea de las celdas A(i), los volmenes Vt(i) y Vs(i) o las partes metlicas xt(i) y xs(i) no son
los mismos de una celda a otra. Estos valores dependen del tipo intercambiador de calor. Para el caso en
estudio, por ejemplo, todas las celdas tienen la misma rea de transferencia que se puede expresar mediante
donde De es el dimetro externo del tubo, es el espaciado entre bafles y n el nmero de paso por tubos. El
volumen de fluido en la celda por el lado de tubos, Vt(i), y del lado de coraza, Vs(i), se pueden estimar, con las
ecuaciones
donde Din y De son los dimetros interno y externo de los tubos, respectivamente, y Ds es el dimetro de la
coraza.
El sistema dado por las Ecuaciones (11) y (13) se cede resolver con las condiciones iniciales
Las condiciones de las entradas para el lado de tubos y el lado de coraza son dadas, respectivamente por
donde: M = densidad del metal, lb / ft3, CpM = calor especifico del metal, Btu / lb.F, VM = volumen del fluido
en la chaqueta, ft3, TM = temperatura del metal, F, ho = coeficiente convectivo de transferencia de calor, Btu
/h.ft2.F, Ao = rea de transferencia de calor entre la chaqueta y el reactor, ft2, Tc = temperatura del fluido en la
chaqueta, F.
El coeficiente convectivo de transferencia de calor, est dado por (Kern, 1965)
donde: Re = nmero de Reynolds del fluido de la chaqueta, adim, Pr = nmero de Prandatl, adim, Cps = calor
especifico del fluido de la chaqueta, Btu / lb.F, G'ps = flujo msico por unidad de rea, lb / h.pie2. El nmero
de Reynolds del fluido de la chaqueta se calcula mediante la siguiente ecuacin
donde: s = viscosidad del fluido de la chaqueta, lb / h.ft, DH = dimetro hidrulico de la chaqueta, pie, DH =
JH/4, JH = distancia entre la pared del reactor y la chaqueta, (JH = 50 mm).
La ecuacin (19) puede ser rearreglada adimensionalmente como sigue
en la cual
donde: Tc(ent) = Temperatura del aceite entrando a la chaqueta, F, Tc(sal) = Temperatura del aceite saliendo
de la chaqueta, F. Esta ecuacin en forma adimensional se puede escribir
donde:
Usando las aproximaciones dadas por las Ecuacin (7) y (8), las ecuaciones (22) y (25) se discretizan resultando
las siguientes expresiones para los balances de energa para la pared y chaqueta del reactor, respectivamente
Las ecuaciones (11) y (12), que describen el comportamiento del intercambiador, en conjunto con las
ecuaciones que representan el comportamiento de la pared del reactor (27) y de la chaqueta (28) permiten
simular el comportamiento dinmico del sistema de enfriamiento.
Funciones de transferencia
Para disear los controladores se determinan las funciones de transferencia para las variaciones de la
temperatura de la pared del reactor y del aceite a la salida del intercambiador con respecto a las variables
manipulables. La funcin de transferencia para la variacin de la temperatura del aceite a la salida del
intercambiador con respecto a la temperatura del agua a la entrada del intercambiador se expresa como un
sistema de primer orden con retardo
del mismo modo, la funcin de transferencia para la variacin de la temperatura de la pared del reactor con
respecto a la temperatura del aceite a la salida del intercambiador es
valores de
Para trabajos futuros referentes a la simulacin de intercambiadores de calor de tubos y coraza se recomienda
modificar la estructura del programa para adaptarlos a intercambiadores ms complejos. Tambin se sugiere
estudiar otros modelos de las funciones de transferencia con el fin de evitar oscilaciones en las respuestas.
REFERENCIAS
1.- Burden, R.L. & J.D. Faires. 1996. Anlisis Numrico. Iberoamrica, Mxico.
2.- Correa, D. and Marchetti, J. (1987). Dynamic simulation of shell and tube heat exchangers. Heat transfer
engineering, vol. 8, n 1, pp 50-59, 1987.
3.- Gaddis, E.S. and Schlnder, E.U. (September, 1979). Temperature distribution and heat exchange in
multipass sell and tube exchangers with baffles. Heat transfer engineering. vol. 1, n 1, pp 43-52, July September, 1979.
4.- Kern, Q. Donald. (Agosto, 1965). Procesos de transferencia de calor. Editorial Mc Graw Hill. 10
edicin. Mexico, 1965.
5.- Seborg, Dale E. and Thomas F, Edgar. . (1989). Proc dynamics and control. John Wiley and Sons, New
Y 1989.
6.- Sundaresan, K. R., and P. R. Krishnaswamy. (1977 Estimation of delay, time constant parameters in
time, frequency and Laplace domains. Can. J. Chem. Eng., 56, 257 (1977).
AGRADECIMIENTOS
Este trabajo se realiz gracias a la colaboracin de la empresa DUPONT DE VENEZUELA C.A y al CDCHT
de la Universidad de Carabobo.