Convection to jacketed vessels (ht.conv_jacket)

ht.conv_jacket.Lehrer(m, Dtank, Djacket, H, Dinlet, rho, Cp, k, mu, muw=None, isobaric_expansion=None, dT=None, inlettype='tangential', inletlocation='auto')[source]

Calculates average heat transfer coefficient for a jacket around a vessel according to [1] as described in [2].

NuS,L=[0.03ReS0.75Pr1+1.74(Pr1)ReS0.125](μμw)0.14Nu_{S,L} = \left[\frac{0.03Re_S^{0.75}Pr}{1 + \frac{1.74(Pr-1)} {Re_S^{0.125}}}\right]\left(\frac{\mu}{\mu_w}\right)^{0.14}
dg=(83)0.5δd_g = \left(\frac{8}{3}\right)^{0.5}\delta
vh=(vSvinlet)0.5+vAv_h = (v_Sv_{inlet})^{0.5} + v_A
vinlet=Qπ4dinlet2v_{inlet} = \frac{Q}{\frac{\pi}{4}d_{inlet}^2}
vs=Qπ4(Djacket2Dtank2)v_s = \frac{Q}{\frac{\pi}{4}(D_{jacket}^2 - D_{tank}^2)}

For Radial inlets:

vA=0.5(2gHβδΔT)0.5v_A = 0.5(2g H \beta\delta \Delta T)^{0.5}

For Tangential inlets:

vA=0v_A = 0
Parameters
mfloat

Mass flow rate of fluid, [kg/s]

Dtankfloat

Outer diameter of tank or vessel surrounded by jacket, [m]

Djacketfloat

Inner diameter of jacket surrounding a vessel or tank, [m]

Hfloat

Height of the vessel or tank, [m]

Dinletfloat

Inner diameter of inlet into the jacket, [m]

rhofloat

Density of the fluid at Tm [kg/m^3]

Cpfloat

Heat capacity of fluid at Tm [J/kg/K]

kfloat

Thermal conductivity of fluid at Tm [W/m/K]

mufloat

Viscosity of fluid at Tm [Pa*s]

muwfloat, optional

Viscosity of fluid at Tw [Pa*s]

isobaric_expansionfloat, optional

Constant pressure expansivity of a fluid, [m^3/mol/K]

dTfloat, optional

Temperature difference of fluid in jacket, [K]

inlettypestr, optional

Either ‘tangential’ or ‘radial’

inletlocationstr, optional

Either ‘top’ or ‘bottom’ or ‘auto’

Returns
hfloat

Average heat transfer coefficient inside the jacket [W/m^2/K]

Notes

If the fluid is heated and enters from the bottom, natural convection assists the heat transfer and the Grashof term is added; if it were to enter from the top, it would be subtracted. The situation is reversed if entry is from the top.

References

1

Lehrer, Isaac H. “Jacket-Side Nusselt Number.” Industrial & Engineering Chemistry Process Design and Development 9, no. 4 (October 1, 1970): 553-58. doi:10.1021/i260036a010.

2(1,2,3)

Gesellschaft, V. D. I., ed. VDI Heat Atlas. 2nd edition. Berlin; New York:: Springer, 2010.

Examples

Example as in [2], matches completely.

>>> Lehrer(m=2.5, Dtank=0.6, Djacket=0.65, H=0.6, Dinlet=0.025, dT=20.,
... rho=995.7, Cp=4178.1, k=0.615, mu=798E-6, muw=355E-6)
2922.128124761829

Examples similar to in [2] but covering the other case:

>>> Lehrer(m=2.5, Dtank=0.6, Djacket=0.65, H=0.6, Dinlet=0.025, dT=20.,
... rho=995.7, Cp=4178.1, k=0.615, mu=798E-6, muw=355E-6,
... inlettype='radial', isobaric_expansion=0.000303)
3269.4389632666557
ht.conv_jacket.Stein_Schmidt(m, Dtank, Djacket, H, Dinlet, rho, Cp, k, mu, muw=None, rhow=None, inlettype='tangential', inletlocation='auto', roughness=0.0)[source]

Calculates average heat transfer coefficient for a jacket around a vessel according to [1] as described in [2].

lch=[(π2)2Dtank2+H2]0.5l_{ch} = \left[\left(\frac{\pi}{2}\right)^2 D_{tank}^2+H^2\right]^{0.5}
dch=2δd_{ch} = 2\delta
Rej=vchdchρμRe_j = \frac{v_{ch}d_{ch}\rho}{\mu}
GrJ=gρ(ρρw)dch3μ2Gr_J = \frac{g\rho(\rho-\rho_w)d_{ch}^3}{\mu^2}
ReJ,eq=[ReJ2±(GrJHdch50)]0.5Re_{J,eq} = \left[Re_J^2\pm \left(\frac{|Gr_J|\frac{H}{d_{ch}}}{50} \right)\right]^{0.5}
NuJ=(NuA3+NuB3+NuC3+NuD3)1/3(μμw)0.14Nu_J = (Nu_A^3 + Nu_B^3 + Nu_C^3 + Nu_D^3)^{1/3}\left(\frac{\mu} {\mu_w}\right)^{0.14}
NuJ=hdchkNu_J = \frac{h d_{ch}}{k}
NuA=3.66Nu_A = 3.66
NuB=1.62Pr1/3ReJ,eq1/3(dchlch)1/3Nu_B = 1.62 Pr^{1/3}Re_{J,eq}^{1/3}\left(\frac{d_{ch}}{l_{ch}} \right)^{1/3}
NuC=0.664Pr1/3(ReJ,eqdchlch)0.5Nu_C = 0.664Pr^{1/3}(Re_{J,eq}\frac{d_{ch}}{l_{ch}})^{0.5}
if ReJ,eq<2300:NuD=0\text{if } Re_{J,eq} < 2300: Nu_D = 0
NuD=0.0115Pr1/3ReJ,eq0.9(1(2300ReJ,eq)2.5)(1+(dchlch)2/3)Nu_D = 0.0115Pr^{1/3}Re_{J,eq}^{0.9}\left(1 - \left(\frac{2300} {Re_{J,eq}}\right)^{2.5}\right)\left(1 + \left(\frac{d_{ch}}{l_{ch}} \right)^{2/3}\right)

