Heat transfer and pressure drop across tube bundles (ht.conv_tube_bank)

ht.conv_tube_bank.ESDU_tube_angle_correction(angle)[source]

Calculates the tube bank inclination correction factor according to [1] for heat transfer across a tube bundle.

F3=NuθNuθ=90=(sin(θ))0.6F_3 = \frac{Nu_{\theta}}{Nu_{\theta=90^{\circ}}} = (\sin(\theta))^{0.6}
Parameters
anglefloat

The angle of inclination of the tuba bank with respect to the longitudinal axis (90° for a straight tube bank)

Returns
F3float

ESDU tube inclination correction factor, [-]

Notes

A curve is given in [1] but it is so close the function, it is likely the function is all that is used. [1] claims this correction is valid for 100<Re<106100 < Re < 10^{6}.

For angles less than 10°, the problem should be considered internal flow, not flow across a tube bank.

References

1(1,2,3)

“Convective Heat Transfer During Crossflow of Fluids Over Plain Tube Banks.” ESDU 73031 (November 1, 1973).

Examples

>>> ESDU_tube_angle_correction(75)
0.9794139080247666
ht.conv_tube_bank.ESDU_tube_row_correction(tube_rows, staggered=True, Re=3000.0, method='Hewitt')[source]

Calculates the tube row correction factor according to [1] as shown in [2] for heat transfer across a tube bundle. This is also used for finned bundles. The correction factors are slightly different for staggered vs. inline configurations.

This method is a tabular lookup, with values of 1 when the tube row count is 10 or more.

Parameters
tube_rowsint

Number of tube rows per bundle, [-]

staggeredbool, optional

Whether in the in-line or staggered configuration, [-]

Refloat, optional

The Reynolds number of flow through the tube bank using the bare tube outer diameter and the minimum flow area through the bundle, [-]

methodstr, optional

‘Hewitt’; this may have another option in the future, [-]

Returns
F2float

ESDU tube row count correction factor, [-]

Notes

In [1], for line data, there are two curves given for different Reynolds number ranges. This is not included in [2] and only an average curve is given. This is not implemented here; Re is an argument but does not impact the result of this function.

For tube counts 1-7, [3] claims the factors from [1] are on average: [0.65, 0.77, 0.84, 0.9, 0.94, 0.97, 0.99].

References

1(1,2,3)

“Convective Heat Transfer During Crossflow of Fluids Over Plain Tube Banks.” ESDU 73031 (November 1, 1973).

2(1,2)

Hewitt, G. L. Shires, T. Reg Bott G. F., George L. Shires, and T. R. Bott. Process Heat Transfer. 1st edition. Boca Raton: CRC Press, 1994.

3

Rabas, T. J., and J. Taborek. “Survey of Turbulent Forced-Convection Heat Transfer and Pressure Drop Characteristics of Low-Finned Tube Banks in Cross Flow.” Heat Transfer Engineering 8, no. 2 (January 1987): 49-62.

Examples

>>> ESDU_tube_row_correction(4, staggered=True)
0.8984
>>> ESDU_tube_row_correction(6, staggered=False)
0.9551
ht.conv_tube_bank.Nu_ESDU_73031(Re, Pr, tube_rows, pitch_parallel, pitch_normal, Pr_wall=None, angle=90.0)[source]

Calculates the Nusselt number for crossflow across a tube bank with a specified number of tube rows, at a specified Re according to [1], also shown in [2].

Nu=aRemPr0.34F1F2\text{Nu} = a \text{Re}^m\text{Pr}^{0.34}F_1 F_2

The constants a and m come from the following tables:

In-line tube banks:

Re

a

m

10-300

0.742

0.431

300-2E5

0.211

0.651

2E5-2E6

0.116

0.700

Staggered tube banks:

Re

a

m

10-300

1.309

0.360

300-2E5

0.273

0.635

2E5-2E6

0.124

0.700

Parameters
Refloat

Reynolds number with respect to average (bulk) fluid properties and tube outside diameter, [-]

Prfloat

Prandtl number with respect to average (bulk) fluid properties, [-]

tube_rowsint

Number of tube rows per bundle, [-]

pitch_parallelfloat

Distance between tube center along a line parallel to the flow; has been called longitudinal pitch, pp, s2, SL, and p2, [m]

pitch_normalfloat

Distance between tube centers in a line 90° to the line of flow; has been called the transverse pitch, pn, s1, ST, and p1, [m]

Pr_wallfloat, optional

Prandtl number at the wall temperature; provide if a correction with the defaults parameters is desired; otherwise apply the correction elsewhere, [-]

anglefloat, optional

The angle of inclination of the tuba bank with respect to the longitudinal axis (90° for a straight tube bank)

Returns
Nufloat

Nusselt number with respect to tube outside diameter, [-]

Notes

The tube-row count correction factor F2 can be disabled by setting tube_rows to 10. The property correction factor F1 can be disabled by not specifying Pr_wall. A Prandtl number exponent of 0.26 is recommended in [1] for heating and cooling for both liquids and gases.

The pitches are used to determine whhether or not to use data for staggered or inline tube banks.

The inline coefficients are valid for a normal pitch to tube diameter ratio from 1.2 to 4; and the staggered ones from 1 to 4. The overall accuracy of this method is claimed to be 15%.

References

1(1,2)

“High-Fin Staggered Tube Banks: Heat Transfer and Pressure Drop for Turbulent Single Phase Gas Flow.” ESDU 86022 (October 1, 1986).

2

Hewitt, G. L. Shires, T. Reg Bott G. F., George L. Shires, and T. R. Bott. Process Heat Transfer. 1st edition. Boca Raton: CRC Press, 1994.

Examples

>>> Nu_ESDU_73031(Re=1.32E4, Pr=0.71, tube_rows=8, pitch_parallel=.09,
... pitch_normal=.05)
98.2563319140594
ht.conv_tube_bank.Nu_Grimison_tube_bank(Re, Pr, Do, tube_rows, pitch_parallel, pitch_normal)[source]

Calculates Nusselt number for crossflow across a tube bank of tube rows at a specified Re, Pr, and D using the Grimison methodology as described in [1].

NuDˉ=1.13C1ReD,maxmPr1/3C2\bar{Nu_D} = 1.13C_1Re_{D,max}^m Pr^{1/3}C_2
Parameters
Refloat

Reynolds number with respect to average (bulk) fluid properties and tube outside diameter, [-]

Prfloat

Prandtl number with respect to average (bulk) fluid properties, [-]

Dofloat

Tube outer diameter, [m]

tube_rowsint

Number of tube rows per bundle, [-]

pitch_parallelfloat

Distance between tube center along a line parallel to the flow; has been called longitudinal pitch, pp, s2, SL, and p2, [m]

pitch_normalfloat

Distance between tube centers in a line 90° to the line of flow; has been called the transverse pitch, pn, s1, ST, and p1, [m]

Returns
Nufloat

Nusselt number with respect to tube outside diameter, [-]

Notes

Tube row correction factors are applied for tube row counts less than 10, also published in [1].

References

1(1,2)

Grimson, E. D. (1937) Correlation and Utilisation of New Data on Flow Resistance and Heat Transfer for Cross Flow of Gases over Tube Banks. Trans. ASME. 59 583-594

Examples

>>> Nu_Grimison_tube_bank(Re=10263.37, Pr=.708, tube_rows=11,
... pitch_normal=.05, pitch_parallel=.05, Do=.025)
79.07883866010
>>> Nu_Grimison_tube_bank(Re=10263.37, Pr=.708, tube_rows=11,
... pitch_normal=.07, pitch_parallel=.05, Do=.025)
79.92721078571
ht.conv_tube_bank.Nu_HEDH_tube_bank(Re, Pr, Do, tube_rows, pitch_parallel, pitch_normal)[source]

Calculates Nusselt number for crossflow across a tube bank of tube rows at a specified Re, Pr, and D using the Heat Exchanger Design Handbook (HEDH) methodology, presented in [1].

