Reactor CSTR

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 88

Facultad de Ingeniería Química

Maestría en Ingeniería Química

REACTORES CONTINUOS

Mg. Pablo Díaz Bravo


Reactores Continuos
El flujo de reactantes y productos se mantienen constantes. Operan en
estado estacionario. Pueden ser: reactor tanque agitado CSTR y reactor
tubular PFR.

REACTOR TANQUE AGITADO (Continuos Stirred Tank Reactor)


Alim.

Reactor de Mezclado completo


o perfectamente agitada

Prod. Su composición en cada


instante es la misma en todo el
tanque y es igual a la salida

Ecuaciones de Diseño

1.- Balance de materia: 𝐹𝑗𝑜 − 𝐹𝑗 + න 𝑟𝑗 𝑑𝑉 = 0 (1)


𝑣

Para mezclado perfecto න 𝑟𝑗 𝑑𝑉 = 𝑟𝑗 𝑉


𝑣
Reemplazando en el B. M. 𝐹𝑗𝑜 − 𝐹𝑗 + 𝑟𝑗 𝑉 = 0

𝐹𝐴𝑜 − 𝐹𝐴 (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.

Velocidad espacial (s)


Es el numero de volumen de alimentación en condiciones
determinadas que pueden tratarse en unidad de tiempo medidos en
volumen de reactor. Así una velocidad espacial de 5 h-1 significa que el
volumen de alimentación que se trata en el reactor cada hora, es 5
veces el volumen del reactor.
1
Luego: 𝑆=
𝜏
𝑚𝑜𝑙𝑒𝑠 𝐴 𝑞𝑢𝑒 𝑖𝑛𝑔𝑟𝑒𝑠𝑎𝑛
1 1 𝐶𝐴𝑜 𝑉 𝑉 𝑉 𝑑𝑒 𝑟𝑒𝑎𝑐𝑡𝑜𝑟
𝑝𝑜𝑟 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑎𝑙𝑖𝑚𝑒𝑛𝑡𝑎𝑐𝑖𝑜𝑛
𝜏= = = = =
𝑆 𝐹𝐴𝑜ൗ 𝐹𝐴𝑜 𝑣𝑜 𝑀𝑜𝑙𝑒𝑠 𝐴 𝑞𝑢𝑒 𝑖𝑛𝑔𝑟𝑒𝑠𝑎𝑛 𝑝𝑜𝑟 𝑢𝑛𝑖𝑑𝑎𝑑 𝑑𝑒 𝑡𝑖𝑒𝑚𝑝𝑜
𝐶𝐴𝑜 𝑉

Tiempo de residencia y tiempo de espacial


CASO I
Para sistemas en fase liquida, el tiempo de permanencia y el tiempo espacial son
idénticos debido a que los líquidos no sufren expansión.

𝑡ҧ = 𝜏

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.

Para reactores CSTR


𝑉
𝜏=
𝑣𝑜
Reemplazando en la ecuación (2)
𝑉 𝐹𝐴𝑜 − 𝐹𝐴
= 𝜏= (3)
𝑣𝑜 −𝑟𝐴 𝑣𝑜
En términos de concentración (volumen constante)
𝐹𝐴𝑜 − 𝐹𝐴 𝐶𝐴𝑜 − 𝐶𝐴
𝜏= = (4)
−𝑟𝐴 𝑣𝑜 −𝑟𝐴 𝐴𝑟𝑒𝑎 = 𝜏
1
𝐶𝐴 = 𝐶𝐴𝑜 1 − 𝑋𝐴 −𝑟𝐴
Representación grafica

𝐶𝐴 𝐶𝐴𝑜
En términos de conversión

𝐹𝐴𝑜 − 𝐹𝐴 𝐶𝐴𝑜 𝑣𝑜 − 𝐶𝐴 𝑣𝑜 𝐶𝐴𝑜 − 𝐶𝐴𝑜 (1 − 𝑋𝐴 )


𝜏= = =
−𝑟𝐴 𝑣𝑜 −𝑟𝐴 𝑣𝑜 −𝑟𝐴

𝐶𝐴𝑜 𝑋𝐴 )
𝜏=
−𝑟𝐴
𝑉 𝑋𝐴
= (5)
𝐹𝐴𝑜 −𝑟𝐴 𝜏
𝐴𝑟𝑒𝑎 =
1 𝐶𝐴𝑜
Representación grafica −𝑟𝐴

Para sistemas isotérmicas


en fase gaseosa 0 𝑋𝐴
𝑛𝐴𝑜 (1 − 𝑋𝐴 ) 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
𝐶𝐴 = = (6)
𝑉𝑜 1 + 𝜀𝐴 𝑋𝐴 1 + 𝜀𝐴 𝑋𝐴

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

Para reactor Adiabático:


𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − =0 (9)
𝑇 𝑎
𝑖=1

Para reactor Adiabático con multireacciones:


𝑚 𝑇𝑜 𝑚

෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 + ෍ 𝑟𝑖𝑗 ∆𝐻𝑅 𝑇 𝑖𝑗 +Q=0


𝑖=1 𝑇 𝑗=1
REACTOR CSTR: Proceso Isotérmico
Ejemplo 1
La conversión de 𝐴՜ 𝑅 se efectúa en fase gaseosa. Las
𝑘
mediciones de laboratorio muestran la velocidad de reacción en
función de la concentración en la tabla adjunta. La temperatura fue
de 406.5 K y la presión total de 10 atm. La carga inicial fue de una
mezcla equimolar de componentes A e Inertes con una alimentación
volumétrica de 5 L/s. Calcule:
a) Los volúmenes necesarios para alcanzar conversiones de 60 y 80
% en un reactor CSTR.
b) Grafique las áreas que darían los volúmenes para las
conversiones solicitadas.
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
Determine la concentración de C en condiciones estacionarias para
una concentración máxima de B en el producto, si : k1 = 5.8x10-2 min-1
y k2 = 1.0x10-2 min-1. CAo = 10 mol/L
Tabla de datos experimentales

