Facultad de

Ingeniera Química
Ingeniería de las
Reacciones II
Diseño de reactores no isotérmicos CSTR y PFR
Integrantes:
Dueñas Ricardo
Mantilla Arianna
Mendoza Sara
Portero María Fernanda
 Diseño de reactores no isotérmicos CSTR y PFR
0 Balance molar
1 𝑑𝑋 −𝑟𝐴
=
𝑑𝑉 𝐹𝐴0

02 Ley de velocidad
−𝑟𝐴 = 𝑘𝐶𝐴 La principal diferencia reside en el
método para evaluar la ecuación de
diseño cuando la temperatura varía a
03 Ecuación de Arrhenius
lo largo de un PFR o cuando se retira
𝐸 1 1 calor de un CSTR.
𝑘 = 𝑘1 exp −
𝑅 𝑇1 𝑇

04 Estequiometria
𝐹𝐴 = 𝐶𝐴 𝑣0
𝐹𝐴0 = 𝐶𝐴 𝑣0

𝐶𝐴 = 𝐶𝐴0 1 − 𝑋
 05 Combinando
𝐸 1 1
−𝑟𝐴 = 𝑘1 𝑒𝑥𝑝 − 𝐶𝐴0 1 − 𝑋
𝑅 𝑇1 𝑇

06 Cancelamos la concentración de entrada 𝐶𝐴0

𝑑𝑋 𝑘 1 − 𝑋
=
𝑑𝑉 𝑣0

07 Como T varia a lo largo del reactor, k también lo hace, algo que no ocurre en
reactores isotérmicos tubulares.

𝑑𝑋 𝐸 1 1 1−𝑋
= 𝑘1 𝑒𝑥𝑝 −
𝑑𝑉 𝑅 𝑇1 𝑇 𝑣0

08 Balance de energía

−∆𝐻𝑅𝑥
𝑇 = 𝑇0 + 𝑋
𝐶𝑃𝐴
 El balance de Energía
Primera ley de la termodinámica
Comenzaremos aplicando la primera ley de la termodinámica primero a un
sistema cerrado y después a un sistema abierto. Un sistema es cualquier porción
del universo limitada, en movimiento o estacionaria, que se elija para aplicar las
diversas ecuaciones termodinámicas. En un sistema cerrado, en el cual la masa no
atraviesa los límites de éste, el cambio de energía total del sistema, dE, es igual al
flujo de calor hacia el interior del mismo, dQ, menos el trabajo realizado por el
sistema sobre los alrededores, dW.
Balance de energía

𝛿 𝐸෠ = 𝛿𝑄 − 𝛿𝑊

La letra δ significa que δQ y δW no son diferenciales exactas de una función de estado.

 Reactor CSTR no Isotérmico
• Ecuación de Energía

𝐸𝑛𝑒𝑟𝑔í𝑎 𝑒𝑛 𝑙𝑎𝑠 𝑐𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒𝑠 𝐸𝑛𝑒𝑟𝑔í𝑎 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑖𝑑𝑎 𝑑𝑒 𝐸𝑛𝑒𝑟𝑔í𝑎 𝑒𝑛 𝑙𝑎𝑠 𝑐𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒𝑠

• 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑎𝑙 + 𝑙𝑜𝑠 𝑎𝑙𝑟𝑒𝑑𝑒𝑑𝑜𝑟𝑒𝑠 𝑎𝑙 = 𝑞𝑢𝑒 𝑠𝑎𝑙𝑒𝑛 𝑑𝑒𝑙 +
𝑒𝑙𝑒𝑚𝑒𝑛𝑡𝑟𝑜 𝑑𝑒 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑒𝑙𝑒𝑚𝑒𝑛𝑡𝑜 𝑑𝑒 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑒𝑙𝑒𝑚𝑒𝑛𝑡𝑟𝑜 𝑑𝑒 𝑣𝑜𝑙𝑢𝑚𝑒𝑛

𝐴𝑐𝑢𝑚𝑢𝑙𝑎𝑐𝑖ó𝑛 𝑑𝑒
𝑒𝑛𝑒𝑟𝑔í𝑎 𝑒𝑛 𝑒𝑙
𝑒𝑙𝑒𝑚𝑒𝑛𝑡𝑜 𝑑𝑒 𝑣𝑜𝑙𝑢𝑚𝑒𝑛

• Ecuación de Diseño

𝑉 𝑋𝐴
=
𝐹𝐴𝑜 −𝑟𝐴
 Si tomamos A como RL
ሶ
𝑄ሶ − ∆𝐻
෡𝑅𝑜 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝐹𝐴0 𝑋𝐴 = 𝐹𝐴0 ෍ 𝑀𝑖 𝐶𝑝𝑖 (𝑇 − 𝑇𝑖𝑜 )

