Diseño de Hornos de Proceso

Método de Lobo y Evans

 Introducción - Justificación

• Estimar dimensiones

• Estimar consumo de combustible

• Evaluar y comparar propuestas de vendedores

• Predecir efecto del cambio de variables

• Predecir el efecto de modificaciones

Tipos básicos de hornos de proceso

Hasta 40 MMkcal/h
Cilíndrico
vertical
Cilíndrico vertical
Tipos cabina y caja
Cabina
Caja doble cámara parte radiante
Partes Chimenea

Arco (Arch)

Material
refractario

Entrada Breeching
Corbel

Sección convectiva

Tubos escudo Sección

Placa radiante
portatubos
Cruce Tabique Envolvente
(casing)
Cabezal

Tubos Plataforma

Salida

Quemador
Pilar
Quemadores
Quemadores
Quemadores
 Material Refractario
• Material aislante dentro del horno
• Propósito:
– reducir pérdidas de calor al ambiente
– Reducir temperatura de la envolvente (ej: 82°C
con 27°C amb)
– Material: ladrillos, colado o fibras cerámicas
 Material Refractario

Colado

Ladrillos

Fibras cerámicas
 Material de Tubos
Límite a las temperaturas de diseño
Material Tipo o Grado Máxima Temperatura
Carbon steel B 538°C
C-1/2Mo T1 or P1 593°C
1 1/4Cr-1/2Mo T11 or P11 593°C
2 1/4Cr-1Mo T22 or P22 649°C
5Cr-1/2Mo T5 or P5 649°C
7Cr-1/2Mo T7 or P7 704°C
9Cr-1Mo T9 or P9 704°C
18Cr-8Ni 304 or 304H 815°C
16Cr-12Ni-2Mo 316 or 316H 815°C
18Cr-10Ni-Ti 321 or 321H 815°C
18Cr-10Ni-Cb 347 or 347H 815°C
Ni-Fe-Cr Alloy 800H 982°C
25Cr-20Ni HK-40 1010°C
 Cuerpo Gris
• Poder emisivo de un cuerpo es la energía que
este emite por unidad de tiempo y área
• Emisividad (𝜖) de un cuerpo: relación de su poder
emisivo al de un cuerpo negro a igual T
• Poder emisivo cuerpo negro (absorbe toda la
radiación incidente) 𝑊𝑏 = 𝜎𝑇 4
• Para un cuerpo gris, la emisividad es constante
para cualquier
– longitud de onda
– temperatura
Factor de visión entre cuerpos negros
Flujo de calor de dA1 a dA2

Integrando considerando que se

absorbe la energía en su totalidad

Energía total que sale de A1

Definición del factor de visión

 Modelo Zona Radiante
• Principal dificultad → Factores de visión
• Simplificación → Plano gris equivalente
 Factor de visión - Eficiencia
• Entre plano y una o dos
filas de tubos
• Para más de dos filas,
αpe=1
• Con o sin refractario

αpe : factor de comparación

con dos planos paralelos
Pitch
Do
 Factor de visión - Eficiencia
Banco de tubos con refractario a T2
con emisor a TE
αpe
Superficies paralelas F=1

Superficie 1 a T2 con eficiencia

αpe es entonces el factor de Pitch

eficiencia del plano equivalente Do
 Gases activos e inertes
• Ciertos gases pueden absorber o irradiar
energía.
• Gases de moléculas iguales N2, O2, H2 son
transparentes
• Gases con moléculas asimétricas CO2, CO, H2O
son capaces de absorber y emitir
• Gráficos de McAdams permiten calcular
emisividad y absorción
 Circuito Equivalente del Horno
Zona Radiante
Factor para
Resistencia a la transferencia de
superficie gris
energía entre el gas y los tubos
WG: Poder emisivo del T2: Temperatura
gas superficie tubos

Resistencia a la transferencia de
energía entre el gas y el refractario Emitido directamente
entre refractario y los
tubos

: Transparencia
Resistencia equivalente Resistencia equivalente
Total zona punteada
 Circuito Equivalente del Horno
Zona Radiante
Cálculo
1 1
=
1 1 1
+ + 𝑅𝑓
𝑈 ℎ𝑖𝑜 ℎ𝑜
Despejando Comparación con cálculo
coeficiente global de
transferencia

