Diseño de Intercambiadores
Diseño de Intercambiadores
Diseño de Intercambiadores
Taller No.1
An old shell and tubes heat exchanger is about to be sold as scrap metal because it is at the
end of its useful life. The following are the specifications of this equipment:
Shell: Tubes:
Inside diameter: 23¼” OD: ¾” – BWG 17 – Pitch: 1” – Layout: 60º
Baffles spacing: 12” Fixed tubes sheet
Baffles cut: 25% (Type: segmental) Passes: 4
Passes: 1 No.: 360
Nozzles diam: 6” Length: 11 ft.
Material: SS 316L Nozzles diam: 4”
Material: Monel
Nonetheless, the process engineer thinks that this measure may be a little extreme and the
heat exchanger can be revamped for being re-used for heating, cooling, condensing, or
boiling some stream of the current process. He meets with other processes engineers for
finding a suitable solution and he asks that they propose three potential alternatives where
this old heat exchanger can be used. In this meeting, three ideas came up:
These alternatives arose because the engineers identified that the heat exchangers
associated to these sub-processes are close to their maintenance dates. They thought that
one of them can be replaced by the revamped heat exchanger to avoid a plant shutdown
with the consequent costs. As design team, the following activities were suggested:
Test every alternative on the old heat exchanger, in order to evaluate how the current
equipment works without modifications.
Modify the setup of the existing heat exchanger (number of tubes, OD, BWG, pitch,
layout, baffles spacing and cut) to improve the heat transfer performance for each
alternative.
Show the results using the standard TEMA sheet, highlighting the most relevant and
useful information to take an adequate decision.
Display, as a ranking, from the best option to the worst one, the evaluation of each
alternative (including before and after modifications).
In a written report, depicts the calculations (handmade or with a spreadsheet) for the best three
options and compare the results obtained with a design software. Explain the reasons for
the discrepancies between results. Finally, justify the best selection from an economic and
technical point of view.
1. Test every alternative on the old heat exchanger, to evaluate how the current equipment
works without modifications.
Para evaluar el uso alternativo del intercambiador en la operación de enfriamiento se uso una hoja
de Excel y se compararon los valores obtenidos en la herramienta de office con la simulación en
CHEMCAD.
TEMPERATURE
T1 220 °F
T2 180 °F
t1 120 °F
t2 145 °F
Table 1. Temperaturas de entrada y salida.
Para los cálculos es importante definir cual seré el fluido frio y el caliente. Los resultados fueron
consignados en las tablas a continuación.
Fluid properties
COLD
Fluid Water
Temperature 328,98 K Fluid properties
Cp 0,9985 BTU/lb*°F HOT
Viscosity 1,2272 lb/ft*h Fluid Toluene
k 3,7183 BTU/h*ft*°F Temperature 366,48 K
Density 60,6408 lb/ft^3 Cp 0,4608 BTU/lb*°F
Viscosity 0,6883 lb/ft*h
W lb/h k 0,6688 BTU/h*ft*°F
T1 120 °F Density 49,8026 lb/ft^3
T2 145 °F
Tprom 132,5 °F W 33000 lb/h
Table 2. Propiedades del Agua T1 220 °F
T2 180 °F
Tprom 200 °F
Table 3. Propiedades del Tolueno
Paso 1. Estimate the heat duty from de condition of the hot fluid.
El calor transferido este dado por la ecuación 1, aclarando que W es el flujo masico del fluido
caliente.
Para este paso es importante identificar la configuración de nuestro equipo, en este caso es
un intercambiador de calor de tubos y coraza. La media logarítmica de temperatura esta
dada por la ecuación 3,
( T 1−T 2 )
S= ( a−1
a−1 )
[ 8]
R= [5]
( t 2−t 1 ) 1−S
Fc=
√ R 2+ 1 ln ( 1−RS ) [9 ]
( t 2−t 1 )
2−S (R+1−√ R 2+1)
P=
( t 2−t 1 )
[6 ]
1/N
(R−1) ln
[ 2−S( R+1+ √ R2 +1) ]
( 1−RP )
a=
( ( 1−P ) ) [7]
Los resultados obtenidos al realizar los cálculos en Excel son los siguientes:
En primer lugar, debemos hallar el área de transferencia pues está hace parte de las
variables de la ecuación de Ud. Tenemos que el área está dada por la ecuación 10 y Ud por
la ecuación 11.
A=NtπLd o [10]
Q
Ud= [11 ]
A ∆ Tc
Los resultados obtenidos en Excel
Para el cálculo del hi debemos conocer Gt, Re, Pr, hi y hio; estas variables están dadas por las
ecuaciones 12, 13, 14 y 17. El hi depende del Re, si este es menor o mayor que 2100, las ecuaciones
15 y 16 nos ayudan a ver un poco la diferencia.
w Cp μ
Gt = [12] Pr= [14]
a' k
d i Gt di
ℜ=
μ
[13] hio =hi ( )
do
[17]
hi d i d hi d i
k [ ( )]
=1.86 RePr i [ 15 ]
L k
0,8 1 /3
=1.86 ( ℜ Pr ) [16]
También es importante conocer el área de cada tubo (a’) y el número de tubos por paso (n’), que se
obtienen por las expresiones 18 y 19, respectivamente
Para calcular la caída de presión, debemos de forma individual en los tubos (Ec. 32) y la coraza (Ec.
33). también debemos conocer la velocidad (Ec. 36).
En los cálculos anteriores el valor de f era supuesto (0,003 h*ft^2°F/BTU) pero para hallar el valor
preciso de caída de presión debemos calcularlo para el caso de los tubos. Es importante tener en
cuenta que el f depende del régimen del fluido (laminar o tubulento) y este lo define el Reynolds, en
el caso de ser laminar Re<2100 (Ec. 30) y de ser turbulento Re>2100 (Ec. 31)
64 0.4137
f= [30] f= [31]
144 ℜ 144 ℜ0.2585
La caída de presión en los tubos va a estar dada por la suma (Ec.32) de la caída de presión en los
tubos (Ec.32a) y en los retornos( Ec . b)
B
[
f =144 f 1−1.25 1−
( ) ]
( f −f ) [34 ]
ds 1 2
f d s G2s (N c +1)
∆ p s= 10
[33]
5.22 x 10 d s ( ρ/62.4)
12 L
N c +1=
B
De estos cálculos en Excel obtenemos
Finalmente, este mismo intercambiador de calor fue simulado en CHEMCAD y los resultados
obtenidos fueron los siguientes
SUMMARY REPORT
--------------
General Data: Heat Transfer Data:
Exch Class/Type R/AEL Effective Transfer Area
758.77
Shell I.D. 1.94 Area Required 110.43
Shell in Series/Parallel 1/1 COR LMTD
64.66
Number of Tubes 360 U (Calc/Service) 47.39/12.39
Tube Length 11.00 Heat Calc
2.33
Tube O.D./I.D. 0.0625/0.0528 Heat Spec
0.61
Excess %
282.39
Tube Pattern TRI60 Foul(S/T) 1.500E-
003/0.000E+000
Tube Pitch 0.08 Del P(S/T) -0.95/0.10
Number of Tube Passes 4 SS Film Coeff
71.55
Number of Baffles 9 SS CS Vel
0.24
Baffle Spacing 1.00 TW Resist
0.000854
Baffle Cut, % Area 25 TS Film Coeff
247.94
Baffle Type SSEG TS Vel
0.56
Baffle space def. Edge-Edge
Thermodynamics:
K: NRTL
H: Latent Heat
D: Library
Number of Components: 2
Engineering Units:
Temperature F
Flow/Hour lb/h)/h
Pressure psia
Enthalpy MMBtu
Diameter/Area ft/ft2
Length/Velocity ft/(ft/sec)
Film Btu/hr-ft2-F
Fouling hr-ft2-F/Btu
Podemos notar variación en los resultados, aunque sea el mismo intercambiador de calor. Esto se
debe a que el programa de simulación tiene en cuenta factores como el corte de los bafles, en
cambien para los cálculos en Excel no
Condensing 4,500 lb/h of saturated formaldehyde (@50 psig)
SUMMARY REPORT
--------------
General Data: Heat Transfer Data:
Exch Class/Type R/AEL Effective Transfer Area
764.29
Shell I.D. 1.94 Area Required
700.83
Shell in Series/Parallel 1/1 COR LMTD
43.23
Number of Tubes 360 U (Calc/Service) 45.39/41.43
Tube Length 11.00 Heat Calc
1.49
Tube O.D./I.D. 0.0625/0.0528 Heat Spec
1.37
Excess %
9.05
Tube Pattern TRI60 Foul(S/T) 1.500E-
003/1.500E+003
Tube Pitch 0.08 Del P(S/T) -0.10/28.33
Number of Tube Passes 4 SS Film Coeff
304.53
Number of Baffles 9 SS CS Vel
2.50
Baffle Spacing 1.00 TW Resist
0.000889
Baffle Cut, % Area 25 TS Film Coeff
76.34
Baffle Type SSEG TS Vel
5.08
Baffle space def. Edge-Edge
Thermodynamics:
K: NRTL
H: Latent Heat
D: Library
Number of Components: 2
Engineering Units:
Temperature F
Flow/Hour lb/h)/h
Pressure psia
Enthalpy MMBtu
Diameter/Area ft/ft2
Length/Velocity ft/(ft/sec)
Film Btu/hr-ft2-F
Fouling hr-ft2-F/Btu
CHEMCAD 7.1.2 25/07/2020
4:02 PM
Para el caso del condensador, nos damos cuenta de que él %de exceso es muy pequeño. Sin
embargo, el flujo requerido de flujo masico del Etilenglicol (Fluidos utilizado para condensar el
formaldehido) es muy alto. Esto puede tener relación con la densidad de los fluidos.
SUMMARY REPORT
--------------
General Data: Heat Transfer Data:
Exch Class/Type R/AEL Effective Transfer Area
905.83
Shell I.D. 1.94 Area Required
311.20
Shell in Series/Parallel 1/1 COR LMTD
88.78
Number of Tubes 360 U (Calc/Service) 141.89/48.75
Tube Length 11.00 Heat Calc
11.41
Tube O.D./I.D. 0.0625/0.0528 Heat Spec
3.92
Excess %
191.08
Tube Pattern TRI60 Foul(S/T) 1.500E-
003/1.500E+003
Tube Pitch 0.08 Del P(S/T) -0.18/3.41
Number of Tube Passes 4 SS Film Coeff
922.65
Number of Baffles 9 SS CS Vel
0.27
Baffle Spacing 1.00 TW Resist
0.000778
Baffle Cut, % Area 25 TS Film Coeff
393.79
Baffle Type SSEG TS Vel
6.84
Baffle space def. Edge-Edge
Thermodynamics:
K: NRTL
H: Latent Heat
D: Library
Number of Components: 2
Engineering Units:
Temperature F
Flow/Hour lb/h)/h
Pressure psia
Enthalpy MMBtu
Diameter/Area ft/ft2
Length/Velocity ft/(ft/sec)
Film Btu/hr-ft2-F
Fouling hr-ft2-F/Btu
2. Modify the setup of the existing heat exchanger (number of tubes, OD, BWG, pitch,
layout, baffles spacing and cut) to improve the heat transfer performance for each
alternative.
Para las modificaciones del intercambiador para el caso del enfriamiento de tolueno se usó una hoja
de Excel y se compararon los valores obtenidos en la herramienta de office con la simulación en
CHEMCAD.
TEMPERATURE
T1 220 °F
T2 180 °F
t1 120 °F
t2 145 °F
Table 12. Temperaturas de entrada y salida.
Paso 1. Estimate the heat duty from de condition of the hot fluid.
SUMMARY REPORT
--------------
General Data: Heat Transfer Data:
Exch Class/Type C/AEL Effective Transfer Area
423.38
Shell I.D. 1.94 Area Required
207.27
Shell in Series/Parallel 1/1 COR LMTD
64.66
Number of Tubes 200 U (Calc/Service) 45.37/22.21
Tube Length 11.00 Heat Calc
1.24
Tube O.D./I.D. 0.0625/0.0519 Heat Spec
0.61
Excess %
104.26
Tube Pattern SQUAR Foul(S/T) 1.000E-003/1.000E-
003
Tube Pitch 0.08 Del P(S/T) -0.97/3.49
Number of Tube Passes 4 SS Film Coeff
62.68
Number of Baffles 8 SS CS Vel
0.27
Baffle Spacing 1.16 TW Resist
0.000943
Baffle Cut, % Area 25 TS Film Coeff
409.28
Baffle Type SSEG TS Vel
1.04
Baffle space def. Edge-Edge
Thermodynamics:
K: NRTL
H: Latent Heat
D: Library
Number of Components: 2
Engineering Units:
Temperature F
Flow/Hour (lb/h)/h
Pressure psia
Enthalpy MMBtu
Diameter/Area ft/ft2
Length/Velocity ft/(ft/sec)
Film Btu/hr-ft2-F
Fouling hr-ft2-F/Btu
CHEMCAD 7.1.2
Condensing 4,500 lb/h of saturated formaldehyde (@50 psig)
Para este caso, al obtener un %exceso tan bajo decidimos no realizar modificaciones .
SUMMARY REPORT
--------------
General Data: Heat Transfer Data:
Exch Class/Type C/AEL Effective Transfer Area
642.43
Shell I.D. 1.94 Area Required
319.68
Shell in Series/Parallel 1/1 COR LMTD
88.31
Number of Tubes 300 U (Calc/Service) 138.37/69.09
Tube Length 11.00 Heat Calc
7.88
Tube O.D./I.D. 0.0625/0.0519 Heat Spec
3.92
Excess %
100.96
Tube Pattern SQUAR Foul(S/T) 1.000E-003/1.000E-
003
Tube Pitch 0.08 Del P(S/T) -0.17/3.51
Number of Tube Passes 1 SS Film Coeff
843.43
Number of Baffles 8 SS CS Vel
0.32
Baffle Spacing 0.00 TW Resist
0.000943
Baffle Cut, % Area 30 TS Film Coeff
387.22
Baffle Type NOBF TS Vel
7.18
Baffle space def. Edge-Edge
Thermodynamics:
K: NRTL
H: Latent Heat
D: Library
Number of Components: 2
Engineering Units:
Temperature F
Flow/Hour (lb/h)/h
Pressure psia
Enthalpy MMBtu
Diameter/Area ft/ft2
Length/Velocity ft/(ft/sec)
Film Btu/hr-ft2-F
Fouling hr-ft2-F/Btu
CHEMCAD 7.1.2
3. Show the results using the standard TEMA sheet, highlighting the most relevant and
useful information to take an adequate decision.
Las hojas TEMA C para cada intercambiador de calor se realizaron en Excel teniendo en cuenta la
información obtenida de la simulación de cada intercambiador en CHEMCAD. El libro de Excel
con cada hoja tema será anexado el momento de enviar el informe. Las variables que tendremos en
cuenta para la comparación entre los intercambiadores son: Calor transferido (Q), Area (A), Area
requerida (Areq), %de exceso (%Exc) y la caída de presión (ΔP).
4. Display, as a ranking, from the best option to the worst one, the evaluation of each
alternative (including before and after modifications).
Las variables anteriores nos permiten evaluar la transferencia de calor entre los fluidos (Q), la
eficiencia térmica (A – Areq), la existencia o ausencia de sobredimensionamiento represento en el
porcentaje de exceso pues afecta la caída de presión; por último, la caída de presión nos ayuda a
reconocer la eficiencia hidrodinámica del intercambiador.
Conjugando las consideraciones anteriores podemos clasificar de la siguiente manera: