Problema 5.

5 (Levenspiel)
Una alimentación acuosa de A y B (400 L/min, 100 mmol/L de A, 200
mmol/L de B) va a ser convertida en producto en un reactor de flujo en
pistón. La cinética de la reacción está representada por:
𝐴+𝐵 →𝑅
𝑚𝑜𝑙
−𝑟𝐴 = 200𝐶𝐴 𝐶𝐵
𝐿 𝑚𝑖𝑛

Halle el volumen del reactor requerido para alcanzar el 99.9% de

conversión de A en producto.
A partir de la ecuación de diseño de un reactor PFR:
𝑋𝑓
𝑑𝑋𝐴
𝑉 = 𝐹𝐴𝑜 ∫
0 −𝑟𝐴
Y se tiene la tasa de reacción del sistema:
𝑚𝑜𝑙
−𝑟𝐴 = 𝑘 𝐶𝐴 𝐶𝐵
𝐿 𝑚𝑖𝑛
Dejar las variables de tasa de reacción en términos del reactivo limitante A
utilizando la estequiometría requerida:
𝐶𝐴 = 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
𝐶𝐵 = 𝐶𝐵𝑜 − 𝐶𝐴𝑜 𝑋𝐴
−𝑟𝐴 = 𝑘 𝐶𝐴𝑜 (1 − 𝑋𝐴 )(𝐶𝐵𝑜 − 𝐶𝐴𝑜 𝑋𝐴 )
Simplificando la ecuación anterior:
𝐶𝐵𝑜
−𝑟𝐴 = 𝑘 𝐶𝐴𝑜 2 (1 − 𝑋𝐴 ) ( − 𝑋𝐴 )
𝐶𝐴𝑜
𝐶𝐵𝑜 200
𝑀𝐵 = = =2
𝐶𝐴𝑜 100
−𝑟𝐴 = 𝑘 𝐶𝐴𝑜 2 (1 − 𝑋𝐴 )(2 − 𝑋𝐴 )
Sustituir la ecuación de tasa de reacción en la ecuación de diseño
𝑋𝑓
𝑑𝑋𝐴
𝑉 = 𝐹𝐴𝑜 ∫ 2
0 𝑘 𝐶𝐴𝑜 (1 − 𝑋𝐴 )(2 − 𝑋𝐴 )

Si 𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜
𝑋𝑓
𝑣𝑜 𝑑𝑋𝐴
𝑉= ∫
𝑘 𝐶𝐴𝑜 0 (1 − 𝑋𝐴 )(2 − 𝑋𝐴 )
𝑚𝑚𝑜𝑙 1𝑚𝑜𝑙 𝑚𝑜𝑙
𝐶𝐴𝑜 = 100 × = 0.1
𝐿 1000 𝑚𝑚𝑜𝑙 𝐿
Siendo k=200 1/min y vo=400L/min, resolver la ecuación del volumen:
0.999
400 𝑑𝑋𝐴
𝑉= ∫ = 124.312 𝐿
200 ∗ 0.1 0 (1 − 𝑋𝐴 )(2 − 𝑋𝐴 )
Problema 5.11
La enzima E cataliza la fermentación del sustrato A (el reactivo) para que se
convierta en producto R. Calcular el tamaño del tanque agitado necesario para una
conversión 95 % del reactivo con una corriente de alimentación (25 L/min) de
reactivo (2 mol/L) y enzima. La cinética de la fermentación a esta
concentración de enzima está dada por:
enzima
𝐴 → 𝑅

0.1 𝐶𝐴 𝑚𝑜𝑙
−𝑟𝐴 =
1 + 0.5𝐶𝐴 𝐿 𝑚𝑖𝑛
La ecuación de diseño de un reactor CSTR se expresa como:
𝐹𝐴𝑜 𝑋𝑓
𝑉=
−𝑟𝐴
En donde, 𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜 y 𝐶𝐴 = 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
Sustituyendo cada variable en la ecuación de diseño se tiene que:
𝐶𝐴𝑜 𝑣𝑜 𝑋𝑓
𝑉=
0.1 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
1 + 0.5 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
Sabiendo que XA=0.95, vo=25L/min y CAo=2mol/L
2 ∗ 25 ∗ 0.95
𝑉= = 4987.5 𝐿
0.1 ∗ 2 ∗ (1 − 0.95)
1 + 0.5 ∗ 2 ∗ (1 − 0.95)
Problema 5.9 (Levenspiel)
Una enzima específica actúa como catalizador en la fermentación de A. Para una
concentración dada de la enzima en la corriente acuosa (25 litros/min), Calcular el
volumen del reactor de flujo en pistón requerido para el 95 % de conversión del
reactivo A (CA0 = 2 mol/L) a una concentración dada de la enzima. La cinética de la
fermentación a esta concentración de enzima es:

enzima
0.1C A mol
A  R rA 
1  0.5C A litro  min

Solución

A partir de la ecuación de diseño de un reactor FPR


FAO  C AO  (1.1)

C A  C AO (1  X A ) (1.2)
XA dx
V  FAO  (1.3)
0 r A

 XA dx
