MaríaFernandaCañon Punto 5
MaríaFernandaCañon Punto 5
MaríaFernandaCañon Punto 5
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.
1. Nomenclatura:
Tiempo de residencia τ ks
Temperatura en el reactor 1 T1 K
Temperatura en el reactor 2 T2 K
Fi ,o + rV
i R = Fi ,1 (2)
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:
- 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)
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:
FR ,o − FR ,1 + rAVR = 0 ( 10 )
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
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:
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
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 .
k 0.27
1
F 13.15
R ,1
−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)
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 )
−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):
1. Reactor 1:
FA ,o − FA ,o (1 − x A ,1 )
VR = ( 21 )
k1C A2 ,o (1 − x A ,1 )
2
20 − 20 (1 − 0.33)
VR ,1 =
( 0.27)(1) (1 − 0.33)
2 2
2. Reactor 2:
FA ,o (1 − x A ,1 ) − FA ,o (1 − x A ,2 )
VR = ( 22 )
k2C A2 ,o (1 − x A ,2 )
2
20 (1 − 0.33) − 20 (1 − 0.7 )
VR =
(1.5)(1) (1 − 0.7)
2 2