Reactor CSTR
Reactor CSTR
Reactor CSTR
REACTORES CONTINUOS
Ecuaciones de Diseño
𝐹𝐴𝑜 − 𝐹𝐴 (2)
Para j = A 𝑉=
−𝑟𝐴
Caudal de alimentación
Donde: 𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜
𝐹𝐴 = 𝐶𝐴 𝑣
Caudal de salida
Tiempo espacial
(𝜏)
Es el tiempo necesario para tratar un volumen de alimentación igual al
volumen del reactor en condiciones determinadas (P, T, estado). Así un
tiempo espacial de 3 minutos significa que cada 3 minutos se trata en
el reactor un volumen de alimentación igual al volumen del reactor.
𝑡ҧ = 𝜏
CASO II
Para reacciones en fase gaseosa, debido a la estequiometria podría sufrir expansión o
contracción, luego. 𝑡ҧ ≠ 𝜏
CASO III
Para reactor PFR (flujo pistón) como la expansión no ocurre inmediatamente a la
entrada, ni a la salida, sino lo hace progresivamente, entonces:
𝑡ҧ ≅ 𝜏
Conclusiones
El tiempo de permanencia se determina de acuerdo a la cinética
particular del sistema.
Para el diseño de reactores continuos, se usara el tiempo espacial,
mientras que para sistemas discontinuos el tiempo de reacción.
𝐶𝐴 𝐶𝐴𝑜
En términos de conversión
𝐶𝐴𝑜 𝑋𝐴 )
𝜏=
−𝑟𝐴
𝑉 𝑋𝐴
= (5)
𝐹𝐴𝑜 −𝑟𝐴 𝜏
𝐴𝑟𝑒𝑎 =
1 𝐶𝐴𝑜
Representación grafica −𝑟𝐴
Para sistemas No
𝑛𝐴𝑜 (1 − 𝑋𝐴 ) 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
isotérmicas en fase 𝐶𝐴 = = (7)
𝑇 𝑇
gaseosa 𝑉𝑜 1 + 𝜀𝐴 𝑋𝐴 1 + 𝜀𝐴 𝑋𝐴
𝑇𝑜 𝑇𝑜
Las expresiones de diseño relacionan XA, -rA, V y FAo. Conociendo 3
de estas variables se obtiene la otra.
Para Reactores Continuos No Isotérmicos
Se requiere el balance de energía:
𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇
𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − + 𝑄ሶ = 0 (8)
𝑇 𝑎
𝑖=1
𝒎𝒐𝒍 𝑳
𝑪𝑨 ( ) −𝒓𝑨 𝒔−𝟏
𝑳 𝒎𝒐𝒍
0.15 0.0053
0.135 0.0052
0.120 0.0050
0.105 0.0045
0.09 0.0040
0.075 0.0033
0.060 0.0025
0.045 0.0018
0.030 0.00125
0.0225 0.0010
Ejercicio 1
Un reactante liquido se convierte al 50% en un reactor continuo de
mezcla perfecta. Se decide cambiar al reactor por otro también de
mezcla perfecta con un volumen 2 veces mayor que el original. Si se
mantiene constante la composición de alimentación, la temperatura,
y el caudal de alimentación. Determine la nueva conversión cuando la
reacción sea:
1ൗ
a) 𝑨՜𝑹 −𝑟𝐴 = 𝑘𝐶𝐴 2
Ejercicio 2
La siguiente reacción consecutiva ocurre en un reactor continuo CSTR
y es de primer orden:
𝑨 𝑩 𝑪 𝑫
𝒌𝟏 𝒌𝟐 𝒌𝟑
900
800
600
500
1/-RA
400
300
200
100
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
X
Ejemplo 2
Una corriente de alimentación CAo ingresa con un caudal Vo a un reactor CSTR
de volumen constante. La sustancia A reacciona consecutivamente mediante
una reacción de primer orden según:
𝐴 𝐵 𝐶
𝑘1 𝑘2
−𝐹𝐵 + 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝐵 𝑉 = 0
−𝐶𝐵 + 𝑘1 𝐶𝐴 𝜏 − 𝑘2 𝐶𝐵 𝜏 = 0 𝐶𝐵 1 + 𝑘2 𝜏 = 𝑘1 𝐶𝐴 𝜏
𝑘1 𝐶𝐴 𝜏 𝑘1 𝐶𝐴𝑜 𝜏
𝐶𝐵 = = (4)
1 + 𝑘2 𝜏 (1 + 𝑘1 𝜏)(1 + 𝑘2 𝜏)
𝜕𝐶𝐵
La CB será máxima cuando, =0
𝜕𝜏
`
Derivando (4): Identidad 𝑈 𝑈` 𝑈 `
= − 2𝑉
𝑉 𝑉 𝑉
Sea: 𝑈=𝜏
𝑉 = (1 + 𝑘1 𝜏)(1 + 𝑘2 𝜏)𝜏
𝜕𝐶𝐵 1 𝜏 𝑘1 (1 + 𝑘2 𝜏 + 𝑘2 (1 + 𝑘1 𝜏
= 𝑘1 𝐶𝐴𝑜 −
𝜕𝜏 (1 + 𝑘1 𝜏)(1 + 𝑘2 𝜏) (1 + 𝑘1 𝜏)2 (1 + 𝑘2 𝜏)2
1 𝜏 𝑘1 (1 + 𝑘2 𝜏 + 𝑘2 (1 + 𝑘1 𝜏
0= −
(1 + 𝑘1 𝜏)(1 + 𝑘2 𝜏) (1 + 𝑘2 𝜏)2 (1 + 𝑘2 𝜏)2
1 + 𝑘1 𝜏 1 + 𝑘2 𝜏 − 𝜏𝑘1 − 𝑘1 𝑘2 𝜏 2 − 𝜏𝑘2 − 𝑘1 𝑘2 𝜏 2
0=
(1 + 𝑘1 𝜏)2 (1 + 𝑘2 𝜏)2
0 = 1 + 𝑘1 𝜏 + 𝑘2 𝜏 + 𝑘1 𝑘2 𝜏 2 − 𝜏𝑘1 − 𝑘1 𝑘2 𝜏 2 − 𝜏𝑘2 − 𝑘1 𝑘2 𝜏 2
𝑘1 𝑘2 𝜏 2 = 1
1 (5)
𝜏=
𝑘1 𝑘2
Luego, 104
𝜏= = 41,523 𝑚𝑖𝑛.
5,8
5,8𝑥10−2 41,523 10
𝐶𝐵 = = 4,9933 𝑚𝑜𝑙/𝐿
(1 + 5,8𝑥10−2 (41,523))(1 + 1,0𝑥10−2 41,523 )
𝐶𝐴𝑜 = 𝐶𝐴 + 𝐶𝐵 + 𝐶𝐶
𝐶𝑐 = 𝐶𝐴𝑜 − 𝐶𝐵 − 𝐶𝐶
1 8.41 0.048
2 4.882 0.0398
3 1.8 0.025
4 0.75 0.0166
5 0.394 0.0125
Solución
Reactor CSTR, fase gaseosa, isotérmico
1.- Ecuación estequiometrica
𝑘 𝜀𝐴 = 𝛿𝐴 𝑦𝐴𝑜
2𝐴 ՜ 𝑅 𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋)
𝑡=0 𝑛𝐴𝑜 𝑦𝐴𝑜 = 1
𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
1 1−2
𝑡>0 𝑛𝐴𝑜 𝑋 𝑛𝐴𝑜 𝑋 𝑛𝐴𝑜 (1 − 𝑋) 𝛿𝐴 = = −0,5
2 𝐶𝐴 = 2
𝑉 (1 + 𝜀 𝑋)
𝐶𝐴𝑜 (1 − 𝑋)
𝐶𝐴 = (1)
(1 − 0,5𝑋)
𝐶𝐴𝑜 −𝐶𝐴 (2)
𝑋=
𝐶𝐴𝑜 −0,5𝐶𝐴
𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 = 0
𝐹𝐴𝑜 − 𝐹𝐴
𝑉=
−𝑟𝐴
1
recta
0.5 Lineal (recta)
ln(-rA)
Lineal (recta)
0
-5 -4 -3 -2 -1 0
-0.5
-1
ln(CA) -1.5
n= 2.00
𝑙𝑛𝑘 = 7.8213 𝑘 = 2493,15 L/mol h-1
Reactor CSTR en estado No estacionario
𝑦𝑒 𝑃 𝑥 𝑑𝑥
= න 𝑄 𝑥 𝑒𝑃 𝑥 𝑑𝑥
𝑑𝑥 + 𝐶 (4)
1 + 𝑘𝜏 Condiciones de contorno
𝑃 𝑥 =
𝜏
𝑡 = 0, 𝐶𝐴 = 𝐶𝐴𝑠𝑜
𝐶𝐴𝑜
𝑄 𝑥 = Concentración de
𝜏 salida en t = 0
Reemplazando
Reemplazando
1+𝑘𝜏 𝐶𝐴𝑜 1+𝑘𝜏
𝐶𝐴𝑜
𝐶𝐴 𝑒 𝜏 𝑡 =න 𝑒 𝜏 𝑡 𝑑𝑡 +𝐶
𝜏 𝐶 = 𝐶𝐴𝑠𝑜 −
1 + 𝑘𝜏
1+𝑘𝜏
𝑡 𝐶𝐴𝑜 1+𝑘𝜏
𝑡 Luego,
𝐶𝐴 𝑒 𝜏 = 𝑒 𝜏 +𝐶
1 + 𝑘𝜏
𝐶𝐴𝑜 𝐶𝐴𝑜 −
1+𝑘𝜏
𝐶𝐴𝑠 = − − 𝐶𝐴𝑠𝑜 𝑒 𝜏 𝑡
𝐶𝐴𝑜 −
1+𝑘𝜏
1 + 𝑘𝜏 1 + 𝑘𝜏
𝐶𝐴 = + 𝐶𝑒 𝜏 𝑡 (5) (6)
1 + 𝑘𝜏
Ecuación de la evolución de la concentración
a la salida del reactor a lo largo del tiempo
Casos: 𝐶𝐴𝑠)𝐸𝐸
1 + 𝑘𝜏
1.- Cuando 𝐶𝐴𝑠𝑜 = 0 arranque 𝑡 = ln
𝜏 𝐶𝐴𝑠)𝐸𝐸 − 𝐶𝐴𝑠
𝐶𝐴𝑜 −
1+𝑘𝜏 Suponiendo que el estado estacionario se
𝐶𝐴𝑠 = 1−𝑒 𝜏 𝑡 (7)
1 + 𝑘𝜏 alcanza cuando 𝐶𝐴𝑠 = 0.99𝐶𝐴𝑠)𝐸𝐸
−
1+𝑘𝜏 𝐶𝐴𝑠 (8)
𝑒 𝜏 𝑡 =1−
𝐶𝐴𝑠)𝐸𝐸
1 + 𝑘𝜏 𝐶𝐴𝑠)𝐸𝐸 − 𝐶𝐴𝑠
− 𝑡 = ln
𝜏 𝐶𝐴𝑠)𝐸𝐸
Para reacciones lentas
𝑡 = 4.605𝜏 (10)
4.605 (11)
𝑡=
𝑘
Para la mayoría de las reacciones de primer orden el estado estacionario se alcanza a
tiempos iguales o menores de 3 a 4 veces el 𝝉 del reactor.
Ejemplo:
En un reactor continuo de mezcla perfecta se lleva a cabo una reacción de primer
orden. 𝐴՜𝑅
Donde 𝒌 = 𝟎, 𝟓 𝒉−𝟏 ; 𝝉 = 𝟐 𝒉.
En estado estacionario, al reactor ingresa una alimentación con CAo = 30 mol/m3.
En un instante dado, la concentración de alimentación disminuye hasta 20 mol/m3 ,
determine:
a) La evolución a la concentración de salida a lo largo del tiempo.
b) La nueva concentración de salida en el estado estacionario.
c) El tiempo en que se alcanza el nuevo estado estacionario.
Solución: 𝑪𝑨𝒐
En el estado estacionario
La concentración de entrada:
𝑡՜∞
𝑪𝑨𝒐 = 𝟑𝟎 mol/m3
𝐶𝐴𝑜
𝐶𝐴𝑠)𝐸𝐸 =
𝑪𝑨 1 + 𝑘𝜏
30
La evolución de la concentración de A en la salida es: 𝐶𝐴𝑠)𝐸𝐸 = 1+0,5(2) = 15 mol/m3
𝐶𝐴𝑜 −
1+𝑘𝜏
𝑡
𝐶𝐴𝑠 = 1−𝑒 𝜏
1 + 𝑘𝜏
Cuando se inicia el cambio de la concentración de entrada.
𝐶𝐴𝑠)𝐸𝐸 = 15 = 𝐶𝐴𝑠𝑜
Luego, 𝐶𝐴𝑜 𝐶𝐴𝑜 1+𝑘𝜏
− 𝑡
𝐶𝐴𝑠 = − − 𝐶𝐴𝑠𝑜 𝑒 𝜏
1 + 𝑘𝜏 1 + 𝑘𝜏
20 20
𝐶𝐴𝑠 = − − 15 𝑒 −𝑡
1 + 0,5(2) 1 + 0,5(2)
20 20
𝐶𝐴𝑠 = − − 15 𝑒 −𝑡
1 + 0,5(2) 1 + 0,5(2)
𝐶𝐴𝑠 = 10 + 5𝑒 −𝑡
b) La nueva concentración de salida en el nuevo estado estacionario. 𝑡 ՜ ∞
𝐶𝐴𝑠)𝐸𝐸 = 10 mol/m3
𝑡 = 0 𝐶𝐴𝑜
1 1 a) Operación isotérmica a 25 ºC
𝑡 > 0 𝐶𝐴𝑜 𝑋 𝐶𝐴𝑜 𝑋 𝐶𝐴𝑜 𝑋
2 2 30500
18 −
𝑡 = 𝑡 𝐶𝐴 𝐶𝑅 𝐶𝑆 𝑘= 4.68𝑥10 𝑒 1.987(298) = 2𝑥10−4
𝐶𝐴 = 𝐶𝐴𝑜 (1 − 𝑋) Reemplazando en (1)
2.- Ecuación cinética
100 𝑋
=
−𝑟𝐴 = 𝑘𝐶𝐴 2 5 2𝑥10−4 (0,5) 1 − 𝑋 2
2
100 30500
18 − 1.987𝑇
4172 − 14𝑇 4172 − 14𝑇 (4)
4.68𝑥10 𝑒 0,5 1 − =
5 687 − 6,5𝑇 687 − 6,5𝑇
𝑇 = 338,63 𝐾
𝑋 = 0,3756
Ejemplo 2
La siguiente reacción: 𝐴 + 𝐵 ՜ 𝑃 ocurre en un reactor CSTR en fase liquida. El
volumen del reactor es de 1 L. La alimentación consiste de CAo = 4,2 y CBO = 21 mol/L
respectivamente a 20 ºC.
El caudal de alimentación es de 100 cc/min. Las capacidades caloríficas de las
soluciones son: CPA = 42,5 y CPB = 2,4 cal/mol ºC
La constante de velocidad esta dada por:
−19500
9
𝑘 = 2,9𝑥10 𝑒 𝑅𝑇 L/mol min-1
R = 1,987 cal/mol K
El calor de reacción es -21 kcal/mol y se supone constante. El tanque se encuentra
sumergido en un baño a 77 ºC, el área de contacto de intercambio es 250 cm2 y el
coeficiente global de transferencia de calor U = 0,1 cal/cm2min ºC.
Determine la conversión y la temperatura en el reactor.
Solución
Reactor CSTR, fase liquida, No isotermico
Reemplazando en (2).
−19500
51,156𝑥109 𝑒 1.987(308,281+94,585𝑋) 1 − 𝑋 5 − 𝑥 = 0,42𝑋
Resolviendo.
𝑋 = 0,825 𝑇 = 386,31 𝐾
C
X B.M Entre B y C el calor producido es superior al
B.E (exotérmica) eliminado por lo que el reactor no puede
funcionar continuamente de forma estable,
de manera que evoluciona hacia el punto C.
B Este punto B se denomina de encebamiento o
ignición ya que es la temperatura mínima a la
A que debe calentarse el reactor durante la
puesta en marcha para que pueda alcanzar
To To T por sus propios medios el punto estable C de
mayor conversión.
Ejemplo 3
Una reacción homogénea en fase liquida de primer orden se lleva a
cabo en un reactor CSTR. La concentración inicial del reactante es 3
mol/m3 y el caudal vo = 60x10-6 m3/s. La densidad y el calor especifico
de la mezcla reaccionante son constantes e iguales a 103 kg/m3 y 4.19
kJ/kg K respectivamente. El volumen del reactor es 18x10-3 m3. La
corriente de alimentación no contiene productos y el reactor opera
adiabáticamente. El calor de reacción es -2.09x108 J/mol y la ecuación
cinética es: −
62800
−𝑟𝐴 = 4.48𝑥106 𝑒 𝑅𝑇 𝐶𝐴 mol/m3 s-1
Solución
Reactor CSTR, fase liquida, Adiabatico
𝑘𝜏 𝑉 18𝑥10−3
𝑋= (1) 𝜏= = −6 = 3𝑥102 𝑠
1 + 𝑘𝜏 𝑣𝑜 60𝑥10
Luego, 62800
4.48𝑥106 𝑒 − 𝑅𝑇 (3𝑥102 )
𝑋= 62800
−
1+ 4.48𝑥10 𝑒 𝑅𝑇 (3𝑥102 )
6
62800
9 −
1,34𝑥10 𝑒 𝑅𝑇
𝑋= 62800 (2) Ecuación del B.M
9 − 𝑅𝑇
1+ 1,34𝑥10 𝑒
4.- Balance de energía
𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇
𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − =0
𝑇 𝑎
𝑖=1
62800
9 −
1.34𝑥10 𝑒 𝑅𝑇
𝑋𝐴 = 62800
B.M
9 −
1+ 1.34𝑥10 𝑒 𝑅𝑇
Ejercicio 1
Se desea producir 1000 kg/h de una solución de acido acético al 40%
(en peso) efectuando la hidrolisis del anhídrido acético en un reactor
CSTR operando a 50 ºC. Calcule:
a) El volumen del reactor para la conversión de 95%
b) La concentración inicial del anhídrido acético para lograr el
porcentaje de acido exigido a la salida.
c) Si el reactor opera adiabáticamente, a que temperatura deberá
alimentarse la solución de anhídrido acético.
Se dispone de los siguientes datos:
Para solución de A. acético al
𝐶𝐻3 𝐶𝑂 2 𝑂 + 𝐻2 𝑂 ՜ 2𝐶𝐻2 𝐶𝑂𝑂𝐻 40% a 50 ºC la densidad es
𝑀: 102 18 60 1,027 kg/L.
Alim. Prod.
T1,x1 T2,x2
Consiste de un tubo cilíndrico que opera en estado
estacionario, flujo altamente turbulento y el campo de
flujo se puede modelar como el flujo pistón.
El modelo de flujo tapón supone que no existe mezclado
ni difusión en la dirección del flujo e ignora los efectos de
entrada y salida al reactor. No existe un perfil radial de la
velocidad lineal, esto es , la velocidad es idéntica en todo
el elemento de volumen que forma el disco o pistón.
Fj(V) Fj(V+ΔV)
Reactantes Productos
Xj(V) Xj(V+ΔV)
𝑑𝐹𝐴
lim 𝑟𝐴 =
∆𝑉՜0 𝑑𝑉
𝑑𝐹𝐴
𝑟𝐴 = (1)
𝑑𝑉
En términos de concentración (volumen constante)
𝑣𝑜 𝑑𝐶𝐴
𝑑𝑉 =
𝑟𝐴
𝐶𝐴
Integrando 𝑣𝑜 𝑑𝐶𝐴
න 𝑑𝑉 = න
𝑉 𝐶𝐴𝑜 𝑟𝐴
𝐶𝐴
𝑉 𝑑𝐶𝐴
= 𝜏 = −න (2)
𝑣𝑜 𝐶𝐴𝑜 (−𝑟𝐴 )
Representación grafica 1
−𝑟𝐴 𝐴𝑟𝑒𝑎 = 𝜏
𝐶𝐴 𝐶𝐴𝑜
En términos de conversión
𝐶𝐴 = 𝐶𝐴𝑜 1 − 𝑋𝐴
𝐶𝐴
Reemplazando en (2) 𝑉 −𝐶𝐴𝑜 𝑑𝑋𝐴
= 𝜏 = −න
𝑣𝑜 𝐶𝐴𝑜 (−𝑟𝐴 )
𝑋𝐴
𝑉 𝑉 𝑑𝑋𝐴
= =න (3)
𝑉
𝑣𝑜 𝐶𝐴𝑜 𝐹𝐴𝑜 0 (−𝑟𝐴 ) 𝐴𝑟𝑒𝑎 =
𝐹𝐴𝑜
1
Representación grafica −𝑟𝐴
0 𝑋𝐴
𝑇
V = 𝑉𝑜 1 + 𝜀𝐴 𝑋𝐴
𝑇𝑜
𝐶𝐴𝑜 (1 − 𝑋𝐴 )
𝐶𝐴 =
𝑇
1 + 𝜀𝐴 𝑋𝐴
𝑇𝑜
REACTOR TUBULAR: Proceso Isotérmico
Ejemplo: 1
Para la reacción en fase gaseosa homogénea a 227 ºC
𝐴 ՜ 3𝑅
𝑘
La ecuación cinética es : −𝑟𝐴 = 0.012𝐶𝐴1/2
Calcule el tiempo espacial para lograr una conversión de 80% a
partir de una alimentación equimolar de reactante A e Inertes en
un reactor de flujo pistón que opera a 10 atm y 227 °C.
solución: 1
Reactor PFR, fase gaseosa, isotérmico
1.- Ecuación estequiometrica
𝐴 ՜ 3𝑅 Reactante Producto
𝑘
𝑡 = 0 𝑛𝐴𝑜
𝐶𝐴𝑜 (1 − 𝑋)
𝑡 > 0 𝑛𝐴𝑜 𝑋 𝐶𝐴 =
(1 + 𝜀𝐴 𝑋)
𝑡 = 𝑡 𝑛𝐴
𝜀𝐴 = 𝛿𝐴 𝑦𝐴𝑜 𝑦𝐴𝑜 = 0,5
𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋) 3−1
𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋) 𝛿𝐴 = =2
1
2- Ecuación cinética
1Τ
1ൗ 1Τ 1−𝑋 2
−𝑟𝐴 = 𝑘𝐶𝐴 2 −𝑟𝐴 = 𝑘𝐶𝐴𝑜 2
1+𝑋
3.- Ecuación de diseño
𝑋 1 1/2
𝜏 𝜕𝑋 𝐶𝐴𝑜 ൗ2 𝑋 1 + 𝑋
=න 𝜏= න 𝜕𝑋 (1)
𝐶𝐴𝑜 𝑜 −𝑟𝐴 𝑘 𝑜 1−𝑋
0.6 4 2 1,7
1.5
0.8 9 3
1
0.5
1/2 X
𝑋
1 + 𝜀𝐴 𝑋
න 𝜕𝑋 = 1,36
𝑜 1−𝑋
Método numérico
Se recurre a la regla de Simpson. Para 5 puntos:
grafico
3.5
𝑋
ℎ 3
න 𝑓(𝑥)𝜕𝑋 = 𝑓 + 4𝑓1 + 2𝑓2 + 4𝑓3 + 𝑓4
𝑜 3 𝑜 2.5
f(x)
1.5
0,8 − 0 1
ℎ= = 0,2
4 0.5
0
𝑋
0,2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
න 𝑓(𝑥)𝜕𝑋 = 1 + 4(1.227) + 2(1.528) + 4(2) + 3 X
𝑜 3
1/2
𝑋 (1 + 𝑋) 1+𝑋
X (1 − 𝑋) 1−𝑋
න 𝑓(𝑥)𝜕𝑋 = 1,331 0 1 1 𝑓𝑜
𝑜 0.2 1.5 1.227 𝑓1
0.4 2.3 1.528 𝑓2
Método analítico 0.6 4 2 𝑓3
0.8 9 3 𝑓4
𝑋 1/2 1/2
1+𝑋 1+𝑋
න 𝜕𝑋
𝑜 1−𝑋 1+𝑋
𝑋
𝑋
1+𝑋 1+𝑋
1/2 (1 + 𝑋)
න 𝑥 𝜕𝑋 න 𝜕𝑋
1−𝑋 1+𝑋 𝑜 1− 𝑋2
𝑜
𝑋 𝑋 𝑋
(1 + 𝑋) 𝜕𝑋 𝑋𝜕𝑋
න 𝜕𝑋 = න +න
𝑜 1 − 𝑋2 𝑜 1 − 𝑋2 𝑜 1 − 𝑋2
𝑋
(1 + 𝑋)
න 𝜕𝑋 = 𝑎𝑟𝑐 𝑠𝑒𝑛𝑋 − 1 − 𝑋 2
𝑜 1 − 𝑋2
Realizando expansión en serie,
1 ∗ 𝑋3 1 ∗ 3 ∗ 𝑋5
𝑎𝑟𝑐 𝑠𝑒𝑛𝑋 = 𝑋 + + +⋯
2(3) 2 4 5
Luego, 𝑋
(1 + 𝑋) 1 ∗ 𝑋3 1 ∗ 3 ∗ 𝑋5
න 𝜕𝑋 = 𝑋 + + − 1 − 𝑋2
𝑜 1 − 𝑋2 2(3) 2 4 5
𝑋
(1 + 𝑋) 1 ∗ 𝑋3 1 ∗ 3 ∗ 𝑋5
න 𝜕𝑋 = 𝑋 + + − 1 − 𝑋2
𝑜 1 − 𝑋2 2(3) 2 4 5
Tomando limites,
𝑋
(1 + 𝑋) 1 ∗ 𝑋3 1 ∗ 3 ∗ 𝑋5
න 𝜕𝑋 = 𝑋 + + + 1 − 1 − 𝑋2
𝑜 1 − 𝑋2 2(3) 2 4 5
Para X = 0,8, 𝑋
(1 + 𝑋) 0,83 3 ∗ 0,85
න 𝜕𝑋 = 0,8 + + + 1 − 1 − 0,82
𝑜 1 − 𝑋2 6 40
𝑋
(1 + 𝑋)
න 𝜕𝑋 = 1,31
𝑜 1− 𝑋2
𝑦𝐴𝑜 𝑃 0,5(10)
𝐶𝐴𝑜 = = = 0,122
𝑅𝑇 0,082(500)
Luego,
0,1220.5
𝜏= ∗ 1,31 = 38,13 𝑠
0,012
Ejemplo: 2
La siguiente reacción en fase gaseosa ocurre en un reactor PFR.
𝐴՜𝑅
𝑘
Para los mismos datos ejemplo 1 (reactor CSTR) determine el volumen
requerido para alcanzar conversiones de 80 y 60 %.
Solución
Reactor PFR, fase gaseosa, isotérmico
1.- Ecuación estequiometrica
𝒌
𝑨՜ 𝑹
𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋)
𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
1−1
𝛿𝐴 = =0
1
𝜀𝐴 = 0
0,6
𝑉 𝜕𝑋
= න
𝐹𝐴𝑜 0 −𝑟𝐴
a) 𝐴՜𝑅
𝑘
Reactante Producto
b) 𝐴 ՜ 𝑅+S
𝑘
c) 𝐴 ՜ 2𝑅 + 𝑆
𝑘
Solución
Reactor PFR, fase gaseosa, isotérmico
1.- Ecuación estequiometrica
𝐴՜𝑅 𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
𝑘
1−𝑋
𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋) 𝐶𝐴 = 𝐶𝐴𝑜
1 + 𝜀𝐴 𝑋
2.- Ecuación cinética
2
1−𝑋
−𝑟𝐴 = 𝑘𝐶𝐴 2 = 𝑘𝐶𝐴𝑜 2
1 + 𝜀𝐴 𝑋
3.- Ecuación de diseño
𝑋 𝑋 𝑋
𝑉 𝜕𝑋 𝜕𝑋 1 1 + 𝜀𝐴 𝑋 2 𝜕𝑋
=න =න 2 = 2න
𝐹𝐴𝑜 −𝑟𝐴 1−𝑋 𝑘𝐶𝐴𝑜 0 1−𝑋 2
0 0
𝑘𝐶𝐴𝑜 2
1 + 𝜀𝐴 𝑋
1 𝑋
(1 + 2𝜀𝐴 𝑋 +𝜀𝐴 2 𝑋 2 )𝜕𝑋
𝜏= න Identidades:
𝑘𝐶𝐴𝑜 0 1−𝑋 2 𝑋
𝜕𝑥 1 1
න 2
=−
𝑋 𝑎𝑥 + 𝑏 𝑎 𝑎𝑥 + 𝑏
𝜕𝑋 1 0
න 2
= 𝑋
𝑥𝜕𝑥 𝑏 1 1
0 1−𝑋 1−𝑋 න = + 𝑙𝑛 𝑎𝑥 + 𝑏
𝑎𝑥 + 𝑏 2 𝑎2 𝑎𝑥 + 𝑏 𝑎2
0
𝑋
𝑋𝜕𝑋 1
2𝜀𝐴 න = + ln(1 − 𝑋) 𝑋
𝑥 2 𝜕𝑥 𝑎𝑥 + 𝑏 𝑏2 1 2𝑏
1−𝑋 2 1−𝑋 න 2
= − − 𝑙𝑛 𝑎𝑥 + 𝑏
0 0 𝑎𝑥 + 𝑏 𝑎3 𝑎3 𝑎𝑥 + 𝑏 𝑎3
𝑋
𝜀𝐴 2 න 𝑋 2 𝜕𝑋 1
= − 1 − 𝑋 + + 2ln(1 − 𝑋)
0 1−𝑋 2 1−𝑋 𝑎 = −1
𝑏=1
Agrupando y factorizando
1 + 2𝜀𝐴 + 𝜀𝐴 2
+ 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 − 𝜀𝐴 2 (1 − 𝑋)
1−𝑋
Luego,
1 (1 + 𝜀𝐴 )2
𝜏= + 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 − 𝜀𝐴 2 (1 − 𝑋)
𝑘𝐶𝐴𝑜 1 − 𝑋
Tomando el limite cuando X = 0,
1 (1 + 𝜀𝐴 )2 𝑋 (1)
𝜏= + 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 + 𝜀𝐴 2 𝑋
𝑘𝐶𝐴𝑜 1−𝑋
𝑉 1 𝑣𝑜 1,5 1000
= 𝐿= 𝐿= 𝑥
𝑣𝑜 𝑘𝐶𝐴𝑜 𝑘𝐶𝐴𝑜 𝐴 0,2 0,75 0,1 1002
(1 + 𝜀𝐴 )2 𝑋
𝐿 = 10 𝑐𝑚 + 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 + 𝜀𝐴 2 𝑋 (2)
1−𝑋
Caso a) 𝐴՜𝑅 1−1
𝑘 𝛿𝐴 =
=0
1
0,8 𝜀𝐴 = 0
𝐿 = 10 𝑐𝑚 +0+0
0,2
𝐿 = 40 𝑐𝑚
Caso b) 𝐴 ՜ 𝑅+S
𝑘
1+1−1 𝑦𝐴𝑜 = 1 𝜀𝐴 = 1
𝛿𝐴 = =1
1
Luego,
(2)2 𝑋
𝐿 = 10 𝑐𝑚 + 2 2 ln 1 − 𝑋 + 𝑋
1−𝑋
4(0,8)
𝐿 = 10 𝑐𝑚 + 4 ln 0,2 + 0,8
0,2
𝐿 = 103,62 𝑐𝑚
Similarmente,
Caso c) 𝐴 ՜ 2𝑅+S
𝑘
𝐿 = 198,87 𝑐𝑚
Ejemplo: 4
En un reactor PFR isotérmico ocurre la siguiente reacción en fase
gaseosa a 250 °C.
𝑘1
𝐴 𝐵 + 3𝐶
𝑘2
𝑃𝐵 𝑃𝐶3
La ecuación cinética es: −𝑟𝐴 = 𝑘 𝑃𝐴 −
𝑚𝑜𝑙 𝑘𝑃
𝑘 = 1.747𝑥10−3
𝐿 − 𝑎𝑡𝑚 − 𝑠
𝑘𝑃 = 3.561 𝑎𝑡𝑚3
Para un caudal de 120 L/s y P = 1 atm. Calcule el volumen del reactor para
una conversión de 90% de la fracción de conversión en el equilibrio en las
siguientes condiciones: 𝑦 =1𝐴𝑜
𝑦𝐴𝑜 = 0.5 ; 𝑦𝐼𝑜 = 0.5
Solución
Reactor PFR, fase gaseosa, isotérmico
1.- Ecuación estequiometrica
𝑘1 𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋)
𝐴 𝐵 + 3𝐶
𝑘2 𝑛𝐵 = 𝑛𝐴𝑜 𝑋)
𝑡 = 0 𝑛𝐴𝑜
𝑛𝐶 = 3𝑛𝐴𝑜 𝑋)
𝑡 > 0 𝑛𝐴𝑜 𝑋 𝑛𝐴𝑜 𝑋 3𝑛𝐴𝑜 𝑋
𝑡 = 𝑡 𝑛𝐴 𝑛𝐵 𝑛𝐶 𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
1−𝑋
𝐶𝐴 = 𝐶𝐴𝑜
1 + 𝜀𝐴 𝑋
𝑋
𝐶𝐵 = 𝐶𝐴𝑜
1 + 𝜀𝐴 𝑋
𝑋
𝐶𝐶 = 3 𝐶𝐴𝑜
1 + 𝜀𝐴 𝑋
2.- Ecuación cinética
𝑃𝐵 𝑃𝐶 3 𝑃𝐴 = 𝐶𝐴 𝑅𝑇
−𝑟𝐴 = 𝑘 𝑃𝐴 −
𝐾𝑃
Reemplazando
𝐶𝐵 𝑅𝑇 (𝐶𝐶 𝑅𝑇) 3
−𝑟𝐴 = 𝑘 𝐶𝐴 𝑅𝑇 −
𝐾𝑃
𝐶𝐵 (𝐶𝐶 ) 3
−𝑟𝐴 = 𝑘𝑅𝑇 𝐶𝐴 −
𝐾𝑃
൘(𝑅𝑇)3
𝐶𝐵 (𝐶𝐶 ) 3
−𝑟𝐴 = 𝑘 ` 𝐶𝐴 − (1)
𝐾𝑐
𝐾𝑃
𝑘 ` = 𝑘𝑅𝑇 𝐾𝑐 = ൘(𝑅𝑇)3
𝑘 ` = 1,747𝑥10−3 0,082 523 = 0,075 𝑠 −1
3
3,561 −5
𝑚𝑜𝑙 1+3−1
𝐾𝑐 = = 4,5𝑥10 𝛿𝐴 = =3
(0,082𝑥523)3 𝐿 1
𝑦𝐴𝑜 =1
En términos de conversión
𝜀𝐴 = 3
1−𝑋 𝑋 3𝐶𝐴𝑜 X 3
−𝑟𝐴 = 𝑘 ` 𝐶𝐴𝑜 −
1 + 𝜀𝐴 𝑋 𝐾𝑐 (1 + 𝜀𝐴 𝑋) 1 + 𝜀𝐴 𝑋 3
1 𝑋 (1 + 3𝑋)𝜕𝑋
𝜏 = `න
𝑘 0 27(0,0233)3 𝑋 4 (3)
1−𝑋−
4,5𝑥10−5 (1 + 3𝑋)3
X F(x)
Resolviendo numéricamente,
0 1
0,20275 2,0250
0,811
ℎ 0,4055 3,8502
න 𝑓 𝑥 𝜕𝑋 = 𝑓𝑜 + 4𝑓1 + 2𝑓2 + 4𝑓3 + 𝑓4
0 3 0,60825 8,1710
0,811 31,8073
0,811
0,20275
න 𝑓 𝑥 𝜕𝑋 = 1 + 4 2,025 + 2 3,8502 + 4(8,171) + 31,8073
0 3
0,811
න 𝑓 𝑥 𝜕𝑋 = 5,494
0
5,494 𝑉 = 120 73,25 8,79 m3
Luego, 𝜏= = 73,25 𝑠
0,075
Resolviendo para : 𝑦𝐴𝑜 = 0,5
𝑋𝑒 = 0,9463
𝑋 = 0,85167
0,85167
න 𝑓 𝑥 𝜕𝑋 = 4,106
0
𝑉 = 6,569 m3
Ejercicio: 1
La descomposición en fase gaseosa del acetaldehído procede
de acuerdo con la siguiente reacción:
𝐶𝐻3 𝐶𝐻𝑂 ՜ 𝐶𝐻4 + 𝐶𝑂
Se desea producir 0.257 kg/día de metano a 518 °C y 1 atm. La
constante cinética es, 𝐿
𝑘 = 1.72
𝑚𝑜𝑙 − 𝑠
Determine el volumen del reactor para una alimentación de
acetaldehído puro y conversión de 80% operando en:
a) Reactor BATCH, con tiempo muerto de 3 horas aproximadamente
b) Reactor CSTR
c) Reactor PFR.
a) 300 L
Respuestas: b) 2214.46 L
c) 354 L.
Ejercicio: 2
Se desea realizar la siguiente reacción irreversible en fase
liquida: 𝐴 + 3𝐵 ՜ 2𝑅
En un reactor batch a 1 atm de presión y 60 °C constantes. La
reacción es de primer orden respecto a cada reactivo con un
valor de la constante cinética,
𝐿 𝑚𝑜𝑙
𝑘 = 0.25 𝐶𝐴𝑜 = 1 𝑦 𝐶𝐵𝑜 = 2 𝑚𝑜𝑙/𝐿
𝑚𝑜𝑙 − 𝑚𝑖𝑛 𝐿
Para una producción de 100 kmol/dia de producto R y
conversión de 90% del reactivo limitante determine:
a) Volumen del reactor batch si el tiempo muerto es de 900
segundos aproximadamente
b) El volumen del reactor su opera en forma continua
c) El volumen de un reactor tubular para el mismo propósito.
a) 1835.18 L
Respuestas: b) 5247.65 L
c) 967.12 L.
Ejercicio: 3
Se desea realizar la siguiente reacción irreversible en fase liquida:
𝐴՜𝑅
El orden de reacción no se conoce, por lo que se ha decidido estudiar la
cinética de reacción en un reactor batch, obteniéndose los siguientes datos
experimentales de concentración y velocidad a 30 ºC.,
CA (mol/L) -rA (mol/L min)
Calcule para una concentración inicial
de A, CAo = 1,3 mol/L. 0,1 0,1
0,3 0,5
a) El tiempo de reacción para 0,5 0,5
disminuir la concentración hasta 0,7 0,1
0,3 mol/L.
0,9 0,06
b) El volumen de un reactor CSTR
para lograr una conversión de 75% 1,1 0,0478
con una alimentación de 1000 1,3 0,045
mola/h. 2,0 0,042
c) El volumen de un reactor tubular
para el mismo propósito de b.
a) 12,72 min.
b) 25 L.
Respuestas:
c) 158,2 L.
REACTOR PFR: Proceso Adiabático
Ejemplo 5
Considere la siguiente reacción en fase gaseosa en un reactor tubular continuo.
𝐴+𝐵՜ 𝑅
𝑘
Se dispone como alimentación una mezcla equimolar de reactivos A y B a
500 °C y 1 atm de presión.
Determine :
a) El tiempo espacial para convertir 10 % de los reactivos, si la operación se
realiza isotérmicamente
b) El tiempo espacial y la temperatura si la operación es adiabática.
27500
Se dispone de los siguientes datos: 𝑘= 107.5 𝑒 − 𝑅𝑇 L/mol s-1
R en cal/mol K
∆𝐻𝑅 25 °𝐶 = −30 000 𝑐𝑎𝑙/𝑚𝑜𝑙
Las capacidades caloríficas en fase gas pueden suponerse constantes
en el rango de temperaturas utilizadas.
Componente Cp (cal/mol K)
A 36,8
B 20,2
R 59,5
Solución
Reactor PFR, fase gaseosa
1.- Ecuación estequiometrica
𝐴 + 𝐵 ՜ 𝑅
𝑘
𝑡 = 0 𝑛𝐴𝑜 𝑛𝐵𝑜
𝑡 > 0 𝑛𝐴𝑜 𝑋 𝑛𝐴𝑜 𝑋 𝑛𝐴𝑜 𝑋
𝑡=𝑡 𝑛𝐴 𝑛𝐵 𝑛𝑅
𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋)
𝑛𝐴𝑜 = 𝑛𝐵𝑜
𝑛𝐵 = 𝑛𝐵𝑜 − 𝑛𝐴𝑜 𝑋
𝑇
𝑛𝐵 = 𝑛𝐴𝑜 (1 − 𝑋) 𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
𝑇𝑜
2.- Ecuación cinética
𝑛𝐴𝑜 (1 − 𝑋) 𝑛𝐴𝑜 (1 − 𝑋)
−𝑟𝐴 = 𝑘𝐶𝐴 𝐶𝐵 = 𝑘
𝑉 𝑉
𝐶𝐴𝑜 2 (1 − 𝑋)2 𝑇𝑜 2
−𝑟𝐴 = 𝑘
(1 + 𝜀𝐴 𝑋) 2 𝑇
𝑋
𝑉 𝜕𝑋
3.- Ecuación de diseño =න
𝐹𝐴𝑜 0 −𝑟𝐴
𝑋
𝑉 𝜕𝑋
=න
𝐹𝐴𝑜 0 𝐶𝐴𝑜 2 (1 − 𝑋)2 𝑇𝑜 2
𝑘
(1 + 𝜀𝐴 𝑋) 2 𝑇
0,1
1 𝑇 2 (1 + 𝜀𝐴 𝑋) 2 𝜕𝑋
(1)
𝜏= 2න 𝑘 (1 − 𝑋) 2
𝐶𝐴𝑜 𝑇𝑜 0
a) Operación isotérmica
𝑇𝑜 = 450 + 273 = 723 𝐾
27500
−
𝑘= 107.5 𝑒 1,987(723) =0,153 L/mol s-1
𝑦𝐴𝑜 𝑃 0,5(1)
𝐶𝐴𝑜 = = = 8,43𝑥10−3 mol/L
𝑅𝑇 0,082(723)
Luego,
0,1
1 (1 + 𝜀𝐴 𝑋) 2 𝜕𝑋
𝜏= න
𝐶𝐴𝑜 𝑘 0 (1 − 𝑋) 2
Pero, el termino integral es:
(1 + 𝜀𝐴 )2 𝑋
+ 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 + 𝜀𝐴 2 𝑋
1−𝑋
(1 + 𝜀𝐴 )2 𝑋
1 (2)
𝜏= + 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 + 𝜀𝐴 2 𝑋
𝐶𝐴𝑜 𝑘 1−𝑋
1−2
𝛿𝐴 = = −1
1
𝜀𝐴 = −0,5
Luego,
1 (0,5)2 0,1
𝜏= −3 − 0,5 ln 0,9 + 0,52 (0,1)
8,43𝑥10 (0,153) 0,9
b) Operación adiabática
𝑇
∆𝐻𝑅 𝑇 = ∆𝐻𝑅 298 + න ∆𝐶𝑃 𝑑𝑇
298
0,1
ℎ
න 𝑓 𝑋 𝑑𝑋 = 𝑓 + 4𝑓1 + 2𝑓2 + 4𝑓3 + 𝑓4
0 3 𝑜
X T exp() f(x)
0 723 205798126 1.07577E+14
0.025 735.677985 147970575 8.21516E+13 Integral= 6.67612E+12
0.05 748.328228 107661699 6.35049E+13
0.075 760.95082 79218393.8 4.96658E+13
0.1 773.545852 58913622.9 3.92781E+13
6.67612𝑥1012
𝜏 = (107.5 )8,43𝑥10−3 (723)2 = 47,909 s
T vs X
780
770
760
T (K)
750
740
730
720
0 0.02 0.04 0.06 0.08 0.1 0.12
X
Ejemplo 6
𝑘1
Para la reacción elemental en fase liquida, 𝐴 𝐵
𝑘2
En el equilibrio: −𝑟 = 0 (1)
𝐴
1
1− 𝑋 =1
𝑘𝑐 𝑒
𝑘𝑐
𝑋𝑒 = (2)
1 + 𝑘𝑐
∆𝐻 1 1
− 𝑅 𝑇−298
Pero, 𝑘𝑐 (𝑇) = 𝑘𝑐 298 𝑒
50 𝑇 − 300 (4)
𝑋𝑒 =
20000
298 1 -0.005
1
300 1 0
320 1 0.05
0.8
340 1 0.1
360 0.99665 0.15 XBM
0.6
380 0.9856 0.2 XBE
Al reactor ingresa una alimentación de una mezcla que contiene 20% molar de A y 80%
molar de inertes a 227 ºC y 2 atm de presión. Para una conversión de 90% calcule:
a) El tiempo espacial para un reactor PFR isotérmico
b) El tiempo espacial y la temperatura de salida para un reactor PFR adiabático
c) El tiempo de reacción y la temperatura final para un reactor batch de volumen
constante. 9000
Solución
Reactor PFR, fase gaseosa:
1.- Ecuación estequiometrica
𝑛𝐴 (1 − 𝑋) 𝑇𝑜
4𝐴 ՜ 𝑅 + 2𝑆 𝐶𝐴 = = 𝐶𝐴𝑜
𝑉 (1 + 𝜀𝐴 𝑋) 𝑇
𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋)
𝑇
𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
𝑇𝑜
2.- Ecuación cinética
−𝑟𝐴 = 𝑘𝐶𝐴 2
Para PFR Isotérmico
3.- Ecuación de diseño
𝑋
𝑉 𝑋
𝜕𝑋 𝑉 1 + 𝜀𝐴 𝑋 2 𝜕𝑋
=න =න 2
𝐹𝐴𝑜 𝐹𝐴𝑜 0 𝑘𝐶𝐴𝑜 1 − 𝑋
2
0 −𝑟𝐴
𝑋
1 1 + 𝜀𝐴 𝑋 2 𝜕𝑋
𝜏= න (1)
𝑘𝐶𝐴𝑜 0 1−𝑋 2
Resolviendo la integración:
1 (1 + 𝜀𝐴 )2 𝑋
𝜏= + 2𝜀𝐴 1 + 𝜀𝐴 ln 1 − 𝑋 + 𝜀𝐴 2 𝑋 (2)
𝑘𝐶𝐴𝑜 1−𝑋
A 500 K,
9000
6 − 500
𝑘 = 5𝑥10 𝑒 = 0,07615 L/mol min-1
𝑦𝐴𝑜 𝑃 0,2(2)
𝐶𝐴𝑜 = = = 9,756𝑥10−3 mol/L
𝑅𝑇 0,082(500)
3−4
𝛿𝐴 = = −1/4
4
𝑦𝐴𝑜 = 0,2 𝜀𝐴 = −0,05
Reemplazando valores:
1 (1 − 0,05)2 0,9 2
𝜏= − 0,1 0,95 ln 0,1 + 0,05 (0,9)
0,07615(9,756𝑥10−3 ) 0,1
𝜏 = 11230,7 min
𝜏 = 187,178 h
Para PFR Adiabático
𝑋
1 𝑇 2 1 + 𝜀𝐴 𝑋 2 𝜕𝑋
𝜏= න
𝐶𝐴𝑜 𝑇𝑜 2 0 𝑘 1−𝑋 2
9000
𝑋 2
1 𝑇 𝑒 𝑇 1 + 𝜀𝐴 𝑋 2 𝜕𝑋 (3)
𝜏= 2 න
𝐶𝐴𝑜 𝑇𝑜 5𝑥10 6
0 1−𝑋 2
4.- Balance de energía:
𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇 (4)
𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − =0
𝑇 𝑎
𝑖=1
𝑇
∆𝐻𝑅 𝑇 = ∆𝐻𝑅 298 + න ∆𝐶𝑃 𝑑𝑇
298
∆𝐶𝑃 = 𝑣𝑖 𝐶𝑃𝑖 = 4 + 4 − 4 3 = −4
𝑖
∆𝐻𝑅 𝑇 = −15000 − 4(𝑇 − 298)
∆𝐻𝑅 𝑇 = −13808 − 4𝑇
∆𝐻𝑅 𝑇
= −3452 − 𝑇
𝑎
Luego, 𝑚 𝑇𝑜
𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 = (𝐹𝐴𝑜 𝐶𝑃𝐴 + 𝐹𝐼𝑜 𝐶𝑃𝐼 )(500 − 𝑇)
𝑖=1 𝑇
integral 1.00203E+13
1,00203𝑥1013
𝜏= = 821,67 min
5𝑥106 (0,07615) 9,756𝑥10−3 5002
𝜏 = 821,67 min
𝜏 = 13,69 h
Balance de energía:
𝑛𝐴𝑜 ∆𝐻𝑅 𝜕𝑋 𝑛𝐴𝑜 𝜕𝑇
+ 𝑛𝑖𝑜 + 𝑣𝑖 𝑋 𝐶𝑃𝑖 =0
𝑎 𝜕𝑡 𝑎 𝜕𝑡
𝑖
Desarrollando cada termino
𝑛𝑖𝑜 𝑣𝑖 𝑋 𝑋
+ 𝑋 𝐶𝑃𝑖 = 1 − 𝑋 𝐶𝑃𝐴 + 𝐶𝑃𝑅 + 2 𝐶𝑃𝑠
𝑛𝐴𝑜 𝑎 4 4
𝑖
𝑋 𝑋
1−𝑋 3+ 4+2 2=3−X
4 4
Reemplazando en la ecuación del balance de energía
(−13808 − 4T)
(3 − X)𝜕𝑇 = 𝜕𝑋
4
3 − X 𝜕𝑇 = (−3452 − 𝑇)𝜕𝑋
𝜕𝑋 𝜕𝑇
= 𝑋
𝜕𝑋 𝑇
𝜕𝑇
3 − 𝑋 −3452 − 𝑇 න =න
0 3 − 𝑋 500 −3452 − 𝑇
Integrando
3 3452 + 𝑇
𝑙𝑛 = 𝑙𝑛
3−𝑋 3952
3 3452 + 𝑇 3(3952)
= 𝑇= − 3452 (6)
3−𝑋 3952 3−𝑋
Resolviendo la integral numéricamente (ecuación 5)
X T exp() f(x)
0 500 65659969.1 65659969.1
0.3 939.111111 14523.6048 29640.0099
0.6 1488 423.429528 2646.43455
0.9 2193.71429 60.4992455 6049.92455
Integral 7398323.82
Reemplazando en (5)
7398323.82
𝑡 = 5𝑥106 (9,756𝑥10−3 ) = 151,667 min
𝑡 = 2,527 h
Ejercicio 4
En un reactor tubular adiabático, se lleva a cabo la reacción en fase
gaseosa:
𝐴 ՜ 2𝑅 + 𝑆
𝑘
Determine :
El volumen del reactor y la temperatura de salida para lograr una
conversión de 70%.
Se dispone de los siguientes datos:
10000
Componente A B R Inertes
Cp (cal/molK) 5 6 7 10