Nu=NumfNNu = Nu_m f_N
Num=0.3+Num,lam2+Num,turb2Nu_m = 0.3 + \sqrt{Nu_{m,lam}^2 + Nu_{m,turb}^2}
Num,turb=0.037Re0.8Pr1+2.443Re0.1(Pr2/31)Nu_{m,turb} = \frac{0.037Re^{0.8} Pr}{1 + 2.443Re^{-0.1}(Pr^{2/3} -1)}
Num,lam=0.664Re0.5Pr1/3Nu_{m,lam} = 0.664Re^{0.5} Pr^{1/3}
ψ=1π4a if b >= 1\psi = 1 - \frac{\pi}{4a} \text{ if b >= 1}
ψ=1π4abif b < 1\psi = 1 - \frac{\pi}{4ab} \text{if b < 1}
fA=1+0.7ψ1.5b/a0.3(b/a)+0.7)2if inlinef_A = 1 + \frac{0.7}{\psi^{1.5}}\frac{b/a-0.3}{(b/a) + 0.7)^2} \text{if inline}
fA=1+23belif partly staggeredf_A = 1 + \frac{2}{3b} \text{elif partly staggered}
fN=1+(n1)fAnf_N = \frac{1 + (n-1)f_A}{n}
Parameters
Refloat

Reynolds number with respect to average (bulk) fluid properties and tube outside diameter, [-]

Prfloat

Prandtl number with respect to average (bulk) fluid properties, [-]

Dofloat

Tube outer diameter, [m]

tube_rowsint

Number of tube rows per bundle, [-]

pitch_parallelfloat

Distance between tube center along a line parallel to the flow; has been called longitudinal pitch, pp, s2, SL, and p2, [m]

pitch_normalfloat

Distance between tube centers in a line 90° to the line of flow; has been called the transverse pitch, pn, s1, ST, and p1, [m]

Returns
Nufloat

Nusselt number with respect to tube outside diameter, [-]

Notes

Prandtl number correction left to an outside function, although a set of coefficients were specified in [1] because they depent on whether heating or cooling is happening, and for gases, use a temperature ratio instaed of Prandtl number.

The claimed range of validity of these expressions is 10<Re<1E510 < Re < 1E5 and 0.6<Pr<10000.6 < Pr < 1000.

References

1(1,2)

Schlunder, Ernst U, and International Center for Heat and Mass Transfer. Heat Exchanger Design Handbook. Washington: Hemisphere Pub. Corp., 1987.

2

Baehr, Hans Dieter, and Karl Stephan. Heat and Mass Transfer. Springer, 2013.

Examples

>>> Nu_HEDH_tube_bank(Re=1E4, Pr=7., tube_rows=10, pitch_normal=.05,
... pitch_parallel=.05, Do=.03)
382.4636554404698

Example 3.11 in [2]:

>>> Nu_HEDH_tube_bank(Re=10263.37, Pr=.708, tube_rows=11, pitch_normal=.05,
... pitch_parallel=.05, Do=.025)
149.18735251017594
ht.conv_tube_bank.Nu_Zukauskas_Bejan(Re, Pr, tube_rows, pitch_parallel, pitch_normal, Pr_wall=None)[source]

Calculates Nusselt number for crossflow across a tube bank of tube number n at a specified Re according to the method of Zukauskas [1]. A fit to graphs from [1] published in [2] is used for the correlation. The tube row correction factor is obtained from digitized graphs from [1], and a lookup table was created and is used for speed.

The formulas are as follows:

Aligned tube banks:

NˉuD=0.9CnReD0.4Pr0.36(PrPrw)0.25 for 1<Re<100\bar Nu_D = 0.9 C_nRe_D^{0.4}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \text{ for } 1 < Re < 100
NˉuD=0.52CnReD0.5Pr0.36(PrPrw)0.25 for 100<Re<1000\bar Nu_D = 0.52 C_nRe_D^{0.5}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \text{ for } 100 < Re < 1000
NˉuD=0.27CnReD0.63Pr0.36(PrPrw)0.25 for 1000<Re<20000\bar Nu_D = 0.27 C_nRe_D^{0.63}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \text{ for } 1000 < Re < 20000
NˉuD=0.033CnReD0.8Pr0.36(PrPrw)0.25 for 20000<Re<200000\bar Nu_D = 0.033 C_nRe_D^{0.8}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \text{ for } 20000 < Re < 200000

