Procesos de Transferencia de Calor

Alneira Cuellar Burgos

([email protected])

Evaluación de un Intercambiador TC 1-2

Presentado Por
Daniel Aguirre López, cód. 0316002 ([email protected])
Camila Andrea Barreto Cubides, cód. 0316007 ([email protected])
Jhony Alexander Galvis Corredor, cód. 0316030 ([email protected])
Manuel Uribe Villa, cód. 0316078 ([email protected]).
Alejandro Velásquez Castro, cód. 0316079 ([email protected])
Fecha
21 de Octubre, 2019

Objetivo
Se tiene un intercambiador TC 1-2, de 17 1/4 in de diámetro interno, con 182 tubos de
diámetro externo 3/4 in, BWG 14 y 12 ft de longitud, organizados en 4 pasos y arreglo
triangular con paso de 1 in, con deflectores espaciados 6 in. Según la terminología TEMA, el
equipo posee una configuración AES. Determinar si éste intercambiador se puede usar para
enfriar 78359 lb/h de Isobutano (118ºAPI) de 203ºF a 180ºF, calentando Butano (111.5ºAPI)
de 154ºF a 177ºF, calculando caídas de presión y exceso de área para el sistema.

Falta: nomenclatura, ubicación de los fluidos, bibliografía,

Nomenclatura
W, flujo másico del fluido caliente. ho, coeficiente de película del fluido en el
T1, temperatura de entrada del fluido caliente ánulo.
al intercambiador. UD, coeficiente global referido al área externa
T2, temperatura de salida del fluido caliente del tubo interno.
del intercambiador. Re, número de Reynolds.
Tavg, temperatura promedio del fluido Pr, número de Prandtl.
caliente. Nu, número de Nusselt.
Tc, temperatura calórica del fluido caliente. Di, diámetro interno del tubo interno.
w, flujo másico del fluido frío. Do, diámetro externo del tubo interno.
t1, temperatura de entrada del fluido frío al Ds, diámetro interno de la coraza.
intercambiador. BWG, espesor de la coraza.

1
t2, temperatura de salida del fluido frío del DE, diámetro equivalente para el ánulo.
intercambiador.
tavg, temperatura promedio del fluido frío. Además, excepto se exprese lo contrario, los
tc, temperatura calórica del fluido frío. subíndices c y f indican que el valor
tw, temperatura de la pared externa del tubo pertenece al fluido caliente o al fluido frío,
interno. respectivamente. Por ejemplo: Nuc significa
tp, temperatura de la pared interna del tubo que es el número de Nusselt para el fluido
interno. caliente. Además, los colores rojo y azul se
q, flujo de calor. usan para diferenciar entre caliente y frío,
Cp, capacidad calorífica. respectivamente. De igual manera, los
µ, viscosidad. subíndices i y o indican que el valor es
k, conductividad térmica. referente a interno o externo.
ρ, densidad.
At, área transversal de los tubos.
As, área transversal de la coraza.
G, gasto másico.
ΔPt, caída de presión en los tubos.
ΔPs, caída de presión en la coraza.
fF, factor de fricción de Fanning.

Contenido
Parte I. Propiedades de los Fluidos 3
Parte II. Descripción del Equipo Disponible 4
Parte III. Evaluación Térmica 11
Parte IV. Perfil de Temperaturas 13
Parte V. Evaluación Hidráulica 14
Referencias 17

Parte I. Propiedades de los Fluidos

Para encontrar las propiedades de los fluidos, es necesario definir la fase en la que éstos se
encuentran. Según el enunciado, los fluidos están reportados en grados API, por lo que
deben encontrarse en fase líquida. Debido a que las temperaturas de los fluidos ya están
definidas y tanto el Butano como el Isobutano son gases a condiciones de temperatura y
presión normales, ambos fluidos deben estar a mayor presión que la atmosférica. Por tanto,
se deben conocer las presiones de saturación, éstas son las presiones mínimas a las que
deben estar ambos fluidos para estar en fase líquida. Las propiedades necesarias son:
conductividad térmica (k), viscosidad (µ), capacidad calorífica (Cp), densidad (ρ) y presión de
saturación (Psat).