𝑚𝑜𝑙𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠 𝑑𝑒 𝑖
𝑀𝑖 =
𝑚𝑜𝑙𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠 𝑑𝑒 𝑅𝐿

𝑣𝑖 𝑣𝑖
∆𝐶𝑝 = ෍ 𝐶𝑝 − ෍ 𝐶𝑝
𝑎 𝑎
𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑟𝑒𝑎𝑐𝑡𝑖𝑣𝑜𝑠
𝑣𝑖 𝑣𝑖
෡𝑅𝑜 =
∆𝐻 ෍ ෡𝑓𝑖
∆𝐻 𝑜
− ෍ ෡𝑓𝑖
∆𝐻 𝑜
𝑎 𝑎
𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑟𝑒𝑎𝑐𝑡𝑖𝑣𝑜𝑠

𝑄ሶ ሶ
෡𝑅 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝑋𝐴 = ෍ 𝑀𝑖 𝐶𝑝𝑖 (𝑇 − 𝑇𝑖𝑜 )
− ∆𝐻 𝑜
𝐹𝐴0
 Generalidades sobre balances de energía
𝑛 𝑛
𝑑 𝐸෠𝑠𝑖𝑠𝑡𝑒𝑚𝑎
𝑄ሶ − 𝑊ሶ 𝑆 + ෍ 𝐹𝑖0 𝐻𝑖0 − ෍ 𝐹𝑖 𝐻𝑖 =
𝑑𝑡
𝑖=1 𝑖=1

Relacionaremos la temperatura con la conversión y los flujos molares, así como