Resistencia eq. total ¿Cálculo? → Lobo y Evans

 Radiación a tubos escudo y efecto de
convección
Tubos escudo Efecto convectivos
• No tienen refractario detrás • Aunque casi la totalidad del
de ellos calor se transfiere por
• La energía que no impacta radiación, existen efectos
directamente, es absorbida convectivos
por los tubos de la zona
convectiva
• αpe = 1 para estos tubos

Calor transferido considerando efectos convectivos

 Flujo de calor – Sección radiante
• Primer paso: selección de flujo promedio de
calor radiante (Q/A2)
– Q: Calor recibido por el fluido
– A2: Área externa de los tubos
• Relación de compromiso, altos Q/A2:
– ↓ tamaño del horno
– ↑ T de gases de combustión ↑ mantenimiento
– ↓ vida útil de los materiales
– ↑ probabilidad de formación de coque
 Flujo de calor – Sección radiante
Valores típicos de densidad promedio de flujo de calor en la zona radiante

Densidad radiante promedio

Tipo de horno
(W/m2)
Calentadores de crudo para destilación atmosférica 31400–44000
Calentadores de crudo para destilación al vacío 25000–31400
Reboilers 31400–37600
Hornos de hot oil 25000–34500
Hornos de refroming catalítico 23500–37600
Hornos de coque retardado 31400–34599
Desasfaltado de propano 25000–28000
Calentadores de carga a procesos de hidrotratamiento 31400
Calentadores de carga a cracking catalítico 31400–34500
Sobrecalentadores de vapor 28200–40000
Hornos de plantas de gasolina 31400–37600
 Flujo de calor máximo

El flujo de calor no es uniforme puede encontrarse

entre 1,5 y 3 veces la media
 Método de Lobo y Evans
Limitaciones
• Considera temperatura del gas de combustión
igual a la de salida de la zona radiante
– Buenos resultados en hornos tipo caja o cabina
– No en cilíndricos o con recorridos largos → Wimpress
sugiere tomar una diferencia entre esas temperaturas
• Asume T uniforme de pared de tubo
• Publicación original daba el estimado de la
longitud semiesférica equivalente para solo 6
casos de relación de dimensiones del hogar
Longitud semiesférica equivalente
• Longitud ficticia característica de la geometría
bajo consideración que permite obtener las
emisividades y factores de absorción del gas
para que la transferencia pueda expresarse
como:
Relación dimensional (Longitud-ancho-alto
L
en cualquier orden)
Hornos rectangulares
1-1-1 to 1-1-3 2/3(volumen del horno)1/3
1-1-4 to 1-1-∞ 1 × dimensión menor
1-2-5 to 1-2-8 1.3 × dimensión menor
1-3-3 to 1-∞-∞ 1.8 × dimensión menor
Hornos cilíndricos
d×d 2/3 × diámetro
d × 2d a d × ∞ 1 × diámetro
 Cálculo factor de visión
• Propuesta de Lobo y Evans

AR: Área de refractario

A1: Área del plano equivalente
 Emisividad de los gases
• Gráficos de McAdams ó
• Gráfico simple cuando
se usan combustibles
de refinería (H2O y CO2
cte.)
• P parcial H2O y CO2
– Gráfico simple para
Comb. de refinería
– Estequiometría
Gráfico para combustibles de refinería
 Emisividad de los gases

Gráfico simple para

utilizar en lugar del
de McAdams

Presiones cercanas a la atmosférica

Emisividad equivalente
de gas gris
 Temperatura de pared de tubo
• La ecuación:

permite obtener Q a una dada T2

• Acople transf. pared tubo - fluido incluyendo:
– Resistencia de la pared del tubo
– Resistencia de ensuciamiento o coque
– Resistencia del fluido (capa límite)

o
A2: Área externa de los tubos
 Procedimiento de cálculo - Diseño
• Adoptar un Q/A2 deseado ⇒ 𝐴2 [Q = 𝑊𝑐 (𝑖𝑒 − 𝑖𝑠 )]
• Definir configuración (𝑁𝑡 , 𝐿, 𝐷0 , 𝑃𝑡 )
• Definir número de pasos → ramas en paralelo (vel. y ΔP)
• Calcular hio (suponer Ø=1ra aprox.)
• Calcular T2 (para el Q/A2 deseado)
• Calcular TG
– iterativo ƐG depende de TG
– asumir hG=11,5 W/m2K (bajo impacto)
– Suponer ℱ1𝐺 y calcular TG de la sig. ecuación

