过冷盐溶液流动沸腾传热预测关联式
顾君苹1,2, 刘琦1, 吴玉新1, 王庆功2, 吕俊复1    
1. 清华大学 能源与动力工程系, 热科学与动力工程教育部重点实验室, 北京 100084;
2. 中国空间技术研究院 钱学森空间技术实验室, 北京 100094
摘要:为保证以含盐溶液为工质的蒸发换热设备内的水动力及传热安全,要求能够对盐溶液在受热管内的传热特性实现准确预测。该文开展了垂直上升受热管内NaCl溶液传热实验,获得了盐溶液管内强迫对流及过冷流动沸腾传热系数实验值。基于实验结果,校验了现有传热预测关联式对盐溶液不同传热阶段的传热预测准确性。结合公开报道中的盐溶液及纯水的传热数据进行分析,获得了可用于表征盐溶液过冷流动沸腾传热变化的主要参数及影响规律,据此建立了适应于盐溶液的过冷流动沸腾传热预测关联式,该关联式的预测准确性高,对指导和优化此类换热设备的设计、提高工业生产安全具有重要价值。
关键词传热关联式    盐溶液    过冷流动沸腾    强迫对流    
Heat transfer correlation for subcooled flow boiling of saline solutions
GU Junping1,2, LIU Qi1, WU Yuxin1, WANG Qinggong2, LYU Junfu1    
1. Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Energy and Power Engineering, Tsinghua University, Beijing 100084, China;
2. Qian Xuesen Laboratory of Space Technology, China Academy of Space Technology, Beijing 100094, China
Abstract: Accurate predictions of the heat transfer characteristics of saline solutions is of great importance for modeling the hydrodynamics and heat transfer in evaporative heat exchangers using saline solutions as the heat transfer medium. The heat transfer characteristics of NaCl solutions in a vertical heated tube were investigated experimentally to determine the heat transfer coefficients of NaCl solutions for single-phase forced convection and subcooled flow boiling conditions. The present experimental data was used to verify the prediction accuracies of existing heat transfer correlations for saline solutions. The experimental heat transfer data for saline solutions and pure water were compared to identify the main factors characterizing the heat transfer coefficient variations. A new heat transfer correlation is presented that accurately predicts the subcooled flow boiling heat transfer coefficient of saline solutions for the design and safe operation of industrial applications.
Key words: heat transfer correlation    saline solution    subcooled flow boiling    forced convection    

两相过冷流动沸腾传热广泛存在于核能、航空、航天等各领域的蒸发换热设备中,其传热过程的影响因素众多、物理过程复杂,传热计算分析与预测充满了挑战。为了满足工业设计安全性需求,学者们通常利用已有的传热实验数据进行经验/半经验关联式的建立,从而对不同工况下的设备受热面温度或加热热流密度等进行预测。

近年来,淡水资源短缺的危机日趋加重,海水(或者污水回用水)直流冷却/循环冷却技术得以加速发展,特别是在火电、核电等行业中应用较广[1-2]。我国自20世纪90年代初开始了对海水冷却技术的实验研究,现有研究主要聚焦于发生核态沸腾时加热表面上的结垢及腐蚀等现象,并在防腐、阻垢、防生等关键技术上取得了重大的突破[2]。然而,考虑到海水中含有大量以氯化钠(NaCl)为主的可溶性盐类,使得工质物性与纯水相比差异较大。过冷流动沸腾发生时,汽泡在壁面处产生,伴随蒸汽产生而遗留下来的盐分由液相带走。与纯水相比,盐溶液物性的不同将直接影响流动沸腾中的汽泡参数及行为变化,进而使换热设备内的工质传热特性发生变化[3-5],因此需建立相应的传热预测关联式,对盐溶液在换热设备中的水动力及壁面传热进行准确预测,以保证工业运行安全。

李修伦等[6]、Prodanovic等[7]和方贤德等[8]对现有的过冷流动沸腾传热预测关联式及其适应范围进行了总结。从形式上看,现有过冷流动沸腾传热关联式主要分为:增强型,例如Shah模型[9]、Hata-Noda关联式[10]以及Baburajan等[11]的关联式;叠加型,例如Chen模型[12]、以及在Chen模型基础上发展而来的Steiner等[13]、Yan等[14]以及Gungor等[15]的模型;渐进型,例如Hua等[16]的关联式、Liu-Winterton关联式[17]以及Ramstorfer等[18]的关联式;qTsatn (q为热流密度,ΔTsat为壁面过热度,n为指数)型式关联式,例如McAdams等[19]、Jens & Lottes[20]以及Thom等[21]的关联式;另外,还有部分学者通过无量纲参数分析的方法,获得对传热系数影响较大的相关无量纲参数,例如沸腾数(Bo)、Prandtl数(Pr)和密度比等,进一步通过实验数据进行拟合得到表达式。

由于纯水是工程应用及实验研究中使用的主要工质,因此大多数过冷流动沸腾传热预测关联式是基于纯水实验数据获得的,仅有少数关联式的建立采用了制冷剂和盐溶液等其他工质的实验数据。

根据对现有过冷流动沸腾传热预测关联式的总结,发现现有研究中仅Najibi[22]针对含盐溶液的流动沸腾进行了传热预测模型的建立。Najibi[22]采用盐溶液进行传热实验,基于所获得的实验数据,校验了现有强迫对流传热和池沸腾传热关联式对盐溶液传热预测的适应性,并根据校验结果分别采用Gorenflo关联式[23]和Gnielinski关联式[24],对Chen模型[12]中的核态沸腾传热系数和单相液体强迫对流传热系数进行计算。采用修正后的传热预测模型对该实验中盐溶液过冷流动沸腾传热数据进行预测,预测值与实验值的平均偏差为14.5%。

Yusufova等[25]对浓度为0~220 g/L过冷NaCl溶液进行流动沸腾实验,实验工况为热流密度q=500~1 600 kW/m2,压力P=1~10 MPa,质量流速G=300~1 100 kg/(m2·s)。在他们的实验研究范围内,发现纯水与NaCl溶液的过冷流动沸腾传热系数h与质量流速无关,与q2/3成正相关;壁面过热度ΔTsatR0.25(R为盐溶液浓度,g/L)以及q1/3成正相关。实验中还发现,不同压力下浓度对传热的影响规律不同,溶液的过冷流动沸腾传热系数与压力也有关。但是文中未能给出NaCl溶液过冷流动沸腾传热预测关联式。

经验/半经验传热预测关联式的建立依赖于所采用的实验数据,其使用具有一定局限性,而由于目前含盐溶液流动沸腾传热实验数据的缺失,现有关联式在盐水传热特性预测方面的适应性有待考究。同时为了更好地服务于工程应用[26-28],应力求对形式更为简单、适应范围更广的传热预测关联式的建立方法进行探索。

因此本文首先对比分析了报道中对强迫对流以及过冷流动沸腾传热预测性能较好(通常预测偏差 < 20%)的关联式在盐溶液工况下的传热预测准确性。随后,基于盐溶液在垂直上升受热管内的传热实验数据,推荐了能够用于盐溶液强迫对流的传热预测关联式,分析获得了可用于表征溶液过冷流动沸腾传热变化的主要参数及影响规律,进而提出了盐溶液过冷流动沸腾传热预测关联式的建立方法,并对关联式的预测性能进行了验证,为相关换热设备的设计及安全运行提供依据。

1 盐溶液管内流动沸腾传热实验数据

由于现有公开报道中可溶性盐溶液流动沸腾传热实验数据匮乏,因此本文工作首先搭建了垂直上升受热管内盐溶液核态沸腾传热可视化实验系统,如图 1所示。

图 1 实验系统示意图

实验加热段为方形截面的透明石英管,采用外壁面单面镀ITO透明导电涂层的方法进行加热。实验工况范围为:常压、质量流速G=200~600 kg/(m2·s)、热流密度q=10~180 kW/m2、工质过冷度ΔTsub=8.5~36 K,实验工质包括去离子水以及质量分数为1%~6%的NaCl溶液(海水平均含盐量为3.5%左右,实验中NaCl溶液取质量分数为1%~6%以覆盖海水浓缩倍数为1.5倍的工况)。实验中测量参数主要包括加热段进出口工质温度、加热侧外壁温度、实验段加热量、流动工质的体积流量和实验段散热量等,后通过传热计算得到其相应工况下的传热系数。根据误差传递公式分别计算去离子水、1%、3%、6%的NaCl溶液在不同工况下传热系数的不确定度,其传热系数平均不确定度分别为12.1%、13.3%、12.8%和10.43%,为了保证各图中数据点清晰,仅以6% NaCl溶液过冷流动沸腾为例标注实验所得传热系数不确定度(见图 4c)。具体实验系统介绍、参数测量、不确定度分析及实验过程已在文[29-30]中进行了详细描述,这里不再赘述。基于本实验系统,获得了盐溶液在垂直上升受热管内的单相液体强迫对流以及过冷流动沸腾传热系数实验值。