los parámetros del sistema, o sea el coeficiente total de transferencia de calor y el
área, Ua también la temperatura ambiente correspondiente Ta y el calor de
reacción ΔHRx.
1. Adiabático (𝑄ሶ = 0) CSTR, PFR, reactor intermitente o reactor empacado. La
ሶ = 0) A Cp constante y
relación entre la conversión XEB, y la temperatura (𝑊
ΔCp=0
σ Θ𝑖 𝐶𝑝𝑖 𝑇 − 𝑇0
𝑋𝐸𝐵 =
−Δ𝐻𝑅𝑋

−Δ𝐻𝑅𝑋 𝑋
𝑇 = 𝑇0 +
σ Θ𝑖 𝐶𝑃𝑖
 2. PFR/ reactor empacado con cambiador de calor

2.1. PFR en términos de conversión

𝑑𝑇 𝑈𝑎 𝑇𝑎 − 𝑇 + 𝑟𝐴 Δ𝐻𝑅𝑋 𝑇
=
𝑑𝑉 𝐹𝐴0 ቀσΘ𝑖 𝐶𝑝 + Δ𝐶𝑝 𝑋൯

3. Para reacciones múltiples en PFR

𝑑𝑇 𝑈𝑎 𝑇𝑎 − 𝑇 + σ 𝑟𝑖𝑗 Δ𝐻𝑅𝑋𝑖𝑗
=
𝑑𝑉 σ 𝐹𝑖 𝐶𝑝
𝑖
 PFR adiabático
𝑋 −Δ𝐻𝑅𝑥 𝑇𝑅 + σ Θ𝑖 𝐶𝑃 𝑖 𝑇0 + 𝑋Δ𝐶𝑝 𝑇𝑅
𝑇=
σ Θ𝑖 𝐶𝑃 𝑖 + 𝑋Δ𝐶𝑃

Con esta ecuación podemos obtener la

temperatura, la conversión y los perfiles de
concentración, junto con la longitud del reactor;
para construir una tabla de T en función de X.

En esta ecuación se acopla con el balance molar diferencial

𝑑𝑋
𝐹𝐴0 = −𝑟𝐴 𝑋 𝑇
𝑑𝑉
Algoritmo para PFR y el reactor
empacado adiabáticos

La reacción elemental reversible es fase gaseosa

𝐴՞𝐵

Se lleva a cabo en un PFR, en el cual la caída de presión es despreciable y

se alimenta A puro al reactor.

Balance molar:
𝑑𝑋 −𝑟𝐴
=
𝑑𝑉 𝐹𝐴0
Ley de velocidad
𝐶𝐵
−𝑟𝐴 = 𝑘 𝐶𝐴 −
𝐾𝐶
 Con 𝐸 1 1
𝑘 = 𝑘1 𝑇1 𝑒𝑥𝑝 −
𝑅 𝑇1 𝑇

Δ𝐻°𝑅𝑥 1 1
𝐾𝑐 = 𝑘𝑐2 𝑇2 𝑒𝑥𝑝 −
𝑅 𝑇2 𝑇

Estequiometria
Gas, 𝜀 = 0, 𝑝 = 𝑃0
𝑇0
𝐶𝐴 = 𝐶𝐴0 1 − 𝑋
𝑇
𝑇0
𝐶𝐵 = 𝐶𝐴0 𝑋
𝑇
Combinar
𝑋 𝑇0
−𝑟𝐴 = 𝑘𝐶𝐴0 1 − 𝑋 −
𝐾𝑐 𝑇
Balance de energía
● Para relacionar la temperatura y la conversión aplicamos el balance de energía a un
PFR adiabático. Si todas las especies entran a la misma temperatura 𝑇𝑖0 = 𝑇0

Resolvemos la siguiente ecuación con 𝑄ሶ = 0, 𝑊ሶ 𝑠 = 0 para obtener T en

función de la conversión.
σ Θ𝑖 𝐶𝑃 𝑖 𝑇 − 𝑇𝑖0
𝑋=
− Δ𝐻°𝑅𝑥 𝑇𝑅 + Δ𝐶𝑃 ൫𝑇 − 𝑇𝑅

Se tiene
𝑋 −Δ𝐻°𝑅𝑥 𝑇𝑅 + σ Θ𝑖 𝐶𝑃 𝑖 𝑇0 + 𝑋Δ𝐶𝑃 𝑇𝑅
𝑇=
σ Θ𝑖 𝐶𝑃 𝑖 + 𝑋Δ𝐶𝑝

Si A es puro y si Δ𝐶𝑝 = 0 entonces

𝑋 −Δ𝐻°𝑅𝑥 𝑇𝑅
𝑇 = 𝑇0 +
𝐶𝑃 𝐴
 Ejercicio CSTR
Se desea estudiar la reacción de un gas que se comporta como ideal:
𝐴 → 2𝑅. La alimentación consiste en una mezcla que contiene 90%
de A y 10% de inerte a 340 K y presión de 2 atm. El flujo total es de
2 L/s.
a. Suponiendo que el reactor opera en forma isotérmica y se desea
lograr una conversión del 50% determine el volumen del reactor
si es de mezcla completa.
b. Estime la cantidad de calor que es necesario retirar o suministrar
para mantener la temperatura en 340K para el reactor.
c. Suponga que se desea obtener una conversión del 90% y el
reactor opera en forma adiabática; determine la temperatura de
salida para el reactor de mezcla completa.
Datos complementarios Solución:

∆𝐻𝑅𝑜 =
1100𝑐𝑎𝑙 Considerando la ecuación de diseño para un reactor de mezcla
𝑔𝑚𝑜𝑙
1165
− 𝑇
completa
𝑘 = 150𝑒 , 𝑠 −1
𝑐𝑎𝑙 𝑉 𝑋𝐴
𝐶𝑝𝐴 = 20 =
𝑔𝑚𝑜𝑙 ∗ 𝐾 𝐹𝐴𝑂 −𝑟𝐴 ∣𝑋𝐴
𝑐𝑎𝑙
𝐶𝑝𝑅 = 15 −𝑟𝐴 = 𝑘 ∗ 𝐶𝐴
𝑔𝑚𝑜𝑙 ∗ 𝐾 𝐶𝐴 = 𝐶𝐴𝑂 ∗ 1 − 𝑋𝐴
𝑐𝑎𝑙
𝐶𝑝𝑖𝑛 = 10
𝑔𝑚𝑜𝑙 ∗ 𝐾 Sustituyendo
−𝒓𝑨 = 𝒌 ∗ 𝑪𝑨𝑶 ∗ (𝟏 − 𝑿𝑨 )
Usando Arrhenius
1165
−
𝑘 = 150𝑒 𝑇
𝑇 = 340𝐾
𝒌 = 𝟒, 𝟖𝟕𝟓𝟐 𝒔−𝟏
 Considerando Gases Ideales
𝑃𝑇 ∗ 𝑦𝐴𝑂 2𝑎𝑡𝑚 ∗ 0,9
𝐶𝐴𝑂 = =
𝑅∗𝑇 0,082 𝑎𝑡𝑚 ∗ 𝐿
∗ 340𝐾
𝑚𝑜𝑙 ∗ 𝐾
𝑪𝑨𝑶 = 𝟎, 𝟎𝟔𝟒𝟓𝒈𝒎𝒐𝒍/𝑳
Se sabe:
𝐿 𝑔𝑚𝑜𝑙
𝐹𝐴𝑂 = 𝐶𝐴𝑂 ∗ 𝑣𝑂 = 2 ∗ 0,0645
𝑠 𝐿
𝒈𝒎𝒐𝒍
𝑭𝑨𝑶 = 𝟎, 𝟏𝟑𝟎𝟖
𝒔
A partir de la ecuación de diseño
𝐹𝐴𝑂 ∗ 𝑋𝐴
𝑉=
−𝑟𝐴 ∣𝑋𝐴
𝐹𝐴𝑂 ∗ 𝑋𝐴
𝑉=
𝑘 ∗ 𝐶𝐴𝑂 ∗ (1 − 𝑋𝐴 )
𝑽 = 𝟎, 𝟒𝟏𝟓𝟗 𝑳

b) Considerando el balance de energía