𝒎𝒐𝒍 𝑳
𝑪𝑨 ( ) −𝒓𝑨 𝒔−𝟏
𝑳 𝒎𝒐𝒍
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

b) 𝑨՜𝑹 −𝑟𝐴 = 𝑘𝐶𝐴


c) 𝑨՜𝑹 −𝑟𝐴 = 𝑘𝐶𝐴 2

Ejercicio 2
La siguiente reacción consecutiva ocurre en un reactor continuo CSTR
y es de primer orden:
𝑨 𝑩 𝑪 𝑫
𝒌𝟏 𝒌𝟐 𝒌𝟑

Determine la concentración de C en condiciones estacionarias para


una concentración máxima de B en el producto, si : k1 = 0,5 min-1 , k2
= 0,1 min-1, y k3 = 0,05 min-1 , y CAo = 10 mol/L
Ejercicio 3
La descomposición reversible en fase gaseosa del reactante A ocurre a
temperatura y presión constantes según:
𝒌𝟏
𝑨 ⇌ 𝟐𝑹
𝒌𝟐
La alimentación consiste de componente A puro a 373 ºC y 5 atm. La constante de
equilibrio a esta temperatura es de 1 mol/dm3.
Calcule la conversión de A para un reactor CSTR de volumen constante con:
𝝉=𝟓𝒔 𝒌𝟏 = 𝟏 𝒔−𝟏
Ejercicio 4
La cinética de la reacción en fase gaseosa: 𝑨 + 𝟑𝑩 ՜ 𝑹 es:
−10000
−𝑟𝐴 = 0,8𝑒 𝑅𝑇 𝑃𝐴 1,25 𝑃𝐵 0,3
El proceso se realizara en un reactor de mezcla completa a la presión de 33.9
atm y 227 ºC. Se dispone de una corriente de alimentación molar de B/A = 5. Se
requiere producir 200 toneladas mol de C cada 250 días, con una operación de
24 horas diarias y conversión de 85%. ¿Cual será el volumen del reactor?
Para el volumen calculado anteriormente, ¿Cual será la temperatura que
deberá operar el reactor para obtener una conversión de 99,9%?.
Ejemplo 1
La conversión de 𝐴 ՜ 𝑘
𝑅 se efectúa en fase gaseosa. Las mediciones de
laboratorio muestran la velocidad de reacción en función de la concentración
en la tabla adjunta. La temperatura fue de 406.5 K y la presión total de 10 atm.
La carga inicial fue de una mezcla equimolar de componentes A e Inertes con
una alimentación volumétrica de 5 L/s. Calcule:
a) Los volúmenes necesarios para alcanzar conversiones de 60 y 80 % en un
reactor CSTR.
b) Grafique las áreas que darían los volúmenes para las conversiones
solicitadas.
Tabla de datos experimentales Solución
𝒎𝒐𝒍 𝑳
𝑪𝑨 ( ) −𝒓𝑨 𝒔−𝟏 Reactor CSTR, fase gaseosa, isotérmico
𝑳 𝒎𝒐𝒍
0.15 0.0053 1.- Ecuación estequiometrica
0.135 0.0052 𝒌
𝑨՜ 𝑹
0.120 0.0050
𝑡=0 𝑛𝐴𝑜
0.105 0.0045 𝑡>0 𝑛𝐴𝑜 𝑋 𝑛𝐴𝑜 𝑋
0.09 0.0040
𝑡=𝑡 𝑛𝐴 𝑛𝑅
0.075 0.0033
0.060 0.0025 𝑛𝐴 = 𝑛𝐴𝑜 (1 − 𝑋)
0.045 0.0018
𝑉 = 𝑉𝑜 (1 + 𝜀𝐴 𝑋)
0.030 0.00125
1−1
𝛿𝐴 = =0 𝜀𝐴 = 0
0.0225 0.0010 1
Luego, volumen constante
𝐶𝐴𝑜 − 𝐶𝐴 𝐶𝐴𝑜 =
𝑦𝐴𝑜 𝑃 0,5(10)
= 0,082(406,5) = 0,15 mol/L
𝑋= 𝑅𝑇
𝐶𝐴𝑜

2.- Ecuación de diseño


𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 = 0
Se requiere de X, y 1/-rA, luego:
𝐶𝐴𝑜 𝑣𝑜 − 𝐶𝐴 𝑣𝑜 = −𝑟𝐴 𝑉
𝑉 𝑋 (1) CA X -rA 1/-rA
=
𝐹𝐴𝑜 −𝑟𝐴 0.15 0 0.0053 188.679245

Para X = 0,8 0.135 0.1 0.0052 192.307692


0.12 0.2 0.005 200
𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜 = 0,15 5 = 0,75 mol/s
0.105 0.3 0.0045 222.222222
0.09 0.4 0.004 250
𝐹𝐴𝑜 𝑋
𝑉= = 0,75 0,8 800 = 480 𝐿
−𝑟𝐴 0.075 0.5 0.0033 303.030303
Para X = 0,6 0.06 0.6 0.0025 400

𝐹𝐴𝑜 𝑋 0.045 0.7 0.0018 555.555556


𝑉= = 0,75 0,6 400 = 180 𝐿 0.03 0.8 0.00125 800
−𝑟𝐴 0.0225 0.85 0.001 1000
Área = 640

900

800

Área = 240 700

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

Determine la concentración de C en condiciones estacionarias para una


concentración máxima de B en el producto, si : k1 = 5.8x10-2 min-1 y k2 = 1.0x10-2
min-1. CAo = 10 mol/L
Solución
Reactor CSTR, fase liquida, isotérmico

1.- Ecuación estequiometrica