For Radial inlets:

vch=vMit(lnbMitbEin1bEinbMit)v_{ch} = v_{Mit}\left(\frac{\ln\frac{b_{Mit}}{b_{Ein}}}{1 - \frac{b_{Ein}}{b_{Mit}}}\right)
bEin=π8Dinlet2δb_{Ein} = \frac{\pi}{8}\frac{D_{inlet}^2}{\delta}
bMit=π2Dtank1+π24Dtank2H2b_{Mit} = \frac{\pi}{2}D_{tank}\sqrt{1 + \frac{\pi^2}{4}\frac {D_{tank}^2}{H^2}}
vMit=Q2δbMitv_{Mit} = \frac{Q}{2\delta b_{Mit}}

For Tangential inlets:

vch=(vx2+vz2)0.5v_{ch} = (v_x^2 + v_z^2)^{0.5}
vx=vinlet(ln[1+fdDtankHDinlet2vx(0)vinlet]fdDtankHDinlet2)v_x = v_{inlet}\left(\frac{\ln[1 + \frac{f_d D_{tank}H}{D_{inlet}^2} \frac{v_x(0)}{v_{inlet}}]}{\frac{f_d D_{tank}H}{D_{inlet}^2}}\right)
vx(0)=K3+(K32+K4)0.5v_x(0) = K_3 + (K_3^2 + K_4)^{0.5}
K3=vinlet4Dinlet2vinlet4fdDtankHK_3 = \frac{v_{inlet}}{4} -\frac{D_{inlet}^2v_{inlet}}{4f_d D_{tank}H}
K4=Dinlet2vinlet22fdDtankHK_4 = \frac{D_{inlet}^2v_{inlet}^2}{2f_d D_{tank} H}
vz=QπDtankδv_z = \frac{Q}{\pi D_{tank}\delta}
vinlet=Qπ4Dinlet2v_{inlet} = \frac{Q}{\frac{\pi}{4}D_{inlet}^2}
Parameters
mfloat

Mass flow rate of fluid, [kg/m^3]

Dtankfloat

Outer diameter of tank or vessel surrounded by jacket, [m]

Djacketfloat

Inner diameter of jacket surrounding a vessel or tank, [m]

Hfloat

Height of the vessel or tank, [m]

Dinletfloat

Inner diameter of inlet into the jacket, [m]

rhofloat

Density of the fluid at Tm [kg/m^3]

Cpfloat

Heat capacity of fluid at Tm [J/kg/K]

kfloat

Thermal conductivity of fluid at Tm [W/m/K]

mufloat

Viscosity of fluid at Tm [Pa*s]

muwfloat, optional

Viscosity of fluid at Tw [Pa*s]

rhowfloat, optional

Density of the fluid at Tw [kg/m^3]

inlettypestr, optional

Either ‘tangential’ or ‘radial’

inletlocationstr, optional

Either ‘top’ or ‘bottom’ or ‘auto’

roughnessfloat, optional

Roughness of the tank walls [m]

Returns
hfloat

Average transfer coefficient inside the jacket [W/m^2/K]

Notes

[1] is in German and has not been reviewed. Multiple other formulations are considered in [1].

If the fluid is heated and enters from the bottom, natural convection assists the heat transfer and the Grashof term is added; if it were to enter from the top, it would be subtracted. The situation is reversed if entry is from the top.

References

1(1,2,3)

Stein, Prof Dr-Ing Werner Alexander, and Dipl-Ing (FH) Wolfgang Schmidt. “Wärmeübergang auf der Wärmeträgerseite eines Rührbehälters mit einem einfachen Mantel.” Forschung im Ingenieurwesen 59, no. 5 (May 1993): 73-90. doi:10.1007/BF02561203.

2(1,2)

Gesellschaft, V. D. I., ed. VDI Heat Atlas. 2nd edition. Berlin; New York:: Springer, 2010.

Examples

Example as in [2], matches in all but friction factor:

>>> Stein_Schmidt(m=2.5, Dtank=0.6, Djacket=0.65, H=0.6, Dinlet=0.025,
... rho=995.7, Cp=4178.1, k=0.615, mu=798E-6, muw=355E-6, rhow=971.8)
5695.2041698088615