Staggered tube banks:

NˉuD=1.04CnReD0.4Pr0.36(PrPrw)0.25 for 1<Re<500\bar Nu_D = 1.04C_nRe_D^{0.4}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \text{ for } 1 < Re < 500
NˉuD=0.71CnReD0.5Pr0.36(PrPrw)0.25 for 500<Re<1000\bar Nu_D = 0.71C_nRe_D^{0.5}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \text{ for } 500 < Re < 1000
NˉuD=0.35CnReD0.6Pr0.36(PrPrw)0.25(XtXl)0.2 for 1000<Re<20000\bar Nu_D = 0.35 C_nRe_D^{0.6}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \left(\frac{X_t}{X_l}\right)^{0.2} \text{ for } 1000 < Re < 20000
NˉuD=0.031CnReD0.8Pr0.36(PrPrw)0.25(XtXl)0.2 for 20000<Re<200000\bar Nu_D = 0.031 C_nRe_D^{0.8}Pr^{0.36}\left(\frac{Pr}{Pr_w}\right)^{0.25} \left(\frac{X_t}{X_l}\right)^{0.2} \text{ for } 20000 < Re < 200000
Parameters
Refloat

Reynolds number with respect to average (bulk) fluid properties and tube outside diameter, [-]

Prfloat

Prandtl number with respect to average (bulk) fluid properties, [-]

tube_rowsint

Number of tube rows per bundle, [-]

pitch_parallelfloat

Distance between tube center along a line parallel to the flow; has been called longitudinal pitch, pp, s2, SL, and p2, [m]

pitch_normalfloat

Distance between tube centers in a line 90° to the line of flow; has been called the transverse pitch, pn, s1, ST, and p1, [m]

Pr_wallfloat, optional

Prandtl number at the wall temperature; provide if a correction with the defaults parameters is desired; otherwise apply the correction elsewhere, [-]

Returns
Nufloat

Nusselt number with respect to tube outside diameter, [-]

Notes

If Pr_wall is not provided, the Prandtl number correction is not used and left to an outside function. A Prandtl number exponent of 0.25 is recommended in [1] for heating and cooling for both liquids and gases.

References

1(1,2,3,4)

Zukauskas, A. Heat transfer from tubes in crossflow. In T.F. Irvine, Jr. and J. P. Hartnett, editors, Advances in Heat Transfer, volume 8, pages 93-160. Academic Press, Inc., New York, 1972.

2

Bejan, Adrian. “Convection Heat Transfer”, 4E. Hoboken, New Jersey: Wiley, 2013.

Examples

>>> Nu_Zukauskas_Bejan(Re=1E4, Pr=7., tube_rows=10, pitch_parallel=.05, pitch_normal=.05)
175.9202277145248
ht.conv_tube_bank.Zukauskas_tube_row_correction(tube_rows, staggered=True, Re=10000.0)[source]

Calculates the tube row correction factor according to a graph digitized from [1] for heat transfer across a tube bundle. The correction factors are slightly different for staggered vs. inline configurations; for the staggered configuration, factors are available separately for Re larger or smaller than 1000.

This method is a tabular lookup, with values of 1 when the tube row count is 20 or more.

Parameters
tube_rowsint

Number of tube rows per bundle, [-]

staggeredbool, optional

Whether in the in-line or staggered configuration, [-]

Refloat, optional

The Reynolds number of flow through the tube bank using the bare tube outer diameter and the minimum flow area through the bundle, [-]

Returns
Ffloat

Tube row count correction factor, [-]

Notes

The basis for this method is that an infinitely long tube bank has a factor of 1; in practice the factor is reached at 20 rows.

References

1

Zukauskas, A. Heat transfer from tubes in crossflow. In T.F. Irvine, Jr. and J. P. Hartnett, editors, Advances in Heat Transfer, volume 8, pages 93-160. Academic Press, Inc., New York, 1972.

