Flow boiling (ht.boiling_flow)

ht.boiling_flow.Chen_Bennett(m, x, D, rhol, rhog, mul, mug, kl, Cpl, Hvap, sigma, dPsat, Te)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is developed in [1] and [2], and reviewed in [3]. This model is one of the most often used, and replaces the Chen_Edelstein correlation. It uses the Dittus-Boelter correlation for turbulent convection and the Forster-Zuber correlation for pool boiling, and combines them with two factors F and S.

htp=Shnb+Fhsp,lh_{tp} = S\cdot h_{nb} + F \cdot h_{sp,l}
hsp,l=0.023Rel0.8Prl0.4kl/Dh_{sp,l} = 0.023 Re_l^{0.8} Pr_l^{0.4} k_l/D
Rel=DG(1x)μlRe_l = \frac{DG(1-x)}{\mu_l}
hnb=0.00122(λl0.79cp,l0.45ρl0.49σ0.5μ0.29Hvap0.24ρg0.24)ΔTsat0.24Δpsat0.75h_{nb} = 0.00122\left( \frac{\lambda_l^{0.79} c_{p,l}^{0.45} \rho_l^{0.49}}{\sigma^{0.5} \mu^{0.29} H_{vap}^{0.24} \rho_g^{0.24}} \right)\Delta T_{sat}^{0.24} \Delta p_{sat}^{0.75}
F=(Pr1+12)0.444(1+Xtt0.5)1.78F = \left(\frac{Pr_1+1}{2}\right)^{0.444}\cdot (1+X_{tt}^{-0.5})^{1.78}
S=1exp(FhconvX0/kl)FhconvX0/klS = \frac{1-\exp(-F\cdot h_{conv} \cdot X_0/k_l)} {F\cdot h_{conv}\cdot X_0/k_l}
Xtt=(1xx)0.9(ρgρl)0.5(μlμg)0.1X_{tt} = \left( \frac{1-x}{x}\right)^{0.9} \left(\frac{\rho_g}{\rho_l} \right)^{0.5}\left( \frac{\mu_l}{\mu_g}\right)^{0.1}
X0=0.041(σg(ρlρv))0.5X_0 = 0.041 \left(\frac{\sigma}{g \cdot (\rho_l-\rho_v)}\right)^{0.5}
Parameters
mfloat

Mass flow rate [kg/s]

xfloat

Quality at the specific tube interval []

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

rhogfloat