3.- Ecuación de diseño
𝐴 𝐵 𝐶
𝑘1 𝑘2 Para A: 𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 = 0
2.- Ecuación cinética 𝐶𝐴𝑜 𝑣𝑜 − 𝐶𝐴 𝑣𝑜 = −𝑟𝐴 𝑉 = 𝑘1 𝐶𝐴 𝑉
−𝑟𝐴 = 𝑘1 𝐶𝐴 (1) 𝑉
𝐶𝐴𝑜 − 𝐶𝐴 = 𝑘1 𝐶𝐴 = 𝑘1 𝐶𝐴 𝜏
𝑣𝑜
Para B: 𝑟𝐵 = 𝑘1 𝐶𝐴
−𝑟𝐵 = 𝑘2 𝐶𝐵 𝐶𝐴𝑜 = 𝐶𝐴 (1 + 𝑘1 𝜏)
𝐶𝐴𝑜
𝑟𝐵 = 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝐵 (2) 𝐶𝐴 = (3)
1 + 𝑘1 𝜏
Para B: 𝐹𝐵𝑜 − 𝐹𝐵 + 𝑟𝐵 𝑉 = 0

−𝐹𝐵 + 𝑘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

Reemplazando en (3) y (4)


𝐶𝐴𝑜 10
𝐶𝐴 = = = 2,934 𝑚𝑜𝑙/𝐿
1 + 𝑘1 𝜏 1 + 5,8𝑥10−2 (41,523)

5,8𝑥10−2 41,523 10
𝐶𝐵 = = 4,9933 𝑚𝑜𝑙/𝐿
(1 + 5,8𝑥10−2 (41,523))(1 + 1,0𝑥10−2 41,523 )

𝐶𝐴𝑜 = 𝐶𝐴 + 𝐶𝐵 + 𝐶𝐶
𝐶𝑐 = 𝐶𝐴𝑜 − 𝐶𝐵 − 𝐶𝐶

𝐶𝑐 = 10 − 2,934 − 4,9933 = 2,073 𝑚𝑜𝑙/𝐿


Ejemplo 3
En un reactor CSTR de 0,1 Litros, ingresa con caudal constante de
alimentación un reactante gaseoso puro con concentración de 0,10 mol/L.
Luego, ocurre la dimerización según:
𝑘
2𝐴 ՜ 𝑅
Calcule la ecuación cinética de la reacción a partir de la siguiente información
obtenida experimentalmente para distintos caudales de alimentación.
Nº experimentos 𝒗𝒐 L/h 𝑪𝑨

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𝐶𝐴

2.- Ecuación cinética

−𝑟𝐴 = 𝑘𝐶𝐴 𝑛 𝐶𝐴𝑜 𝑣𝑜 𝑋


−𝑟𝐴 = (4)
Aplicando logaritmo 𝑉
ln(−𝑟𝐴 ) = 𝑙𝑛𝑘 + 𝑛𝑙𝑛𝐶𝐴 (3)
Se prepara la siguiente tabla:
Y = a + mX

2.- Ecuación de diseño

𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 = 0
𝐹𝐴𝑜 − 𝐹𝐴
𝑉=
−𝑟𝐴

𝐶𝐴𝑜 𝑣𝑜 − 𝐶𝐴𝑜 (1 − 𝑋)𝑣𝑜


−𝑟𝐴 =
𝑉
vo Ca X rA ln(rA) ln(A
8.41 0.048 0.68421053 5.75421053 1.74993185 -3.03655427
4.882 0.0398 0.75156055 3.6691186 1.29995147 -3.22388837
1.8 0.025 0.85714286 1.54285714 0.43363599 -3.68887945
0.75 0.0166 0.90948746 0.68211559 -0.38255614 -4.09835258
0.394 0.0125 0.93333333 0.36773333 -1.00039724 -4.38202663
recta
𝑦 = 2.0079𝑋 + 7.8213
2
y = 2.0079x + 7.8213
R² = 0.9991
1.5

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

El estado transiente o no estacionario en condiciones isotérmicas no es muy


frecuente pero si ocurre ocasionalmente en operaciones en marcha (arranque)
o bien cuando se cambia o interrumpe súbitamente el caudal y/o la
concentración de entrada al reactor.

La respuesta del reactor provoca que la concentración a la salida evoluciona


través de sucesivos estados transitorios hasta alcanzar un nuevo estado
estacionario.
El balance de moles resulta:
𝜕𝑛𝐴 (1)
𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 =
𝜕𝑡

Para reacciones de primer orden y densidad constante (líquidos).


𝑉𝜕𝐶𝐴
𝐶𝐴𝑜 𝑣𝑜 − 𝐶𝐴 𝑣𝑜 − (−𝑟𝐴 )𝑉 = 𝜕𝐶𝐴 1 + 𝑘𝜏 𝐶𝐴𝑜
𝜕𝑡 + 𝐶 = (2)
𝜏𝜕𝐶𝐴 𝜕𝑡 𝜏 𝐴
𝜏
𝐶𝐴𝑜 − 𝐶𝐴 − (𝑘𝐶𝐴 )𝜏 =
𝜕𝑡
𝜏𝜕𝐶𝐴
𝐶𝐴𝑜 − 𝐶𝐴 (1 + 𝑘𝜏) = Ec. Diferencial de primer orden,
𝜕𝑡 cuya solución es.
𝑦ƴ + 𝑃 𝑥 𝑦 = 𝑄(𝑥) (3)

𝑦𝑒 ‫𝑃 ׬‬ 𝑥 𝑑𝑥
= න 𝑄 𝑥 𝑒‫𝑃 ׬‬ 𝑥 𝑑𝑥
𝑑𝑥 + 𝐶 (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𝐶𝐴𝑠)𝐸𝐸

2.- Para alcanzar el estado estacionario Entonces, 1 + 𝑘𝜏 1


𝑡 = ln
𝐶𝐴𝑜 𝜏 0.01
𝐶𝐴𝑠)𝐸𝐸 = 𝜏
1 + 𝑘𝜏 𝑡 = ln100
1 + 𝑘𝜏
Luego,

1+𝑘𝜏 4.605𝜏
𝐶𝐴𝑠 = 𝐶𝐴𝑠)𝐸𝐸 1 − 𝑒 𝜏 𝑡 𝑡= (9)
1 + 𝑘𝜏
𝐶𝐴𝑠 −
1+𝑘𝜏 Tiempo requerido para alcanzar el estado
= 1−𝑒 𝜏 𝑡
𝐶𝐴𝑠)𝐸𝐸 estacionario.


