PROCEDIMIENTO DE DISEÑO DE UN INTERCAMBIADOR DE CALOR

Diseño de un intercambiador de calor BES para calentar 500000 lb/h de un

Crudo desde 90°F hasta 149°F, utilizando un gasoil liviano que se enfría desde
330°F hasta 241°
Tabla 1. Datos de los fluidos

Fluido Crudo Gasoil Ligero

Fluido Frio Fluido Caliente
Propiedades
Entrada Salida Entrada Salida
Temperatura °F 90 149 330 241
Densidad relativa Sg(60°F /60°F) 0.876 0.876 0.742 0.742
Densidad lb/pie³
Viscosidad cP 5,923 3,046 0,135 0,207
Calor específico BTU/lb°F
Conductividad BTU/h pie°F
Flujo Másico lb/h 500000 500000
P Psia 421 415 105 101
∆P Permitido Psi 10 5
Masa Molecular
Factor de ensuciamiento BTU/hpie2ºF

Gasoil Liviano Gasoil Liviano

330°F 241°F

Crudo 90°F

Crudo 149°F

1. Determinación de las propiedades de los fluidos.

Cálculo de los ° API:
141.5
° 131.5
60/60
1
 °API Crudo=30.02 °API Gasoil = 59.2

Cálculo de la temperatura promedio de los fluidos:

TpromCrudo=

= = 119,5°F TpromGasoil=

= =285,5°F

Lectura de las propiedades de los fluidos:

3
Densidad: ρ(lb/pie ):

Se determina con la Tpromedio y con los API del fluido, su densidad relativa
en Tablas y Graficas Unidad I Figura 1 Gravedades específicas de
hidrocarburos

Densidad relativa= ρFluidoTProm /ρAgua60°F

ρFluidoTpromedio=Densidad relativa* ρAgua60°F ρAgua60°F = 62,4lb/pie3

Con Tprom Crudo=119,5 y 30 API Sg=0,85 ρCrudo 119,5° F = 53,04 lb/pie3

Con TpromGasoil=285,5 y 59,20 API Sg=0,634 ρGas 285,5°F = 39,56 lb/pie3

Conductividad térmica: K (BTU / Pie lb °F):

Se determina con la Tpromedio y con los API del fluido, su K en Tablas y

Graficas Unidad I Figura 2 Conductividades térmica de hidrocarburos líquidos.

Calor específico: Cp (BTU/lb°F):

Se determina con la Tpromedio y con los API del fluido, su Cp en Tablas y

Graficas Unidad I Figura 3 Calores específicas de hidrocarburos

Viiscosidad: µ(Cp)

En la Tabla 1 se muestran las viscosidades del crudo y del gasoil a la

temperatura de entrada y de salida, por lo tanto se realiza una interpolación de

2
 la viscosidad a la temperatura promedio de cada fluido, aplicando logaritmo
Neperiano (LN) tanto a las temperaturas como a las viscosidades. A
continuación se muestra el procedimiento para el gasoil.

Te=330 °F µTe = 0,135 cP

TProm=285,5 µTProm=?

Ts=241 °F µ Ts= 0,207cP

Aplicando LN tanto a las temperaturas como a las viscosidades e interpolando:

LN 330 = 5,799 ---------- LN 0,135 = - 2,002

LN 285,5 = 5,654 ---------- LNµTProm = X

LN 241= 5,484 ---------- LN 0,207= -1,575

5,799 − 5,484 = 0,315.---------−2,002 − −1,575 = −0,427

5,799 − 5,654 = 0,145 --------−2,002 − LNμ'()*+ = X

0,145 ∗ −0,427
−2,002 − LNμ'()*+ = = −0,197
0,315

−2,002 + 0,197 = LNμ'()*+ LNμ' /0,0 = −1,805

μ123*45' /0,0°6 = AntiLN −1,805 = 0,164 cP

Para el crudo se utiliza el mismo procedimiento obteniéndose una

μ=)>?* ,0°6 = 4,080 cP

En la Tabla 2 se muestran todas las propiedades del crudo y gasoil a la temperatura