2
1. Conductividad térmica (Yaws):
a. Butano: En W/m·K si t en K.

𝑘 = 0.23840 − 3.7626 × 10−4 𝑡 − 2.5291 × 10−7 𝑡 2 (1)

b. Isobutano: En W/m·K si T en K.

𝑘 = 0.25340 − 2.2734 × 10−4 𝑇 − 5.7145 × 10−7 𝑇 2 (2)

2. Viscosidad (Yaws):
a. Butano: En cP si t en K.

674.1
Log 𝜇 = −6.859 + + 0.022𝑡 − 3.07 × 10−5 𝑡 2 (3)
𝑡

b. Isobutano: En cP si T en K.

1310
Log 𝜇 = −13.421 + + 0.0443𝑇 − 5.58 × 10−5 𝑇 2 (4)
𝑇

3. Capacidad calorífica (Kern, 1950):

a. Butano: En BTU/lb·ºF si t en ºF.

𝐶𝑝 = 0.0009𝑡 + 0.4923 (5)

b. Isobutano: En BTU/lb·ºF si T en ºF.

𝐶𝑝 = 0.001𝑇 + 0.5037 (6)

4. Densidad (Yaws):
a. Butano: En g/cm3 si t en K.

0.2863
𝑇
𝜌 = 0.2283 × 0.2724
−(1−
425.18
) (7)

b. Isobutano: En g/cm3 si T en K.

0.273
𝑇
𝜌 = 0.2228 × 0.2729 −(1−
408.14
) (8)

5. Presión de saturación (Yaws): La presión se halla por medio de la ecuación de Antoine

para cada fluido:
a. Butano: En mmHg si t en K.

3
 1904.9
Log 𝑃𝑠𝑎𝑡 = 27.044 − − 7.1805 Log 𝑡 − 6.6845 × 10−11 𝑡 + 4.2190 × 10−6 𝑡 2 (9)
𝑡

b. Isobutano: En mmHg si T en K.

1953.2
Log 𝑃𝑠𝑎𝑡 = 31.254 − − 8.806 Log 𝑇 + 8.9246 × 10−11 𝑇 + 5.7501 × 10−6 𝑇 2 (10)
𝑇

Parte II. Descripción del Equipo Disponible

Para el proceso se tiene disponible un intercambiador de tubo y coraza 1-2 con ciertas
características. Se debe determinar si éste equipo es adecuado para que los fluidos logren
llegar a la temperatura que se requiere. Las características del equipo son:

Tabla 1. Características de los tubos.

Tubos
Cantidad, Nt 182
Diámetro externo, Do 0.75 in
BWG 14; 0.083 in
Diámetro interno, Di 0.667 in
Longitud, L 12 ft
Número de pasos, n 4
Arreglo Triangular
Paso entre tubos 1 in
Coraza
Diámetro interno, Ds 17 1/4 in
Distancia entre deflectores, B 6 in

Además, el intercambiador está dispuesto en forma AES. Esto significa que el cabezal frontal
es tipo A:

Figura 1. Cabezal frontal tipo A.

El cabezal tipo A es fácil de reparar y reemplazar, y permite buen acceso al arreglo de tubos
para realizar la limpieza necesaria sin necesidad de desarmar el mismo, pero tiene
problemas de fugas al tener dos sellos. La coraza es tipo E:

4
 Figura 2. Coraza tipo E.

Ésta coraza es la más común entre todos los intercambiadores, los otros tipos son más
especializados para otros procesos. Y el cabezal posterior es tipo S:

Figura 3. Cabezal posterior tipo S.