本文对盐溶液传热预测关联式的分析及建立基于上述实验结果[29-30]以及Najibi[22]报道的NaCl溶液传热系数实验值进行。另外,为了丰富实验数据,获得更为普适性的、工质种类适应范围更广的参数影响规律及传热预测关联式建立方法,下文部分分析中也采用了Najibi[22]、Ahmadi等[31-32]、Okawa等[33]以及Phillips[34]公布的垂直上升受热管内纯水的过冷流动沸腾传热实验数据。

一般情况下,海水冷却技术中工质为过冷状态,工质在受热管内将经历强迫对流和过冷流动沸腾的不同传热阶段,因此下文分别针对这2个传热阶段进行盐溶液管内传热预测关联式的研究。

2 盐溶液单相强迫对流传热预测

针对垂直上升受热管内盐溶液强迫对流流动传热,本文校验了Petukhov-Popov关联式[35]、Gnielinski关联式[24]、Kays-Crawford关联式[36]、Whitaker关联式[37]以及Dittus-Boelter预测关联式对实验中纯水及NaCl溶液强迫对流传热系数的预测性能,各传热预测关联式的具体表达式列于表 1

表 1 单相强迫对流传热预测关联式
关联式 表达式
Petukhov-Popov[35] $ \begin{align} & Nu=\frac{\frac{F}{8}{Re}Pr}{1.07+12.7\sqrt{\frac{F}{8}\left( {{Pr }^{\frac{2}{3}}}-1 \right)}}\text{, } \\ & F={{(1.82\lg Re-1.64)}^{-2}} \\ \end{align}$
Gnielinski [24] $\begin{align} & Nu=\frac{\frac{F}{8}({Re}-1000)Pr }{1+12.7\sqrt{\frac{F}{8}\left( Pr \frac{2}{3}-1 \right)}}\left[ 1+{{\left( \frac{D}{L} \right)}^{\frac{2}{3}}} \right], \\ & F={{(1.82\lg Re-1.64)}^{-2}} \\ \end{align} $
Kays-Crawford[36] $\begin{gathered} N u=5+0.015 {Re}^{a} P r^{b}, \\ a=0.88-\frac{0.24}{4.0+P r}, \\ b=0.333+\exp (-0.6 P r) \end{gathered} $
Whitaker[37] $ N u=0.015 R e^{0.83} P r^{0.42}\left(\frac{\mu_{1}}{\mu_{\mathrm{w}}}\right)^{0.14}$
Dittus-Boelter Nu=0.023Re0.8Pr0.4

表 1Nu表示Nusselt数;Re表示Reynolds数;Pr表示Prandtl数;μ为动力黏度,N·s/m2,下角标l、w代表液体和壁面。

上述强迫对流传热预测关联式对本文实验不同工况下NaCl溶液强迫对流传热系数的预测值与实验值对比如图 2所示。可以看出,所选取的关联式都可以在可接受误差范围内对NaCl溶液的强迫对流传热系数进行预测。Petukhov-Popov关联式[35]、Dittus-Boelter预测关联式、Gnielinski关联式[24]、Kays-Crawford关联式[36]以及Whitaker关联式[37]对本实验中盐溶液强迫对流传热系数预测的平均误差分别为12.0%、8.8%、10.1%、24.3%以及8.1%。Dittus-Boelter预测关联式以及Whitaker关联式[37]的预测性能尤为突出。

图 2 盐溶液强迫对流传热预测值与实验值对比

考虑到在Najibi[22]的研究中,发现Whitaker关联式[37]对他们实验中NaCl溶液强迫对流传热系数的预测值偏低。因此,进一步采用Dittus-Boelter预测关联式对Najibi[22]的NaCl溶液强迫对流实验工况进行预测,其预测值与实验值的对比如图 3所示。可以看出,Dittus-Boelter关联式可以准确地对Najibi[22]实验中的NaCl溶液强迫对流传热系数进行预测。

图 3 Dittus-Boelter关联式预测值与Najibi实验中NaCl溶液强迫对流传热系数实验值对比

基于图 2图 3中预测值与实验值的对比结果,本文推荐可采用Dittus-Boelter关联式对垂直上升受热管内NaCl盐溶液强迫对流传热进行预测。

3 盐溶液过冷流动沸腾传热预测 3.1 过冷流动沸腾传热预测关联式的比较分析

