19 junio del 2020

Universidad nacional de Colombia

Ingeniería de las reacciones químicas

[email protected]

Baterías de CSTRs, reactores con recirculación, caída de presión en reactores,

múltiples reacciones en reactores químicos y estado dinámico de reactores CSTR.
María Fernanda Cañón Ochoa[1]
Estudiante de ingeniería química, código. 316516

Resumen: En el presente escrito se hace el desarrollo para una batería de dos CSTRs que operan a condiciones
diferentes de operación. El calculo se divide en 3 parte, primero el cálculo de las conversiones para cada reactor,
segundo el cálculo de la temperatura en el reactor 1 y finalmente el cálculo del volumen para cada reactor. El
planteamiento matemático parte del balance molar por componente de un CSTR y se extrapola al análisis de
baterías; además, se hace uso de la ecuación de Arrhenius para conocer la condición de temperatura en el reactor
1.

PROBLEMA 5
• Planteamiento:

La reacción irreversible, A →2R, se lleva a cabo en fase líquida en dos reactores CSTR en serie. La temperatura
del segundo reactor es de 120 °C. Se alimenta al primer reactor 20 kmol de A/ks y el caudal de salida de R del
último reactor es 28 kmol/ks. La constante específica de velocidad de reacción es de 1.5 m3/(kmol.ks) a 120 °C
y la energía de activación es 84 kJ/mol. El tiempo de residencia promedio en el primer y segundo reactor es de
2.75 ks. La concentración de reactivo es de 1 kmol/m3. Los dos reactores funcionan de forma isotérmica, pero
con diferentes temperaturas. Calcule las conversiones, la temperatura del primer reactor, y los volúmenes de
cada reactor que permiten alcanzar el desempeño especificado.

Figura 1. Arreglo de n CSTRs en serie

 • Planteamiento matemático(Cálculo de conversiones):

1. Nomenclatura:

Variable Símbolo Unidades

Flujo inicial de A FA,o kmol/ks

Flujo de A final en el reactor 1 FA,1 kmol/ks

Flujo de A final en el reactor 2 FA,2 kmol/ks

Flujo de B final en el reactor 1 FB,1 kmol/ks

Flujo final de B FB,2 kmol/ks

Tiempo de residencia τ ks

Constante especifica de velocidad kA,1 m3/(kmol.ks)

Conversión en el reactor 1 XA,1 Adimensional

Conversión en el reactor 2 XA,2 Adimesional

Ley de velocidad rA Kmol/(m3.ks)

Flujo volumétrico ν m3/ks

Constate especifica de velocidad k1 m3/(kmol.ks)

en el reactor 1

Constate especifica de velocidad k2 m3/(kmol.ks)

en el reactor 2

Energía de activación E Kj/mol

Constante de gases R Kj/mol.K

Temperatura en el reactor 1 T1 K

Temperatura en el reactor 2 T2 K

2. Balance molar para un CSTR en estado estable.

 E +G = S + A (1)

Considerando un estado estable entonces el balance sería:

Fi ,o + rV
i R = Fi ,1 (2)

3. Balance para reactor 1.

Siendo el análisis para una batería de dos CSTRs se analiza cada uno de los reactores y cada uno de los
compuestos por separado como se va a desarrollar a continuación. Además, se tiene en cuenta que las
condiciones de temperatura son diferentes en los dos reactores por lo que se diferencia entre k 1 y k2.

- Ley de velocidad

Considerando las unidades reportadas para la constante especifica de velocidad, se infiere que la ley de
velocidad es de segundo orden y vendría expresada por la siguiente ecuación:

−rA = kC A2 (3)

- Tabla estequiométrica:

Tabla 1. Tabla estequiométrica.

Especie Inicio Cambio Remanente
A FA,o -FA,oxA FA,o(1-XA)
B 0 2FA,oXA 2FA,oXA
Total FA,o FA,oXA FA,o(1+XA)

- Componente A:

Partiendo de la ecuación (2), conociendo la estequiometria y considerando la ley de velocidad planteada, siendo
A un reactivo se tiene:

Fi ,o − Fi ,o (1 − X A ,1 ) − rAVR = 0 (4)

Redistribuyendo la ecuación (4) se obtiene:

Fi ,o XA,1 − rAVR = 0 (5)

Haciendo el reemplazo de la educación (3) en la ecuación (5) se tiene:

 Fi ,o X A,1 − ( k1C A2 ,1 )VR = 0 (6)
F
Sabiendo que C = i
y que de la estequiometria FA ,1 = FA ,o (1 − X A ,1 )
ν
i

 F 2 (1 − X )2 
Fi ,o X A,1 −  k1 VR = 0
A ,o A ,1
(7)
 ν 2

 
V
Reorganizando la ecuación de forma tal que se incluya el termino de tiempo de residencia ( τ = ):
R

 F 2 (1 − X )2 
FA ,o X A ,1 −  k1 τ = 0
A ,o A ,1
(8)
 ν 
 

