Caida de Presion en Reactores Part1

Diseño de reactores.
Caída de presión
Ing. Daissy Lorena Restrepo Serna
Caída de presión en reactores
Fase liquida Fase gaseosa
• No hay cambios en las • Las concentraciones de los

concentraciones de los reactivos cambian en función
reactivos por efecto de la de la presión
presión • Cambios en la presión pueden
• Se pueden ignorar los ocasionar cambios en el
cambios de presión (P=Po) volumen del reactor
• P ≠ Po
𝜃𝑖 + 𝑣𝑖 𝑋 𝑃 𝑇𝑜
𝐶𝑖 = 𝐶𝐴𝑜 1
1 + 𝜀𝑋 𝑃𝑜 𝑇
𝐹𝑖𝑜
𝜃𝑖 = 𝜀 = 𝑦𝐴𝑜 𝛿 𝑣𝑖 = 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑒𝑠𝑡𝑒𝑞𝑢𝑖𝑜𝑚𝑒𝑡𝑟𝑖𝑐𝑜
𝐹𝐴𝑜
𝑑 𝑐 𝑏 𝑁𝐴𝑜
𝜀= + − −1
𝑎 𝑎 𝑎 𝑁𝑇𝑜
𝑐𝑎𝑚𝑏𝑖𝑜 𝑒𝑛 𝑒𝑙 𝑛ú𝑚𝑒𝑟𝑜 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑚𝑜𝑙𝑒𝑠 𝑝𝑎𝑟𝑎 𝑐𝑜𝑛𝑣𝑒𝑟𝑠𝑖ó𝑛 𝑡𝑜𝑡𝑎𝑙

𝜀 =
𝑀𝑜𝑙𝑒𝑠 𝑡𝑜𝑡𝑎𝑙𝑒𝑠 𝑎𝑙𝑖𝑚𝑒𝑛𝑡𝑎𝑑𝑎𝑠
𝐶𝑎𝑚𝑏𝑖𝑜 𝑒𝑛 𝑒𝑙 𝑛𝑢𝑚𝑒𝑟𝑜 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑚𝑜𝑙𝑒𝑠
𝛿=
𝑀𝑜𝑙𝑒𝑠 𝑑𝑒 𝑎 𝑞𝑢𝑒 𝑟𝑒𝑎𝑐𝑐𝑖𝑜𝑛𝑎𝑟𝑜𝑛
Determinación de la relación P/Po en función de la concentración de

catalizador para una reacción de segundo orden
Combinación de (5) con
2𝐴 → 𝐵 + 𝐶 (2) y asumiendo T=To
𝑑𝑋 2 2
𝐹𝐴𝑜 = −𝑟𝐴 2 𝑑𝑋 𝐶𝐴𝑜 1 − 𝑋 𝑃
𝑑𝑊 𝐹𝐴𝑜 =𝑘 6
𝑑𝑊 1 − 𝜀𝑋 𝑃𝑜
−𝑟𝐴 = 𝑘𝐶𝐴2 3 2 2
𝑑𝑋 𝑘𝐶𝐴𝑜 1 − 𝑋 𝑃
= (7)
𝐶𝐴𝑜 1 − 𝑋 𝑃 𝑇𝑜 𝑑𝑊 𝑣𝑜 1 + 𝜀𝑋 𝑃𝑜
𝐶𝐴 = 4
1 + 𝜀𝑋 𝑃𝑜 𝑇 Dividiendo entre 𝐹𝐴𝑜
(𝑣𝑜 CAo ) 𝑑𝑋
2
= 𝐹1 𝑋, 𝑃 (8)
𝐶𝐴𝑜 1 − 𝑋 𝑃 𝑇𝑜 𝑑𝑊
−𝑟𝐴 = 𝑘 5 Falta una ecuación que correlacione la
1 − 𝜀𝑋 𝑃𝑜 𝑇 presión con el peso del catalizador
(P=f(w))
Caída de presión en reactores de lecho empacado
Ecuación de Ergun
𝑑𝑃 𝐺 1−𝜑 150 1 − 𝜑 𝜇
=− + 1.75𝐺 (9)
𝑑𝑧 𝜌𝑔𝑐 𝐷𝑝 𝜑 3 𝐷𝑝
𝜇= viscosidad del gas que pasa por el
P = presión, lbf /ft2 (kPa)
lecho, lbm/pies h (kg/m s)
𝜑 = porosidad =
𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒 𝑣𝑎𝑐í𝑜 z = longitud o largo del lecho empacado de
ൗ𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑔𝑙𝑜𝑏𝑎𝑙 𝑑𝑒𝑙 𝑙𝑒𝑐ℎ𝑜 tubería, ft (m)
𝑔𝑐 = 32.174 lbm ft/s2 lbf u = velocidad superficial ft/h (m/s)
𝑔𝑐 = 4.17 x 108 lbm ft/ h2 lbf 𝜌 = densidad del gas, lbm/ft3 (kg/m3)
(recuerde que en el sistema métrico gc = 1.0)
G = pu = densidad de flujo másico,
Dp = diámetro de partículas en el lecho, ft (m) lbm/ft2 h (kg/m2s)
Ecuación de Ergun
Por la ecuación de continuidad y la ecuación de gases ideales
𝑚ሶ 𝑜 = 𝑚ሶ
𝜌𝑜 𝑣𝑜 = 𝜌𝑣
𝑃𝑜 𝑇 𝐹𝑇
𝑣 = 𝑣𝑜 10
𝑃 𝑇𝑜 𝐹𝑇𝑜
𝑣𝑜 𝑃 𝑇𝑜 𝐹𝑇𝑜
𝜌 = 𝜌𝑜 = 𝜌𝑜 (11)
𝑣 𝑃𝑜 𝑇 𝐹𝑇
Ecuación de Ergun
Corrigiendo la densidad del gas en la Ecuación de Ergun
𝑑𝑃 𝐺 1−𝜑 150 1 − 𝜑 𝜇 𝑃𝑜 𝑇 𝐹𝑇
=− + 1.75𝐺 (12)
𝑑𝑧 𝜌𝑜 𝑔𝑐 𝐷𝑝 𝜑 3 𝐷𝑝 𝑃 𝑇𝑜 𝐹𝑇𝑜
Simplificando
𝑑𝑃 𝑃𝑜 𝑇 𝐹𝑇
= −𝐵𝑜 13
𝑑𝑧 𝑃 𝑇𝑜 𝐹𝑇𝑜
𝐺 1−𝜑 150 1 − 𝜑 𝜇
𝐵𝑜 = + 1.75𝐺 (14)
𝜌𝑜 𝑔𝑐 𝐷𝑝 𝜑 3 𝐷𝑝
Ecuación de Ergun
Con el propósito de expresar la caída del presión en función del catalizador se
introducen las ecuaciones 15 y 16
𝑊 = 1 − 𝜑 𝐴𝑐 𝑧 ∗ 𝜌𝑐 15 𝜌𝑏 = 𝜌𝑐 1 − 𝜑 16
Reemplazando 15 y 16 en 13 y haciendo y = P/Po:
𝑑𝑦 𝛼 𝑇 𝐹𝑇 2𝐵𝑜
=− 17 𝛼= (18)
𝑑𝑊 2𝑦 𝑇𝑜 𝐹𝑇𝑜 𝐴𝑐 𝜌𝑐 1 − 𝜑 𝑃𝑜
La ecuación 17 se aplica para caída de presión cuando ocurren