– con pCO2 + pH2O y long. semiesférica equiv. → εg

– Calcular 𝐹1𝐺 y verificar ℱ1𝐺
• Verificar ΔP
 Conclusiones
• El método de Lobo y Evans permite el cálculo
sencillo de hornos de procesos con un ajuste
adecuado a los resultados de campo
• Las simplificaciones en el cálculo de ƐG cuando
se utilizan combustibles comunes en refinería
aceleran aún más los cálculos facilitando las
iteraciones
 Bibliografía Recomendada
• Cao, E. (2006). Transferencia de calor en ingeniería de
procesos. Buenos Aires: Nueva librería.
• Cao, E. (2010). Heat transfer in process engineering.
New York: McGraw-Hill. Retrieved from
http://accessengineeringlibrary.com/browse/heat-
transfer-in-process-engineering
• Kern, D. Q. (1965). Procesos de transferencia de calor.
México, D.F.: Compañía Editorial Continental.
• Al-Haj Ibrahim, H., & Mourhaf Al-Qassimi, M. (2013).
Calculation of Radiant Section Temperatures in Fired Matlab
Process Heaters. Chemical Engineering and Science,
1(4), 55–61. https://doi.org/10.12691/ces-1-4-2
• Ibrahim, H. A.-H., & Al-Qassimi, M. (2010). Simulation of
heat transfer in the convection section of fired process Matlab
heaters. Periodica Polytechnica Chemical Engineering,
54(1), 33. https://doi.org/10.3311/pp.ch.2010-1.05
 Otras lecturas
• Baukal, C. E., & John Zink Hamworthy Combustion. (2013). The John
Zink Hamworthy combustion handbook.
• Mullinger, P., & Jenkins, B. (2008). An introduction to heat transfer
in furnaces. In Industrial and Process Furnaces (pp. 89–140).
Elsevier. Retrieved from http://linkinghub.elsevier.com/retrieve/pii
/B9780750686921000041
• API Standard 560 - Fired Heaters for General Refinery Service.
• API Standard 530 - Calculation of Heater Tube thickness in
Petroleum Refineries
• API RP 536- Post Combustion NOx Control in Fired Equipment in
General Refinery Services
• API RP 535- Burners for Fired Heaters in General Refinery Service
 Problema de ejemplo
• Estimar el tamaño de un horno tipo cabina para calentar 45
kg/s de un crudo de 26°API desde 250°C a 375°C y consumo
de gas combustible.
• Crudo Cp=2978 J/kg K
• Gas combustible PCI: 60 MJ/Nm3, Densidad: 0,88 kg/m3
• 25% exceso de aire
• Tubos DN: 4”(Do:0,1143m Di:0,103m) LT: 11,5m Pt:8”
• hi: 757 W/m2K Rfio: 0,0003 m2K/W
• Flujo medio de calor deseado: 35500 W/m2
 Problema de ejemplo
• Calor intercambiado
𝑄 = 𝑊𝑐 𝐶𝑐 𝑡2 − 𝑡1 = 45 2978 357 − 250 = 14,34 𝑀𝑊

• Estimación del área de transferencia

𝑄 14,34𝑀𝑊
𝐴2 = = = 403,9𝑚2
𝑞 35500 𝑊
𝑚2
• Número de tubos
𝐴2 403,9
𝑁= = = 95,7 ∴ 𝑁 = 96
𝜋𝐷𝑜 𝐿 𝑇 3,14 0,1143 11,75

• Por lo que el área es A2=404,8m2

 Problema de ejemplo
• Área de la zona radiante Trapecio
Tubos escudo
𝐴𝑇
= 2 6,96 ⋅ 5,08 + 0,5 6,96 + 1,22 2,87
+ 2 4,06 + 5,08)12,2 + 6,96 ⋅ 12,2 = 402,1
Piso
• Área del plano equivalente
(tubos radiante + escudo)
𝐴1𝑅 = 𝐿 𝑇 𝑁 𝑅 𝑃𝑡 = 11,75 ⋅ 90 ⋅ 0,203 = 214,8𝑚2
𝐴1𝐸 = 𝐿 𝑇 𝑁 𝐸 𝑃𝑡 = 6 ⋅ 11,75 ⋅ 0,203 = 14,3𝑚2

