This worksheet allows calculation of the Overall Heat Transfer Coefficient and pressure drops in a TEMA type E shell
and tubes heat exchanger using Bell-Delaware method.
Data input corresponds to yellow cells and the tube array must be selected in the list box. Results are shown in cells K48 to K50
Process Data
(see drawing at right)
tube side shell side N= Number of tubes per shell
Ntw= Number of tubes in window
Flowrate (kg/s) 38.88 33.33 AUX NCL=Number of tubes in central row
Inlet Temp. (C) 40 104 44 NBT= Number of tubes crossing
Outlet Temp. (C) 60 85 45 the baffle
Inlet Pressure (bar abs) 4 3 NC= Number of tubes rows between
Q= Heat Duty (W) 1664064 1671833 consecutives baffles edges
Thermal Balance sigma=separation between rows
check OK- error <2% DS= Shell Diameter
DB=Baffle Diameter
Density (kg/m3) 716 578 BC= Baffle cut
Viscosity (cp) 0.62 0.16 A= Baffle central angle in degrees
Specific Heat. (J/kg.K) 2140 2640 Do= Tubes external diameter
Thermal Conduct. (W/m.K) 0.129 0.0917 Di= Tubes internal diameter
Visc at wall temperature (cp) 0.59 0.16 DBT=Baffle holes diameter
Pr 10.29 4.61 Pt= Tubes separation (Pitch)
c= clearance between tubes
Tube side Calculations n= number of tube passes
NB= Number of baffles
at= Flow Area 0.0401 m2 NS= Number of sealing strips pairs
Gt= mass velocity 970.2 kg/m2/s L= tubes length
Ret= Reynolds Number 24568 B= Baffles spacing
hi= internal coeff. 1327 W/m2K
hio= 1096 W/m2K
friction factor 0.00758 Internal fouling ressistance(m2.K/W)
Pressure drop 16099 N/m2 External fouling ressistance
Shell side calculations- Bell method
effective number of tubes rows 334 Select the tubes array
Total number of tubes rows correction 1
Sm=Cross flow area 0.0491 m2
Gm= Mass velocity 678.3 kg/s.m2
Rem 80544
j=Colburn Coefficient 3.76E-03
FBP=By pass fraction 0.390
alfa 1.35
by pass correction for h 0.8426 Ns/0.5*Nc
r 0.1812 0.307692
Sw 0.0356 m2
window correction factor 1.1251
STB=tube-baffle leakage area 0.0083 m2
SSB=shell baffle leakage area 0.0030 m2
SL= total leakage area 0.0112 m2 U= Overall heat transfer coefficient calculated
SL/Sm= 0.2288 Tubeside pressure drop per shell
(1-hL/hNL)o 0.2029 Shellside pressure drop per shell
1-hL/hNL 0.2567
hL/hNL 0.7433 Wall temperature calculation
ho 1714.8 W/m2K
shell side friction factor 0.0908 1/(1/hio+Rfi)=
alfa2 4 1/(1/ho+Rfo)=
by pass correction for deltap 0.602
(deltap)BP 1131.3 N/m2 Wall mean temp
vZ 1.38 m/s
Nw 5.54
(deltap)w 2921 N/m2 Number of shells
(1-deltapL/deltapNL)o 0.398 in series 1
(1-deltapL/deltapNL) 0.503
deltapL/deltapNL 0.497 Heat transfer
deltap shell 36886 N/m2 Total area 105.4
LMTD 44.5
Q 1894475
% Overdesign 14%
This worksheet allows calculation of the Overall Heat Transfer Coefficient and pressure drops in a TEMA type E shell and tubes heat exchanger using Bell-Delaware method.
Data input corresponds to yellow cells and the tube array must be selected in the list box. Results are shown in cells K48 to K50
GEOMETRIC DATA
(see drawing at right)
N= Number of tubes per shell 414
Ntw= Number of tubes in window 75
NCL=Number of tubes in central row 21
NBT= Number of tubes crossing
339
NC= Number of tubes rows between
consecutives baffles edges 13
sigma=separation between rows 0.022 m
0.609 m
0.6043 m
0.15225 m
A= Baffle central angle in degrees 120
Do= Tubes external diameter 0.019 m
de tubos 0.0157 m
DBT=Baffle holes diameter 0.0198 m
Pt= Tubes separation (Pitch) 0.0254 m
c= clearance between tubes 0.0064
n= number of tube passes 2
17
NS= Number of sealing strips pairs 2
4.267 m
0.234 m
Internal fouling ressistance(m2.K/W) 0.0005
External fouling ressistance 0.0004
U= Overall heat transfer coefficient calculated 417.5 W/m2K
Tubeside pressure drop per shell 16099 N/m2
Shellside pressure drop per shell 36886 N/m2
Wall temperature calculation
708 W/m2K Ts=shell side fluid mean temp. 94.5 C
1017 W/m2K Tt=Tubeside fluid mean temp. 50 C
76.24 C
DPt 16099 Pa
0.217 DPc 36886 Pa
m2
C Ft 0.97 MTD 43.0 C
W
LMTD Correction factor Ft Ft
T1 104
T2 85.0
t1 40
t2 60.0
R 1.0
S 0.31
Number of shells Numero de carcasas 1
Tube passes 2
FT 0.967152
Rem<100
j=a*Rem^b
a
array
1 triangle (30) Pt/D= 1.25 1.81
2 Square Pt/D=1.25 0.97
3 Rotated Square Pt/D=1.25 1.81
4 triangle (30) Pt/D= 1.5 1.34
5 Square Pt/D=1.5 0.88
3 Selected array 1.81
For Rem = 80544.11 j =
Coefficients for f calculation - ln(f)=a+blnRem+c*(lnRem)^2+d*(lnRem)^3
a
arreglo
1 tringulo (60) Pt/D= 1.25 4.376
2 Cuadrado Pt/D=1.25 3.923
3 Cuadro rotado Pt/D=1.25 3.923
4 tringulo (60) Pt/D= 1.5 3.196
5 Cuadrado Pt/D=1.5 3.04
3 Arreglo seleccionado 3.923
Coefficients for Rem 80544.11
Calculation of friction factor and J factor
coefficients for J calculation
100<Rem<3000 Rem>3000
j=a*Rem^b lnj=a+b*lnRem+c(lnRem)^2 j=aRem^b
b a b c a b
-0.72 1.7 -1.25 0.065 0.275 -0.38
-0.62 3.29 -1.9 0.121 0.275 -0.38
-0.72 0.32 -0.76 0.025 0.275 -0.38
-0.68 0.37 -0.56 0.0087 0.275 -0.38
-0.64 1.61 -1.34 0.082 0.275 -0.38
-0.72 0.32 -0.76 0.025 0.275 -0.38
3.76E-03
Coefficients for f calculation - ln(f)=a+blnRem+c*(lnRem)^2+d*(lnRem)^3
Rem<100 100<Rem<1300 1300<Rem<3000
b c d a b c d a b
-1.025 0 0 5.293 -1.864 0.1584 -0.00472 5.293 -1.864
-0.984 0 0 3.871 -0.498 -0.2052 0.0221 3.871 -0.498
-0.984 0 0 6.3 -2.403 0.226 -0.00721 6.3 -2.403
-1.176 0.0557 0 3.472 -1.547 0.1425 -0.00454 3.472 -1.547
-1.13 0.031 0 -8.015 5.152 -1.166 0.0771 -6.666 1.152
-0.984 0 0 6.3 -2.403 0.226 -0.00721 6.3 -2.403
a = 6.3
b = -2.403 f = 0.090804
c = 0.226
d = -0.00721
1300<Rem<3000 Rem>3000
c d a b c d
0.1584 -0.00472 5.293 -1.864 0.1584 -0.00472
-0.2052 0.00221 -7.907 1.774 -0.16 0.00407
0.226 -0.00721 6.3 -2.403 0.226 -0.00721
0.1425 -0.00454 3.472 -1.547 0.1425 -0.00454
-0.0981 0.00237 -6.666 1.152 -0.0981 0.00237
0.226 -0.00721 6.3 -2.403 0.226 -0.00721