múltiples reacciones o se esta trabajando un reactor de membrana
𝑑𝑦 𝛼 𝑇 La ecuación 19 para reacciones únicas en

= − (1 + 𝜀𝑋) 19 lechos empacados
𝑑𝑊 2𝑦 𝑇𝑜
Ecuación de Ergun – solución analítica
𝑑𝑦 𝛼 𝑇
= − (1 + 𝜀𝑋) 19
𝑑𝑊 2𝑦 𝑇𝑜
𝜀 = 0 ∴ 1 ≫ 𝜀𝑋 𝑇 = 𝑇𝑜
Reacomodando 2𝑦
𝑑𝑦
= −𝛼 Integrando y=1 (P=Po) en W=0
𝑑𝑊
𝑦 2 = 1 − 𝛼𝑊
Introduciendo a la derivada
𝑑𝑦 2 Sacando raíz a cada lado

= −𝛼
𝑑𝑊
𝑃 1/2
Expresando la caída de presión en términos de la 𝑦 = = 1 − 𝛼𝑊 (21)
𝑃𝑜
longitud del reactor (z)
1/2
𝑃 2𝐵𝑜 𝑧
𝑦= = 1− (22)
𝑃𝑜 𝑃𝑜
Ecuación de Ergun – solución analítica – ecuación de segundo orden
Balance molar Combinando
𝑑𝑋 2
𝐹𝐴𝑜 = −𝑟𝐴 2 𝑑𝑋 𝑘𝐶𝐴𝑜 2 1/2
2
𝑑𝑊 = 1−𝑋 1 − 𝛼𝑊
𝑑𝑊 𝐹𝐴𝑜
Ley de velocidad
Separando variables
−𝑟𝐴 = 𝑘𝐶𝐴2 3
𝐹𝐴𝑜 𝑑𝑋
2 1−𝑋 2
= 1 − 𝛼𝑊 𝑑𝑊
𝑘𝐶𝐴𝑜
Estequiometria 𝑣 = 𝑣𝑜 /𝑦 , T=To, 𝜀 = 0
𝐹𝐴 Integrando con X=0 cuando W=0 y