Éste cabezal es especial para permitir la expansión térmica, pues es desplazable axialmente
y no tiene límites en cuánto a qué tanto puede expandirse. Sin embargo, es muy complejo
desmontar el cabezal para limpiar el arreglo de tubos. El intercambiador se ve de la siguiente
manera:

Figura 4. Vista 3D del intercambiador.

Se puede observar la forma AES del intercambiador, ubicándose el cabezal A en la parte

izquierda de la figura 4, de forma plana, y el cabezal S ubicado a la derecha con forma
redondeada. El arreglo en 2D se puede observar en la figura 5:

Figura 5. Arreglo en 2D del intercambiador.

5
Parte III. Evaluación Térmica
Al intercambiador entrarán Butano e Isobutano, fluidos entre los cuales se dará el
intercambio de calor, con el fin de enfriar 78359 lb/h del segundo. El flujo del primero es
desconocido, pero se puede calcular a partir del balance de energía de ambos fluidos,
conociendo la cantidad de calor que se transferirá:

𝑞 = 𝑊𝐶𝑝𝑐 𝑑𝑇
203º𝐹
lb
𝑞 = 78359 ∫ (0.001 + 0.5037𝑇)𝑑𝑇 (11)
h
180º𝐹
BTU
𝑞 = 1.2529 × 106
h

Con el valor del calor, se halla el flujo másico del Butano:

𝑞 = 𝑤𝐶𝑝𝑓 𝑑𝑡
𝑞 1.2529 × 106 BTU⁄h
𝑤= =
𝐶𝑝𝑓 𝑑𝑡 ∫177º𝐹(0.0009𝑡 + 0.4923𝑇)𝑑𝑡 (12)
154º𝐹
lb lb
𝑤 = 84951.543 ≅ 84952
h h

Adicional a conocer éste flujo másico, es necesario conocer las presiones de ambos fluidos
(usando las ecuaciones de Antoine) para luego ubicarlos dentro del intercambiador. Las
presiones se evalúan a la temperatura media de ambos fluidos, y se observan en las tabla 2
y 3:

Tabla 2. Propiedades del Butano (fluido frío) en general.

Propiedad Butano Valor
Flujo másico, w 84952 lb/h
Temperatura de entrada, t1 154ºF
Temperatura de salida, t2 177ºF
Temperatura promedio, tavg 165.5ºF
Presión de saturación, Psat, a t1 111.58 psi
Presión de saturación a t2 148.64 psi

Tabla 3. Propiedades del Isobutano (fluido caliente) en general.

Propiedad Isobutano Valor
Flujo másico, W 78359 lb/h
Temperatura de entrada, T1 203ºF
Temperatura de salida, T2 180ºF
Temperatura promedio, Tavg 191.5ºF
Presión de saturación, Psat, a T1 261.73 psi

6
Presión de saturación a T2 203.71 psi

Como se puede observar, las presiones de saturación de ambos fluidos son altas a las
diferentes temperaturas, pero, si se trabaja sobre éstas presiones para el diseño del
intercambiador, el Butano, fluido que se está calentando, cambiará de fase en cuanto cambie
su temperatura, por lo que debe ingresarse a mayor presión y garantizar que, definiendo una
caída de presión para el diseño, salga a una presión de igual manera mayor que la presión
de saturación calculada a t2. Para éste caso, se ha definido una caída de presión máxima de
15 psi. El siguiente paso es el cálculo de la media logarítmica de la diferencia de
temperaturas, usando la ecuación de flujo a contracorriente:

(𝑇1 − 𝑡2 ) − (𝑇2 − 𝑡1 )
𝑴𝑳𝑫𝑻 =
𝑇 −𝑡
Ln (𝑇1 − 𝑡2 )
2 1
(203º𝐹 − 177º𝐹) − (180º𝐹 − 154º𝐹)
𝑴𝑳𝑫𝑻 = (13)
203º𝐹 − 177º𝐹
Ln ( )
180º𝐹 − 154º𝐹
0
𝑴𝑳𝑫𝑻 =
0