promedio
Tabla 2. Propiedades de los Fluidos
Propiedades Crudo Gasoil Ligero
Temperatura (°F) 119,5 285,5
3
ρ(lb/pie ) 53,04 39,56
Viscosidad ( cP) 4.080 0,164
Cp (BTU/lbºF) 0,475 0,62
Conductividad (BTU/lbºF) 0,0752 0,0845
3
2. Factor de ensuciamiento

En Tablas y Graficas Unidad I Tabla 3 Factores de ensuciamiento PDVSA,

donde dice Miscelaneos: Calentadores de aceite combustible 5x10-3°Fpie2h/BTU y
para Líquidos orgánicos 1x10-3°Fpie2h/BTU.

3. Localización de los fluidos a través de los tubos y la carcasa.

Para determinar la localización de los fluidos a través de los tubos y la carcasa, se
consideró que aquel fluido con mayor ensuciamiento irá por los tubos.
Quedando los fluidos a través de los tubos y la carcasa de la siguiente forma:

Lado Tubo: Crudo ri=5x10-3 °Fpie2h/BTU

Lado Carcasa: Gasoil Ligero ro = 1x10-3 °Fpie2h/BTU

4. Calor Transferido:
El calor transferido en los hidrocarburos presentes en el intercambiador se evalúo
en función de las propiedades físicas del Crudo, haciendo un balance de energía
tenemos que:

Calor recibido por el fluido que entra a menor temperatura:

@ABCDEDFG HIJDFG HKDG @LKJFG mNLKJFG ∗ OPLKJFG ∗ ∆RLKJFG
5T V'W
Q=)>?* 500000 U ∗ 0,475 5TU
∗ 149 °F − 90°F QCrudo= 14012500 Btu/h

Flujo másico del fluido que entra a mayor temperatura:

