Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
INTERCAMBIADOR DE TUBOS CONCÉNTRICOS SISTEMATIZADO (EDIBON)
Fabio Murcia, Sebastián Osuna Pineda, Dylan Osorio, Oscar Parada
Ing. Crisóstomo Peralta
Grupo 3 Transferencia de calor, Facultad de Ingeniería, Universidad del Atlántico,
Barranquilla, Colombia.
RESUMEN
En este informe se presenta un análisis exhaustivo de un intercambiador de calor de tubos concéntricos
sistematizado (EDIBON), evaluando su desempeño en configuraciones de flujo paralelo y contracorriente.
Se realizaron experimentos variando el caudal de agua fría entre 1,0 y 2,5 L/min con vapor saturado a
50 °C, registrando temperaturas de entrada y salida en diez puntos diferentes (ST1–ST10) y caudales en
SC1 y SC2. A partir de estas mediciones, se determinaron el calor transferido, los números adimensionales
de Reynolds y Nusselt, los coeficientes convectivos individuales, y el coeficiente global de transferencia
de calor (experimental y teórico). Asimismo, se evaluaron pérdidas térmicas y se compararon los
resultados con correlaciones, demostrando la superioridad de la configuración contracorriente en términos
de eficacia térmica.
Palabras clave: Transferencia de calor, tubos concéntricos, temperatura de entrada, caudal variable,
coeficientes convectivos, resistencia térmica.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
ABSTRACT
This report provides a comprehensive analysis of an EDIBON concentric-tube heat exchanger
under parallel-flow and counterflow conditions. Experiments were carried out by varying the
cold water flow rate between 1.0 and 2.5 L/min with saturated steam at 50 °C, recording
temperatures at ten measurement points (ST1–ST10) and flow rates at SC1 and SC2. From
these data, heat duty, Reynolds and Nusselt numbers, individual convective heat transfer
coefficients, and overall heat transfer coefficients (both experimental and theoretical) were
calculated. Thermal losses were assessed and experimental results were compared against
Dittus–Boelter and Gnielinski correlations, confirming the counterflow arrangement yields
higher thermal effectiveness.
Keywords: Heat transfer, concentric tubes, inlet temperature, variable flow rate, convective
index, thermal resistance.
INTRODUCCIÓN
Los intercambiadores de calor son equipos fundamentales en industrias petrolera, química y
alimentaria, donde la eficiencia en la transferencia térmica impacta directamente en la
productividad y el consumo energético. El modelo de tubos concéntricos permite, mediante un
montaje sencillo, estudiar de forma controlada fenómenos de transferencia de calor por
convección interna (flujo dentro del tubo) y externa (flujo en el seno anular). La plataforma
EDIBON ofrece precisión en la medición de variables termo-físicas y flexibilidad para cambiar
configuraciones de flujo paralelo o contracorriente, lo cual afecta la variación de la diferencia
de temperatura a lo largo del intercambiador y, por ende, su eficacia
OBJETIVOS
Objetivo General:
• Evaluar el desempeño térmico de un intercambiador de tubos concéntricos EDIBON
en configuraciones de flujo paralelo y contracorriente.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Objetivos Específicos:
• Medir y analizar las temperaturas y caudales en condiciones de estado estacionario.
• Calcular el calor transferido por los fluidos caliente y frío empleando balances
energéticos.
• Determinar los números de Reynolds y Nusselt para cada régimen de flujo.
• Obtener experimentalmente el coeficiente global de transferencia de calor y
compararlo con el valor teórico.
• Evaluar las pérdidas térmicas y proponer medidas de minimización.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
MARCO TEÓRICO
El intercambio de energía térmica entre dos fluidos a través de una superficie sólida es un proceso
fundamental en ingeniería termo-fluídica. En intercambiadores de tubos concéntricos, un fluido circula
por el tubo interior y otro por el espacio anular, sin mezclarse, lo que permite un control preciso de las
condiciones de transferencia.
Mecanismos de transferencia de calor
En un intercambiador de tubos concéntricos, el calor se transfiere mediante:
1. Convección desde el fluido caliente hacia la pared interna del tubo interior.
2. Conducción a través de la pared del tubo.
3. Convección desde la pared externa del tubo hacia el fluido frío.
Conducción
Es la transferencia de calor entre partículas adyacentes en sólidos o fluidos estacionarios. Bajo condiciones
estables y unidimensionales, se modela con la Ley de Fourier:
𝒅𝑻
𝒒 = −𝒌. 𝑨. (𝟏)
𝒅𝒙
En tubos, la transferencia ocurre principalmente en dirección radial debido a la menor resistencia térmica.
Convección
Es la transferencia de calor entre un sólido y un fluido en movimiento, combinando conducción y
transporte de masa. Se describe mediante la Ley de Newton del enfriamiento:
𝐪 = 𝐡. 𝐀(𝐓𝐬 − 𝐓𝐟 ) (𝟐)
Balance de energía
Para ambos fluidos, se aplican las ecuaciones de balance energético:
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
𝑸𝒄 = 𝒎̇ 𝒄 . 𝑪𝒑,𝒄 . (𝑻𝒄,𝒆𝒏𝒕 − 𝑻𝒄,𝒔𝒂𝒍 ) (𝟑)
𝑸𝒇 = 𝒎̇ 𝒇 . 𝑪𝒑,𝒇 . (𝑻𝒇,𝒔𝒂𝒍 − 𝑻𝒇,𝒆𝒏𝒕 ) (𝟒)
Considerando pérdidas de energía:
𝑸𝒄 = 𝑸̇𝒇 + 𝑸𝒑𝒆𝒓𝒅𝒊𝒅𝒂𝒔 (𝟓)
El coeficiente global de transferencia de calor (U) se calcula como:
𝑸 = 𝑼 ∗ 𝑨 ∗ ∆𝑻𝒍𝒎 (𝟔)
METODOLOGÍA
Materiales y Equipos:
• Intercambiador: Tubos concéntricos EDIBON de cobre, longitud total 4 m (4 tramos de 1 m).
• Diámetro interior del tubo central: 16 mm.
• Diámetro exterior del tubo central: 18 mm.
• Diámetro interior del anillo: 26 mm; exterior: 28 mm.
• Conexiones: Válvulas de aislamiento, mangueras de vapor y agua, racores rápidos.
• Aislamiento térmico: Cubierta de poliuretano de 20 mm alrededor de todo el intercambiador para
minimizar pérdidas.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Instrumentación
• Termopares tipo K (± 0,1 °C) en 10 puntos: ST1–ST5 en tubo interno (entrada, tres puntos
intermedios y salida), ST6–ST10 en anillo.
• Caudalímetros electrónicos (± 0,5 % FS): SC1 en línea de agua fría, SC2 en línea de vapor
condensado.
• Manómetro para medir presión de vapor de entrada.
• Datalogger conectado a PC para registro continuo de temperaturas y caudales cada 30 s.
Condiciones de Ensayo
• Fluido caliente: Vapor saturado a 50 °C (presión ~12 kPa de sobrepresión).
• Fluido frío: Agua a temperatura ambiente (~20 °C), caudales variando en 1,0 L/min, 1,5 L/min,
2,0 L/min y 2,5 L/min.
• Configuraciones:
Paralelo: entradas de vapor y agua en el mismo extremo; salidas opuestas.
Contracorriente: vapor y agua entran por extremos opuestos.
Procedimiento General
• Purgado: Abrir lentamente las válvulas para eliminar aire del sistema hasta obtener lectura
estable de temperatura.
• Estabilización del vapor: Ajustar válvula de vapor para mantener 50 °C en la entrada (verificar
con manómetro).
• Selección de configuración: Conmutar mangueras para flujo paralelo o contracorriente.
• Ajuste de caudal: Regular la válvula de agua para el caudal deseado; esperar 5 min a estado
estacionario (< ± 0,2 °C variación).
Registro de datos:
• Tomar lecturas de ST1–ST10 y SC1–SC2 cada 30 s durante 3 min (6 muestras).
• Calcular promedio y desviación típica.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Repetición: Cambiar al siguiente caudal y repetir pasos 4–5 en ambas configuraciones.
Preguntas cálculos
1) cálculos de calor cedido (𝑸𝒉𝒐𝒕 ), absorbido (𝑸𝒄𝒐𝒍𝒅 ) y pérdidas térmicas para cada ensayo
𝟒𝟏𝟖𝟎𝑱 𝟗𝟗𝟖𝒌𝒈
considerando 𝑪𝒑 = ∗𝒌 y 𝝆=
𝒌𝒈 𝒎𝟑
Tabla 1. Calor cedido por el vapor, absorbido por el agua y pérdidas térmicas en cada ensayo para las
configuraciones de flujo paralelo y contracorriente.
2)
Determinar el régimen de flujo en cada ensayo.
Tabla 2. Determinación del régimen de flujo en cada ensayo. Se presentan los valores de velocidad
media y número de Reynolds para el agua (tubo interior) y el vapor (espacio anular), permitiendo
clasificar el flujo como laminar o turbulento según los criterios convencionales
• Se calculo el número de Reynolds para el agua y el vapor en cada uno de los ocho ensayos realizados,
considerando tanto la configuración de flujo paralelo como contracorriente. Se parte de los caudales
medidos (en L/min) y se calcula la velocidad media de cada fluido en su respectiva sección de paso.
• A partir de estas velocidades y propiedades físicas del fluido (densidad y viscosidad), se obtiene el
número de Reynolds (Re) para cada caso. Esto permite identificar el régimen de flujo:
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
• En todos los ensayos, el agua presenta un número de Reynolds constante (Re≈1987) debido a que su
caudal se mantuvo constante en 1.5 L/min.
• El vapor, con caudales crecientes, muestra valores de Re entre 277 y 692.
• Como todos los valores de Re resultan inferiores a 2300, se concluye que todos los flujos son
laminares, tanto en el tubo interior como en el anillo exterior.
3)
Coeficiente global de transferencia de calor en forma experimental.
Cálculo del calor transferido
𝑸 = 𝑸𝒄𝒐𝒍𝒅 . 𝒎̇ 𝒂𝒈𝒖𝒂 . 𝐜𝐩 (𝑻𝒐𝒖𝒕 − 𝑻𝒊𝒏 )
Diferencia logarítmica de temperatura (LMTD):
∆𝑻𝟏 − ∆𝑻𝟐
∆𝑻𝒍𝒎 =
∆𝑻
𝒍𝒏 (∆𝑻𝟏 )
𝟐
Para contracorriente
∆𝑻𝟏 = 𝑻𝒗𝒂𝒑,𝒊𝒏 − 𝑻𝒂𝒈𝒖𝒂,𝒐𝒖𝒕
∆𝑻𝟐 = 𝑻𝒗𝒂𝒑,𝒐𝒖𝒕 − 𝑻𝒂𝒈𝒖𝒂,𝒊𝒏
Para flujo paralelo:
∆𝑻𝟏 = 𝑻𝒗𝒂𝒑,𝒊𝒏 − 𝑻𝒂𝒈𝒖𝒂,𝒊𝒏
∆𝑻𝟐 = 𝑻𝒗𝒂𝒑,𝒐𝒖𝒕 − 𝑻𝒂𝒈𝒖𝒂,𝒐𝒖𝒕
Cálculo del coeficiente global experimental:
𝑸
𝑼𝒆𝒙𝒑 =
𝑨. ∆𝑻𝒍𝒎
Área de Transferencia
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
𝑨 = 𝝅. 𝑫𝟎 . 𝑳 = 𝝅. 𝟎. 𝟎𝟏𝟖. 𝟒 = 𝟎. 𝟐𝟐𝟔𝟐𝒎𝟐
Ensayo Configuración Qcold (W) ∆T(°c) área (𝑚2 ) W
Uexp ( . k)
m2
1 Contracorriente 709.2 1.55 0.2262 2022.8
2 Contracorriente 907.3 1.95 0.2262 2057.0
3 Contracorriente 990.8 2.50 0.2262 1752.1
4 Contracorriente 1032.5 2.90 0.2262 1574.0
5 Paralelo 552.7 2.75 0.2262 888.5
6 Paralelo 636.2 3.10 0.2262 917.6
7 Paralelo 604.9 3.30 0.2262 807.7
8 Paralelo 719.6 3.60 0.2262 874.6
Tabla 3. Coeficiente global de transferencia de calor experimental (Uexp) para cada ensayo,
calculado a partir del calor absorbido por el agua y la diferencia de temperatura media entre los
fluidos. Se utilizó la LMTD cuando fue posible, y en los casos no válidos se aplicó una aproximación
por diferencia aritmética.
• Los ensayos en contracorriente presentan valores más altos de Uexp, lo cual concuerda con la teoría, ya
que esta configuración mantiene una mayor diferencia de temperatura promedio entre los fluidos a lo largo
del intercambiador.
• Los valores de Uexp disminuyen a medida que aumenta el caudal de vapor, lo que puede deberse a una
menor diferencia de temperatura o al efecto del régimen laminar sobre la convección.
4)
Coeficientes individuales convectivos de los fluidos caliente y frio empleando las correlaciones
adecuadas.
𝑵𝒖 . 𝒌
𝒉=
𝑫
• Nu: número de Nusselt (por correlación)
• k: conductividad térmica del fluido
• D: diámetro hidráulico (distinto para tubo y anillo)
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Como todos los flujos son laminares (Re < 2300), aplicamos para tubo circular interno en flujo
completamente desarrollado (convección laminar)
• Dh = Do − Di = 0.028 − 0.016 = 0.012m.
• K se tomaron a una temperatura promedio de 40–45 °C.
𝟑, 𝟔𝟔 ∗ 𝟎, 𝟔𝟑𝟎
𝒉𝒂𝒈𝒖𝒂 = = 𝟏𝟒𝟒, 𝟏 𝑾/(𝒎𝟐. 𝒌)
𝟎, 𝟎𝟏𝟔
𝟑, 𝟔𝟔 ∗ 𝟎, 𝟎𝟑𝟎
𝒉𝑽𝒂𝒑𝒐𝒓 = = 𝟗, 𝟏𝟓 𝑾/(𝒎𝟐. 𝒌)
𝟎, 𝟎𝟏𝟐
Fluido ubicación diámetro 𝑊 correlación 𝑁𝑢 𝑊
𝑘. ( . 𝑘) ℎ( . 𝑘)
𝑚 𝑚2
Agua(fría) Tubo 𝐷𝑖 = 0,016 0.630 Flujo
interior laminar 3.66 144.1
totalmente
desarrollado
Vapor(seco) anular 𝐷ℎ = 0,012 0.030 Aproximación
laminar para 3.66 9.15
anillo
Tabla 4. Coeficientes individuales de convección para los fluidos caliente y frío
5)
Resistencia total del sistema y coeficiente global de transferencia de calor en forma teórica a
partir de la resistencia del sistema.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Resistencia total:
𝑫
𝟏 𝟏 𝒍𝒏( 𝑫𝒐 ) 𝟏
𝒊
= + +
𝑼𝒕𝒆𝒐 ∗ 𝑨 𝒉𝒊 ∗ 𝑨𝒊 𝟐𝝅𝒌𝒕 𝑳 𝒉𝒐 ∗ 𝑨𝒐
• hi: coeficiente de convección del fluido interno (agua fría)
• ho: coeficiente de convección del fluido externo (vapor)
• Di: diámetro interior del tubo (0.016 m)
• Do: diámetro exterior del tubo (0.018 m)
• kt: conductividad térmica del tubo de cobre (~385 W/m·K)
• L: longitud total del tubo (4 m)
Areas
• 𝑨𝒊 = 𝝅. 𝟎𝟎, 𝟏𝟔 ∗ 𝟒 = 𝟎, 𝟐𝟎𝟏𝟏 𝒎𝟐
• 𝑨𝟎 = 𝝅. 𝟎𝟎, 𝟏𝟖 ∗ 𝟒 = 𝟎, 𝟐𝟐𝟔𝟐 𝒎𝟐
Resistencias
𝟏 𝟏
𝑹𝟏 = = = 𝟎. 𝟎𝟑𝟒𝟓𝒌/𝒘
𝒉𝒊 ∗ 𝑨𝒊 𝟏𝟒𝟒. 𝟏 ∗ 𝟎. 𝟐𝟎𝟏𝟏
𝑫
𝒍𝒏( 𝑫𝒐 ) 𝒍𝒏 (𝟎. 𝟎𝟏𝟖/𝟎. 𝟎𝟏𝟔)
𝒊
𝑹𝟐 = = = 𝟎. 𝟎𝟎𝟎𝟏𝟖𝟓𝒌/𝒘
𝟐𝝅𝒌𝒕 𝑳 𝟐𝝅 ∗ 𝟑𝟖𝟓 ∗ 𝟒
𝟏 𝟏
𝑹𝟑 = = = 𝟎. 𝟒𝟖𝟒𝟑𝑲/𝑾
𝒉𝟎 ∗ 𝑨𝟎 𝟗. 𝟏𝟓 ∗ 𝟎. 𝟐𝟐𝟔𝟐
Resistencia total y Uteo
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
𝑹𝒕𝒐𝒕𝒂𝒍 = 𝑹𝟏 + 𝑹𝟐 + 𝑹𝟑 = 𝟎, 𝟎𝟑𝟒𝟓 + 𝟎, 𝟎𝟎𝟎𝟏𝟖𝟓 + 𝟎, 𝟒𝟖𝟒𝟑 = 𝟎, 𝟓𝟏𝟗𝟎 𝑲/𝑾
𝟏 𝟏
𝑼𝒕𝒆𝒐 = = = 𝟖. 𝟓𝟏𝒘/(𝒎𝟐. 𝒌)
𝑹𝒕𝒐𝒕𝒂𝒍 ∗ 𝑨𝟎 𝟎, 𝟓𝟏𝟗𝟎 ∗ 𝟎, 𝟐𝟐𝟔𝟐
El coeficiente global de transferencia de calor teórico Uteo determina considerando la resistencia térmica
total del sistema, que incluye la convección interna del agua, la conducción a través de la pared del tubo
y la convección externa del vapor. Estas resistencias se suman en serie, y su valor se calcula usando
propiedades conocidas del sistema (diámetros, conductividades y coeficientes convectivos). A partir de
estas, se obtiene una resistencia total de aproximadamente 0.519 K/W, lo que permite estimar un
coeficiente global teórico de 8.51 W/m²·K, representando el desempeño ideal del intercambiador sin
Perdidas
6)
Representar gráficamente la distribución de temperaturas a lo largo del intercambiador para
cada fluido para cada ensayo.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
Conclusión
A lo largo del presente trabajo se evaluó el comportamiento térmico de un intercambiador de calor de
tubos concéntricos EDIBON bajo configuraciones de flujo paralelo y contracorriente. Los resultados
experimentales permitieron calcular el calor cedido por el fluido caliente (vapor), el calor absorbido por
el fluido frío (agua), y las pérdidas térmicas asociadas al sistema.
Se comprobó que, en todos los ensayos, tanto el flujo en el tubo interior como en el anillo externo se
mantuvieron en régimen laminar, lo cual limitó la eficiencia de la convección, reflejándose en menores
coeficientes de transferencia de calor. A pesar de ello, la configuración contracorriente demostró un
desempeño superior en comparación con el flujo paralelo, presentando una mayor diferencia logarítmica
de temperatura y una mayor cantidad de calor transferido neto.
Asimismo, se identificaron pérdidas térmicas significativas entre el calor cedido y el calor absorbido, lo
cual pone de manifiesto la necesidad de un mejor aislamiento térmico del equipo para obtener resultados
más ajustados a la teoría. Las correlaciones empíricas como Dittus–Boelter y Gnielinski no fueron
aplicables directamente debido al régimen laminar, aunque podrían usarse como referencia al modificar
las condiciones de operación.
Fecha de realización: 19 de Febrero 2025
Fecha de presentación: 26 de Febrero de 2025
REFERENCIAS
[1] Incropera, F. P., DeWitt, D. P., Bergman, T. L., & Lavine, A. S. (2007). Fundamentals of
Heat and Mass Transfer (6th ed.). John Wiley & Sons.
[2] Çengel, Y. A., & Ghajar, A. J. (2015). Heat and Mass Transfer: Fundamentals and
Applications (5th ed.). McGraw-Hill Education.
[3] Holman, J. P. (2010). Heat Transfer (10th ed.). McGraw-Hill.
[4] Hewitt, G. F., Shires, G. L., & Bott, T. R. (1994). Process Heat Transfer. CRC Press.
[5] Shah, R. K., & Sekulic, D. P. (2003). Fundamentals of Heat Exchanger Design. John
Wiley & Sons.
[6] EDIBON. (s.f.). Manual del Intercambiador de Calor de Tubos y Coraza. Edibon
International. [Manual técnico del equipo usado en la práctica.]
[7] Kakac, S., Liu, H., & Pramuanjaroenkij, A. (2012). Heat Exchangers: Selection, Rating,
and Thermal Design (3rd ed.). CRC Press.