Como se puede ver, el valor de MLDT es indeterminado; lo que se da porque las diferencias
de temperaturas son iguales en ambas terminales. Esto significa que la diferencia de
temperaturas a través de todo el intercambiador es constante, por lo que el perfil de
temperaturas es lineal. Esto se puede demostrar aplicando la regla de L’Hôpital en la
ecuación (13):

(𝑇1 − 𝑡2 ) − (𝑇2 − 𝑡1 ) 𝛥𝑡2 − 𝛥𝑡1

𝑴𝑳𝑫𝑻 = =
𝑇 −𝑡 𝛥𝑡
Ln (𝑇1 − 𝑡2 ) Ln (𝛥𝑡2 )
2 1 1
𝛥𝑡2
Definiendo: 𝒙 = ⟶ 𝜟𝒕𝟐 = Δ𝑡1 𝑥
𝛥𝑡1
𝛥𝑡2 − 𝛥𝑡1 Δ𝑡1 𝑥 − 𝛥𝑡1 𝛥𝑡1 (𝑥 − 1) (14)
𝑴𝑳𝑫𝑻 = = =
Ln(𝑥) Ln(𝑥) Ln(𝑥)
𝛥𝑡1
Derivando numerador y denominador con respecto a 𝑥: 𝑴𝑳𝑫𝑻 = = Δ𝑡1 𝑥
1⁄𝑥
𝛥𝑡2
𝑴𝑳𝑫𝑻 = Δ𝑡1 = Δ𝑡2 = 𝑇1 − 𝑡2
𝛥𝑡1

El valor del MLDT es, entonces, igual al valor de esa diferencia:

𝑴𝑳𝑫𝑻 = 𝑇1 − 𝑡2 = 𝑇2 − 𝑡1
𝑴𝑳𝑫𝑻 = 203º𝐹 − 177º𝐹 = 180º𝐹 − 154º𝐹 (15)
𝑴𝑳𝑫𝑻 = 26º𝐹

7
Éste valor se usa en el balance de energía dentro del intercambiador:

𝒒 = 𝑈𝐷 𝐴𝐷 𝐹𝑇 𝑀𝐿𝐷𝑇 (16)

En ésta ecuación hay un factor de corrección de la media logarítmica de temperaturas, pues,

dentro de los intercambiadores tubo y coraza de tipo 1-2, el flujo interno de los fluidos se
presenta de dos maneras: en paralelo y en contracorriente. En éste caso, en dos pasos de
tubos hay flujo paralelo, y en los otros dos pasos, el flujo es contracorriente. Por ésta razón,
el valor de MLDT, tomado como contracorriente, debe corregirse por un factor, expresado
como FT. Éste se calcula de la siguiente manera:

1−𝑆
(𝑟 2 + 1)0.5 Ln
𝑭𝑻 = 1 − 𝑟𝑆
2 − 𝑆[𝑟 + 1 − (𝑟 2 + 1)0.5 ] (17)
(𝑟 − 1) Ln
2 − 𝑆[𝑟 + 1 + (𝑟 2 + 1)0.5 ]

Donde r y S son, respectivamente:

𝑇1 − 𝑇2
𝒓= (18)
𝑡2 − 𝑡1
1−𝛼
𝑺= (19)
𝑟−𝛼

En la ecuación (19), alfa es:

1
𝑟𝑃 − 1 𝑁 (20)
𝜶=( )
𝑃−1

Donde N es el número de pasos por coraza, y P es:

𝑡2 − 𝑡1
𝑷= (21)
𝑇1 − 𝑡1

En el caso de que r tome el valor de 1, el cálculo de S cambia a:

𝑃
𝑺= (22)
𝑁 − 𝑃(𝑁 − 1)

Y el cálculo de FT cambia de igual manera a:

8
 𝑆√2
𝑭𝑻 =
2 − 𝑆(2 − √2) (23)
(1 − 𝑆) Ln
2 − 𝑆(2 + √2)

