INSTITUTO POLITÉCNICO

NACIONAL
Escuela Superior De Ingeniería Química E Industrias
Extractivas

DEPARTAMENTO DE INGENIERÍA QUÍMICA INDUSTRIAL

DISEÑO DE EQUIPOS INDUSTRIALES

“DISEÑO TERMODINÁMICO DE UN INTERCAMBIADOR DE CALOR

PARA EL SISTEMA QUEROSENO 42° API DESTILADO 35° API”

INTEGRANTES:
 GUEVARA SAMANIEGO CYNTHIA
 LECONA ROMÁN ANGÉLICA

PROFESOR:
 INOCENCIO CASTILLO TERAN

Grupo: 3IM82
Periodo 22-2
Fecha de entrega: 12/03/2022
Realice el diseño termodinámico de un intercambiador de calor de tubos y coraza que resulte eficiente y
económico justificando las suposiciones que realice. Se permite una caída de presión de 10 PSI para ambas
corrientes. El sistema ha realizar es queroseno 42° API de 375°F a 220°F cuyo flujo másico es de 600,000
lb/h, con destilado 35° API que entra a 90°F y sale a 180°F.
1. Determinación de propiedades.
Condición Queroseno Destilado
Tentrada (°F) 375 90
Tsalida (°F) 220 180
Cp (BTU/lb °F) 0.57 0.5
S 0.77 0.82
m (Cps) 0.5 2.3
k(BTU/h ft2) 0.079 0.077
Rd (h ft2 BTU) 0.003 0.003

Q=Cp Gm ∆ T
lb
Q=(0.57 BTU /lb ° F)(600000 )(375 ° F−220 ° F )
h
Q=53,010,000 BTU /h

Q
G mdestilado=
Cp ∆ T
53,010,000 BTU /h
G mdestilado=
( 0.5 BTU /lb° F )(180 ° F−90 ℉ )
G mdestilado=1,178,000 lb/h

2. Localización de los fluidos

Realizando el análisis correspondiente de cada uno de los parámetros de los fluidos mediante la siguiente
tabla, se toma la decisión del fluido que circulara por los tubos y el fluido por la coraza.

Factor Tubos Coraza

Gm Destilado Queroseno
Pop
Corrosión
Temperatura Queroseno Destilado
Rd Queroseno Destilado
μ Queroseno Destilado
Toxicidad
ΔP
3. A continuación, se calcula la ∆ T ML y determinación de la dirección del flujo.

EN PARALELO:
∆ T 1=T 1−t 1

∆ T 1=375 ℉−90℉
 CONTRACORRIENTE:
∆ T 1=T 1−t 2

∆ T 1=375 ℉−180 ℉

∆ T 1=195 ℉

∆ T 2=T 2−t 1

∆ T 2=220 ℉−90℉

∆ T 2=130 ℉

∆ T 1−∆T 2
∆ T ML = C
∆T1
ln ln
∆T2

195 ℉−130 ℉
∆ T ML =
C
195℉
ln ln
130℉
∆ T ML =160.3097 ℉
C

Mayor eficiencia térmica en flujo de contracorriente, por lo tanto, se utiliza ese arreglo.

4. Cálculo de la ∆ T M L por correcciones de R y S.

V

T 1−T 2
R=
t 2−t 1
375 ℉−220 ℉
R=
180 ℉−90 ℉
 R=1.72

t 2−t 1
S=
T 1−t 1
180 ℉−90℉
S=
375 ℉−90℉
S=0.32

Con el gráfico 18, se encuentra el factor de corrección de diferencia de temperatura (FT) y se calcula ∆ T M L .
V

F T =0.91

∆ T M L =∆ T M L ∗FT
V c

∆ T M L =(160.3097℉ )(0.91)
V

∆ T M L =145.882℉
V

5. Determinación del área de transferencia de calor.

De la ecuación de Fourier:

Q
A=
UD ∆TML ¿ V

Fluido Fluido frío

caliente
Ligero Pesado
Queroseno Destilado
μ= 0.5 μ= 2.3
10 - 40

Con el análisis de viscosidades, se determina la transferencia de calor de la tabla 8.

S.O.P 10 – 40