Density of the gas [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

mugfloat

Viscosity of gas [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

Cplfloat

Heat capacity of liquid [J/kg/K]

Hvapfloat

Heat of vaporization of liquid [J/kg]

sigmafloat

Surface tension of liquid [N/m]

dPsatfloat

Difference in Saturation pressure of fluid at Te and T, [Pa]

Tefloat

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

See also

Chen_Edelstein
turbulent_Dittus_Boelter
Forster_Zuber

Notes

[1] and [2] have been reviewed, but the model is only put together in the review of [3]. Many other forms of this equation exist with different functions for F and S.

References

1(1,2)

Bennett, Douglas L., and John C. Chen. “Forced Convective Boiling in Vertical Tubes for Saturated Pure Components and Binary Mixtures.” AIChE Journal 26, no. 3 (May 1, 1980): 454-61. doi:10.1002/aic.690260317.

2(1,2)

Bennett, Douglas L., M.W. Davies and B.L. Hertzler, The Suppression of Saturated Nucleate Boiling by Forced Convective Flow, American Institute of Chemical Engineers Symposium Series, vol. 76, no. 199. 91-103, 1980.

3(1,2)

Bertsch, Stefan S., Eckhard A. Groll, and Suresh V. Garimella. “Review and Comparative Analysis of Studies on Saturated Flow Boiling in Small Channels.” Nanoscale and Microscale Thermophysical Engineering 12, no. 3 (September 4, 2008): 187-227. doi:10.1080/15567260802317357.

Examples

>>> Chen_Bennett(m=0.106, x=0.2, D=0.0212, rhol=567, rhog=18.09,
... mul=156E-6, mug=7.11E-6, kl=0.086, Cpl=2730, Hvap=2E5, sigma=0.02,
... dPsat=1E5, Te=3)
4938.275351219369
ht.boiling_flow.Chen_Edelstein(m, x, D, rhol, rhog, mul, mug, kl, Cpl, Hvap, sigma, dPsat, Te)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is developed in [1] and [2], and reviewed in [3]. This model is one of the most often used. It uses the Dittus-Boelter correlation for turbulent convection and the Forster-Zuber correlation for pool boiling, and combines them with two factors F and S.

htp=Shnb+Fhsp,lh_{tp} = S\cdot h_{nb} + F \cdot h_{sp,l}
hsp,l=0.023Rel0.8Prl0.4kl/Dh_{sp,l} = 0.023 Re_l^{0.8} Pr_l^{0.4} k_l/D
Rel=DG(1x)μlRe_l = \frac{DG(1-x)}{\mu_l}
hnb=0.00122(λl0.79cp,l0.45ρl0.49σ0.5μ0.29Hvap0.24ρg0.24)ΔTsat0.24Δpsat0.75h_{nb} = 0.00122\left( \frac{\lambda_l^{0.79} c_{p,l}^{0.45} \rho_l^{0.49}}{\sigma^{0.5} \mu^{0.29} H_{vap}^{0.24} \rho_g^{0.24}} \right)\Delta T_{sat}^{0.24} \Delta p_{sat}^{0.75}
F=(1+Xtt0.5)1.78F = (1 + X_{tt}^{-0.5})^{1.78}
Xtt=(1xx)0.9(ρgρl)0.5(μlμg)0.1X_{tt} = \left( \frac{1-x}{x}\right)^{0.9} \left(\frac{\rho_g}{\rho_l} \right)^{0.5}\left( \frac{\mu_l}{\mu_g}\right)^{0.1}
S=0.96220.5822(tan1(ReLF1.256.18104))S = 0.9622 - 0.5822\left(\tan^{-1}\left(\frac{Re_L\cdot F^{1.25}} {6.18\cdot 10^4}\right)\right)
Parameters
mfloat

Mass flow rate [kg/s]

xfloat

Quality at the specific tube interval []

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

rhogfloat

Density of the gas [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

mugfloat

Viscosity of gas [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

Cplfloat

Heat capacity of liquid [J/kg/K]

Hvapfloat

Heat of vaporization of liquid [J/kg]

sigmafloat

Surface tension of liquid [N/m]

dPsatfloat

Difference in Saturation pressure of fluid at Te and T, [Pa]

Tefloat

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

See also

turbulent_Dittus_Boelter
Forster_Zuber

Notes

[1] and [2] have been reviewed, but the model is only put together in the review of [3]. Many other forms of this equation exist with different functions for F and S.

References

1(1,2)

Chen, J. C. “Correlation for Boiling Heat Transfer to Saturated Fluids in Convective Flow.” Industrial & Engineering Chemistry Process Design and Development 5, no. 3 (July 1, 1966): 322-29. doi:10.1021/i260019a023.

2(1,2)

Edelstein, Sergio, A. J. Pérez, and J. C. Chen. “Analytic Representation of Convective Boiling Functions.” AIChE Journal 30, no. 5 (September 1, 1984): 840-41. doi:10.1002/aic.690300528.

3(1,2)

Bertsch, Stefan S., Eckhard A. Groll, and Suresh V. Garimella. “Review and Comparative Analysis of Studies on Saturated Flow Boiling in Small Channels.” Nanoscale and Microscale Thermophysical Engineering 12, no. 3 (September 4, 2008): 187-227. doi:10.1080/15567260802317357.

Examples

>>> Chen_Edelstein(m=0.106, x=0.2, D=0.0212, rhol=567, rhog=18.09,
... mul=156E-6, mug=7.11E-6, kl=0.086, Cpl=2730, Hvap=2E5, sigma=0.02,
... dPsat=1E5, Te=3)
3289.058731974052
ht.boiling_flow.Lazarek_Black(m, D, mul, kl, Hvap, q=None, Te=None)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in vertical tubes for either upward or downward flow. Correlation is as shown in [1], and also reviewed in [2] and [3].

Either the heat flux or excess temperature is required for the calculation of heat transfer coefficient.

Quality independent. Requires no properties of the gas. Uses a Reynolds number assuming all the flow is liquid.

htp=30Relo0.857Bg0.714klDh_{tp} = 30 Re_{lo}^{0.857} Bg^{0.714} \frac{k_l}{D}
Relo=GtpDμlRe_{lo} = \frac{G_{tp}D}{\mu_l}
Parameters
mfloat

Mass flow rate [kg/s]

Dfloat

Diameter of the channel [m]

mulfloat

Viscosity of liquid [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

Hvapfloat

Heat of vaporization of liquid [J/kg]

qfloat, optional

Heat flux to wall [W/m^2]

Tefloat, optional

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

Notes

[1] has been reviewed.

[2] claims it was developed for a range of quality 0-0.6, Relo 860-5500, mass flux 125-750 kg/m^2/s, q of 1.4-38 W/cm^2, and with a pipe diameter of 3.1 mm. Developed with data for R113 only.

References

1(1,2)

Lazarek, G. M., and S. H. Black. “Evaporative Heat Transfer, Pressure Drop and Critical Heat Flux in a Small Vertical Tube with R-113.” International Journal of Heat and Mass Transfer 25, no. 7 (July 1982): 945-60. doi:10.1016/0017-9310(82)90070-9.

2(1,2)

Fang, Xiande, Zhanru Zhou, and Dingkun Li. “Review of Correlations of Flow Boiling Heat Transfer Coefficients for Carbon Dioxide.” International Journal of Refrigeration 36, no. 8 (December 2013): 2017-39. doi:10.1016/j.ijrefrig.2013.05.015.

3

Bertsch, Stefan S., Eckhard A. Groll, and Suresh V. Garimella. “Review and Comparative Analysis of Studies on Saturated Flow Boiling in Small Channels.” Nanoscale and Microscale Thermophysical Engineering 12, no. 3 (September 4, 2008): 187-227. doi:10.1080/15567260802317357.

Examples

>>> Lazarek_Black(m=10, D=0.3, mul=1E-3, kl=0.6, Hvap=2E6, Te=100)
9501.932636079293
ht.boiling_flow.Li_Wu(m, x, D, rhol, rhog, mul, kl, Hvap, sigma, q=None, Te=None)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is as shown in [1], and also reviewed in [2] and [3].

Either the heat flux or excess temperature is required for the calculation of heat transfer coefficient. Uses liquid Reynolds number, Bond number, and Boiling number.

htp=334Bg0.3(BoRel0.36)0.4klDh_{tp} = 334 Bg^{0.3}(Bo\cdot Re_l^{0.36})^{0.4}\frac{k_l}{D}
Rel=G(1x)DμlRe_{l} = \frac{G(1-x)D}{\mu_l}
Parameters
mfloat

Mass flow rate [kg/s]

xfloat

Quality at the specific tube interval []

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

rhogfloat

Density of the gas [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

Hvapfloat

Heat of vaporization of liquid [J/kg]

sigmafloat

Surface tension of liquid [N/m]

qfloat, optional

Heat flux to wall [W/m^2]

Tefloat, optional

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

Notes

[1] has been reviewed.

[1] used 18 sets of experimental data to derive the results, covering hydraulic diameters from 0.19 to 3.1 mm and 12 different fluids.

References

1(1,2,3)

Li, Wei, and Zan Wu. “A General Correlation for Evaporative Heat Transfer in Micro/mini-Channels.” International Journal of Heat and Mass Transfer 53, no. 9-10 (April 2010): 1778-87. doi:10.1016/j.ijheatmasstransfer.2010.01.012.

2

Fang, Xiande, Zhanru Zhou, and Dingkun Li. “Review of Correlations of Flow Boiling Heat Transfer Coefficients for Carbon Dioxide.” International Journal of Refrigeration 36, no. 8 (December 2013): 2017-39. doi:10.1016/j.ijrefrig.2013.05.015.

3

Kim, Sung-Min, and Issam Mudawar. “Review of Databases and Predictive Methods for Pressure Drop in Adiabatic, Condensing and Boiling Mini/micro-Channel Flows.” International Journal of Heat and Mass Transfer 77 (October 2014): 74-97. doi:10.1016/j.ijheatmasstransfer.2014.04.035.

Examples

>>> Li_Wu(m=1, x=0.2, D=0.3, rhol=567., rhog=18.09, kl=0.086, mul=156E-6, sigma=0.02, Hvap=9E5, q=1E5)
5345.409399239492
ht.boiling_flow.Liu_Winterton(m, x, D, rhol, rhog, mul, kl, Cpl, MW, P, Pc, Te)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is as developed in [1], also reviewed in [2] and [3].

Excess wall temperature is required to use this correlation.

htp=(Fhl)2+(Shnb)2h_{tp} = \sqrt{ (F\cdot h_l)^2 + (S\cdot h_{nb})^2}
S=(1+0.055F0.1ReL0.16)1S = \left( 1+0.055F^{0.1} Re_{L}^{0.16}\right)^{-1}
hl=0.023ReL0.8Prl0.4kl/Dh_{l} = 0.023 Re_L^{0.8} Pr_l^{0.4} k_l/D
ReL=GDμlRe_L = \frac{GD}{\mu_l}
F=[1+xPrl(ρl/ρg1)]0.35F = \left[ 1+ xPr_{l}(\rho_l/\rho_g-1)\right]^{0.35}
hnb=(55ΔTe0.67PPc(0.120.2log10Rp)(log10PPc)0.55MW0.5)1/0.33h_{nb} = \left(55\Delta Te^{0.67} \frac{P}{P_c}^{(0.12 - 0.2\log_{10} R_p)}(-\log_{10} \frac{P}{P_c})^{-0.55} MW^{-0.5}\right)^{1/0.33}
Parameters
mfloat

Mass flow rate [kg/s]

xfloat

Quality at the specific tube interval []

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

rhogfloat

Density of the gas [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

Cplfloat

Heat capacity of liquid [J/kg/K]

MWfloat

Molecular weight of the fluid, [g/mol]

Pfloat

Pressure of fluid, [Pa]

Pcfloat

Critical pressure of fluid, [Pa]

Tefloat, optional

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

Notes

[1] has been reviewed, and is accurately reproduced in [3].

Uses the Cooper and turbulent_Dittus_Boelter correlations.

A correction for horizontal flow at low Froude numbers is available in [1] but has not been implemented and is not recommended in several sources.

References

1(1,2,3)

Liu, Z., and R. H. S. Winterton. “A General Correlation for Saturated and Subcooled Flow Boiling in Tubes and Annuli, Based on a Nucleate Pool Boiling Equation.” International Journal of Heat and Mass Transfer 34, no. 11 (November 1991): 2759-66. doi:10.1016/0017-9310(91)90234-6.

2

Fang, Xiande, Zhanru Zhou, and Dingkun Li. “Review of Correlations of Flow Boiling Heat Transfer Coefficients for Carbon Dioxide.” International Journal of Refrigeration 36, no. 8 (December 2013): 2017-39. doi:10.1016/j.ijrefrig.2013.05.015.

3(1,2)

Bertsch, Stefan S., Eckhard A. Groll, and Suresh V. Garimella. “Review and Comparative Analysis of Studies on Saturated Flow Boiling in Small Channels.” Nanoscale and Microscale Thermophysical Engineering 12, no. 3 (September 4, 2008): 187-227. doi:10.1080/15567260802317357.

Examples

>>> Liu_Winterton(m=1, x=0.4, D=0.3, rhol=567., rhog=18.09, kl=0.086,
... mul=156E-6, Cpl=2300, P=1E6, Pc=22E6, MW=44.02, Te=7)
4747.749477190532
ht.boiling_flow.Sun_Mishima(m, D, rhol, rhog, mul, kl, Hvap, sigma, q=None, Te=None)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is as shown in [1], and also reviewed in [2].

Either the heat flux or excess temperature is required for the calculation of heat transfer coefficient. Uses liquid-only Reynolds number, Weber number, and Boiling number. Weber number is defined in terms of the velocity if all fluid were liquid.

htp=6Relo1.05Bg0.54Wel0.191(ρl/ρg)0.142klDh_{tp} = \frac{ 6 Re_{lo}^{1.05} Bg^{0.54}} {We_l^{0.191}(\rho_l/\rho_g)^{0.142}}\frac{k_l}{D}
Relo=GtpDμlRe_{lo} = \frac{G_{tp}D}{\mu_l}
Parameters
mfloat

Mass flow rate [kg/s]

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

rhogfloat

Density of the gas [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

Hvapfloat

Heat of vaporization of liquid [J/kg]

sigmafloat

Surface tension of liquid [N/m]

qfloat, optional

Heat flux to wall [W/m^2]

Tefloat, optional

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

Notes

[1] has been reviewed.

[1] used 2501 data points to derive the results, covering hydraulic diameters from 0.21 to 6.05 mm and 11 different fluids.

References

1(1,2,3)

Sun, Licheng, and Kaichiro Mishima. “An Evaluation of Prediction Methods for Saturated Flow Boiling Heat Transfer in Mini-Channels.” International Journal of Heat and Mass Transfer 52, no. 23-24 (November 2009): 5323-29. doi:10.1016/j.ijheatmasstransfer.2009.06.041.

2

Fang, Xiande, Zhanru Zhou, and Dingkun Li. “Review of Correlations of Flow Boiling Heat Transfer Coefficients for Carbon Dioxide.” International Journal of Refrigeration 36, no. 8 (December 2013): 2017-39. doi:10.1016/j.ijrefrig.2013.05.015.

Examples

>>> Sun_Mishima(m=1, D=0.3, rhol=567., rhog=18.09, kl=0.086, mul=156E-6, sigma=0.02, Hvap=9E5, Te=10)
507.6709168372167
ht.boiling_flow.Thome(m, x, D, rhol, rhog, mul, mug, kl, kg, Cpl, Cpg, Hvap, sigma, Psat, Pc, q=None, Te=None)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is as developed in [1] and [2], and also reviewed [3]. This is a complicated model, but expected to have more accuracy as a result.

Either the heat flux or excess temperature is required for the calculation of heat transfer coefficient. The solution for a specified excess temperature is solved numerically, making it slow.

h(z)=tlτhl(z)+tfilmτhfilm(z)+tdryτhg(z)h(z) = \frac{t_l}{\tau} h_l(z) +\frac{t_{film}}{\tau} h_{film}(z) + \frac{t_{dry}}{\tau} h_{g}(z)
hl/g(z)=(Nulam4+Nutrans4)1/4k/Dh_{l/g}(z) = (Nu_{lam}^4 + Nu_{trans}^4)^{1/4} k/D
Nulaminar=0.91Pr1/3ReD/L(z)Nu_{laminar} = 0.91 {Pr}^{1/3} \sqrt{ReD/L(z)}
Nutrans=(f/8)(Re1000)Pr1+12.7(f/8)1/2(Pr2/31)[1+(DL(z))2/3]Nu_{trans} = \frac{ (f/8) (Re-1000)Pr}{1+12.7 (f/8)^{1/2} (Pr^{2/3}-1)} \left[ 1 + \left( \frac{D}{L(z)}\right)^{2/3}\right]
f=(1.82log10Re1.64)2f = (1.82 \log_{10} Re - 1.64 )^{-2}
Ll=τGtpρl(1x)L_l = \frac{\tau G_{tp}}{\rho_l}(1-x)
Ldry=vptdryL_{dry} = v_p t_{dry}
tl=τ1+ρlρgx1xt_l = \frac{\tau}{1 + \frac{\rho_l}{\rho_g}\frac{x}{1-x}}
tv=τ1+ρgρl1xxt_v = \frac{\tau}{1 + \frac{\rho_g}{\rho_l}\frac{1-x}{x}}
τ=1fopt\tau = \frac{1}{f_{opt}}
fopt=(qqref)nff_{opt} = \left(\frac{q}{q_{ref}}\right)^{n_f}
qref=3328(PsatPc)0.5q_{ref} = 3328\left(\frac{P_{sat}}{P_c}\right)^{-0.5}
tdry,film=ρlΔHvapq[δ0(z)δmin]t_{dry,film} = \frac{\rho_l \Delta H_{vap}}{q}[\delta_0(z) - \delta_{min}]
δ0D=Cδ0(3νlvpD)0.84[(0.07Bo0.41)8+0.18]1/8\frac{\delta_0}{D} = C_{\delta 0}\left(3\sqrt{\frac{\nu_l}{v_p D}} \right)^{0.84}\left[(0.07Bo^{0.41})^{-8} + 0.1^{-8}\right]^{-1/8}
Bo=ρlDσvp2Bo = \frac{\rho_l D}{\sigma} v_p^2
vp=Gtp[xρg+1xρl]v_p = G_{tp} \left[\frac{x}{\rho_g} + \frac{1-x}{\rho_l}\right]
hfilm(z)=2klδ0(z)+δmin(z)h_{film}(z) = \frac{2 k_l}{\delta_0(z) + \delta_{min}(z)}
δmin=0.3106m\delta_{min} = 0.3\cdot 10^{-6} \text{m}
Cδ,0=0.29C_{\delta,0} = 0.29
nf=1.74n_f = 1.74

if t dry film > tv:

δend(x)=δ(z,tv)\delta_{end}(x) = \delta(z, t_v)
tfilm=tvt_{film} = t_v
tdry=0t_{dry} = 0

Otherwise:

δend(z)=δmin\delta_{end}(z) = \delta_{min}
tfilm=tdry,filmt_{film} = t_{dry,film}
tdry=tvtfilmt_{dry} = t_v - t_{film}
Parameters
mfloat

Mass flow rate [kg/s]

xfloat

Quality at the specific tube interval []

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

rhogfloat

Density of the gas [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

mugfloat

Viscosity of gas [Pa*s]

klfloat

Thermal conductivity of liquid [W/m/K]

kgfloat

Thermal conductivity of gas [W/m/K]

Cplfloat

Heat capacity of liquid [J/kg/K]

Cpgfloat

Heat capacity of gas [J/kg/K]

Hvapfloat

Heat of vaporization of liquid [J/kg]

sigmafloat

Surface tension of liquid [N/m]

Psatfloat

Vapor pressure of fluid, [Pa]

Pcfloat

Critical pressure of fluid, [Pa]

qfloat, optional

Heat flux to wall [W/m^2]

Tefloat, optional

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

Notes

[1] and [2] have been reviewed, and are accurately reproduced in [3].

[1] used data from 7 studies, covering 7 fluids and Dh from 0.7-3.1 mm, heat flux from 0.5-17.8 W/cm^2, x from 0.01-0.99, and G from 50-564 kg/m^2/s.

Liquid and/or gas slugs are both considered, and are hydrodynamically developing. Ll is the calculated length of liquid slugs, and L_dry is the same for vapor slugs.

Because of the complexity of the model and that there is some logic in this function, Te as an input may lead to a different solution that the calculated q will in return.

References

1(1,2,3)

Thome, J. R., V. Dupont, and A. M. Jacobi. “Heat Transfer Model for Evaporation in Microchannels. Part I: Presentation of the Model.” International Journal of Heat and Mass Transfer 47, no. 14-16 (July 2004): 3375-85. doi:10.1016/j.ijheatmasstransfer.2004.01.006.

2(1,2)

Dupont, V., J. R. Thome, and A. M. Jacobi. “Heat Transfer Model for Evaporation in Microchannels. Part II: Comparison with the Database.” International Journal of Heat and Mass Transfer 47, no. 14-16 (July 2004): 3387-3401. doi:10.1016/j.ijheatmasstransfer.2004.01.007.

3(1,2)

Bertsch, Stefan S., Eckhard A. Groll, and Suresh V. Garimella. “Review and Comparative Analysis of Studies on Saturated Flow Boiling in Small Channels.” Nanoscale and Microscale Thermophysical Engineering 12, no. 3 (September 4, 2008): 187-227. doi:10.1080/15567260802317357.

Examples

>>> Thome(m=1, x=0.4, D=0.3, rhol=567., rhog=18.09, kl=0.086, kg=0.2,
... mul=156E-6, mug=1E-5, Cpl=2300, Cpg=1400, sigma=0.02, Hvap=9E5,
... Psat=1E5, Pc=22E6, q=1E5)
1633.008836502032
ht.boiling_flow.Yun_Heo_Kim(m, x, D, rhol, mul, Hvap, sigma, q=None, Te=None)[source]

Calculates heat transfer coefficient for film boiling of saturated fluid in any orientation of flow. Correlation is as shown in [1] and [2], and also reviewed in [3].

Either the heat flux or excess temperature is required for the calculation of heat transfer coefficient. Uses liquid Reynolds number, Weber number, and Boiling number. Weber number is defined in terms of the velocity if all fluid were liquid.

htp=136876(BgWel)0.1993Rel0.1626h_{tp} = 136876(Bg\cdot We_l)^{0.1993} Re_l^{-0.1626}
Rel=GD(1x)μlRe_l = \frac{G D (1-x)}{\mu_l}
Wel=G2DρlσWe_l = \frac{G^2 D}{\rho_l \sigma}
Parameters
mfloat

Mass flow rate [kg/s]

xfloat

Quality at the specific tube interval []

Dfloat

Diameter of the tube [m]

rholfloat

Density of the liquid [kg/m^3]

mulfloat

Viscosity of liquid [Pa*s]

Hvapfloat

Heat of vaporization of liquid [J/kg]

sigmafloat

Surface tension of liquid [N/m]

qfloat, optional

Heat flux to wall [W/m^2]

Tefloat, optional

Excess temperature of wall, [K]

Returns
hfloat

Heat transfer coefficient [W/m^2/K]

Notes

[1] has been reviewed.

References

1(1,2)

Yun, Rin, Jae Hyeok Heo, and Yongchan Kim. “Evaporative Heat Transfer and Pressure Drop of R410A in Microchannels.” International Journal of Refrigeration 29, no. 1 (January 2006): 92-100. doi:10.1016/j.ijrefrig.2005.08.005.

2

Yun, Rin, Jae Hyeok Heo, and Yongchan Kim. “Erratum to ‘Evaporative Heat Transfer and Pressure Drop of R410A in Microchannels; [Int. J. Refrigeration 29 (2006) 92-100].” International Journal of Refrigeration 30, no. 8 (December 2007): 1468. doi:10.1016/j.ijrefrig.2007.08.003.

3

Bertsch, Stefan S., Eckhard A. Groll, and Suresh V. Garimella. “Review and Comparative Analysis of Studies on Saturated Flow Boiling in Small Channels.” Nanoscale and Microscale Thermophysical Engineering 12, no. 3 (September 4, 2008): 187-227. doi:10.1080/15567260802317357.

Examples

>>> Yun_Heo_Kim(m=1, x=0.4, D=0.3, rhol=567., mul=156E-6, sigma=0.02, Hvap=9E5, q=1E4)
9479.313988550184