6.6.
A p l i c a c i ó n del B a l a n c e d e E n e r g í a 313
nivel intermedio de complejidad pudiera emplearse la Ecuación 6.18 para cuando fuera
aceptable suponer que la temperatura de la mezcla reaccionante permanece "constante"
durante el tiempo de residencia del fluido en el serpentín.
Por otro lado, en este ejemplo se trató un análisis en que el flujo volumétrico del agua de
enfriamiento se mantiene constante durante todo el tiempo de retención. Sin embargo, en
un grado de considerable mayor dificultad se podría optimizar un perfil de dicho flujo vo-
lumétrico contra el tiempo de retención. Entonces, el flujo volumétrico sería relativamente
bajo al inicio para permitir cierto aumento en la temperatura; posteriormente se aumen-
taría dicho flujo gradualmente para lograr remover el calor generado que aumenta debido
al aumento en las velocidades de reacción; y finalmente, al irse agotando los reactivos la
velocidad de reacción decaería y entonces el flujo volumétrico del agua de enfriamiento se
disminuiría gradualmente para evitar que la temperatura caiga demasiado. Por supuesto
que la función objetivo deberá reflejar la economía de toda la planta.
6.6.2. Reactor de Tanque Agitado No-Isotérmico
Existen casos muy simples en que es posible desarrollar una solución analítica para las
ecuaciones no-lineales que describen a los reactores continuos de t a n q u e agitado, pero
para un gran número de problemas deberemos recurrir a m é t o d o s numéricos. A u n q u e es
posible plantear todas las ecuaciones involucradas e i n t e n t a r resolverlas s i m u l t á n e a m e n t e ,
se requerirán de excelentes aproximaciones iniciales para lograr convergencia. Siempre que
sea factible, se recomienda un procedimiento secuencial c o m o el que utilizaremos en el
primer ejemplo. Para el primer ejemplo se requiere resolver t a n solo una ecuación implícita
porque se t r a t a de una sola reacción independiente. Para varias reacciones, aunque esto no
se necesitará para nuestro segundo ejemplo, t a m b i é n es posible emplear el procedimiento
secuencial en el que se resuelven s i m u l t á n e a m e n t e las ecuaciones no-lineales para el primer
reactor, los resultados se incluyen en las ecuaciones del segundo reactor que a su vez se
resuelven solas, y así sucesivamente hasta resolver aisladamente las ecuaciones del ú l t i m o
reactor del t r e n .
E j e m p l o 6.3. Se requiere diseñar un tren de tres reactores de tanque agitado y lograr
3
un fracción conversión de 0.9 del reactivo limitante para 1.5 de alimentación a 28°C
con CAO = 0.9 M y CBO — 1-8 M . La reacción en fase liquida es
2 A + 3 B -* P r o d u c t o s
314 C a p í t u l o 6. R e a c t o r e s Ideales No-isotérmicos
La velocidad de reacción está dada por:
-22,300 ^ a l
.i 5 mol
r=kC C
A
l 5
B = 1.3 x 1 0 1 4
-^ e 1 9 8 7
^ T
C C\f
A
mol ' - min
1 5
La reacción es endotérmica y el calor de reacción es 350 — ¡ ¿^reacción • ^ ¡ tar los a r a ca en
reactores se dispone de vapor saturado a 1.34 b a r . Debido al tipo de intercambiador que a
se seleccionó, el área de transferencia dependerá del volumen de reactor de acuerdo a:
A = 0.85 V R , donde el área estará en m y el volumen en m . El coeficiente global de 2 3
c a l
transferencia de calor se estimó en 41,000 . „ . La densidad de la solución es 0.95 2 Q
min m C
y su capacidad calorífica es 0.89 ^ J j -
Si los reactores deben tener las mismas dimensiones, determine el volumen de cada reactor
y el consumo esperado de vapor. Verifique si las soluciones encontradas para los tres
reactores son estables.
Solución: Vamos a comenzar planteando los tres balances de energía independientes. Acla-
ramos que en esta solución emplearemos la temperatura en °C y siempre que utilicemos
ésta en la ecuación de Arrhenius sumaremos 273.2 para obtener temperatura absoluta.
Esto lo haremos porque normalmente pensamos en °C y tenemos un buen sentido de
qué representa su magnitud mientras que si se trata de unidades absolutas debemos ha-
cer la conversión antes de tener una buena idea de que tan caliente o frío está el sistema
en cuestión. Puesto que en el balance de energía se implican incrementos de temperatura,
mientras seamos consistentes, podemos utilizar °C o K. Se recomienda analizar lo anterior
en base a la ecuación que aparece justo antes de la Ecuación 6.23 en la página 293. De
s a t
tablas de vapor obtenemos: T = 108°C y X = 2,235.4 | . Con las concentraciones vap
iniciales y los coeficientes estequiométricos determinamos que A es el rl.
La temperatura a la salida del primer reactor la obtenemos mediante la Ecuación 6.21,
T = UAT +
C1 VQPCPT 0 - AH A V 0 C A o fA?
1
VopC P + UA
_ 41,000 • 0.85V R • 1 0 8 ° C + 1.5 x 1 0 • 0.95 • 0.89 • 2 8 ° C -
6
• 1.5 x 1 0 • 0.9 •
3
fAl
~ 1.5 x 1 0 • 0.95 • 0.89 + 4 1 . 0 0 0 • 0 . 8 5 V
6
ñ
6 3 3
En la expresión anterior el factor 1 0 es para convertir de m a c m mientras que el factor
3 3
10 es para transformar de m a lt. Simplificando,
6.6. A p l i c a c i ó n del B a l a n c e d e E n e r g í a 315
donde T está en °C y VR es el volumen de un reactor en m . Con la Ecuación 6.23 3
calculamos las temperaturas a las salidas del segundo y tercer reactor,
UATc + VppCp T - AH
t A V C 0 Ao (f A2 - f A l )
-Í2 =
VopCp + UA
UAT c + V pC T
0 P 2 - AH A V C 0 AQ (f A3 - f ) A2
T z =
VopCp + UA
Nuevamente simplificando,
= 108 VR + l, 0 1 9 T - 6,779 ( f A 2 - f A l )
2
36.39 + V R
[
'
T = 108 VR + 1, 019 T 2 - 6, 779 (f A3 - f A 2 )
3
36.39 + V R
1
'
Antes de acondicionar las ecuaciones de diseño simplifiquemos un poco la expresión de
velocidad de reacción:
— 11 2 2 3
l 5 1 4 1 5
(-r )
A = 2r = 2k C A C B = 2.6 x I 0 e ^ + ^ (C Ao - C Ao f)
A (C Bo - 1.5 CAo f)-
A
— 11 2 2 3
= 5.651 x 1 0 1 4
e ™ n ¡ (i - f) A (i - 0.75/A) 1 5
Sustituyendo en las ecuaciones de diseño para los tres reactores, la Ecuación 4.11,
VR _ C A O fAl
-11,223 (D)
V
Q 5.651 x i0 1 4
e i - " . 2 (i
T
3 _ j A ¿ ^ _ o 7 5 / A l )i-5
VR C A O (f A2 - f A i) (E)
— 11 2 2 3
V
o 5.651 x 10 e 2+273.2 (x _ f
14 r
A 2 ) (! _ o.75 f A 2 ) l b
(F)
VR CAO (jUg ~ ÍA ) 2
223
V
o 5.651 x 1 0 e 3 + " 3 ( i _ f ) ( i _ 0.75 f ) ' °
1 4
T 2
A z A i
l
Puesto que f Az = 0.9 la Ecuación F se puede simplificar un poco. Se tiene un sistema
de 6 ecuaciones, la Ecuaciones A a F, y seis incógnitas: VR, f , Al f , A2 T\,T y 2 T . Este
3
sistema se puede resolver simultáneamente pero lo podemos reducir a un sistema de 3 x 3
si sustituimos T i de la Ecuación A en las Ecuaciones B y D, la ecuación resultante para T 2
(ya sin T i ) a su vez la sustituimos en las Ecuaciones C y E. Con el nuevo balance de energía
obtenemos una ecuación para T 3 sin T\ ni T 2 que a su vez utilizaremos para eliminarla de
la Ecuación F. El resultado serán 3 ecuaciones (muy largas) con tres incógnitas: V , A f Al
y f , A2 las cuales podemos resolver simultáneamente con aproximaciones iniciales VR > 0
y 0 < f Al < f A2 < 0.9.
316 C a p í t u l o 6. Reactores Ideales No-isotérmicos
En lugar de tratar de resolver las ecuaciones para los tres reactores a la vez, utilizaremos
un procedimiento secuencial. Sustituimos las Ecuaciones A a C en las Ecuaciones D a F,
pero a diferencia del procedimiento del párrafo anterior, ahora no intentaremos eliminar
T\ ni T de las Ecuaciones B a C. Las ecuaciones resultantes son:
2
(G)
Vn
TOS V f í +1,019-6.780 S \ 1
5.651 x 1 0 e V 14 3 6 3 9 + V
H (l-/xÍ)(l-0.75/ A l ) - 1 8
VR C (f -f i)
Ao A2 A
(H)
V 0 "
108 V + 3 6 39 Ti - 6 , 7 8 0 ( / . A 2 ~ / . A l J
f í
5.651 x 1 0 e V 14
3639+vs •-+»»...; (1 _ / A ¡ 2 ) ( i _ 0 . 7 5 / A a ) - 1
Vfl _ C Ao (f A3 - f )
A2 (0
108 V ' + 3 6 - 3 9 T
R 2 - 6 , 7 8 0 {f A 3 ~ f A 2)
5.651 x 1 0 e V 14
*-+273. V (1_/ A 3 )(I_0.75/A'3) - 1 5
El procedimiento es:
1. Suponer VR;
2. Resolver en forma implícita la Ecuación G para f;
Al
3. Evaluar T\ de la Ecuación A con el valor calculado para f;
Al
4. Resolver en forma implícita la Ecuación H para f A2 incluyendo los valores calculados
para f A l y 2\;
5. Evaluar T 2 de la Ecuación B con el valor calculado para f A l y f ;
A2
6. Resolver en forma implícita la Ecuación I para f A3 incluyendo los valores calculados
para f A2 y T; 2
7. Evaluar T 3 de la Ecuación B con el valor calculado para f A2 y f¿ ; 3 y
8. Finalmente, si [f ] aicuiado A3 C 0.9, suponer otro valor para VR hasta lograr conver-
gencia.
Para lograr convergencia podemos emplear un procedimiento gráfico en que suponemos
valores para VR a intervalos constantes, construimos un gráfico y de éste determinamos
el valor de VR correspondiente a f A3 = 0.9. Puesto que podemos esperar tiempos de
residencia de varios minutos vamos a suponer valores de V R = 10,15 . . . 45 m . Si el valor 3
deseado para la fracción conversión no se encuentra en este rango habrá que ampliarlo de
acuerdo a la tendencia observada. La Figura 6.4 presenta la manera gráfica para encontrar
el valor de los volúmenes de reactor. Los puntos marcados con x muestran los resultados
calculados con el procedimiento anteriormente descrito. La curva une los puntos en forma
suavizada. De la curva observamos que una fracción conversión de 0.9 corresponde a
aproximadamente V R = 33.5 m . 3
6.6. A p l i c a c i ó n del B a l a n c e d e E n e r g í a 317
10 15 20 25 30 35 40 45
3
V o l u m e n de c a d a r e a c t o r , m
Figura 6.4: Solución gráfica para encontrar el v o l u m e n de cada reactor para el Ejemplo
6.3.
Los resultados intermedios para VR = 33.5 m , en orden de cálculo, son: f
3
Al = 0.3850,
T i = 29.01°C, f A2 = 0.6878, T 2 = 37.50°C, f Az = 0.9031 y T 3 = 50.41°C.
Flujo másico de vapor:
El flujo másico de vapor utilizado en los tres reactores lo calculamos mediante:
_ UA{T C - T i ) + UA(T C - T ) + UA{T
2 C - T)
^vap
_ 41,000 • 0.85 • 33.5 (3 • 108 - 29.01 - 37.50 - 50.41)
• 452.5
xfcgj
=
2,235.4 1000 g min
Estabilidad de las soluciones:
Como se trata de una sola reacción independiente, para verificar la estabilidad de las
soluciones vamos a emplear el criterio dado por la Ecuación 6.27 y para evaluarlo necesi-
taremos un resultado de la Ecuación 6.28.
La solución en estado estacionario encontrada para el primer reactor fue f Al — 0.385
y T = 29.0°C. Para verificar su estabilidad debemos perturbar esta temperatura y uti-
lizaremos + 1 ° C . El tiempo espacial para cada reactor, VRI/VQ,
1 4
es 22.33 min. Con
14
La perturbación debe ser relativamente pequeña, pero lo suficientemente grande para que las cifras
significativas de los cálculos permitan evaluar el criterio de estabilidad con cierto grado de contabilidad.
318 C a p í t u l o 6. Reactores Ideales No-isotérmicos
= 30.0°C evaluamos la constante de velocidad y la sustituimos en la ecuación de
diseño escrita como función implícita, modificando la Ecuación G obtenemos:
F ( / + ) = 0 = 22.333
A 0.9 f + A
( — 11 223 \
14
5.651 x I 0 e ^ o . o + 2 7 3 . J ( i - f +) 2 A ( i _ 0.75 ¡AI) 15
La solución numérica es f f A = 0.406. Sustituyendo en el criterio de estabilidad, la
Ecuación 6.27, se tiene:
UA(T C - T + ) - V pCp 0 (T+ - T ) < AH F 0 A Ao (f +A - 0)
7 7
¿8.85 x 1 0 < 9.56 x 10 ?
Por lo tanto, la solución encontrada para el primer reactor en estado estacionario es
estable.
Las soluciones previamente encontradas para los otros dos reactores son: f A2 = 0.
T 2 = 37.5°C, f A3 = 0.903 y T 3 = 50.4°C. Considerando nuevamente un cambio de
+
+ 1 ° C , las temperaturas perturbadas son T 2 = 38.5°C y T + = 51.4°C. Los valores
3
perturbados: j A \ = 0.701 y f $ A = 0.909, resultan de resolver respectivamente:
F C / t f ) = 0 = 22.333 0.8 0.385)
5.651 x I0 e^8.5+273.2 J ( l - f +)
14
( l - 0.75 A f t) -
A
l b
0.9 (f ¡ - 0.688) A
F(f +)
A = 0 = 22.333 / -11,223 \ , .
5.651 x I 0 e ^ i . 4 + 2 7 3 . J ( i _ f +)
14
( i _ 0.75/^3 )
2 A
L 5
Aplicando el criterio de estabilidad a ambos reactores:
UAiTc ~ T+) - V pC 0 P ( T + - T i ) < AH F A Ao (f +A - f A l )
7 7
¿6.91 x 1 0 < 7.46 x 10 ?
UA(T C - T +) - V o p C p (T + - T ) < AH F
3 3 2 A AQ {f +A - f )
A2
7 7
¿4.85 x 1 0 < 5.22 x 10 ?
concluimos que ambas soluciones en estado estacionario también son estables.
E j e m p l o 6.4. Se cuenta con un reactor continuo de tanque agitado de 28 m 3
para
producir D por polimerización de A. El reactor está enfriado por agua y el flujo máximo que
se puede introducir de este liquido es 1,500 ^ . El agua entra a 25°C y su temperatura de