10+40
UD= ¿
2
U D =25 ¿

53,010,000 BTU /h
A=
( 25 ) (145.882° F)
2
A=14,535.035 ft

POR EL LADO DE LOS TUBOS.

 1. Determinación de las condiciones y cálculo de gasto masa.

De tabla 10:
DE= 1 ¼” A
NT= '
L=16 (más comercial) a l∗L
PT= 19/16 (tabla 9) 2
Arreglo en cuadrado. 14,535.035 ft
NT= 2
16 BWG (0.3271 ft )(16)
a’l= 0.3271 ft2
N T =2777.2537 tubos=2777 tubos
a’f=0.985 in2
DI=1.12 in NT∗a ' f
Suponemos que n=1 aF=
144∗n
2
( 2777)(0.985¿ )
aF=
(144 )(1)
2
aF=18.9954 ft
Gm
GT =
aF
1,178,000 lb/ h
GT = 2
18.9954 ft No cumple, esta fuera
del rango 1X106 a
lb
GT =62,015.0141 2
2X106 , 10%
h ft
2. De acuerdo con la tabla 9 se propone el número de tubos para determinar el número de
cuerpos y con ello las nuevas condiciones.
a) CORRIDA 1.
NT
De tabla 10: N c= calculado

NT
DE= ¾” tabla

L=16 (más comercial) 2777

N c=
PT= 1 (tabla 9) 766
Arreglo en cuadrado. N c =3.62=4 cuerpos
16 BWG A R=N T ,tabla∗a' l∗L
a’l= 0.1963 ft2 A R=(766)(0.1963 ft 2 )(16)
a’f= 0.302 in2 2
DI= 0.62 in A R=2,450.853 f t
DS= 35 in
N=6 N T ∗a ' f
Nt=766 aF=
144∗n
2
(766)(0.302¿ )
aF=
144∗6
2
aF=0.2677 f t

Q / Nc
U DR=
A R¿ ∆ T M L V

BTU
(53,010,000 )/(4)
h
U DR= 2
(2405.853 f t )(145.882℉ )
 BTU Se encuentra dentro del rango
U DR=37.7595
f t 2℉ UD de 10 a 40 +- 10%

Gm /Nc
GT =
aF
1,178,000 lb/h/4
GT =
0.2677 f t 2
lb Cumple, dentro del rango
GT =1,100,112.066 2 1X106 a 2X106 , 10%
ht

GT ∗DI
N ReT =
μ∗12∗2.42
lb
N ReT =(1,100,112.066 2 )¿ ¿
ht
N ReT =10,211.844

Con el N ReT en la gráfica 24.

J H =43

Con el N ReT en la gráfica 26.

f =0.000255
TRAMO RECTO.
2
f ∗G T ∗n∗L∗12
∆ PTR = 10
5.22 x 10 ∗DI∗∅ t∗S

∆ PTR =(0.000255)¿ ¿
lb
∆ PTR =21.0069
¿2

RETORNO

( )
2
4∗n 62.4∗v
∆ PTr =
S 2 g∗144

Con GT en la gráfica 27.

( 62.4∗v 2
2 g∗144 )
=0.12
 ( 4 )( 4 )
∆ PTr = (0.12)
( 0.82 )
lb
∆ PTr =2.3414
¿2
Por tanto
∆ PT =∆ PTR +∆ PTr

lb lb
∆ PT =34.3285 2
+2.3414 2
¿ ¿
lb No cumple 10 PSI
∆ PT =36.6699 2
¿
b) CORRIDA 5.
NT
De tabla 10: N c= calculado

NT
DE= ¾”
tabla

2777
L=16 (más comercial) N c=
982
PT= 1 (tabla 9)
Arreglo en cuadrado. N c =3 cuerpos
13 BWG
a’l= 0.1963 ft2
a’f= 0.247 in2 '
DI= 0.56 in A R=N T ,tabla∗a l∗L
2
DS= 39 in A R=(982)(0.1963 ft )(16)
N=4 A R=3084.2656 f t
2

Nt=982

N T ∗a ' f
aF=
144∗n
( 982)(0.247 ¿2)
aF=
144∗4

aF=0.4211 f t 2
Q / Nc
U DR=
A R¿ ∆ T M L V