1+𝑘𝜏 𝐶𝐴𝑠 (8)
𝑒 𝜏 𝑡 =1−
𝐶𝐴𝑠)𝐸𝐸
1 + 𝑘𝜏 𝐶𝐴𝑠)𝐸𝐸 − 𝐶𝐴𝑠
− 𝑡 = ln
𝜏 𝐶𝐴𝑠)𝐸𝐸
Para reacciones lentas
𝑡 = 4.605𝜏 (10)

Para reacciones rápidas

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

c) El tiempo en que se alcanza el nuevo estado estacionario.


4.605𝜏 4.605(2)
𝑡= = = 4.605 ℎ
1 + 𝑘𝜏 1 + 0,5(2)
REACTOR CSTR: Proceso Adiabático
Ejemplo 1
En un reactor de mezcla completa de 100 litros se llevara a cabo la
reacción en fase liquida según:
2𝐴 ՜ 𝑅 + 𝑆
𝑘
La concentración de alimentación de A puro a 25 °C es CAo = 0.5
mol/L y el caudal vo = 5 L/min. La constante cinética de la
reacción es:
30500
𝑘= 4.68𝑥1018 𝑒 − 𝑅𝑇 L/mol min-1
La ecuación cinética es, −𝑟𝐴 = 𝑘𝐶𝐴2
El calor de reacción a 25 °C es -2500 cal/mol. Las capacidades caloríficas
independientes de la temperatura son:
Componente Cp (cal/mol K)
A 14
R 10
S 5
Calcule:
a) La conversión para una operación isotérmica.
b) La conversión y la temperatura de salida para una operación adiabática
Solución Reactantes

Reactor CSTR, fase liquida.


1.- Ecuación estequiometrica
𝑘
2𝐴 ՜ 𝑅 + 𝑆 Productos

𝑡 = 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

3.- Ecuación de diseño 𝑋


= 0,002
𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 = 0 1−𝑋 2
𝐹𝐴𝑜 − 𝐹𝐴 𝑋 = 0,02
𝑉=
−𝑟𝐴
𝑉 𝑋 (1)
= 2
𝑣𝑜 𝑘𝐶𝐴𝑜 1 − 𝑋
b) Operación adiabática Reactantes

4.- Balance de energía


𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − + 𝑄ሶ = 0
𝑇 𝑎 Productos
𝑖=1
𝑚 𝑇𝑜 30500
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇 (2) 𝑘= 18
4.68𝑥10 𝑒

1.987(298) = 2𝑥10−4
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − =0
𝑇 𝑎
𝑖=1
Analizando los términos: 𝑚 𝑇𝑜
𝑇 ෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 = 𝐹𝐴𝑜 𝐶𝑃𝐴 (𝑇𝑜 − 𝑇)
𝑇
∆𝐻𝑅 𝑇 = ∆𝐻𝑅 298 + න ∆𝐶𝑃 𝜕𝑇 𝑖=𝐴
298
Reemplazando en (2)
∆𝐶𝑃 = ෍ 𝑣𝑖 𝐶𝑃𝑖
𝑖
𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 = (687 − 6,5𝑇)𝐹𝐴𝑜 𝑋

∆𝐶𝑃 = 𝐶𝑃𝑅 + 𝐶𝑃𝑠 − 2𝐶𝑃𝐴 14 298 − 𝑇 = (687 − 6,5𝑇) 𝑋


∆𝐶𝑃 = 10 + 5 − 2 14 = −13 4172 − 14𝑇 (3)
𝑇 𝑋=
∆𝐻𝑅 𝑇 = −2500 − න 13𝜕𝑇 687 − 6,5𝑇
298
∆𝐻𝑅 𝑇 = 1374 − 13𝑇 Reemplazando en (1)
∆𝐻𝑅 𝑇
= 687 − 6,5𝑇
𝑎
𝑉 𝑋
= 2
𝑣𝑜 𝑘𝐶𝐴𝑜 1 − 𝑋

2
100 30500
18 − 1.987𝑇
4172 − 14𝑇 4172 − 14𝑇 (4)
4.68𝑥10 𝑒 0,5 1 − =
5 687 − 6,5𝑇 687 − 6,5𝑇

Usando Excel, resulta

𝑇 = 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

1.- Ecuación estequiometrica


𝐴+𝐵 ՜𝑃
𝑡 = 0 𝐶𝐴𝑜 𝐶𝐵𝑜
𝑡 > 0 𝐶𝐴𝑜 𝑋 𝐶𝐴𝑜 𝑋 𝐶𝐴𝑜 𝑋 𝐶𝐵𝑜 21
𝜃𝐵 = = =5
𝑡 = 𝑡 𝐶𝐴 𝐶𝐵 𝐶𝑃 𝐶𝐴𝑜 4,2
𝐶𝐴 = 𝐶𝐴𝑜 (1 − 𝑋)
𝐶𝐵 = 𝐶𝐴𝑜 (𝜃𝐵 − 𝑋)
2.- Ecuación cinética
−𝑟𝐴 = 𝑘 𝐶𝐴 𝐶𝐵
3.- Ecuación de diseño
𝑉 𝑋
= (1)
𝐹𝐴𝑜 −𝑟𝐴
4,2 100
𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜 = = 0,42 𝑚𝑜𝑙/𝑚𝑖𝑛
1000
Luego,
0,42𝑋 (2)
𝑉= −19500
2,9𝑥109 𝑒 𝑅𝑇 4,2 2 (1 − 𝑋)(5 − 𝑥)

Se requiere el balance de energía


4.- Balance de energía:
𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇 𝑄ሶ = 𝑈𝐴(𝑇𝑎 − 𝑇)
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − + 𝑄ሶ = 0
𝑇 𝑎
𝑖=1
Desarrollando el termino sumatoria:
𝑚 𝑇𝑜
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 = 𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 + 𝐹𝐵𝑜 𝐶𝑃𝐵 𝑇𝑜 − 𝑇
𝑖=1 𝑇

Reemplazando en la ecuación del B.E