为了对比已有过冷流动沸腾传热预测关联式对过冷NaCl溶液流动沸腾传热的预测准确性,采用本实验中盐溶液过冷流动沸腾的传热实验数据对文献中报道的预测偏差较小的过冷流动沸腾传热预测关联式进行校验,包括Shah模型[9]、Gungor-Winterton模型[15]、Liu-Winterton关联式[17]以及Najibi模型[22],各关联式的具体表达式列于表 2中。Gungor-Winterton模型[15]、Najibi模型[22]以及Liu-Winterton关联式[17]在计算中均需采用质量流含汽率x进行计算,在过冷度较高的情况下可按照x≈0近似计算。

图 4所示为各关联式的预测值与实验值对比,可以看出Shah模型[9]以及Liu-Winterton关联式[17]对本实验中过冷NaCl溶液流动沸腾传热系数预测较好,其他2个关联式预测值明显偏高。Shah模型[9]、Gungor-Winterton模型[15]、Liu-Winterton关联式[17]以及Najibi模型[22]对本实验中NaCl溶液过冷沸腾传热系数的预测平均偏差分别为18.8%、116.4%、12.3%以及47.9%。

图 4 盐溶液过冷流动沸腾传热预测值与实验值对比

Liu-Winterton关联式[17]对本文实验结果的预测准确度最高,该关联式计算中包含质量流含汽率x,对于流体过冷度较低的工况,沸腾加强,汽泡生成增多,x≈0的假设将会引起较大的计算偏差导致计算结果偏低,因此在过冷度较低时需要采用模型对过冷流动沸腾中的实际质量流含汽率进行计算。经参数敏感度测试发现,Liu-Winterton关联式[17]对质量流含汽率的敏感度较高,x较小的变化将会引起计算结果变化较大,在这种情况下,该关联式的预测准确度将大大依赖于实际质量流含汽率模型的计算精度。同时,采用子模型进行计算使得关联式的形式更加复杂,涉及参数更多。

表 2 过冷流动沸腾传热预测关联式
关联式 表达式
Shah[9] $q=h_{\mathrm{tp}} \Delta T_{\mathrm{sat}}=\varphi h_{\mathrm{sp}, 1} \Delta T_{\mathrm{sat}} \text {, 其中, } \varphi= \begin{cases}\varphi_{0} & \text { 低过冷度情况 } \\ \varphi_{0}+\Delta T_{\mathrm{sub}} / \Delta T_{\mathrm{sat}} & \text { 高过冷度情况 }\end{cases},\\ \varphi_{0}=\left\{\begin{array}{ll} 230 {Bo}^{0.5} & \text { 当 } B o>0.3 \times 10^{-4} \\ 1+46 {Bo}^{0.5} & \text { 当 } B o <0.3 \times 10^{-4} \end{array}, \begin{cases}\Delta T_{\mathrm{sub}} / \Delta T_{\mathrm{sat}} \geqslant 2 & \text { 高过冷度情况 } \\ \Delta T_{\mathrm{sub}} / \Delta T_{\mathrm{sat}} <2 & \begin{cases}\Delta T_{\mathrm{sub}} / \Delta T_{\mathrm{sut}}>6.3 \times 10^{4} B o^{1.25} & \text { 高过冷度情况 } \\ \Delta T_{\mathrm{su}} / \Delta T_{\mathrm{st}} <6.3 \times 10^{4} B o^{1.25} & \text { 低过冷度情况 }\end{cases} \end{cases}\right. $
Gungor-Winterton[15] $\begin{aligned} &h_{\mathrm{tp}}=F h_{\mathrm{sp}, 1}+S h_{\mathrm{pol}} \text {, 其中, } F=1+24000 B o^{1.16}+1.37\left(\frac{1}{X_{\mathrm{t}}}\right)^{0.86} \text {, } \\ &X_{\mathrm{t}}=\left(\frac{1-x}{x}\right)^{0.9}\left(\frac{\rho_{\mathrm{v}}}{\rho_{\mathrm{l}}}\right)^{0.5}\left(\frac{\mu_{\mathrm{l}}}{\mu_{\mathrm{v}}}\right)^{0.1}, S=1 /\left(1+1.15 \times 10^{-6} F^{2} R e_{1}^{1.17}\right) \end{aligned} $
Najibi[22] $\begin{gathered} h_{\mathrm{tp}}=F h_{\mathrm{sp}, 1}+S h_{\mathrm{pool}}, \text { 其中, }\left\{\begin{array}{ll} F=1 & \text { 当 } 1 / X_{\mathrm{tt}} \leqslant 0.1 \\ F=2.35\left(\frac{1}{X_{\mathrm{tt}}}+0.213\right)^{0.736} & \text { 当 } 1 / X_{\mathrm{tt}}>0.1 \end{array}, \right. \\ S=1 /\left(1+2.53 \times 10^{-6} R e_{\mathrm{tp}}^{1.17}\right), R e_{\mathrm{tp}}=R e(1-x) F^{1.25} \end{gathered} $
Liu-Winterton[17] $ \begin{aligned} q^{2} &=\left[F h_{\mathrm{sp}, 1}\left(T_{\mathrm{iw}}-T_{\mathrm{f}}\right)\right]^{2}+\left[S h_{\mathrm{pool}}\left(T_{\mathrm{iw}}-T_{\mathrm{sat}}\right)\right]^{2} \\ \text { 其中, } F &=\left[1+x P r_{1}\left(\frac{\rho_{\mathrm{l}}}{\rho_{\mathrm{v}}}-1\right)\right]^{0.35}, S=\frac{1}{1+0.055 F^{0.1} R e_{1}^{0.16}} \end{aligned}$
注:htphsp, lhpool分别表示两相沸腾传热系数、单相液体强迫对流传热系数、池沸腾传热系数,W/(m2·K);ΔTsat、ΔTsub分别为壁面过热度、工质过冷度,K;Bo表示沸腾数;x表示管内质量流含汽率;Xtt为Martinelli参数;ρ表示密度,kg/m3TiwTf分别表示内壁面温度和流体温度,K;下角标tp、l、v分别表示两相流动、液相、汽相;Shah模型中hsp, l采用Dittus-Boelter公式计算;Gungor-Winterton模型以及Liu-Winterton关联式中hsp, l采用Dittus-Boelter公式计算,hpool采用Cooper关联式[38]计算;Najibi模型中,hsp, lhpool分别采用Gorenflo关联式[23]和Gnielinski关联式[24]计算。

3.2 关联式的建立与验证

本文提出一种过冷流动沸腾传热预测关联式的建立方法,其基本形式借鉴增强型的沸腾传热预测关联式:

$ h_{\mathrm{tp}}=\psi h_{\mathrm{sp}, \text{l}}. $ (1)

其中,ψ为增强因子,基于本文第2部分对强迫对流传热关联式的分析讨论,hsp, l采用Dittus-Boelter关联式进行计算。

下面对增强因子ψ的建立依据进行分析。如图 5所示为增强因子ψ与沸腾数Bo间的关系,ψ通过过冷流动沸腾传热系数实验值与Dittus-Boelter关联式计算的相应工况下的预测值相除得到。由图 5中本文实验数据(质量分数1%~6%的盐溶液)可以看出,ψBo间存在明显的线性关系,ψBo的增大而增加,且在6%质量分数范围内浓度对这一关系的影响不大。除了Boψ产生影响外,流体过冷度ΔTsub也与增强因子ψ存在明显的关系,ψ随ΔTsub的增加而减小。如图 5所示,在Najibi[22]的NaCl溶液过冷流动沸腾传热实验中,ψ同样呈现出随着Bo的增加而增大的趋势。另外,Najibi的数据在ψ-Bo图中分布较为分散,经分析发现,这是由于其不同实验工况下流体过冷度不同导致的,Najibi[22]的实验数据规律与本文实验一致。

图 5 增强因子ψBo及流体过冷度的关系

利用本文中的NaCl溶液实验数据进行关联式的整理可得到如下增强因子的表达式:

$ \begin{gathered} \psi=a \cdot B o+b, \\ a=6.081 \Delta T_{\text {sub }}^{2}-429.51 \Delta T_{\text {sub }}+11066, \\ b=0.0004 \Delta T_{\text {sub }}^{2}-0.0236 \Delta T_{\text {sub }}+0.8566 . \end{gathered} $ (2)

采用Najibi[22]实验中的NaCl溶液过冷流动沸腾传热系数实验数据(压力范围为P=0.101~0.12 MPa)对关联式(1)-(2)进行验证,如图 6所示,关联式对加热管内NaCl溶液过冷流动沸腾传热系数预测准确度较高,平均偏差为9.6%。

图 6 新传热关联式对盐溶液传热预测准确度校验

由于盐溶液传热实验数据有限,上述讨论中未能对不同压力下的NaCl溶液过冷流动沸腾传热情况进行分析,现利用公开报道中纯水的实验数据进一步进行压力对过冷流动沸腾传热系数影响的分析。纯水实验数据采用Najibi[22]、Ahmadi等[31-32]、Okawa等[33]以及Phillips[34]等报道的垂直上升受热管内纯水过冷流动沸腾传热实验数据,压力范围为0.1~0.8 MPa。

图 7a所示,Najibi[22]、Ahmadi等[31-32]、Okawa等[33]以及Phillips[34]的纯水过冷流动沸腾传热实验数据均符合上述讨论的ψBo以及ψ与△Tsub间的关系。Phillips[34]的实验数据显示,当压力从0.15 MPa升高到0.20 MPa时,ψ随之增大,压力在0.105~0.15 MPa范围内变化时,其对ψ的影响不大。

图 7 压力对增强因子的影响分析

为了显示清楚,将Najibi[22]、Ahmadi等[31-32]以及Okawa等[33]的实验数据另作于图 7b中。Ahmadi等[31-32]的实验压力范围为0.1~0.8 MPa,经进一步数据分析发现,若保持流体过冷度一定,Ahmadi等[31-32]的实验数据显示压力对ψ的影响不大。基于以上分析,说明压力可能也是对ψ产生影响的重要参数,但需要进一步利用不同压力下的实验数据进行其影响规律总结。由图 7可知,本节提出的过冷流动沸腾传热预测关联式的建立方法同样适应于纯水工况。工质不同,式(2) 中的系数ab需重新进行参数拟合。综上,本文基于现有实验数据发展的NaCl溶液过冷流动沸腾传热预测关联式(1)-(2) 的适应范围为:P=0.1~0.15 MPa,Bo=3×10-5~0.4×10-3,ΔTsub < 40 K。

4 结论

本文针对含盐溶液在蒸发换热设备内的传热预测关联式进行研究,分析了现有单相强迫对流及过冷流动沸腾传热预测关联式对盐溶液管内传热的预测准确性。推荐采用Dittus-Boelter关联式对管内盐溶液强迫对流传热进行预测。过冷流动沸腾的工况可采用Liu-Winterton关联式进行预测,但过冷度较低时,该关联式的预测精度受真实质量流含汽率计算的影响较大。为此,本文对NaCl溶液过冷流动沸腾传热进行实验研究,并通过实验数据分析,获得了与其传热系数变化密切相关的参数,讨论了沸腾数(Bo)、流体过冷度以及压力对盐溶液过冷流动沸腾传热的影响,建立了适用于盐溶液的过冷流动沸腾传热预测关联式,传热预测准确性得以提高。

参考文献
[1]
李亚红. 我国海水冷却技术的应用现状及发展应对策略[J]. 应用化工, 2017, 46(12): 2431-2434, 2440.
LI Y H. Application-status and developing countermeasures of seawater cooling technology in China[J]. Applied Chemical Industry, 2017, 46(12): 2431-2434, 2440. DOI:10.3969/j.issn.1671-3206.2017.12.036 (in Chinese)
[2]
张玉忠, 彭晓敏. 浅谈海水循环冷却处理技术[J]. 工业水处理, 2004, 24(8): 14-17.
ZHANG Y Z, PENG X M. Treatment technology of seawater as circulating cooling water[J]. Industrial Water Treatment, 2004, 24(8): 14-17. DOI:10.3969/j.issn.1005-829X.2004.08.004 (in Chinese)
[3]
FAZEL S A A, SHAFAEE S B. Bubble dynamics for nucleate pool boiling of electrolyte solutions[J]. Journal of Heat Transfer, 2010, 132(8): 081502. DOI:10.1115/1.4001315
[4]
HAMZEKHANI S, FALAHIEH M M, KAMALIZADEH M R, et al. Experimental study on bubble departure frequency for pool boiling of water/NaCl solutions[J]. Heat and Mass Transfer, 2015, 51(9): 1313-1320. DOI:10.1007/s00231-015-1502-x
[5]
PEYGHAMBARZADEH S M, HATAMI A, EBRAHIMI A, et al. Photographic study of bubble departure diameter in saturated pool boiling to electrolyte solutions[J]. Chemical Industry and Chemical Engineering Quarterly, 2014, 20(1): 143-153. DOI:10.2298/CICEQ120707120P
[6]
李修伦, 黄鸿鼎, 冯亚云. 两相流沸腾传热研究的进展[J]. 化学工程, 1981(2): 3-10.
LI X L, HUANG H D, FENG Y Y. Research progress of two-phase flow boiling heat transfer[J]. Chemical Engineering, 1981(2): 3-10. (in Chinese)
[7]
PRODANOVIC V, FRASER D, SALCUDEAN M. On the transition from partial to fully developed subcooled flow boiling[J]. International Journal of Heat and Mass Transfer, 2002, 45(24): 4727-4738. DOI:10.1016/S0017-9310(02)00197-7
[8]
FANG X D, YUAN Y L, XU A Y, et al. Review of correlations for subcooled flow boiling heat transfer and assessment of their applicability to water[J]. Fusion Engineering and Design, 2017, 122: 52-63. DOI:10.1016/j.fusengdes.2017.09.008
[9]
SHAH M M. A general correlation for heat transfer during subcooled boiling in pipes and annuli[J]. ASHRAE Transactions, 1977, 83: 205-217.
[10]
HATA K, NODA N. Turbulent heat transfer for heating of water in a short vertical tube[J]. Journal of Power and Energy Systems, 2008, 2(1): 318-329. DOI:10.1299/jpes.2.318
[11]
BABURAJAN P K, BISHT G S, GUPTA S K, et al. Measurement of subcooled boiling pressure drop and local heat transfer coefficient in horizontal tube under LPLF conditions[J]. Nuclear Engineering and Design, 2013, 255: 169-179. DOI:10.1016/j.nucengdes.2012.10.012
[12]
CHEN J C. A correlation for boiling heat transfer to saturated fluids in convective flow[J]. Industrial & Engineering Chemistry Process Design and Development, 1966, 5(3): 322-329.
[13]
STEINER H, KOBOR A, GEBHARD L. A wall heat transfer model for subcooled boiling flow[J]. International Journal of Heat and Mass Transfer, 2005, 48(19-20): 4161-4173. DOI:10.1016/j.ijheatmasstransfer.2005.03.032
[14]
YAN J G, BI Q C, LIU Z H, et al. Subcooled flow boiling heat transfer of water in a circular tube under high heat fluxes and high mass fluxes[J]. Fusion Engineering and Design, 2015, 100: 406-418. DOI:10.1016/j.fusengdes.2015.07.007
[15]
GUNGOR K E, WINTERTON R H S. General correlation for flow boiling in tubes and annuli[J]. International Journal of Heat and Mass Transfer, 1986, 29(3): 351-358. DOI:10.1016/0017-9310(86)90205-X
[16]
HUA S Y, HUANG R H, LI Z, et al. Experimental study on the heat transfer characteristics of subcooled flow boiling with cast iron heating surface[J]. Applied Thermal Engineering, 2015, 77: 180-191. DOI:10.1016/j.applthermaleng.2014.11.082
[17]
LIU Z, WINTERTON R H S. A general correlation for saturated and subcooled flow boiling in tubes and annuli, based on a nucleate pool boiling equation[J]. International Journal of Heat and Mass Transfer, 1991, 34(11): 2759-2766. DOI:10.1016/0017-9310(91)90234-6
[18]
RAMSTORFER F, STEINER H, BRENN G, et al. Subcooled boiling flow heat transfer from plain and enhanced surfaces in automotive applications[J]. Journal of Heat Transfer, 2008, 130(1): 011501. DOI:10.1115/1.2780178
[19]
MCADAMS W H, KENNEL W E, MINDEN C S, et al. Heat transfer at high rates to water with surface boiling[J]. Industrial and Engineering Chemistry, 1949, 41(9): 1945-1953. DOI:10.1021/ie50477a027
[20]
JENS W H, LOTTES P A. Analysis of heat transfer, burnout, pressure drop and density date for high-pressure water (No. ANL-4627)[R]. US: Argonne National Lab, 1951.
[21]
THOM J R S, WALKER W M, FALLON T A, et al. Boiling in sub-cooled water during flow up heated tubes or annuli[J]. Proceedings of the Institution of Mechanical Engineers, Conference Proceedings, 1965, 180(3): 226-246. DOI:10.1243/PIME_CONF_1965_180_117_02
[22]
NAJIBI S H. Heat transfer and heat transfer fouling during subcooled flow boiling for electrolyte solutions[D]. Guildford, UK: The University of Surrey, 1997.
[23]
GORENFLO D, SOKOL P. Prediction method of pool boiling heat transfer with cryogenic liquids[J]. International Journal of Refrigeration, 1988, 11(5): 315-320. DOI:10.1016/0140-7007(88)90095-3
[24]
GNIELINSKI V. Forced convection in ducts, section 2.5.1[M]//SCHVNDER E U. Heat Exchanger Design Handbook. Washington: Hemisphere Publishing Corporation, 1983.
[25]
YUSUFOVA V D, UGREKELIDZE G P, BRONSHTEIN A I. The rate of heat transfer of seawater and NaCl solutions under conditions of one-and two-phase flow in tubes[J]. Desalination, 1978, 26(2): 175-180. DOI:10.1016/S0011-9164(00)82199-2
[26]
黄文仕, 吴玉新, 冯乐乐, 等. 高速圆射流中典型非球形颗粒的扩散特性[J]. 清华大学学报(自然科学版), 2020, 60(6): 485-492.
HUANG W S, WU Y X, FENG L L, et al. Dispersion characteristics of typical non-spherical particles in a high-speed round jet[J]. Journal of Tsinghua University (Science and Technology), 2020, 60(6): 485-492. (in Chinese)
[27]
冯乐乐, 吴玉新, 张海, 等. 转轮分离器风量和转速对叶片流道涡的影响[J]. 清华大学学报(自然科学版), 2020, 60(6): 493-499.
FENG L L, WU Y X, ZHANG H, et al. Effect of air flow rate and rotational speed on vortices between neighboring blades in turbo air classifiers[J]. Journal of Tsinghua University (Science and Technology), 2020, 60(6): 493-499. (in Chinese)
[28]
唐国力, 吴玉新, 顾君苹, 等. 垂直上升光管中气液两相摩擦因子分析[J]. 清华大学学报(自然科学版), 2020, 60(6): 500-506.
TANG G L, WU Y X, GU J P, et al. Comparison of two-phase empirical multiplier correlations for high pressure steam-water mixtures flowing upward in a vertical smooth tube[J]. Journal of Tsinghua University (Science and Technology), 2020, 60(6): 500-506. (in Chinese)
[29]
GU J P, WU Y X, TANG G L, et al. Experimental study of heat transfer and bubble behaviors of NaCl solutions during nucleate flow boiling[J]. Experimental Thermal and Fluid Science, 2019, 109: 109907. DOI:10.1016/j.expthermflusci.2019.109907
[30]
GU J P, TANG G L, WU Y X, el al. Heat transfer characteristics and bubble behaviors during nucleate flow boiling for sodium chloride solution[C]//Proceedings of the ASME 2019 Heat Transfer Summer Conference Collocated with the ASME 2019 13th International Conference on Energy Sustainability. Bellevue, USA: ASME, 2019.
[31]
AHMADI R, UENO T, OKAWA T. Bubble dynamics at boiling incipience in subcooled upward flow boiling[J]. International Journal of Heat and Mass Transfer, 2012, 55(1-3): 488-497. DOI:10.1016/j.ijheatmasstransfer.2011.09.050
[32]
AHMADI R, UENO T, OKAWA T. Visualization study on the mechanisms of net vapor generation in water subcooled flow boiling under moderate pressure conditions[J]. International Journal of Heat and Mass Transfer, 2014, 70: 137-151. DOI:10.1016/j.ijheatmasstransfer.2013.10.073
[33]
OKAWA T, KUBOTA H, ISHIDA T. Simultaneous measurement of void fraction and fundamental bubble parameters in subcooled flow boiling[J]. Nuclear Engineering and Design, 2007, 237(10): 1016-1024. DOI:10.1016/j.nucengdes.2006.12.010
[34]
PHILLIPS B A. Experimental investigation of subcooled flow boiling using synchronized high speed video, infrared thermography, and particle image velocimetry[D]. Massachusetts, USA: Massachusetts Institute of Technology, 2014.
[35]
PETUKHOV B S. Heat transfer and friction in turbulent pipe flow with variable physical properties[J]. Advances in Heat Transfer, 1970, 6: 503-564.
[36]
KAYS W M, CRAWFORD M E. Convective heat and mass transfer[M]. 2nd ed. New York: McGraw-Hill Book Company, 1980.
[37]
WHITAKER S. Forced convection heat transfer correlations for flow in pipes, past flat plates, single cylinders, single spheres, and for flow in packed beds and tube bundles[J]. AIChE Journal, 1972, 18(2): 361-371. DOI:10.1002/aic.690180219
[38]
COOPER M G. Heat flow rates in saturated nucleate pool boiling-a wide-ranging examination using reduced properties[J]. Advances in Heat Transfer, 1984, 16: 157-239.