BTU
(53,010,000 )/(3)
h
U DR= 2
Se encuentra dentro del
(3084.2656 f t )(145.882℉ ) rango UD de 10 a 40 +-10%
BTU
U DR=39.272 2
ft ℉

Gm /Nc
GT =
aF
1,178,000 lb/h/( 3)
GT =
0.4211 f t 2
lb Cumple, dentro del rango
GT =932,478.43 2 1X106 a 2X106 , 10%
ht
 GT ∗DI
N ReT =
μ∗12∗2.42
lb lb
N ReT =(932,478.43 2 2 ) ¿ ¿
ht h t
N ReT =7818.1207

Con el N ReT en la gráfica 24.

J H =48

Con el N ReT en la gráfica 26.

f =0.00024
TRAMO RECTO
2
f ∗GT ∗n∗L∗12
∆ PTR = 10
5.22 x 10 ∗DI∗∅ t∗S

∆ PTR =(0.00024) ¿ ¿
lb
∆ PTR =6.6862
¿2

RETORNO

∆ PTr = (
4∗n 62.4∗v 2
S 2 g∗144 )
Con GT en la gráfica 27.

( )
2
62.4∗v
=0.12
2 g∗144

( 4 )( 4 )
∆ PTr = (0.12)
( 0.82 )
lb
∆ PTr =2.3415 2
¿

Por tanto
∆ PT =∆ PTR +∆ PTr

lb lb
∆ PT =6.6862 2
+ 2.3415 2
¿ ¿
Cumple con 10 PSI
 lb
∆ PT =9.028 2
¿
Cálculo de hio.

( ) ( k∗12
DI )
1
Cp∗μ∗2.42 0.14
hio =JH 3
ϕ s
k

( ) ( ( 0.077
0.56 )
1
( 0.5 ) ( 2.3 ) ( 2.42 ) ) 12
hio =48 3
(1)
0.077

BTU
hio =258.7454
h ft 2 ° F

POR EL LADO DE LA CORAZA. 

a) CORRIDA 1.
Distancia entre mamparas, se supone B=29
C=PT −DE
3
C=1− =0.25
4

B . C . DS
a s=
144 PT
( 29 ) ( 0.25 )( 39 )
a s=
144 ( 1 )
2
a s=1.9635 ft

Gm
Nc
Gs =
as

lb
600000 /(3)
h
Gs =
1.9635
lb
Gs =108045.4816 2
h ft
De la figura 28, con DE=0.75 in, PT=1 in y arreglo en cuadro.
Deq=0.95
G s Deq
N ReT =
μ(12)(2.42)

N ReT =
( 108045.4816
h ft )
lb
( 0.95 ) 2

( 0.5 )( 12 ) (2.42 )
 N ReT =7069.0914

De la figura 28:
JH´= 48
De la figura 29:
f’=0.0022
12 L
N +1=
B
12 ( 16 )
N +1=
29
N +1=6.621

f Gs2 ( N +1 ) DS
∆ Ps=
5.22 x 1010 Deq(φt )(S)

( )
2
( 0.0022 ) 108045.4816 lb 2 ( 6.621 ) ( 39 )
h ft
∆ Ps= 10
5.22 x 10 ( 0.95 ) ( 1 )( 0.82 )
∆ Ps=0.1631 PSI
Cálculo de ho.

[ ][ ]
1 /3
Cp μ (2.42) K (12) 0.14
ho=JH ´ φt
K Deq

[ ][ ]
1
( 0.57 )( 0.5 ) ( 2.42 ) 3 ( 0.079 ) ( 12 )
ho=48 ( 1)
0.079 0.95
BTU
ho=97.919 2
h ft ° F
Cálculo del coeficiente global de transferencia de calor.
(hio)(ho)
Uc=
hio+ ho
BTU BTU
(258.7454 2
)(97.919 2 )
h ft ° F h ft ° F
Uc=
BTU BTU
258.7454 2
+97.919 2
h ft ° F h ft ° F
BTU
Uc=71.0362 2
h ft ° F
Uc>U D

Cálculo del factor de incrustación.

Uc−U D
Rd=
Uc U D
 BTU BTU