293 − 𝑇 𝐹𝐴𝑜 𝐶𝑃𝐴 + 𝜃𝐵 𝐹𝐴𝑜 𝐶𝑃𝐵 − 𝐹𝐴𝑜 𝑋 −21 + 𝑈𝐴 350 − 𝑇 = 0

𝐹𝐴𝑜 293 − 𝑇 42,5 + 5(24) + 21000𝐹𝐴𝑜 𝑋 + 0,1(250)(350 − 𝑇) = 0


0,42 293 − 𝑇 42,5 + 5(24) + 21000(0,42)𝑋 + 0,1(250)(350 − 𝑇) = 0
19997,25 − 68,25𝑇 + 8820𝑋 + 8750 − 25𝑇 = 0
−93,25𝑇 + 8820𝑋 + 28747,25 = 0

𝑇 = 308,281 + 94,585𝑋 (3)

Reemplazando en (2).
−19500
51,156𝑥109 𝑒 1.987(308,281+94,585𝑋) 1 − 𝑋 5 − 𝑥 = 0,42𝑋
Resolviendo.
𝑋 = 0,825 𝑇 = 386,31 𝐾

Si fuese isotérmico: 𝑋 = 0,00172


Estabilidad estática de un reactor CSTR adiabático

Los reactores químicos son unidades que presentan frecuentemente grandes


dificultades de control, en particular si las reacciones son rápidas y exotérmicas,
pudiendo ser peligrosas y potencialmente explosivas.

El que el fenómeno sea controlable depende de la velocidad de producción y de


eliminación de calor. Si la velocidad de producción de calor es superior a la de
eliminación, el reactor se irá calentando, pudiendo llegar a fundir y/o producirse
una explosión.
En consecuencia, se procederá a estudiar como es el calor producido sobre el
eliminado en cada uno de los posibles puntos estacionarios para determinar su
estabilidad (diagrama de Van Heerden)
Se observa que el balance de energía es una relación lineal entre T y X (representa
el calor eliminado o consumido). Esta relación lineal da los infinitos puntos de
trabajo posibles del reactor CSTR, siendo el punto real de trabajo la solución
simultanea de la ecuación (3) y el balance de materia (2).

Si incorporamos la ecuación de Arrhenius en el Balance de materia, la relación X


y T, toma una forma sigmoidal (representa el calor generado)

La representación grafica de ambas ecuaciones y de las soluciones al sistema se


representan en la figura.
Para reacciones exotérmicas puede
C
X B.M predecirse que serán mas probables
B.E (exotérmica) las intersecciones para bajas o altas
conversiones que para valores
intermedios. Mientras para reacciones
B endotérmicas la pendiente de la curva
B.E (endotérmica)
de consumo cambia de signo y es
única y la intersección se hace mas
A
factible en la zona central como se
To To T aprecia en la figura.
El punto A es un punto de funcionamiento estable pero de conversión baja. Entre
A y B el calor producido es inferior al eliminado por lo que el reactor no puede
funcionar continuamente de forma estable, de manera que evoluciona hacia el
punto A. En el punto B el calor eliminado es igual al producido, pero sin embargo
es inestable ya que una pequeña perturbación hará que el sistema evolucione
hacia A (enfriamiento) o hacia C (calentamiento). Si aumenta la temperatura se
tiene que Qp > Qe y el sistema se calentará más desplazándose hacia C a través
del correspondiente estado no estacionario, mientras que si disminuye
evolucionarà hacia A ya que se enfriarà más al ser Qp <Qe.

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

La energía de activación esta en SI. Si la corriente de alimentación


esta a 298 K ¿Cuáles son las conversiones y temperaturas en estado
estable en la corriente de productos?

Solución
Reactor CSTR, fase liquida, Adiabatico

1.- Ecuación estequiometrica


𝐴՜𝑅
𝐶𝐴 = 𝐶𝐴𝑜 (1 − 𝑋)
2- Ecuación cinética
−𝑟𝐴 = 𝑘𝐶𝐴
3.- Ecuación de diseño
𝑉 𝑋 𝜏 𝑋 1 − 𝑋 𝑘𝜏 = 𝑋
= =
𝐹𝐴𝑜 −𝑟𝐴 𝐶𝐴𝑜 𝑘𝐶𝐴𝑜 (1 − 𝑋)

𝑘𝜏 𝑉 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

𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 = 𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇


𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 = 𝐹𝐴𝑜 𝑋∆𝐻𝑅 𝑇
𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 𝑤𝐴𝑜 𝐶𝑠𝐴 𝑇𝑜 − 𝑇
𝑋= =
𝐹𝐴𝑜 ∆𝐻𝑅 𝑇 𝐹𝐴𝑜 ∆𝐻𝑅 𝑇
𝜌𝑣𝑜 𝐶𝑠𝐴 𝑇𝑜 − 𝑇
𝑋=
𝐹𝐴𝑜 ∆𝐻𝑅 𝑇
Reemplazando datos:
103 60𝑥10−6 4,19 𝑇𝑜 − 𝑇 1000
𝑋=
3 60𝑥10−6 (−2,09𝑥108 )
4,19106 𝑇𝑜 − 𝑇
𝑋=
−6,27𝑥108 )
𝑇 − 𝑇𝑜 (3)
𝑋= Ecuación del B.E
150
𝑇 − 𝑇𝑜
𝑋𝐴 = B.E
150

62800
9 −
1.34𝑥10 𝑒 𝑅𝑇
𝑋𝐴 = 62800
B.M
9 −
1+ 1.34𝑥10 𝑒 𝑅𝑇

T XBM XBE BM BE Polinómica (BM) Lineal (BE)

298 0.013 0 1.2

313 0.0424 0.1 1


328 0.1178 0.2
0.8
343 0.2677 0.3
358 0.4792 0.4 0.6
373 0.6825 0.5
0.4
380 0.8246 0.56
403 0.9066 0.7 0.2

418 0.95 0.8


0
433 0.9726 0.9 200 250 300 350 400 450 500

445 0.9827 0.98 -0.2