• Eficiencia de plano equival.

𝑅
𝑃𝑡 𝐸 =1
𝛼𝑝𝑒 = 0,92 = 1,78 ; 𝛼𝑝𝑒
𝐷𝑜
Vale 1 porque la radiación es absorbida por los tubos de
choque o por los convectivos
 Problema de ejemplo
• Entonces:
𝐴1 𝛼𝑝𝑒 = 214,8 0,92 + 14,3 = 211,9𝑚2
𝐴𝑅 = 𝐴 𝑇 − 𝐴1 𝛼𝑝𝑒 = 402,1 − 211,9 = 190,2𝑚2
𝐴𝑅 190,2
= = 0,897 ∴ 𝐹𝑅1 =211,9/402,1=0,526
𝐴1 𝛼𝑝𝑒 211,9
• Longitud semiesférica equivalente
6,96 7,95 12,2 2 1/3
Relaciones 6,96 : 6,96 : 6,96 = 1: 1,15: 1,76 ∴ 𝐿 = 3 𝑣𝑜𝑙 = 5,54m
𝑣𝑜𝑙 = 6,96 5,08 + 0,5 6,96 + 1,22 2,87 12,2 = 574,5𝑚3
• Temperatura de pared del tubo
𝑡𝑓𝑙𝑢𝑖𝑑𝑜 = 0,5 250 + 357 = 303,5°𝐶
𝑇2 − 𝑡
𝑄= 𝐴2 1
−1
0,114
−1
𝑊
1
+ 𝑅𝑓 + 𝑅𝑓 = + 0,0003 = 567 2
ℎ𝑖𝑜 ℎ𝑖𝑜 0,103 757 𝑚 𝐾

𝑇2 = 𝑡 + 𝑄/(𝐴2 567) = 303,5 + 14339000/(404,8 ⋅ 567) = 365°𝐶

 Problema de ejemplo
• Iterando con la temperatura de los gases (TG=880°C)

𝑄 − 𝐴2 ℎ𝐺 𝑇𝐺 − 𝑇2 14339000 − 404,8 11,5 (880 − 365)

ℱ1𝐺 = 4 4 = 8 4 4
= 0,62
𝐴1 𝛼𝑝𝑒 𝜎 𝑇𝐺 − 𝑇2 211,9 5,67 10 [ 880 + 273 − 365 + 273 ]

𝑝𝐶𝑂2 + 𝑝𝐻2𝑂 = 0,23𝑎𝑡𝑚

• Verificación de ℱ1𝐺 .
𝑝𝐶𝑂2 + 𝑝𝐻2𝑂 𝐿 = 0,23 5,54 = 1,27𝑎𝑡𝑚. 𝑚

𝜖𝐺 = 0,53
𝐴𝑅 1
𝐹1𝐺 = 𝜖𝑔 1 + 𝜖𝑔
𝐴1 𝛼𝑝𝑒 1 +
𝐹𝑅1 1 − 𝜖𝑔
190,2 1
= 0,53 1 + = 0,681
211,9 1 + 0,53
0,526 1 − 0,53
 Problema de ejemplo
Superficie gris
• Verificando
1 1 1 − 𝜖2 1 1 − 0,9
= + = + = 1,58 ∴ ℱ1𝐺 = 0,633 ≈ 0,62
ℱ1𝐺 𝐹1𝐺 𝜖2 0,681 0,9
• Cálculo de consumo de combustible (NHV=PCI=57,2MJ/kg)

• Reescribiendo para usar gráfico simplificado

Pérdidas al amb
𝑄
• = 𝑁𝐻𝑉 1 − 0,02 − 0,45 =30,3 MJ/kg
𝑊𝐹
𝑄 14339000 𝑘𝑔 𝑘𝑔
𝑊𝐹 = = = 0,473 = 1704
𝑄/𝑊𝐹 57,2 106 0,53 𝑠 ℎ
Relación caudal de combustible / calor entregado al proceso
Camino seguido durante el diseño