𝐶𝐴 = = 𝐶𝐴𝑜 1 − 𝑋 𝑦 23 sustituyendo 𝐹𝐴𝑜 = 𝐶𝐴𝑜 𝑣𝑜
𝑣
𝑦 = 1 − 𝛼𝑊 1/2 (21) 𝑣𝑜 𝑋 𝛼𝑊
=𝑊 1− (25)
𝑘𝐶𝐴𝑜 1 − 𝑋 2
1/2
𝐶𝐴 = 𝐶𝐴𝑜 1 − 𝑋 1 − 𝛼𝑊 (24)
Ecuación de Ergun – solución analítica – ecuación de segundo orden
𝑣𝑜 𝑋 𝛼𝑊
=𝑊 1− (25)
𝑘𝐶𝐴𝑜 1 − 𝑋 2
Despejando conversión Despejando peso del catalizador
𝑘𝐶𝐴𝑜 𝑊 𝛼𝑊 2𝑣𝑜 𝛼 𝑋 1/2

1− 1− 1−
𝑣𝑜 2 𝑘𝐶𝐴𝑜 1 − 𝑋
𝑋= (26) W= (27)
kC W 𝛼𝑊 𝛼
1 + Ao 1−
vo 2
Ejercicio 1 – Caída de presión solución analítica
Se tiene un reactor tubular de 130 m con un área de sección

transversal de 0.0045 m2, empacado con partículas de catalizador
cuyo diámetro es de 0.004 m. El flujo de alimentación del gas al
lecho empacado es de 0.038 kg/h. Las condiciones iniciales de
operación son 260°C y 10 atm. La fracción de vacío es de 35%, la
viscosidad es de 0.1002 kg/(m*h), la densidad es de 6.6156 kg/m3,
la densidad del catalizador solido es de 1765 kg/m3, k = 12
m6/kmol/kgcat/h. Calcular la caída de presión en el reactor.
Ejercicio 2 – solución analítica
Para el mismo reactor tubular con las características presentadas en

la tabla se tiene que el flujo volumétrico a la entrada es de 14 m3/h.
Determine:
1. La conversión en ausencia de caída de presión

2. La conversión considerando la caída de presión
𝐴𝑐 = 0.0045 𝑚2 𝐿 = 130 𝑚 ∅ = 0.35 𝐶𝐴𝑜 = 0.1 𝑘𝑚𝑜𝑙/𝑚3

𝐵𝑜 = 3.33 𝑘𝑃𝑎/𝑚 𝑃𝑜 = 10 𝑎𝑡𝑚 𝜌𝑏 = 1765 𝑘𝑔/𝑚3
Taller 2. Caída de presión
Un reactor tubular empacado opera con una viscosidad de 0.100153 kg/(m h),
una densidad de 6.61563 kg/m3 y un área transversal de 0.0045 m2. El reactor
tiene una longitud de 15 m, el diámetro de las partículas de catalizador es de ¼
de pulgada, G = 56 lbm/h/ft2, la fracción de vacío es de 0.48, Po = 15 atm, la
densidad del catalizador es de 1856 kg/m3. Graficar el cambio en el peso del
catalizador en función de la conversión.
vo 14 m3/h
Adicionalmente analice que sucede si: Cao 0.1 kmol/m3
k 12 m3/(kmol h kg cat)
1. El diámetro de partícula se recude a la tercera parte
2. El diámetro de partícula aumenta 1.5 veces
3. La presión inicial aumenta hasta 13 atm.

Caida de Presion en Reactores Part1

Cargado por

Copyright:

Formatos disponibles

Caida de Presion en Reactores Part1

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Caida de Presion en Reactores Part1

Cargado por

Copyright:

Formatos disponibles

Diseño de reactores.

Fase liquida Fase gaseosa

• No hay cambios en las • Las concentraciones de los

𝑐𝑎𝑚𝑏𝑖𝑜 𝑒𝑛 𝑒𝑙 𝑛ú𝑚𝑒𝑟𝑜 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑚𝑜𝑙𝑒𝑠 𝑝𝑎𝑟𝑎 𝑐𝑜𝑛𝑣𝑒𝑟𝑠𝑖ó𝑛 𝑡𝑜𝑡𝑎𝑙

Determinación de la relación P/Po en función de la concentración de

Por la ecuación de continuidad y la ecuación de gases ideales

Corrigiendo la densidad del gas en la Ecuación de Ergun

Reemplazando 15 y 16 en 13 y haciendo y = P/Po:

La ecuación 17 se aplica para caída de presión cuando ocurren

𝑑𝑦 𝛼 𝑇 La ecuación 19 para reacciones únicas en

𝑑𝑦 2 Sacando raíz a cada lado

𝐹𝐴 Integrando con X=0 cuando W=0 y

Despejando conversión Despejando peso del catalizador

𝑘𝐶𝐴𝑜 𝑊 𝛼𝑊 2𝑣𝑜 𝛼 𝑋 1/2

Se tiene un reactor tubular de 130 m con un área de sección

Para el mismo reactor tubular con las características presentadas en

1. La conversión en ausencia de caída de presión

𝐴𝑐 = 0.0045 𝑚2 𝐿 = 130 𝑚 ∅ = 0.35 𝐶𝐴𝑜 = 0.1 𝑘𝑚𝑜𝑙/𝑚3

También podría gustarte