Para este cálculo, conociendo la concentración inicial y el flujo molar inicial, teniendo en cuenta que la reacción
es en fase liquida y no hay cambio del volumen, el flujo volumétrico se calcula mediante la siguiente expresión:

FA ,o
ν = (9)
C A ,o
En este caso ν = 20 m3/ks y el tiempo de residencia se considera de 2.75 ks.

- Componente R:

Se considera el mismo balance molar en este caso (ecuación (2)).

FR ,o − FR ,1 + rAVR = 0 ( 10 )

Conociendo la estequiometria de R en la tabla 1 y sabiendo que no hay flujo de R en el alimento:

−2FA,o X A,1 + 2rAVR = 0 ( 11 )

Reemplazando la ley de velocidad, teniendo en cuenta el sigo positivo de formación de R:

−2FA,o X A,1 + k1C A2 ,1VR = 0 ( 12 )

Haciendo la misma analogía para la concentración que en el caso del componente A (la definición de
concentración y la estequiometria del componente A):

 F 2 (1 − x )2 
−2FA,o X A,1 + 2k1  τ = 0
A ,o A ,2
( 13 )
 ν 
 
nótese que en la expresión la ley de velocidad es multiplicada por dos, esto de realizar ley de proporciones
múltiples como se muestra a continuación:

rA rR
= ( 14 )
−1 2

4. Balance para reactor 2:

La estequiometria para este caso es igual que en la table 1 pero en este caso se corre la posición de los flujos
como se muestra a continuación.

- Componente A:
 F 2 (1 − X )2 
FA,1 − FA,2 −  k2 τ = 0
A ,o A ,2
( 15 )
 ν 
 
El termino FA,2 por estequiometria va referido a la conversión en el alimento fresco, es decir, respecto a XA,2.
Reemplazando en la ecuación (13)

 F 2 (1 − X )2 
FA ,o ( X A ,2 − X A ,1 ) −  k2 τ = 0
A ,o A ,2
( 16 )
 ν 
 

- Componente B:

Haciendo la misma analogía que en el reactor 1 se tendría:

 F 2 (1 − x )2 
FR ,1 − FR ,2 + 2k2  τ = 0
A ,o A ,2
( 17 )
 ν 
 
5. Muestra de cálculo

 ( 20 ) (1 − X A,1 ) 
2 2

(20 ) X A,1 −  (1.5)  (2.75) = 0

 20 
 

6. Desarrollo

El desarrollo se hizo haciendo uso del programa Matlab. Allí se plateo el sistema de ecuaciones (ecuaciones (8), (13), (16)
y (17) ). Con la ayuda del comando “fsolve” se le dio solución. Nótese que en este caso se tiene un sistema de 4 ecuaciones
con 4 incógnitas, siendo las incógnitas: X A ,1 , X A ,2 , k1 , FR ,1 .

Dando solución, el sistema tiene las siguientes respuestas:

 Incógnita Respuesta
X A ,1
0.33
X A ,2
0.70

k 0.27
1

F 13.15
R ,1

• Planteamiento matemático (cálculo de la temperatura en el reactor 1):

Para este inciso se hizo uso de la ecuación de Arrhenius.

 −E 
k (T ) =A exp   ( 18 )
 RT 
1. Cálculo del factor preexponencial:

El primer cálculo está dirigido a la determinación del factor preexponencial presente en la ecuación de Arrhenius
con las condiciones del reactor 2 que conocemos.

k2
=A
 −E  ( 19 )
exp  
 RT2 
Reemplazando los valores:

1.5
=A
 −84 
exp  
 ( 8.314 x10 −3 ) ( 393.15) 
 

El factor preexponencial tiene un valor de 2.17 x 1011 m3/(kmol.ks).

2. Cálculo de la temperatura:

Del cálculo de las conversiones fue posible conocer la constante especifica de velocidad del reactor 1. Teniendo
esta información entonces:

−E
T1 =
R (ln( k1 ) − ln( A ) )
( 20 )

Reemplazando los valores:

−84
T1 =
(
8.314 x103 ln( 0.27 ) − ln(2.17 x1011 ) )
La temperatura en el reactor 1 es 95.17 °C.
 • Planteamiento matemático (Cálculo del volumen):

Haciendo uso de la ecuación (2) y reemplazando la ley de velocidad:

1. Reactor 1:
FA ,o − FA ,o (1 − x A ,1 )
VR = ( 21 )
k1C A2 ,o (1 − x A ,1 )
2

Reemplazando los valores:

20 − 20 (1 − 0.33)
VR ,1 =
( 0.27)(1) (1 − 0.33)
2 2

El volumen del reactor 1 es de 55 m3.

2. Reactor 2:

De manera análoga como en el reactor 1:

FA ,o (1 − x A ,1 ) − FA ,o (1 − x A ,2 )
VR = ( 22 )
k2C A2 ,o (1 − x A ,2 )
2

Reemplazando los valores conocidos:

20 (1 − 0.33) − 20 (1 − 0.7 )
VR =
(1.5)(1) (1 − 0.7)
2 2

El volumen del reactor 2 es de 55 m3.