Examples

>>> Zukauskas_tube_row_correction(4, staggered=True)
0.8942
>>> Zukauskas_tube_row_correction(6, staggered=False)
0.9465
ht.conv_tube_bank.baffle_correction_Bell(crossflow_tube_fraction, method='spline')[source]

Calculate the baffle correction factor Jc which accounts for the fact that all tubes are not in crossflow to the fluid - some have fluid flowing parallel to them because they are situated in the “window”, where the baffle is cut, instead of between the tips of adjacent baffles.

Equal to 1 for no tubes in the window, increases to 1.15 when the windows are small and velocity there is high; decreases to about 0.52 for very large baffle cuts. Well designed exchangers should typically have a value near 1.0.

Cubic spline interpolation is the default method of retrieving a value from the graph, which was digitized with Engauge-Digitizer.

The interpolation can be slightly slow, so a Chebyshev polynomial was fit to a maximum error of 0.142%, average error 0.04% - well within the margin of error of the digitization of the graph; this is approximately 10 times faster, accessible via the ‘chebyshev’ method.

The Heat Exchanger Design Handbook [4], [5] provides the linear curve fit, which covers the “practical” range of baffle cuts 15-45% but not the last dip in the graph. This method is not recommended, but can be used via the method “HEDH”.

Jc=0.55+0.72FcJ_c = 0.55 + 0.72Fc
Parameters
crossflow_tube_fractionfloat

Fraction of tubes which are between baffle tips and not in the window, [-]

methodstr, optional

One of ‘chebyshev’, ‘spline’, or ‘HEDH’

Returns
Jcfloat

Baffle correction factor in the Bell-Delaware method, [-]

Notes

max: ~1.1536 at ~0.9066 min: ~0.5328 at 0 value at 1: ~1.0314

For the ‘spline’ method, this function takes ~13 us per call. The other two methods are approximately 10x faster.

References

1

Bell, Kenneth J. Final Report of the Cooperative Research Program on Shell and Tube Heat Exchangers. University of Delaware, Engineering Experimental Station, 1963.

2

Bell, Kenneth J. Delaware Method for Shell-Side Design. In Heat Transfer Equipment Design, by Shah, R. K., Eleswarapu Chinna Subbarao, and R. A. Mashelkar. CRC Press, 1988.

3

Green, Don, and Robert Perry. Perry’s Chemical Engineers’ Handbook, Eighth Edition. McGraw-Hill Professional, 2007.

4

Schlünder, Ernst U, and International Center for Heat and Mass Transfer. Heat Exchanger Design Handbook. Washington: Hemisphere Pub. Corp., 1987.

5

Serth, R. W., Process Heat Transfer: Principles, Applications and Rules of Thumb. 2E. Amsterdam: Academic Press, 2014.

Examples

For a HX with four groups of tube bundles; the top and bottom being 9 tubes each, in the window, and the two middle bundles having 41 tubes each, for a total of 100 tubes, the fraction between baffle tubes and not in the window is 0.82. The correction factor is then:

>>> baffle_correction_Bell(0.82)
1.1258554691854046
ht.conv_tube_bank.baffle_leakage_Bell(Ssb, Stb, Sm, method='spline')[source]

Calculate the baffle leakage factor Jl which accounts for leakage between each baffle. Cubic spline interpolation is the default method of retrieving a value from the graph, which was digitized with Engauge-Digitizer.

The Heat Exchanger Design Handbook [4], [5] provides a curve fit as well. This method is not recommended, but can be used via the method “HEDH”.

JL=0.44(1rs)+[10.44(1rs)]exp(2.2rlm)J_L = 0.44(1-r_s) + [1 - 0.44(1-r_s)]\exp(-2.2r_{lm})
rs=SsbSsb+Stbr_s = \frac{S_{sb}}{S_{sb} + S_{tb}}
rlm=Ssb+StbSmr_{lm} = \frac{S_{sb} + S_{tb}}{S_m}
Parameters
Ssbfloat

Shell to baffle leakage area, [m^2]

Stbfloat

Total baffle leakage area, [m^2]

Smfloat

Crossflow area, [m^2]

methodstr, optional

One of ‘spline’, or ‘HEDH’

Returns
Jlfloat

Baffle leakage factor in the Bell-Delaware method, [-]

Notes

Takes ~5 us per call. If the x parameter is larger than 0.743614, it is clipped to it.

The HEDH curve fits are rather poor and only 6x faster to evaluate. The HEDH example in [6]’s spreadsheet has an error and uses 0.044 instead of 0.44 in the equation.

References

1

Bell, Kenneth J. Final Report of the Cooperative Research Program on Shell and Tube Heat Exchangers. University of Delaware, Engineering Experimental Station, 1963.

2

Bell, Kenneth J. Delaware Method for Shell-Side Design. In Heat Transfer Equipment Design, by Shah, R. K., Eleswarapu Chinna Subbarao, and R. A. Mashelkar. CRC Press, 1988.

3

Green, Don, and Robert Perry. Perry’s Chemical Engineers’ Handbook, Eighth Edition. McGraw-Hill Professional, 2007.

4

Schlünder, Ernst U, and International Center for Heat and Mass Transfer. Heat Exchanger Design Handbook. Washington: Hemisphere Pub. Corp., 1987.

5

Serth, R. W., Process Heat Transfer: Principles, Applications and Rules of Thumb. 2E. Amsterdam: Academic Press, 2014.

6

Hall, Stephen. Rules of Thumb for Chemical Engineers, Fifth Edition. 5th edition. Oxford ; Waltham , MA: Butterworth-Heinemann, 2012.

Examples

>>> baffle_leakage_Bell(1, 3, 8)
0.5906621282470
>>> baffle_leakage_Bell(1, 3, 8, 'HEDH')
0.5530236260777
ht.conv_tube_bank.bundle_bypassing_Bell(bypass_area_fraction, seal_strips, crossflow_rows, laminar=False, method='spline')[source]

Calculate the bundle bypassing effect Jb according to the Bell-Delaware method for heat exchanger design. Cubic spline interpolation is the default method of retrieving a value from the graph, which was digitized with Engauge-Digitizer.

The Heat Exchanger Design Handbook [4] provides a curve fit as well. This method is not recommended, but can be used via the method “HEDH”:

Jb=exp[1.25Fsbp(12rss1/3)]J_b = \exp\left[-1.25 F_{sbp} (1 - {2r_{ss}}^{1/3} )\right]

For laminar flows, replace 1.25 with 1.35.

Parameters
bypass_area_fractionfloat

Fraction of the crossflow area which is not blocked by a baffle or anything else and available for bypassing, [-]

seal_stripsint

Number of seal strips per side of a baffle added to prevent bypassing, [-]

crossflow_rowsint

The number of tube rows in the crosslfow of the baffle, [-]

laminarbool

Whether to use the turbulent correction values or the laminar ones; the Bell-Delaware method uses a Re criteria of 100 for this, [-]

methodstr, optional

One of ‘spline’, or ‘HEDH’

Returns
Jbfloat

Bundle bypassing effect correction factor in the Bell-Delaware method, [-]

Notes

Takes ~5 us per call. If the bypass_area_fraction parameter is larger than 0.695, it is clipped to it.

References

1

Bell, Kenneth J. Final Report of the Cooperative Research Program on Shell and Tube Heat Exchangers. University of Delaware, Engineering Experimental Station, 1963.

2

Bell, Kenneth J. Delaware Method for Shell-Side Design. In Heat Transfer Equipment Design, by Shah, R. K., Eleswarapu Chinna Subbarao, and R. A. Mashelkar. CRC Press, 1988.

3

Green, Don, and Robert Perry. Perry’s Chemical Engineers’ Handbook, Eighth Edition. McGraw-Hill Professional, 2007.

4

Schlünder, Ernst U, and International Center for Heat and Mass Transfer. Heat Exchanger Design Handbook. Washington: Hemisphere Pub. Corp., 1987.

Examples