71.0362 2
−39.272 2
h ft ° F ft ℉
Rd=
BTU BTU
(71.0362 2 )(39.272 2 )
h ft ° F ft ℉
Rd=0.01139
b) CORRIDA 3.
Distancia entre mamparas, se supone B=39
C=PT −DE
3
C=1− =0.25
4

B . C . DS
a s=
144 PT
( 39 ) ( 0.25 )( 39 )
a s=
144 ( 1 )
2
a s=2.6406 ft

G mQueroseno
Nc
Gs =
as
lb
600,000 /(3)
h
Gs =
2.6406 ft 2
lb
Gs =80341.512 2
h ft

De la figura 28, con DE=0.75 in, PT=1 in y arreglo en cuadro.

Deq=0.95
Gs Deq
N ℜs =
μ (12)(2.42)

N ℜs =
( 80341.512
h ft )
lb
( 0.95 )
2

( 0.5 )( 12 ) ( 2.42 )
N ReT =5256.5

De la figura 28:
JH´= 27
De la figura 29:
F=0.0031
12 L
N +1=
B
12 ( 16 )
N +1=
31
 N +1=4.923

2
f Gs ( N +1 ) DS
∆ Ps= 10
5.22 x 10 Deq(φt )(S)

( )
2
( 0.0031 ) 80341.512 lb 2 ( 4.923 ) ( 39 )
h ft
∆ Ps= 10
5.22 x 10 ( 0.95 )( 1 ) ( 0.82 )
∆ Ps=0.09448 PSI
Cálculo de ho.

[ ][ ]
1 /3
Cp μ (2.42) K (12) 0.14
ho=JH ´ φt
K Deq

[ ][ ]
1
( 0.57 ) ( 0.5 ) ( 2.42 ) 3 ( 0.079 ) (12 )
ho=27 ( 1)
0.079 0.95
BTU
ho=55.07943 2
h ft ° F
Cálculo del coeficiente global de transferencia de calor.
(hio)(ho)
Uc=
hio+ ho
BTU BTU
(258.7454 2
)(55.07943 2 )
h ft ° F h ft ° F
Uc=
BTU BTU
258.7454 2 +55.07943 2
h ft ° F h ft ° F
BTU
Uc=45.4124
h ft 2 ° F
Uc>U D

Cálculo del factor de incrustación.

Uc−U D
Rd=
Uc U D
BTU BTU
45.4124 2
−39.272 2
h ft ° F ft ℉
Rd=
BTU BTU
(45.4124 2
)(39.272 2 )
h ft ° F ft ℉
Rd=0.0034 Cumple con 0.001≤ Rd ≤ 0.005

h∗ft 2∗° F
No cumple con R D óptimo de 0.0025 Pero esta entre el rango de 0.001 a 0.005
BTU

Corrección de la viscosidad por el lado de los tubos.

a) CORRIDA 1.
Calculando la temperatura de la película.

T 1 +T 2
Tm=
2
375℉ +220 ℉
Tm=
2
Tm=297.5 ℉

t 1+ t 2
tm=
2
90 ℉ +180 ℉
tm=
2
tm=135℉

hi 0
tw=tm+ ∗(Tm−tm)
hi 0 +h 0
BTU
258.7454
h ft 2 ° F
t wsupuesto=135 ℉ + ∗(297.5℉ −135℉ )
BTU BTU
258.7454 2
+55.07943
h ft ° F h ft 2 ° F
tw=268.9796 ° F

De la figura 14, página 928 del libro “Procesos T.C. Kern” la μw =0.75
μ>1 ctp
Destilado, por el lado de los tubos con viscosidad=2.3, por lo tanto:

μ@ tm
ϕt =
μ @tw
2.3
ϕt =
0.75
ϕt =3.0666

( )
0.14
μ@ tm
ϕt 0.14 =
μ@ tw

( )
0.14
0.14 2.3
ϕt =
0.75

ϕt 0.14 =1.1698
Recalculando hi0

( ) ( k∗12
Deq )
1
Cp∗μ∗2.42 3 0.14
hi 0 =JH ϕ s
k