Los puntos de intercepto son A, B y C. Los puntos A y C son estables. A
otras temperatura intermedias no puede existir la operación estable. Si
se opera a temperaturas bajas aproximadamente a 300 K la conversión
seria de 1,5%. En cambio a Temperaturas superiores a 346 K se obtendrá
mayor conversión y a 445 K se obtiene una conversión de 98%.

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.

2595,109 𝐶𝑒𝑠𝑝 = 0,86 𝑐𝑎𝑙/𝑔º𝐶


𝑙𝑜𝑔𝑘 𝑚𝑖𝑛−1= 7,551 −
𝑇(𝐾)
𝑐𝑎𝑙
∆𝐻𝑅 (𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒) = −8950
𝑚𝑜𝑙
Ejercicio 2
En un reactor de mezcla completa de 10 litros de capacidad que
opera adiabáticamente se lleva a cabo una reacción en fase liquida a
densidad constante cuya estequiometria es:
2𝐴 ՜ 𝑅
𝑣𝑜 = 0,2 𝐿/𝑠 , 𝐶𝐴𝑜 = 2 𝑚𝑜𝑙/𝐿
Las condiciones de operación son: 𝑐𝑎𝑙
𝑇𝑜 = 25 º𝐶 , 𝐶 𝑃𝐴 = 14
𝑚𝑜𝑙 − 𝐾

El calor de reacción es constante ∆H = - 4600 cal/mol.


5 −12000/𝑅𝑇 −1
La constante cinética es: 𝑘 = 1,8𝑥10 𝑒 𝑠

Determine la temperatura y conversión a la salida del reactor.

Rpta: 0,933 y T = 451,2 K


REACTOR TUBULAR (Plug Flow Reactor)
Conocido como reactor de flujo pistón o tapón, flujo tubular ideal.

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.

En este reactor, la composición del fluido y la temperatura


varían en la coordenada de posición en la dirección del
flujo.
Para obtener la ecuación de diseño del reactor se dividirá en
varios sub volúmenes ∆𝑽, La velocidad de reacción se puede
considerar espacialmente uniforme.

Fj(V) Fj(V+ΔV)
Reactantes Productos
Xj(V) Xj(V+ΔV)

Realizando el balance de materia en el subvolumen


𝐹𝑗 𝑉 − 𝐹𝑗 𝑉 + ∆𝑉 + 𝑟𝑗 ∆𝑉 = 0 En un subvolumen ∆𝑽,
en estado estacionario,
𝐹𝑗 𝑉+∆𝑉 −𝐹𝑗 (𝑉)
𝑟𝑗 = ∆𝑉
න 𝑟𝑗 𝑑𝑉 = 𝑟𝑗 ∆𝑉
Para j = A 𝐹𝐴 𝑉+∆𝑉 −𝐹𝐴 (𝑉) ∆𝑣
𝑟𝐴 = ∆𝑉

𝑑𝐹𝐴
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 𝑋𝐴

La ecuación de diseño permite calcular el tamaño del reactor


conociendo el caudal, la conversión y la cinética de reacción.

Si la alimentación ingresa parcialmente convertida ,


entonces. 𝑋𝐴
𝜏 𝑉 𝑑𝑋𝐴
= =න (4)
𝐶𝐴𝑜 𝐹𝐴𝑜 𝑋𝐴𝑜 (−𝑟𝐴 )
Para sistemas isotérmicas en fase gaseosa
𝑛𝐴𝑜 (1 − 𝑋𝐴 ) 𝐶𝐴𝑜 (1 − 𝑋𝐴 )
𝐶𝐴 = =
𝑉𝑜 1 + 𝜀𝐴 𝑋𝐴 1 + 𝜀𝐴 𝑋𝐴

Para sistemas No isotérmicas en fase gaseosa: se requiere el


Balance de Energía

𝑇
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−𝑋 2
−𝑟𝐴 = 𝑘𝐶𝐴 2 −𝑟𝐴 = 𝑘𝐶𝐴𝑜 2
1+𝑋
3.- Ecuación de diseño
𝑋 1 1/2
𝜏 𝜕𝑋 𝐶𝐴𝑜 ൗ2 𝑋 1 + 𝑋
=න 𝜏= න 𝜕𝑋 (1)
𝐶𝐴𝑜 𝑜 −𝑟𝐴 𝑘 𝑜 1−𝑋

La integral puede evaluarse gráficamente, numéricamente y/o analíticamente.


Método grafico
Se confecciona la siguiente tabla:
grafico
1/2
(1 + 𝑋) 1+𝑋 3.5
X (1 − 𝑋) 1−𝑋
3
0 1 1
0.2 1.5 1.227 2.5
0.4 2.3 1.528 2
f(x)

0.6 4 2 1,7
1.5
0.8 9 3
1

0.5

Área = 1,7(0,8) = 1,36 0


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

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.15 5 = 0.75 𝑚𝑜𝑙/𝐿


Luego, volumen constante
𝐶𝐴𝑜 − 𝐶𝐴 𝐶𝐴𝑜 =
𝑦𝐴𝑜 𝑃 0,5(10)
= 0,082(406,5) = 0,15 mol/L
𝑋= 𝑅𝑇
𝐶𝐴𝑜

2.- Ecuación de diseño


0,8 Se requiere de X, y 1/-rA, luego:
𝑉 𝜕𝑋 (1)
=න CA X -rA 1/-rA
𝐹𝐴𝑜 0 −𝑟𝐴
0.15 0 0.0053 188.679245 𝑓𝑜
Para X = 0,8
𝑉 ℎ 0.135 0.1 0.0052 192.307692
= 𝑓 + 4𝑓1 + 2𝑓2 + 4𝑓3 + 𝑓4 0.12 0.2 0.005 200 𝑓1
𝐹𝐴𝑜 3 𝑜
0.105 0.3 0.0045 222.222222
𝑉 0,2 0.09 0.4 0.004 250 𝑓2
= 188,67 + 4 200 + 2 250 + 4 400 + 800
𝐹𝐴𝑜 3
0.075 0.5 0.0033 303.030303
𝑉 0.06 0.6 0.0025 400 𝑓3
= 259,245
𝐹𝐴𝑜
0.045 0.7 0.0018 555.555556
0.03 0.8 0.00125 800 𝑓4
𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜 = 0,15 5 = 0,75 mol/s
0.0225 0.85 0.001 1000
𝑉 = 𝐹𝐴𝑜 259,245 = 194,43 𝐿.
Para X = 0,6