>>> bundle_bypassing_Bell(0.5, 5, 25)
0.8469611760884599
>>> bundle_bypassing_Bell(0.5, 5, 25, method='HEDH')
0.8483210970579099
ht.conv_tube_bank.dP_Kern(m, rho, mu, DShell, LSpacing, pitch, Do, NBaffles, mu_w=None)[source]

Calculates pressure drop for crossflow across a tube bank according to the equivalent-diameter method developed by Kern [1], presented in [2].

ΔP=f(m/Ss)2Ds(NB+1)2ρDe(μ/μw)0.14\Delta P = \frac{f (m/S_s)^2 D_s(N_B+1)}{2\rho D_e(\mu/\mu_w)^{0.14}}
SS=DS(PTDo)LBPTS_S = \frac{D_S (P_T-D_o) L_B}{P_T}
De=4(PT2πDo2/4)πDoD_e = \frac{4(P_T^2 - \pi D_o^2/4)}{\pi D_o}
Parameters
mfloat

Mass flow rate, [kg/s]

rhofloat

Fluid density, [kg/m^3]

mufloat

Fluid viscosity, [Pa*s]

DShellfloat

Diameter of exchanger shell, [m]

LSpacingfloat

Baffle spacing, [m]

pitchfloat

Tube pitch, [m]

Dofloat

Tube outer diameter, [m]

NBafflesfloat

Baffle count, []

mu_wfloat

Fluid viscosity at wall temperature, [Pa*s]

Returns
dPfloat

Pressure drop across bundle, [Pa]

Notes

Adjustment for viscosity left out of this function. Example is from [2]. Roughly 10% difference due to reading of graph. Graph scanned from [1], and interpolation is used to read it.

References

1(1,2)

Kern, Donald Quentin. Process Heat Transfer. McGraw-Hill, 1950.

2(1,2)

Peters, Max, Klaus Timmerhaus, and Ronald West. Plant Design and Economics for Chemical Engineers. 5E. New York: McGraw-Hill, 2002.

Examples

>>> dP_Kern(m=11., rho=995., mu=0.000803, mu_w=0.000657, DShell=0.584,
... LSpacing=0.1524, pitch=0.0254, Do=.019, NBaffles=22)
18980.58768759033
ht.conv_tube_bank.dP_Zukauskas(Re, n, ST, SL, D, rho, Vmax)[source]

Calculates pressure drop for crossflow across a tube bank of tube number n at a specified Re. Method presented in [1]. Also presented in [2].

ΔP=NLχ(ρVmax22)f\Delta P = N_L \chi \left(\frac{\rho V_{max}^2}{2}\right)f
Parameters
Refloat

Reynolds number, [-]

nfloat

Number of tube rows, [-]

STfloat

Transverse pitch, used only by some conditions, [m]

SLfloat

Longitudal pitch, used only by some conditions, [m]

Dfloat

Tube outer diameter, [m]

rhofloat

Fluid density, [kg/m^3]

Vmaxfloat

Maximum velocity, [m/s]

Returns
dPfloat

Pressure drop, [Pa]

Notes

Does not account for effects in a heat exchanger. Example 2 is from [2]. Matches to 0.3%; figures are very approximate. Interpolation used with 4 graphs to obtain friction factor and a correction factor.

References

1

Zukauskas, A. Heat transfer from tubes in crossflow. In T.F. Irvine, Jr. and J. P. Hartnett, editors, Advances in Heat Transfer, volume 8, pages 93-160. Academic Press, Inc., New York, 1972.

2(1,2)

Bergman, Theodore L., Adrienne S. Lavine, Frank P. Incropera, and David P. DeWitt. Introduction to Heat Transfer. 6E. Hoboken, NJ: Wiley, 2011.

Examples

>>> dP_Zukauskas(Re=13943., n=7, ST=0.0313, SL=0.0343, D=0.0164, rho=1.217, Vmax=12.6)
235.22916169
>>> dP_Zukauskas(Re=13943., n=7, ST=0.0313, SL=0.0313, D=0.0164, rho=1.217, Vmax=12.6)
217.0750033
ht.conv_tube_bank.laminar_correction_Bell(Re, total_row_passes)[source]

Calculate the correction factor for adverse temperature gradient built up in laminar flow Jr.

