INTERCAMBIADORES DE CALOR: CONSIDERACIONES DE DISEÑO

Tipos de Intercambiadores de Calor

Los intercambiadores de calor permiten el intercambio de calor entre dos fluidos que están
a diferentes temperaturas y separados por una superficie sólida Se clasifican de acuerdo al
arreglo del flujo y el tipo de construcción.

Intercambiadores de calor de tubo concéntrico.

.

Flujo paralelo
Contraflujo
Figura 1. Intercambiadores de calor de tubos concéntricos

Esta es la configuración mas simple y el intercambiador de contraflujo presenta un

desempeño mejor.

Intercambiadores de Calor de Flujo Cruzado.

Figura 2. Intercambiadores de calor de flujo cruzado

Los de flujo cruzado difieren según el fluido que se mueve sobre los tubos esté mezclado o
no mezclado. Para el primer caso se dice que el fluido no está mezclado porque las aletas
impiden el movimiento en una dirección (y) que es transversal a la dirección del flujo
principal (x). En el intercambiador de calor con aletas se dice que ambos fluidos están sin
mezclar mientras que en el intercambiador sin aletas un fluido está mezclado y el otro sin
mezclar.

1
 Intercambiadores de Calor de Coraza y tubos.

Figura 3. Intercambiador de calor de coraza y tubos.

Los bafles son usados para establecer un flujo cruzado e inducir la turbulencia, para
aumentar el coeficiente de convección del fluido del lado de la coraza. El número de tubos
y de pasos en la coraza puede ser variado.

(a)

2
 (b)
Figura 4. Intercambiadores de calor de coraza y tubos. (a) Un paso por la coraza y un
paso por los tubos. (b). Dos pasos por la coraza y cuatro pasos por los tubos.

Intercambiadores de Calor Compactos

Figura 5. Cubiertas de intercambiadores de calor de tipo compacto

a) Tubo con aletas (tubos planos, aletas de placas continuas

b) Tubo con aletas (tubos circulares, aletas de placa continuas)
c) Tubos con aletas (tubos circulares, aletas circulares)
d) Aletas de placa (un solo paso)
e) Aletas de placa (multipaso)

3
Figura 11.6a. Condensador Figura 11.6b Intercambiador de calor

Figura 11.7. Calentador de gas

4
 Figura 11.8 Intercambiadores de calor de coraza y tubo

Coeficiente Global de Transferencia de Calor

1 1 1 1 R "f ,c R "f ,h 1
= = = + + Rw + + (11.1)
UA U c Ac U h Ah (η o hA)c (η o A)c (η o A)h (η o hA)h

Las letras c y h se refieren al fluido frío y caliente respectivamente y R "f es el factor de

impureza debido a las incrustaciones en los tubos cuyos valores representativos se dan en la
siguiente tabla:

5
Valores representativos del factor de impureza
Fluido R "f (m2.K/W)
Agua de mar y agua tratada para alimentación de una 0.0001
caldera (por debajo de 50 C)
Agua de mar y agua tratada para alimentación de una 0.0002
caldera (por arriba de 50 C)
Agua de rìo (por debajo de 50 C) 0.0002 – 0.001
Aceite de motor 0.0009
Lìquidos refrigerantes 0.0002
Vapor (no aceitoso) 0.0001

La cantidad ηo es la eficiencia global de la superficie con aletas:

Af
ηo = 1 − (1 − η f ) (11.2)
At

La transferencia de calor es:

q = η o hA(Tb − T∞ ) (11.3)

Af es el área superficial de la aleta, At es el área total de la superficie con aletas y ηf es la

eficiencia de una aleta.

Para intercambiadores de calor tubulares, sin aletas, la ecuación para obtener U es:

1 1 1 1 R "f ,i ln( Do / Di) R "f ,o 1

= = = + + + + (11.4)
UA U i Ai U o Ao hi Ai Ai 2πkl Ao ho Ao

Donde los subíndices i y o se refieren a las superficies interna y externa del tubo
( Ai = πDi L , Ao = πDo L ), que se pueden exponer al fluido caliente o al frío.
Algunos valores representativos del coeficiente de transferencia de calor global U se
muestran en la siguiente tabla:

Valores representativos del coeficiente de transferencia de calor global U

Combinación de fluidos U (W/m2 K)
Agua con agua 850 - 1700
Agua con aceite 110 - 350
Condensador de vapor (agua en tubos) 1000 - 6000
Condensador de amoniaco (agua en tubos) 800 - 1400
Condensador de alcohol (agua en tubos) 250 - 700
Intercambiador de calor de tubos con aletas (agua en
tubos, aire en flujo cruzado) 25 50