La ecuación (23) sale de aplicar la regla de L’Hôpital a la ecuación (17). Para el caso de
estudio:

𝑇1 − 𝑇2
𝒓=
𝑡2 − 𝑡1
203º𝐹 − 180º𝐹
𝒓=
177º𝐹 − 154º𝐹
𝒓=1

Dado que el valor de r es igual a 1, se usan las ecuaciones (22) y (23), en lugar de usar las
ecuaciones (17), (19) y (20). Primero, se calcula P:

𝑡2 − 𝑡1
𝑷=
𝑇1 − 𝑡1
177º𝐹 − 154º𝐹
𝑷=
203º𝐹 − 154º𝐹
𝑷 = 0.4694

Luego, se calculan S y FT. Sabiendo que sólo hay un paso por coraza:

𝑃
𝑺=
𝑁 − 𝑃(𝑁 − 1)
0.4694
𝑺=
1 − 0.4694(1 − 1)
𝑺 = 0.4694

𝑆√2
𝑭𝑻 =
2 − 𝑆(2 − √2)
(1 − 𝑆) Ln
2 − 𝑆(2 + √2)
0.4694 × √2
𝑭𝑻 =
2 − 0.4694(2 − √2)
(1 − 0.4694) Ln
2 − 0.4694(2 + √2)
𝑭𝑻 = 0.8522

De ésta forma, el valor de FT es 0.8522. Ya que éste debe ser mayor a 0.75, el equipo es
aceptable para continuar con la evaluación térmica. Sabiendo esto, se procede a calcular los

9
coeficientes de película para luego calcular el valor del coeficiente global de diseño, UD.
Primero, se calculan las áreas de flujo para los tubos y para la coraza, respectivamente:

𝜋𝑁𝑡 2
𝑨𝒕 = 𝐷 (24)
4𝑛 𝑖
𝐷𝑠 𝐶𝐵
𝑨𝒔 = (25)
𝑃𝑇

En las cuales Nt es el número de tubos, n es el número de pasos por tubos, C es el claro en

los tubos, B es el espacio entre bafles, y PT es la distancia entre dos tubos, medida desde el
centro de ambos:

Figura 6. Claro y PT en tubos.

Con las áreas calculadas, el paso siguiente es calcular los gastos másicos de ambos fluidos:

𝑚̇
𝑮𝒕 = (26)
𝐴𝑡
𝑚̇
𝑮𝒔 = (27)
𝐴𝑠

Para calcular éstos, se debe determinar primero la ubicación de los fluidos dentro del
intercambiador. Esto se decide mediante una serie de 8 condiciones establecidas por
(Mehra, D.K., 1983):

1. Viscosidad:
2. Toxicidad:
3. Tasa de flujo:
4. Corrosión:
5. Ensuciamiento:
6. Temperatura:
7. Presión:
8. Caída de presión:

El siguiente paso es calcular los valores de los números adimensionales: Reynolds, Prandtl y
Nusselt.

𝐷𝑖 𝐺𝑡
𝑹𝒆 = (28)
𝜇
𝜇𝐶𝑝
𝑷𝒓 = (29)
𝑘

10
Para Nu, se tiene una ecuación diferente para tubo y coraza, asumiendo que el regimen de
flujo de ambos es turbulento:

1 𝜇𝑏 0.14
𝑵𝒖 = 0.027𝑅𝑒 0.8 𝑃𝑟 3 ( ) (30)
𝜇𝑤
1 𝜇𝑏 0.14
0.55
𝑵𝒖 = 0.36𝑅𝑒 𝑃𝑟 3 ( ) (31)
𝜇𝑤

Como se puede observar, ambas ecuaciones tienen un factor de corrección por

viscosidades, siendo µb la viscosidad en el seno del fluido, y µw la viscosidad a temperatura
de pared. Éste factor se asume como 1 inicialmente, y se debe corregir luego de calcular las
temperaturas de pared a partir de los coeficientes de película:

11