V  C AO   (1.4)
0 rA
 XA dX A
Sustituyendo la velocidad de reacción en (1.4) V  C AO  
0 0.1C A
1  0.5C A
En términos de la conversión a partir de la ecuación (1.2)

 X     
 XA dX A dX A 0.5dX A  ln(1  X A )  10
V  C AO   C AO       C AO    5X A 
A

0 0.1C AO (1  X A ) 0
 0.1C AO (1  X A ) 0.1   C AO 
1  0.5C AO (1  X A )

  ln(1  0.95)  10 
V  2 litro
mol
 25 litro
mol
 mol
 5  0.95   986.43litros
 2 litro 

V Reactor  986.43litros
Problema 5.24 (Levenspiel)
Se alimenta continuamente un hidrocarburo gaseoso de alto peso molecular A a
un reactor de tanque agitado que se mantiene a alta temperatura, donde se
craquea térmicamente (reacción homogénea en fase gaseosa) produciendo
materiales de menor peso molecular, llamados conjunto R, con una estequiometria
A→5R. Al cambiar el flujo de alimentación, se obtienen diferentes grados de
rompimiento:

Tabla 4.1. Datos experimentales

FA0, mmol/h 300 1000 3000 5000

CA,OUT, 16 30 50 60
mmol/L

El volumen vacío del reactor es V=0.1 L, y a la temperatura del reactor la

concentración de la alimentación es CA0=100mmol/L. Encontrar una expresión de
velocidad que represente la reacción de craqueo.
Consideraciones:
 Se considera un reactor CSTR
 sistema de densidad variable porque está en estado gaseoso por lo tanto
varia Ftotal
Solución:
Comenzamos realizando la tabla estequiométrica:
Tabla 4.2. Tabla estequiométrica

Especie Fi0 Reacciona Fi 𝑭 𝑭

Mi =𝑭 𝒊𝟎 Ci = 𝑽𝒊 ,V=V0(1+ƐX)
𝑨𝟎
A FA0 - FA0 XA FA= FA0- FA0 XA FA= FA0(1- XA) 𝐹𝐴0 (1 − 𝑋𝐴 )
𝐶𝐴 =
𝑉0 (1 + ƐX𝐴 )
R FR0 𝑟
FA0 XA
𝑟
FR=FR0+𝑎FA0 FR=FA0(MR+𝑎XA)
𝑟 𝑟
𝑎
𝐹𝐴0 (𝑀𝑅 + 𝑎 𝑋𝐴 )
XA 𝐶𝑅 =
𝑉0 (1 + ƐX𝐴 )
𝑟
FT=FT0+FA0(𝑎 1)

Plantemos la ecuación de diseño de un reactor CSTR:

𝐹𝐴0 𝑋𝑓
𝑉= 𝐸𝑐. 4.1
−𝑟𝐴
Despejamos para la velocidad de reacción en función de A, sustituimos el valor de
V=0.1L
−𝑟𝐴 = 10𝐹𝐴0 𝑋𝑓 𝐸𝑐. 4.2

Ahora despejamos para la conversión:

𝐹𝐴0 (1 − 𝑋𝐴 )
𝐶𝐴 = 𝐸𝑐. 4.3
𝑉0 (1 + Ɛ𝑋𝐴 )
En donde
∈𝐴 = 𝑦𝐴0 𝛿 = (5 − 1)(1) = 4 𝐸𝑐. 4.4
𝐹𝐴0 − 𝐶𝐴0 𝑉0
𝑋𝐴 = 𝐸𝑐. 4.5
𝐹𝐴0 + 𝐶𝐴0 𝑉0 ∈
𝐶𝐴0 − 𝐶𝐴
𝑋𝐴 = 𝐸𝑐. 4.6
𝐶𝐴0 + 𝐶𝐴 ∈

Tabla 4.3. Conversiones y avances de reacción calculados

FA0, mmol/h 300 1000 3000 5000

CA,OUT, 16 30 50 60
mmol/L
XA 0.512 0.318 0.167 0.118
-rA 1536.6 3180 5000 5882.4

Linealizamos la ecuación de la velocidad de reacción

−𝑟𝐴 = 𝑘𝐶𝐴 ∝ 𝐸𝑐. 4.7

ln(−𝑟𝐴 ) = ln(𝑘) +∝ ln(𝐶𝐴 ) 𝐸𝑐. 4.8
A continuación, realizamos una regresión para encontrar los parámetros k y ∝:
 y = 1.0114x + 4.5644
Regresión lineal R² = 0.9951

8.8
8.6
8.4
8.2
ln(-rA)

8
7.8
7.6
7.4
7.2
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5
ln(CA)

Figura 4.1. Regresión lineal

∝= 1
𝑘 = 𝑒 4.5644 = 95.58ℎ−1
Problema 2-2 (Fogler)
¿Qué haría si:
a) necesitara estimar el volumen de reactor necesario para alcanzar una
conversión de 98% usando los datos de la tabla 2-1.
b) le pidieran utilizar los datos de la tabla 2-1 para calcular el volumen de un
reactor CSTR a una temperatura 100°C más alta que la temperatura a la
que se obtuvieron los datos?