𝑸ሗ
෡ 𝒐𝑹 + ∆𝑪𝒑 𝑻 − 𝑻𝑹 𝑿𝑨 = ෍ 𝑴𝒊 𝑪𝒑𝒊(𝑻 − 𝑻𝒊𝒐)
− ∆𝑯
𝑭𝑨𝑶

Se sabe que la temperatura de operación (T) es igual a la temperatura de alimentación

(Tio), Entonces:
𝑄ሗ
෡𝑅𝑜 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝑋𝐴 = 0
− ∆𝐻
𝐹𝐴𝑂
𝑸ሗ
= ∆𝑯 ෡ 𝒐𝑹 + ∆𝑪𝒑 𝑻 − 𝑻𝑹 𝑿𝑨
𝑭𝑨𝑶
 Tr=25 C=298,15 K

𝑐𝑎𝑙
∆𝐻𝑅𝑜 = 1100
𝑔𝑚𝑜𝑙
𝑣𝑖 𝑣𝑖
∆𝐶𝑝 = ෍ ∗ 𝑐𝑝 − ෍ ∗ 𝑐𝑝
𝑎 𝑎
𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑟𝑒𝑎𝑐𝑡𝑖𝑣𝑜𝑠
2 1 𝑐𝑎𝑙 𝑐𝑎𝑙
∆𝐶𝑝 = 𝐶𝑝𝑅 − 𝐶𝑝𝐴 = 2 ∗ 15 − 20
1 1 𝑔𝑚𝑜𝑙 ∗ 𝐾 𝑔𝑚𝑜𝑙 ∗ 𝐾
𝒄𝒂𝒍
∆𝑪𝒑 = 𝟏𝟎
𝒈𝒎𝒐𝒍 ∗ 𝑲

La energía que se suministra o retira del reactor es:

𝑄ሶ = 𝐹𝐴𝑂 ∗ ∆𝐻 ෡𝑅𝑜 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝑋𝐴
𝒄𝒂𝒍
𝑸ሶ = 𝟗𝟗, 𝟒𝟎𝟖
𝒔
c ) Considerando el balance de energía

𝑄ሗ
෡𝑅𝑜 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝑋𝐴 = ෍ 𝑀𝑖 𝐶𝑝𝑖 (𝑇 − 𝑇𝑖𝑜 )
− ∆𝐻
𝐹𝐴𝑂
Como se considera un sistema adiabático :
෡𝑅𝑜 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝑋𝐴 = ෍ 𝑀𝑖 𝐶𝑝𝑖 𝑇 − 𝑇𝑖𝑜
− ∆𝐻

Se sabe:
𝑐𝑎𝑙
∆𝐻𝑅𝑜 = 1100
𝑔𝑚𝑜𝑙
𝑐𝑎𝑙
∆𝐶𝑝 = 10
𝑔𝑚𝑜𝑙 ∗ 𝐾
𝑇𝑅 = 2998,15 𝐾
𝑋𝐴 = 0,9
𝑇𝑖𝑜 = 340 𝐾
Calculo de σ 𝑀𝑖 𝐶𝑝𝑖
෍ 𝑀𝑖 𝐶𝑝𝑖 = 𝑀𝐴 ∗ 𝐶𝑝𝐴 + 𝑀𝑖 ∗ 𝐶𝑝𝑖
𝑦𝐴𝑂
𝑀𝐴 = =1
𝑦𝐴𝑂
𝑦𝑖𝑂 0,1
𝑀𝐴 = = = 0,1111
𝑦𝐴𝑂 0,9
𝑐𝑎𝑙 𝑐𝑎𝑙
෍ 𝑀𝑖 𝐶𝑝𝑖 = 1 20 + 0,111 ∗ 10
𝑔𝑚𝑜𝑙 ∗ 𝐾 𝑔𝑚𝑜𝑙 ∗ 𝐾
𝑐𝑎𝑙
෍ 𝑀𝑖 𝐶𝑝𝑖 = 21,111
𝑔𝑚𝑜𝑙 ∗ 𝐾

෡𝑅𝑜 + ∆𝐶𝑝 𝑇 − 𝑇𝑅 𝑋𝐴 = ෍ 𝑀𝑖 𝐶𝑝𝑖 𝑇 − 𝑇𝑖𝑜

− ∆𝐻

− 1100 + 10 𝑇 − 298,15 0,9 = 21,111 𝑇 − 340

𝑻 = 𝟐𝟗𝟒, 𝟔𝟏 𝑲