\R]
@ KYZ[HBKDFG HIJDFG CYIDBZ B = @ABCDEDFG HIJDFG HKDG = 14012500
ℎ
@ KYZ[HBKDFG HIJDFG CYIDBZ B = @_Y[GDI = mN_Y[GDI ∗ OP_Y[GDI ∗ ∆R_Y[GDI

BTU BTU N
14012500 = mN ∗ d0,62 g ∗ (330°F − 241°F
h lbh

mN_Y[GDI =253941,65 lb/h

4
5. Coeficiente Global de transferencia Asumido:

Se define un intervalo de Uo para el flujo caliente: sustancia orgánica ligera y fluido

frío: sustancia orgánica pesada. Debido a estas condiciones se asume un valor de
Uo que estará entre los intervalos 10-40 BTU/hpie2°F. En Tablas y Graficas Unidad
I Tabla 4 Donald Kern Coeficientes globales típicos

Uoasumido = 30 BTU/hpie2°F

6. Temperatura Media Efectiva

6.1 Perfil de Temperatura:

Gasoil
T1=330 liviano
°F

t2=149 T2=241°F

Crudo
t1 =90 °F

6.2 Diferencia de Temperatura Media Logarítmica:

∆R − ∆R
hRij ∆
jk
∆

Observando el perfil de temperatura:

Terminal Caliente (Ambos fluidos a su temperatura mayor)= ∆R = R − l

Terminal Frio (Ambos fluidos a su temperatura menor) = ∆R = R − l

5
 T1 = Temperatura entrada del fluido caliente Gasoil (330°F).
T2 = Temperatura salida del fluido caliente Gasoil (241°F).
t1 = Temperatura entrada del fluido frio Crudo (90°F).
t2 = Temperatura salida del fluido frio Crudo (149°F).

DTML=165,55°°F.

6.3 Factor de corrección de la Temperatura Media Logarítmica:

Este factor de corrección se aplica para obtener la diferencia de temperaturas

media efectiva (DTM).
Los valores de este factor menores a 0.8 son en general inaceptables, debido a
que la configuración escogida del intercambiador se hace ineficiente, pues
termodinámicamente no es funcional.
Con los valores de R y S, entramos a la Figura 5 Factores de corrección MLTD
intercambiadores 1 – 2 (1 paso por carcasa- 2 o más por tubos) en Tablas y
Graficas Unidad I

' n' n o no n
R 1, 51 S 0,246
o no n ' no n

Ft=0.965

Se presenta un factor mayor a 0.8 por lo que es aceptable este valor para los
cálculos posteriores.

6.4 Temperatura Media Efectiva:

MTD=DTML*Ft

MTD=159,76°°F.

6
7. Área de transferencia de calor necesaria ( AO):
Este parámetro representa el área superficial que se requiere para transferir una
cantidad de calor (Q) a través de la superficie del metal.

Q U* ∗ A* ∗ DTM
Q
A*
U* DTM
Dónde:
Q = Calor generado en el intercambiador =14012500 BTU/h.
Uo= Coeficiente de transferencia de calor asumido de acuerdo a los fluidos que
circulan por el intercambiador=30 BTU/h pie2°F.

AO=2923,66 pie2

8. Características de los tubos :

Tabla 4 Características de los tubos

Diámetro Exterior (pulg) 1 (0.083 pie)
Calibre (BWG) 14
Espesor (pulg) 0.083*
Área de flujo por los 0.546*
2
tubos(pulg )
Longitud (pie) 20
Diámetro Interior (pulg) 0.834* (0.0695 pie)
Pitch(Pulg) 11/4
Claro (Pulg) = 1,25 pulg – 0.25 pulg
1 pulg
L/Dinterio=r 20 pie/0,0695 pie 287.77
Arreglo triangular(ro menor a
0.002pie2 h °F / Btu)
Material Acero al carbono
*Tabla 5 Información sobre los tubos Tablas y Graficas Unidad I

7
9. Número de Tubos:
A*
No
π ∗ D' ∗ L

Dónde:
AO = área de transferencia necesaria =2923,66 pie2.
L =Longitud de los tubos.=20 pies.
DT= Diámetro externo de los tubos= 0,083pie.

Ntubos=561

10. Diámetro interno de la carcasa:

En en Tablas y Gráficas Unidad I Tabla 6 Diámetro interno de la carcasa arreglo

triangular, se encontró que para 4 pasos por los tubos, Ntubos=562, Tubos de 1 pulg y
un arreglo triangular de Pt=11/4:
DIcarcaza= 35 pulg.

Recalculando con Ntubos=562 Ao k ∗ u ∗ h ∗ j 2930,86 pie2 y

Uo =29.93 BTU/h pie2°F

11. Selección de las dimensiones físicas del intercambiador de calor.

• El crudo circulará por los tubos.

• Carcasa Tipo E (un solo paso por la carcasa), establecida en el
enunciado
• Cabezal fijo extremidad anterior tipo B (Fluido limpio por la carcasa).
• Cabezal flotante extremidad posterior Tipo S (anillo divisorio).
• Diámetro interno de la carcasa 35 pulg.
• 4 Pasos por los tubos
• Bafles segmentados de 25%
• Espaciado de los Bafles 35pulg ( 35 pulg - 7 pulg)

8
12. Cálculo del área del fluido:

a.- Lado Carcasa:

h O ∗ O′ ∗ \
Avw
144
Dónde:
As= área de flujo lado carcasa, (pie2)
DIC= diámetro interno de la carcasa=35 pulg
C’= espaciado entre tubos =0.25 pulg.
B= espaciado de deflectores=35 pulg
Pt= distancia de centro a centro (Pitch)= 1.25 pulg.

As=1,701 pie2 por carcasa

b.- Lado Tubo:

k ∗ y′
=
z ∗ 144
Dónde:
At= área de flujo lado tubos, (pie2).
a’t= área de flujo por tubos= 0.546 pulg2.
n= número de pasos= 4.
Nt= número de tubos =562

At=0,533pie2 por tubo

13. Velocidad Másica:

a.- Lado Carcasa:

m3
G3 =
A3

9
 Dónde:

m= Flujo másico =253941,65 Lb/h.

As= Área de flujo por cada paso =1,701 pie2.

Gs=149289,62 Lb/pie2h

b.- Lado tubos:

mo
Go
Ao

Dónde:

m= Flujo másico =500000 lb/h.

At= Área de flujo lado tubo =0,533pie2.

Gt=938086, 304 Lb/pie2h

14. Velocidad lineal :

a.- Lado carcasa:

G3
V3
3600 ∗ ρ3

Dónde:

Gs = velocidad másica en la carcasa =149289,62 Lb/pie2h

ρs= densidad del gasoil =39,56 lb/pie3

VS=1,048 pie/s

10
 b.- Lado tubo:
Go
Vo
3600 ∗ ρo

Dónde:

Gt = velocidad másica en los tubos =938086,304 lb/pie2h

ρt= densidad del crudo =53,04 lb/pie3)

Vt= 4,913 pie/s

Sg= ρgasoil a T=119,5ºF/ρagua T=60ºF=0.85

€ , ,0
Velocidad máxima de diseño permitida= 10 × • • = 10,85 pie/s (Unidad I Pag
0 ,

33, Tabla 9.I Límites de erosión Velocidades máximas de diseño para fluidos dentro
de los tubos)

15. Diámetro equivalente de la carcasa:

En Tablas y Graficas Unidad I Tabla 7 Diámetro equivalentes lado carcasa para

varios diametros y arreglos de tubos

Dequivalente=De= 0.72pulg= 0.06 pie

16. Número de Reynolds :

a.- Lado carcasa:

†[ ∗ hB
…B[ =
‡[ 2.42

11
Dónde:
De = Diámetro equivalente =0.06 pie.
µs = viscosidad del fluido lado carcasa =0.164 cp
Gs= velocidad másica lado carcasa=149289,62 lb/pie2h

Res= 22569,48

b.- Lado tubos:

† ∗h
…B
‡ 2.42

Dónde:
Dt = Diámetro interior de los tubos= 0.0695 pie.
Gt = velocidad másica=938086,304 lb/pie2h
µt = viscosidad del fluido lado tubo= 4,080 cp

Ret=6603,16

17. Cálculos de los factores Jh

a.- Lado carcasa:

Para determinar este valor se utiliza la Figura 9 curva de transferencia de calor (jh)
lado carcasa con haz de tubos segmentados, en Tablas y Graficas Unidad I

Jh=89

b.- Lado Tubos:

Para determinar este valor se utiliza la Figura 10 curva de transferencia de calor (jh)
lado tubo , en Tablas y Graficas Unidad I
Jh=25
12
18. Número de Prandtl:
, ∗‰∗=Š
Prandtl =ˆ Œ
‹

a.- Lado carcasa:

Dónde:
Cps= capacidad calorífica lado carcasa =0.62 BTU/lb °F
µs= viscosidad del fluido por la carcasa =0.164Cp)
Ks= conductividad térmica del fluido =0,0845BTU/h pie °F

, ∗ , € LŽ∗ ,€ • •“IE‘’
Prandtl(carcasa)= • •= 2,912
, / 0 V'W/U Š4” °6

[Prandtl(carcasa)]1/3=1,428

b.- Lado Tubos:

Dónde:
Cpt= capacidad calorífica lado Tubos =0.475 BTU/lb °F
µt= viscosidad del fluido por los tubos =4,080 cP
Kt= conductividad térmica del fluido =0,0752 BTU/h pie °F

, –0• •“IE‘’∗ , / LŽ∗ ,

Prandtl(tubos) = • •= 62,366
, –0 V'W/U) Š4” °6

[Prandtl(tubo)]1/3=3,966

19. Cálculo del coeficiente promedio de trasferencia de calor por convección.

a.-Lado Carcasa.
ℎ— šCYKCY[Y • •
= ™ℎ ∗ •yzžlŸ
∅ , h›œ
13
 89 ∗ •0,0845 \R]“P¡› ℉/P¡›• ∗
ŽJI£
ℎ— ∗ 1,428
ŽDB
∅ , 0,72

ℎ—
= 178,988 \R]“
∅ , P¡› ℉ℎ

b.- Lado Tubos.

ℎD š JEG[ • •
= ™ℎ ∗ ∗ •yzžl JEG[
∅ , hD

ℎ¡ 25 0,0752 \R]⁄ℎ℉¤l ∗ 12 3,966

= = 107,282
∅ , 0,834

Corregimos el coeficiente de transferencia al diámetro exterior de la siguiente manera:

= 107,282 \R]“
¦§ ¨ ¬D
ℎ℉P¡› ¬B
*
∅©,ª«

h¡ 0,834 P-Ÿ
= = 0,834
h› 1 P-Ÿ

ℎ¡—
= 89,473 \R]“
∅ , ℎ℉P¡›

Temperatura de la pared
•’G“∅•
tw= tprom + ∗ ®RŽKG‘ − lŽKG‘ ¯
•’DG“∅ ’G“∅•

Dónde:
tw=Temperatura de la pared del tubo=°F
tprom= temperatura promedio del fluido frío(crudo) =119,5ºF.
14
Tprom= temperatura promedio del fluido caliente(gasoil)=285,5ºF.

178,988 \R]“
’G

∅°, ± P¡› ℉ℎ

ℎ¡—
= 89,473 \R]“
∅ , ℎ℉P¡›

d –/, // • •“ g
H ℉
tw= 119,5℉+ ∗ 285,5 − 119,5 ℉
d/ , – • •“ –/, // • •“ g
’℉H H ℉

tw= 230,175 ℉
μs ,
∅² ©,ª«
=d g
μw

µw=0,221 cP (Se busca con el mismo procedimiento anterior extrapolando)

µs=0.164 cP

©, ª´«µ¶ ©,ª«
∅² ©,ª«
=d g
©, ··ª µ¶
∅²©,ª« = ©, ¸¹¸

Con este valor conseguimos ho:

ho=178,988BTU“ ∗ 0,959 = 171,667 BTU“

h℉pie h℉pie

μt ,
∅º ©,ª«
=d g
μw

µw=1,714 cP (Se busca con el mismo procedimiento anterior extrapolando)

µt=4.080 cP

15
 «, ©»©µ¶ ©,ª«
∅º ©,ª«
d g
ª, ¼ª« µ¶
∅º©,ª« ª, ª·¸

Con este valor conseguimos hio:

hio=89,473BTU“ ∗ 1,129 = 101,024 BTU“

h℉pie h℉pie

20. Resistencias al ensuciamiento:

Este parámetro se refiere a cualquier capa o depósito de materias extrañas en

una superficie de transferencia de calor. Normalmente, esos materiales tienen
baja conductividad térmica, lo que provoca una resistencia mayor a la
transferencia de calor. Las características del ensuciamiento como la textura,
espesor, naturaleza, conductividad térmica, dureza, etc., depende del fluido
involucrado.

’ŽDB °½ ’ŽDB °½
ro = 0,001 ri =0,005
• • • •

’ŽDB °½ ¾¿
rio= 0,005
• • ¾§

·
’ŽDB °½ ª ¶ÀÁ -3¦Ä∗¶§¿ ∗℉
rio= 0,005 = 5,995x10
• • ©.»Ã« ¶ÀÁ ÅÆÇ

21. Calculo de la resistencia de la pared del tubo

22.
Ê*nÊ4
rÉ =
ËÌ

16
Dónde:
U)∗Š4” ∗℉
rw= resistencia ofrecida por el material del tubo ( V'W

Do=Diámetro externo del tubo=0,083 pie (1 pulg)

Di=Diámetro interno del tubo=0,0695 pie (0,834 pulg)
Kw=Conductividad térmica del acero al carbono=26BTU“h℉ pie Tabla 8

Conductividades térmicas En Tablas y Graficas Unidad I (Tabla 1.1 Pag. 8

Unidad I guía de teoría)

U)∗Š4” ∗℉
rw= 0,0005192
V'W

23. Cálculo del coeficiente de calor total Uo:

Constituye la suma de los coeficientes individuales de transferencia de todos los

elementos y fluidos colocados en la dirección en la cual fluye el calor:

1
U*
+ r* + rÉ + r4* +
UÍ UÎÍ

1
U* =
+ 0,001 + 0,0005192 + 0,005995 +
– ,€€– ,

1
U* =
0,00583 + 0,001 + 0,0005192 + 0,005995 + 0,00990

Uo= 1/0,02324=43,021 ÏÐÑ“

Ò℉ÓÔÕ·

Ö 0 • •/’
AoCalculado= = = 2038,87P¡›
¬ ×•Ø 0 ,–0℉∗ , • •/’℉ŽDB

UoReq= 29,93 BTU/h pie2°F AOReq=2930,86 pie2

17
 Uactual > Urequerido Sobrediseñado
Uactual < U requerido Subdiseñado

ÇÙÚºÀÙÁnÇÄ¿ÛÀ¿Ä§Ü¨
% Sobre/sub diseño= × ª©©
ÇÄ¿ÛÀ¿Ä§Ü¨

, n ,
% Sobrediseño ∗ 100 = 43,74%
,

,/€ n /,/–
% Sobrediseño= ∗ 100 = 43,75%
/,/–

24. Caída de presión en la carcasa:

2
fs * Gs * Ds ( Nc + 1)
∆Ps =
5.22 E 10 * De * Sg * φ s

En Tablas y Graficas Unidad I Figura 12 Factores de fricción por carcasa para haz
de tubo con deflectores, con Res= 22569,48 y un 25% de corte en los bafles se
obtiene un fs=0.0023 pie2/pulg2.

Gs= 149289,62 Lb/pie2h

Ds=35pulg= 2,917 pie
Nc+1=12*L/B=12*20 pie/35pulg=6,857=7
De= 0.72pulg= 0.06 pie
Densidad relativa en la Figura 1 Gravedades específicas Tablas y Gráficas Unidad
I (en orden), con la con Tp=285,5 ºF y 59,2 API

Sg=ρgasoil a T=285,5ºF/ρagua T=60ºF = 0.634

18
 φs = 0,959

Sustituyendo los valores en la ecuación anterior tenemos:

∆Ps=0,550 psi

Menor ∆Ps permitido de 10 psi (Tabla 10.I Caídas de presión típicas Pag 60
Unidad I)

26. Caída de presión en los tubos:

2
f * Gt * L * N
∆Pt =
5.22 E10 * Di * Sg * φt

Con Ret=6603,16 y tubos para intercambiadores se lee ft=0.0003 pie2/pulg2 en la

Figura 13 Factores de fricción lado tubo, en Tablas y Graficas Unidad I

Gt= 938086,304 lb/hpie2

L=20 pie
n=4
Di=0.0695 pie

Densidad relativa en la Figura 20 (Kern) Gravedades específicas Tablas y Gráficas

Unidad I Página 35, con la Tp=119,5 ºF y 30 API

Sg= ρCrudo a T=119,5ºF/ρagua T=60ºF=0.85

φt =1.129
19
 ∆Pt=6,066 psi.

25. Caída de presión en los cabezales del intercambiador:

 4 nV 2   62.5 
∆ Pr =  * 
 Sg 2 g   144 
 

Con Gt=938086,304 Lb/pie2h en la Figura 14 Pérdida de presión por retorno lado

tubo, en Tablas y Graficas Unidad I

V2  62 .5 
 
 2g  144  = 0.12
 

n=número de pases=4 y Sg=0,85

∆Pr=2,259 psi

26. Caída de presión total:

∆Ptotalt=∆
∆Pt + ∆Pr=8,325 psi

Menor ∆Ps permitido de 10 psi (Tabla 10.I Caídas de presión típicas Pag 60 Unidad
I)

27. Clasificación del intercambiador de calor de tubo y carcasa:

Según normas TEMA:

20
TEMA Clase R (petróleo y aplicaciones relacionadas)
Es un intercambiador BES (de cabezal fijo anterior B, un paso por la carcasa E
y cabezal posterior flotante con anillos divisorio S).

Según su servicio:
Calentador. Es un intercambiador de calor que aumenta la entalpia de una
corriente sin que ocurra un cambio de fase.

21