This correction begins at Re = 100, and is interpolated between the value of the formula until Re = 20, when it is the value of the formula. It is 1 for Re >= 100. The value of the formula is not allowed to be less than 0.4.

Jr=(10Nrow,passes,tot)0.18Jr^* = \left(\frac{10}{N_{row,passes,tot}}\right)^{0.18}
Parameters
Refloat

Shell Reynolds number in the Bell-Delaware method, [-]

total_row_passesint

The total number of rows passed by the fluid, including those in windows and counting repeat passes of tube rows, [-]

Returns
Jrfloat

Correction factor for adverse temperature gradient built up in laminar flow, [-]

Notes

[5] incorrectly uses the number of tube rows per crosslfow section, not total.

References

1

Bell, Kenneth J. Final Report of the Cooperative Research Program on Shell and Tube Heat Exchangers. University of Delaware, Engineering Experimental Station, 1963.

2

Bell, Kenneth J. Delaware Method for Shell-Side Design. In Heat Transfer Equipment Design, by Shah, R. K., Eleswarapu Chinna Subbarao, and R. A. Mashelkar. CRC Press, 1988.

3

Schlünder, Ernst U, and International Center for Heat and Mass Transfer. Heat Exchanger Design Handbook. Washington: Hemisphere Pub. Corp., 1987.

4

Serth, R. W., Process Heat Transfer: Principles, Applications and Rules of Thumb. 2E. Amsterdam: Academic Press, 2014.

5

Hall, Stephen. Rules of Thumb for Chemical Engineers, Fifth Edition. 5th edition. Oxford ; Waltham , MA: Butterworth-Heinemann, 2012.

Examples

>>> laminar_correction_Bell(30, 80)
0.7267995454361379
ht.conv_tube_bank.unequal_baffle_spacing_Bell(baffles, baffle_spacing, baffle_spacing_in=None, baffle_spacing_out=None, laminar=False)[source]

Calculate the correction factor for unequal baffle spacing Js, which accounts for higher velocity of fluid flow and greater heat transfer coefficients when the in and/or out baffle spacing is less than the standard spacing.

Js=(nb1)+(Bin/B)(1nb)+(Bout/B)(1nb)(nb1)+(Bin/B)+(Bout/B)J_s = \frac{(n_b - 1) + (B_{in}/B)^{(1-n_b)} + (B_{out}/B)^{(1-n_b)}} {(n_b - 1) + (B_{in}/B) + (B_{out}/B)}
Parameters
bafflesint

Number of baffles, [-]

baffle_spacingfloat

Average spacing between one end of one baffle to the start of the next baffle for non-exit baffles, [m]

baffle_spacing_infloat, optional

Spacing between entrace to first baffle, [m]

baffle_spacing_outfloat, optional

Spacing between last baffle and exit, [m]

laminarbool, optional

Whether to use the turbulent exponent or the laminar one; the Bell-Delaware method uses a Re criteria of 100 for this, [-]

Returns
Jsfloat

Unequal baffle spacing correction factor, [-]

References

1

Bell, Kenneth J. Final Report of the Cooperative Research Program on Shell and Tube Heat Exchangers. University of Delaware, Engineering Experimental Station, 1963.

2

Bell, Kenneth J. Delaware Method for Shell-Side Design. In Heat Transfer Equipment Design, by Shah, R. K., Eleswarapu Chinna Subbarao, and R. A. Mashelkar. CRC Press, 1988.

3

Schlünder, Ernst U, and International Center for Heat and Mass Transfer. Heat Exchanger Design Handbook. Washington: Hemisphere Pub. Corp., 1987.

4

Serth, R. W., Process Heat Transfer: Principles, Applications and Rules of Thumb. 2E. Amsterdam: Academic Press, 2014.

5

Hall, Stephen. Rules of Thumb for Chemical Engineers, Fifth Edition. 5th edition. Oxford ; Waltham , MA: Butterworth-Heinemann, 2012.

Examples

>>> unequal_baffle_spacing_Bell(16, .1, .15, 0.15)
0.9640087802805195