Tabla de datos 2-1. Datos de velocidad

X -rA
(mol/dm3*s)
0 0.0053
0.1 0.0052
0.2 0.005
0.3 0.0045
0.4 0.004
0.5 0.0033
0.6 0.0025
0.7 0.0018
0.8 0.00125
0.85 0.001

Reacción
𝐴→𝐵+𝐶
−𝑟𝐴 = 𝑘𝐶𝐴 [1]

a) Se propone un tipo de reactor, en este caso se encontrará el volumen

utilizando las ecuaciones de diseño de un reactor CSTR y un PFR.

Primero, utilizando la estequiometría, se calculan las concentraciones finales

del reactivo, mediante la ecuación [2], obteniendo como resultados la tabla
1:

𝐶𝐴 = 𝐶𝐴0 (1 − 𝑥𝐴 ) [2]

Tabla de datos 1.
Concentración final
de A.
CA
0.1442
 0.1298
0.1154
0.1009
0.0865
0.0721
0.0577
0.0433
0.0288
0.0216

Mediante las concentraciones finales, se procede a hacer una regresión lineal para
encontrar el parámetro k, con la ecuación cinética [1], donde las variables
dependiente e independiente son las siguientes:

−𝑟𝐴 = 𝑘𝐶𝐴 [3]

y x
a1

(∑ 𝐶𝐴 ∗ −𝑟𝐴) [4]
𝑎1 =
∑ 𝐶𝐴2

Con la ecuación [2] se encuentra el valor de la constante de velocidad, que es el

siguiente:

𝑘 = 0.04151 1/𝑠

A continuación, para el cálculo del volumen del reactor, se utiliza la ecuación de

diseño del CSTR, [3]

𝐹𝐴0 𝑥𝐴 [5]
𝑣=
−𝑟𝐴
Donde:

𝐹𝐴0 = 𝐶𝐴0 𝑣0 [6]

Las mediciones que se encuentran en la tabla 1, se llevaron a cabo con una
Temperatura de 422.5 K y una presión total de 10 atm y que hay una carga inicial
equimolar del reactivo A e inertes, por lo que se puede calcular la concentración
inicial de A con la ecuación [7]:

𝑃𝐴 [7]
𝐶𝐴0 =
𝑅𝑇

Y sabemos que la presión parcial de una especie es igual a:

𝑃𝐴 = 𝑃𝑇 𝑦𝑖 [8]

Entonces, sustituyendo la ecuación [8] en la [7]:

𝑃𝑇 𝑦𝑖
𝐶𝐴0 =
𝑅𝑇

(10 𝑎𝑡𝑚)(0.5)
𝐶𝐴0 =
𝑑𝑚3 𝑎𝑡𝑚
(0.08206 )(422.5𝑘)
𝑚𝑜𝑙 𝑘

𝑚𝑜𝑙
𝐶𝐴0 = 0.1442
𝑑𝑚3

Las ecuaciones [2], [3] y [6] se sustituyen en la ecuación [5], para obtener el
volumen del reactor CSTR:

𝑣0 𝑥𝐴
𝑉=
𝑘(1 − 𝑥𝐴 )

𝑉𝐶𝑆𝑇𝑅 = 7082.63 𝑑𝑚3

Reactor PFR
Para el cálculo del volumen del reactor PFR, se utiliza su ecuación de diseño, la
cual es la siguiente:

[9]
 𝑥𝐴
𝑑𝑥𝐴
𝑉 = 𝐹𝐴0 ∫
0 −𝑟𝐴

Se sustituyen las ecuaciones [3] y [6] en la [9], y se integra con respecto a la

conversión de A, obteniendo la siguiente expresión para obtener el volumen:
𝑣0
𝑉=− ln(1 − 𝑥𝐴 )
𝑘
6𝑑𝑚3 /𝑠
𝑉𝑃𝐹𝑅 = (− ) (ln(1 − 0.98))
0.04151 1/𝑠

𝑉𝑃𝐹𝑅 = 565.457 𝑑𝑚3

1200

1000

800
1/-rA

600

400

200

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.85
X

Figura 1. Gráfico de representación de ecuaciones de diseño para reactor PFR.

b) Si se pidiera calcular el volumen de un reactor CSTR, con los mismos datos

de la tabla 1, pero ahora a 100°C más alta que la temperatura a la cual se
obtuvieron éstos, se podría utilizar la ecuación de Arrhenius([9])en la
ecuación cinética, para involucrar la nueva temperatura, y obtendríamos la
siguiente expresión [11]:
𝐸𝐴
𝑘 = 𝑘0 𝑒 −𝑅𝑇 [9]
 −𝑟𝐴 = 𝑘𝐶𝐴 [10]

𝐸𝐴
−𝑟𝐴 = 𝑘0 𝑒 −𝑅𝑇 𝐶𝐴 [11]

Sin embargo, la energía de activación (𝐸𝐴 ) y el factor preexponencial (𝑘0 ), serían

dos incógnitas, por lo cual se tendría que buscar otra ecuación ó, los valores de
éstas, para poder utilizar los mismos datos a una temperatura diferente.