6
 Análisis del Intercambiador de Calor: uso de la diferencia de temperatura media
logarítmica

Balance de energía global

Figura 11.9. Balance de energía global

Asumiendo que la transferencia de calor hacia los alrededores es despreciable y que los
cambios de energía cinética y potencial para cada fluido es despreciable.

q = m h (i h , i − i h , o ) (11.5)

q = m c (i c , o − i c , i ) (11.6)

i es la entalpía del fluido

Asumiendo que no hay cambio de fase y que los calores específicos de los fluidos son
constantes

q = m h c p , h (Th ,i − Th ,o ) = C h (Th ,i − Th ,o ) (11.7)

q = mc c p ,c (Tc ,o − Tc ,i ) = C c ((Tc ,o − Tc ,i ) (11.8)
Ch y Cc son las velocidades de capacidad de calor

Otra expresión útil se puede obtener relacionando la velocidad de transferencia de calor

total con la diferencia de temperatura ∆T entre los fluidos caliente y frío, donde:

∆T = Th − Tc (11.9)

Tal expresión debería ser una extensión de la ley de enfriamiento de Newton, usando el
coeficiente global de transferencia de calor U en lugar del coeficiente de transferencia de
calor h. Sin embargo, como ∆T varía con la posición en el intercambiador de calor, es
necesario trabajar con una ecuación de la forma:

7
 q = UA∆Tm (11.10)

Donde ∆Tm es una diferencia de temperatura media apropiada. La ecuación 11.10 puede ser
usada conjuntamente con las 11.7 y 11.8 para realizar un análisis del intercambiador de
calor, estableciendo la forma específica de ∆Tm.

Figura 11.10. Distribución de temperatura en un intercambiador de calor de flujo

paralelo

Par este tipo de intercambiadores se tiene que Th,i = Th,1, Th,o = Th,2, Tc,i = Tc,1 y Tc,o = Tc,2.
La forma de ∆Tm se puede determinar aplicando un balance de energía a los elementos
diferenciales en los fluidos caliente y frío.. Cada elemento es de longitud dx y área de
transferencia de calor dA, como se muestra en la figura 11.10. Para hacer el balance de
energía se consideran las siguientes suposiciones:
1. El intercambiador de calor está aislado térmicamente de los alrededores, en cuyo
caso el único intercambio de energía ocurre entre el fluido caliente y el frío.
2. La conducción axial a lo largo del tubo es despreciable.
3. Los cambios de energía cinética y potencial son despreciables.
4. Los calores específicos de los fluidos son constantes.
5. El coeficiente global de transferencia de calor es constante.

Aplicando un balance de energía a cada uno de los elementos diferenciales de la figura

11.10 se tiene:

*
dq = − m c p ,h dTh = C h dTh (11.11)

8
 *
dq = m c p ,c dTc = C c dTc (11.12)

La transferencia de calor a través del área dA se puede expresar como:

dq = U∆TdA (11.13)

donde ∆T = Th- Tc es la diferencia de temperatura local entre el fluido caliente y el frío.

Para determinar la forma integrada de la ecuación 11.13 se sustituyen las
ecuaciones11.11 y 11.12 en la forma diferencial de la ecuación 11.9

d (∆T ) = dTh − dTc

Para obtener

 1 1 
d (∆T ) = −dq − 
 Ch Cc 

Sustituyendo dq de la ecuación 11.13 e integrando a través del intercambiador se

obtiene

2 d (∆T )  1 1  2
∫1 ∆T
= −U  −  ∫ dA
 Ch Cc 
1

 ∆T   1 1 
ln 2  = −UA −  (11.14)
 ∆T1   Ch Cc 

Sustituyendo Ch y Cc de las ecuaciones 11.7 y 11.8, respectivamente, se sigue que:

 ∆T   Th ,i − Th,o Tc ,o − Tc ,i 
ln 2  = −UA +  = −
UA
[(Th,i − Tc,i ) − (Th,o − Tc,o )]
∆T
 1  q q  q

Reconociendo que para el intercambiador de flujo paralelo de la figura 11.10, ∆T1 =

(Th,i - Tc,i) y ∆T2 = (Th,o – Tc,o), se obtiene

∆T2 − ∆T1
q = UA
 ∆T 
ln 2 
 ∆T1 

q = UA∆Tml (11.15)

9
 ∆T1 − ∆T2 ∆T2 − ∆T1
Donde ∆Tml = = (11.16)
 ∆T1   ∆T 
ln  ln 2 
 ∆T2   ∆T1 

Es importante recordar que para un intercambiador de calor de flujo paralelo

∆T1 = Th ,1 − Tc ,1 = Th ,i − Tc ,i
(11.17)
∆T2 = Th , 2 − Tc , 2 = Th ,o − Tc ,o

Además, para este intercambiador Tc,o no puede exceder a Th,o. Para valores iguales de UA
y de temperatura

Intercambiador de calor en contraflujo:

Figura 11.11. Distribución de temperatura en un intercambiador de calor de

contraflujo
La distribución de temperatura asociada con este intercambiador se muestra en la figura
11.11. Note que la temperatura de salida del fluido frío puede exceder la temperatura de
salida del fluido caliente.
Se puede hacer el mismo análisis que se hizo para el intercambiador de flujo
paralelo y se puede demostrar que también aplican las ecuaciones 11.15 y 11.16. Sin
embargo, para el intercambiador de calor en contraflujo se tiene que:

∆T1 = Th ,1 − Tc ,1 = Th ,i − Tc ,o
(11.18)
∆T2 = Th , 2 − Tc , 2 = Th ,o − Tc ,i

10
Note que, para las mismas temperaturas de entrada y salida, la diferencia de temperatura
media logarítmica para contraflujo es mayor que para flujo paralelo,

∆Tml ,CF > ∆Tml , FP

Intercambiadores de calor de coraza y tubos y flujo cruzado

∆Tml = F∆Tml , FC (11.19)

F se obtiene de las figura 11.12.

Condiciones especiales de operación

Ch>>Cc Ch<<Cc
O vapor que se Lìquido que se
Condensa (Ch → ∞) Evapora (Cc → ∞)

Caso (a): Ch >>Cc o el fluido caliente se esta evaporando (Ch →∞) → cambio despreciable
en Th
(Th,o = Th,i)

Caso (b): Cc >> Ch o el fluido frio se està evaporando (Cc → ∞) → cambio despreciable en
Tc (Tc,o = Tc,i )

Caso ( c ): Ch = Cc → ∆T1 = ∆T2 = ∆Tml

11
 METODO DE LA EFECTIVIDAD – NUT

El método de la Diferencia de Temperatura Media Logarítmica se aplica en

problemas de diseño para los cuales se conocen los flujos másicos y las temperaturas de
entrada de los fluidos, así como una temperatura de salida uno de los fluidos. Con estos
valores se puede determinar el tamaño requerido del intercambiador (área superficial) y la
temperatura de salida que no se conoce.
Cuando se usa el MDTL para determinar las temperaturas de salida de los fluidos,
conociendo las temperaturas de entrada de dichos fluidos, el proceso de solución es
iterativo.
Par realizar los cálculos, tanto de diseño como de desempeño de intercambiadores
de calor, usando el Método de la efectividad – NUT no es necesario realizar iteraciones.

Definiciones

Efectividad del intercambiador de calor, e:

q
ε= (11.20)
q max

0 ≤ ε ≤1

Flujo de calor máximo posible:

q max = C min (Th ,i − Tc ,i ) (11.21)

C h → C h < C c
C min =  (11.22)
C c → C c < C h

Número de Unidades de Transferencia, NUT:

UA
NUT = (11.23)
C min

Es un parámetro adimensional cuya magnitud influencia el desempeño del intercambiador

de calor:

q ↑ con ↑ NUT

12
 Relaciones para intercambiadores de calor

 .

 q = m h h,i −ih, o )