0,6
𝑉 𝜕𝑋
= න
𝐹𝐴𝑜 0 −𝑟𝐴

𝑉 3ℎ Se requiere de X, y 1/-rA, luego:


= 𝑓 + 3𝑓1 + 3𝑓2 + 𝑓3
𝐹𝐴𝑜 8 𝑜
CA X -rA 1/-rA
𝑉 3(0,2) 0.15 0 0.0053 188.679245 𝑓𝑜
= 188,67 + 3 200 + 3 250 + 400
𝐹𝐴𝑜 8
0.135 0.1 0.0052 192.307692
0.12 0.2 0.005 200 𝑓1
𝑉
= 145,4
𝐹𝐴𝑜 0.105 0.3 0.0045 222.222222
0.09 0.4 0.004 250 𝑓2

𝑉 = 𝐹𝐴𝑜 145,4 = 109,1 𝐿. 0.075 0.5 0.0033 303.030303


0.06 0.6 0.0025 400 𝑓3
0.045 0.7 0.0018 555.555556
0.03 0.8 0.00125 800
0.0225 0.85 0.001 1000
Ejemplo: 3
En un reactor PFR se alimenta el reactante A en fase gas con un
caudal de 1.5 L/s. Si la reacción de descomposición es de 2do
orden con una constante .
𝐿 −1
𝑘 = 0.2 𝑠
𝑚𝑜𝑙
La concentración inicial es de 0.75 mol/L. Para una sección
transversal del tubo de 0.1 m2 determine la longitud del reactor para
lograr una conversión de 80% en los siguientes casos:

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−𝑋 27𝐶𝐴𝑜 3 𝑋 4 (2)


−𝑟𝐴 = 𝑘 ` 𝐶𝐴𝑜 −
1 + 3𝑋 𝐾𝑐 (1 + 3𝑋)4

En equilibrio: −𝑟𝐴 = 0 𝑦𝐴𝑜 𝑃 1


𝐶𝐴𝑜 = = = 0,0233
𝑅𝑇 0,082(523)
27𝐶𝐴𝑜 3 𝑋 4
1−𝑋−
4,5𝑥10−5 (1 + 3𝑋)3
𝑋𝑒 = 0,9015
Para X = 0,9 Xe 𝑋 = 0,811
3.- Ecuación de diseño 𝑋
𝑉 𝜕𝑋
=න
𝐹𝐴𝑜 0 −𝑟𝐴
𝑋
𝑉 1 𝜕𝑋
= න
𝐹𝐴𝑜 𝑘 ` 𝐶𝐴𝑜 0 1 − 𝑋 27𝐶𝐴𝑜 3 𝑋 4
1 + 3𝑋 − 𝐾𝑐 (1 + 3𝑋)4

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

𝜏 = 775,19 0,10545 = 81,76 s


𝜏 = 81,76 s

b) Operación adiabática

4.- Balance de energía:


𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − + 𝑄ሶ = 0 (3)
𝑇 𝑎
𝑖=1
Evaluando los términos:
𝑚 𝑇𝑜
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 = 𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 + 𝐹𝐵𝑜 𝐶𝑃𝐵 𝑇𝑜 − 𝑇
𝑖=1 𝑇
= 𝐹𝐴𝑜 (𝐶𝑃𝐴 +𝐶𝑃𝐵 ) 𝑇𝑜 − 𝑇
= 𝐹𝐴𝑜 (36,8 + 20,2) 723 − 𝑇 (2)

𝑇
∆𝐻𝑅 𝑇 = ∆𝐻𝑅 298 + න ∆𝐶𝑃 𝑑𝑇
298

∆𝐶𝑃 = ෍ 𝑣𝑖 𝐶𝑃𝑖 = 59,5 − 36,8 − 20,2 = 2,5


𝑖
𝑇
∆𝐻𝑅 𝑇 = −30000 + න 2,5𝑑𝑇
298
∆𝐻𝑅 𝑇 = −30000 + 2,5(𝑇 − 298)

∆𝐻𝑅 𝑇 = −30745 + 2,5𝑇


Reemplazando en el balance de energía (3)

𝐹𝐴𝑜 57 723 − 𝑇 = −30745 + 2,5𝑇 𝑋𝐹𝐴𝑜


41211 − 57𝑇 = −30745𝑋 + 2,5𝑇𝑋
41211 = −30745𝑋 + 57 + 2,5𝑋 𝑇

41211 + 30745𝑋 (4)


𝑇=
57 + 2,5𝑋
Reemplazando en el balance de materia (1)
27500
0,1 1,987𝑇 𝑇 2 𝜕𝑋
1 𝑒
𝜏= 2න 2
𝐶𝐴𝑜 107.5 𝑇𝑜 0 1−𝑋
1 + 𝜀𝐴 𝑋
27500
0,1 1,987𝑇 𝑇 2 𝜕𝑋
1 𝑒 (5)
𝜏= න 2
(107.5 )8,43𝑥10−3 (723)2 0 1−𝑋
1 − 0,5𝑋

Integrando numéricamente mediante la regla de Simpson

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

a) Grafique la conversión de equilibrio en función de la temperatura.


b) Determine la temperatura adiabática de equilibrio y la conversión cuando se
alimenta A puro a un reactor PFR a 300 K.
𝐻ሶ 𝐴 298 𝐾 = −40 000 𝑐𝑎𝑙/𝑚𝑜𝑙
Datos:
𝐻ሶ 𝐵 298 𝐾 = −60 000 𝑐𝑎𝑙/𝑚𝑜𝑙

𝐶𝑃𝐴 = 50 𝑐𝑎𝑙/𝑚𝑜𝑙 K 𝐶𝑃𝐵 = 50 𝑐𝑎𝑙/𝑚𝑜𝑙 K

𝐾𝑐 (298) = 100 000