( ) ( ( 0.077
0.56 )
1
( 0.5 ) ( 2.3 ) ( 2.42 ) ) 12
hi 0 =48 3
(1.1698)
0.077
BTU
hi 0 =306.3216
h ft 2 ° F

Recalculando tw
hi 0
tw=tm+ ∗( Tm−tm )
hi 0 +h 0
BTU
306.3216 2
h ft ° F
tw c =135 ℉ + ∗(297.5 ℉−135 ℉)
BTU BTU
306.3216 2
+55.07943 2
h ft ° F h ft ° F
tw c =272.7341 ° F

Recalculando, hasta que tw sea constante

b) CORRIDA 2.

De la figura 14, página 928 del libro “Procesos T.C. Kern” la μw @ 272.7341=¿
μ@ tm
ϕt =
μ @tw
2.3
ϕt =
0.73
ϕt =3.1506

( )
0.14
0.14 μ
ϕt = @ tm
μ@ tw
0.14
ϕt =1.1743

Recalculando hi0

( ) ( k∗12
Deq )
1
Cp∗μ∗2.42 0.14
hi 0 =JH 3
ϕ s
k

( ) ( ( 0.077
0.56 )
1
( 0.5 ) ( 2.3 ) ( 2.42 ) ) 12
hi 0 =48 3
(1.1743)
0.077

BTU
hi 0 =307.4999
h ft 2 ° F
Recalculando tw
hi 0
tw=tm+ ∗( Tm−tm )
hi 0 +h 0
BTU
307.4999 2
h ft ° F
tw c =135 ℉ + ∗( 297.5 ℉−135 ℉ )
BTU BTU
307.4999 2
+55.07943
h ft ° F h ft 2 ° F
tw c =272.8146 ° F

tw c =constantes

Recalculando el coeficiente de transferencia de calor y factor de incrustación

hi 0∗h 0
U C=
h i0 + h0

BTU BTU
(307.49992
)(55.07943 )
h ft ° F h ft 2 ° F
U C=
BTU BTU
307.4999 2 +55.07943 2
h ft ° F h ft ° F
BTU
U C =46.7123
h ft 2 ° F
Uc>U D

U C −U DR
Rd=
U C ¿ U DR
BTU BTU
46.7123 2
−39.272 2
h ft ° F ft ℉
Rd=
BTU BTU
(46.7123 2
)(39.272 2 )
h ft ° F ft ℉
Cumple con 0.001≤ Rd ≤ 0.005
Rd=0.004
2
h∗ft ∗° F
No cumple con Rd óptimo de 0.0025 pero se encuentra en el rango.
BTU

La temperatura de la película se considera constante.

 LADO DE LOS TUBOS:
ARREGLO DE DI PT BWG L A´f A´l NT n Nc DS AR aF UDR GT NRET f Jh ΔPR Fig ΔPr ΔPT
27
1 CUADRADO 0.75 0.62 1 16 16 0.302 0.1963 766 6 4 35 2450.853 0.2677 37.7595 1100112 10211.84 0.00025 43 21.0069 0.12 2.3414 36.6699
2 CUADRADO 0.7 0.62 1 16 16 0.302 0.1963 676 6 5 33 2123.181 0.2362 34.2294 997032.4 9255.559 0.00031 38 13.3785 0.12 3.5121 16.8907
3 CUADRADO 0.75 0.62 1 16 16 0.302 0.1963 480 4 6 29 1507.584 0.2516 40.1720 780132.4 7241.6175 0.00031 29 5.4598 0.12 2.3414 7.8013
4 CUADRADO 0.75 0.56 1 16 16 0.247 0.1963 480 4 3 29 3084.266 0.4211 39.272 932478.4 7818.1207 0.00024 48 6.6862 0.12 2.3415 9.028

LADO DE LA CORAZA:
B C Deq As GS NRET F´ JH´ ho Uc N+1 ΔPS Rd
18 0.25 0.95 1.2187 492327.890 32211.5355 0.0017 110 226.024 120.6401 11 4.3471 0.0171
4
29 0.25 0.95 1.9635 108045.48 7069.091 0.0022 48 97.919 71.0362 7 0.1724 0.01139
39 0.25 0.95 2.6406 80341.512 5256.5 0.0031 27 55.0794 45.4124 5 0.0959 0.0034