(i

o
q = C (T − T )
 h h ,i h ,o

 .

q = mc c, o −ic,i )
(i


o

q =C (T −T )
 c c, o c, i


Cálculos del desempeño:

NUT = f (ε , C min / C max )

Para todos los intercambiadores:

ε ↑ si ↓ C r

Para Cr = 0 se usa la siguiente relación para todos los intercambiadores de calor:

ε = 1 − exp(− NUT )
o

NUT = − ln(1 − ε )

13
Tabla 9.1 Ecuaciones de eficiencia de intercambiadores de calor

14
Tabla 9.2 Ecuaciones del número de unidades de transferencia (NTU) de intercambiadores
de calor

15
Figura 11.12 Factor de corrección F

16
 Figura 11.13 Efectividad vs NUT

PROBLEMAS

Problema 11.5 del libro de Incropera.

Un dispositivo de recuperación de calor implica la transferencia de energía de los gases de
escape calientes por una región anular para presurizar agua que fluye por el tubo interno del
anillo. El tubo interno tiene diámetros interno y externo de 24 y 30 mm y està conectado
por ocho puntales a un tubo exterior aislado de 60 mm de diámetro. Cada puntal tiene 3 mm
de espesor y está construido de manera integral con el tubo interno de acero al carbón (k =
50 W/m.K)

17
Solución: Determinación de la transferencia por unidad de longitud de un dispositivo de
recuperación cuyos fluidos son agua y gases calientes.
Se conocen: la geometría del intercambiador de calor anular con aletas, la temperatura de
los gases y el coeficiente de convección, la masa del flujo de agua y la temperatura de la
misma.

Encontrar la velocidad de transferencia de calor por unidad de longitud

Do = 60 mm
Di,2 = 30 mm
Di,1 = 24 mm
t = 3 mm
L = (60-30)/2 mm

Suposiciones: (1) Condiciones de estado estable, (2) propiedades constantes,

(3)Conducción unidimensional, (4)superficie exterior adiabática, (5) radiación del lado del
gas despreciable, (6) flujo interno completamente desarrollado.

Propiedades (tabla A6) Agua a 300 K W

Pr := 5.83 kt := 50
m⋅ K
− 6 N⋅ s
µ := 855⋅ 10 D2 := 0.03m Do := 0.06m
2 D1 := 0.024m
m
(Do − D2) kg W
L1 := Long := 1m ma := 0.161 kagua := 0.613
2 s m⋅ K
W
h := 100
w := Long L1 = 0.015m t := 0.003m 2
m ⋅K

Anàlisis el flujo de calor es q (UAc)⋅(Tm, h − Tm, c)

18
 donde

1 1 1
+ Rw +
( U⋅ A )
c
( h⋅ A )
c ( η o⋅h⋅A )h

W 2 4⋅ ma
kagua = 0.613 Ac = 0.075m ReD :=
m⋅ K µ ⋅ D1⋅ π

0.5
2⋅ h
m1 :=  
3 1
ReD = 9.99 × 10  m1 = 36.515 D1 = 0.024m
 kt⋅ t  m

el flujo interno es turbulento por lo tanto

0.8 0.4
NuD := 0.023⋅ ReD ⋅ Pr NuD = 73.73 Af := 8⋅ 2⋅ ( L1⋅ w)

tanh ( m1⋅ L1)

2 D1 = 0.024m η f :=
Af = 0.24 m m1⋅ L1

η f = 0.911
(
At := Af + π⋅ D2 − 8⋅ t ⋅ w ) Rw = 7.103 × 10
−4 K
W
2
At = 0.31 m
 D2 
ln 
 D1  La eficiencia global es
Rw :=
( )
Af
η o := 1 − ⋅ 1 − ηf
2⋅ π⋅ kt⋅ Long At
3 W
hc = 1.883 × 10
η o = 0.931 m ⋅K
2

1
Uc⋅ Ac = 23.598 W
NuD⋅ kagua K
hc :=
D1
Ac := π⋅ D1⋅ Long 1
W Uc :=
Uc = 312.972
Ac ⋅  
1 1
+ Rw +
η o ⋅ h ⋅ At 
2
m ⋅K hc ⋅ Ac
 4 
q = 1.18 × 10 W
q := Uc⋅ Ac ⋅ ( 800K − 300K)

19
Problema 11.47:
Diseño de un intercambiador de calor de coraza y tubos, un paso por la coraza y dos por los
tubos para suplir vapor para la turbina de un sistema de conversión de energìa basado en un
ciclo de potencia estandar (Rankine). La planta genera 2 MW con una eficiencia de 3%. El
agua del océano entra al intercambiador a 300 K y se desea que salga del mismo a 292 K.
El fluido de trabajo del ciclo de potencia se evapora en los tubos a una temperatura de 290
K y el coeficiente global de transferencia de calor es conocido.
Encontrar: (a) el àrea del evaporador, (b) el flujo màsico de agua.

Soluciòn:

Suposiciones:
a) Pérdidas de calor del intercambiador de calor al medio ambiente despreciables.
b) Cambios de energía cinética y potencial despreciables
c) Propiedades de los fluidos constantes.

Propiedades para el agua (Tabla A-6) a Tm = 296 K : cp = 4181 J/kg.K

(a) La eficiencia es:

W& 2 MW
η= = = 0.03
q q

Por lo tanto el flujo de calor requerido es:

2MW
q= = 66.7 MW
0.03
También:

(300 − 290) − (292 − 290)

∆Tml ,CF = = 5C
300 − 290
ln
292 − 290

Y con P = 0 y R = ∞, de la figura 11.10, se obtiene F = 1. Entonces:

20
 q 6.67 x10 7 W
A= = = 11.1m 2
UF∆Tml ,CF 1200W / m 2 K .1.5C

b) La velocidad de flujo màsico a travès del evaporador es:

q 6.67 x10 7 W
mh = = = 1994kg / s
c p , h (Th ,i − Th ,o ) 4181J / kg.K (300 − 292) K

21