Solución
Reactor PFR, fase gaseosa
1.- Ecuación estequiometrica
𝑘1
𝐴 𝐵 2.- Ecuación cinética
𝑘2
−𝑟𝐴 = 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝐵
𝐶𝐴 = 𝐶𝐴𝑜 (1 − 𝑋)
𝑋
𝐶𝐵 = 𝐶𝐴𝑜 𝑋 −𝑟𝐴 = 𝑘1 𝐶𝐴𝑜 1 − 𝑋 −
𝑘𝑐
1
−𝑟𝐴 = 𝑘1 𝐶𝐴𝑜 1 − 1 − 𝑋
𝑘𝑐

En el equilibrio: −𝑟 = 0 (1)
𝐴

1
1− 𝑋 =1
𝑘𝑐 𝑒

𝑘𝑐
𝑋𝑒 = (2)
1 + 𝑘𝑐
∆𝐻 1 1
− 𝑅 𝑇−298
Pero, 𝑘𝑐 (𝑇) = 𝑘𝑐 298 𝑒

∆𝐻𝑅 = σ𝑖 𝑣𝑖 𝐻ሶ 𝑓(𝑖) = −60 − −40 = −20 kcal/ml


20000 1 1
Luego, −
𝑘𝑐 (𝑇) = 100000𝑒 𝑅 𝑇 298
Reemplazando en (2),
20000 1 1

100000𝑒 𝑅 𝑇 298
𝑋𝑒 = (3)
20000 1 1

1+ 100000𝑒 𝑅 𝑇 298
3.- Balance de energía
𝑚 𝑇𝑜
𝐹𝐴𝑜 𝑋𝐴 ∆𝐻𝑅 𝑇
෍ 𝐹𝑖𝑜 න 𝐶𝑝𝑖 𝑑𝑇 − =0
𝑇 𝑎
𝑖=1

𝐹𝐴𝑜 𝐶𝑃𝐴 𝑇𝑜 − 𝑇 = 𝐹𝐴𝑜 𝑋(−20000)


50 𝑇 − 300 = 20000𝑋

50 𝑇 − 300 (4)
𝑋𝑒 =
20000

Se elabora la siguiente tabla


XA = 0.41 y T = 465 K

T XBM XBE 1.2

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

400 0.9478 0.25 Polinómica (XBM)


0.4
Lineal (XBE)
420 0.8457 0.3
440 0.6484 0.35 0.2
460 0.4055 0.4
480 0.21517 0.45 0
280 330 380 430 480
500 0.1059 0.5
-0.2
Ejemplo 7
Se desea realizar la siguiente reacción en fase gaseosa: 4𝐴 ՜ 𝑅 + 2𝑆

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

Datos: 𝑘 = 5𝑥106 𝑒 − 𝑇 L/mol min-1


∆𝐻(298) = −15000 cal/mol
Las capacidades caloríficas de los gases se suponen constantes en el rango de
temperaturas usadas:
Componente A Inertes R S
CP (cal/molK) 3 5 4 2

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 𝑇

𝐹𝐴𝑜 𝐶𝑃𝐴 + 𝜃𝐼𝑜 𝐶𝑃𝐼 500 − 𝑇 = −3452 − 𝑇 𝑋𝐹𝐴𝑜

3 + 4(5) 500 − 𝑇 = −3452 − 𝑇 𝑋

23 500 + 3452𝑋 (4)


𝑇=
23 − 𝑋

Integrando numéricamente usando la regla de Simpson


X T exp() f(x)
0 500 65659969.1 1.6415E+13
0.3 552.2290749 11966055.4 7.22547E+12
0.6 605.8571429 2827736.23 6.10384E+12
0.9 660.9411765 819906.454 3.26661E+13

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

Para reactor batch adiabático: Volumen constante


Ecuación cinética
𝜕𝐶𝐴 𝜕𝑋
−𝑟𝐴 = − = 𝑘𝐶𝐴 2 𝑑𝑡 =
𝑘𝐶𝐴𝑜 (1 − 𝑋)2
𝜕𝑡
0,9 9000
1 𝜕𝑋
𝑡= න 𝑒 𝑇 (5)
5𝑥106 𝐶𝐴𝑜 0 (1 − 𝑋)2

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𝑅 + 𝑆
𝑘

Se alimenta una corriente de 120 L/min que contiene 20% molar de


A y 80% molar de un gas inerte a 300 °C y 3 atm de presión.

Determine :
El volumen del reactor y la temperatura de salida para lograr una
conversión de 70%.
Se dispone de los siguientes datos:
10000

R en cal/mol K 𝑘= 7𝑥106 𝑒 − 𝑇 L/mol min-1

∆𝐻𝑅 300 °𝐶 = −13500 𝑐𝑎𝑙/𝑚𝑜𝑙


Las capacidades caloríficas en fase gas pueden suponerse constantes
en el rango de temperaturas utilizadas.
Componente Cp (cal/mol K)
A 3,2
R 4,2
S 1,5
I 2,3
Ejercicio 5
Si en el ejercicio 2 de reactor CSTR se reemplazara por un
reactor tubular. Determine el volumen del reactor para
obtener la misma conversión y temperatura de salida.

Respuesta: Aproximadamente 8,97 Ls


Taller 5
En un reactor ideal adiabático se llevara a cabo la reacción irreversible en
fase gaseosa: 𝐴 + 2𝐵 ՜ 𝑅
La alimentación consiste de una mezcla que contiene 10% molar de A,
30% molar de B y 60% molar de inertes a 227 ºC, 10 atm de presión y
caudal de 2,5 L/min. Para una conversión de 90% calcule:
a) El volumen y temperatura de salida para un reactor CSTR de volumen constante
b) El volumen y temperatura de salida para un reactor PFR

Se dispone de los siguientes datos:

El calor de reacción a 25 ºC es -19500 cal/mol


𝑘(400𝐾) = 0,01 L/mol min-1 𝐸𝑎 = 10 000 cal/mol

Las capacidades caloríficas se suponen constantes en el rango de temperaturas utilizadas.

Componente A B R Inertes
Cp (cal/molK) 5 6 7 10

También podría gustarte