Numerical and Experimental Investigation of Boiling Heat Transfer for Subcooled Water Flowing in a Small-Diameter Tube
Numerical and Experimental Investigation of Boiling Heat Transfer for Subcooled Water Flowing in a Small-Diameter Tube
复制标题
DOI：
10.1115/mnhmt2019-4163
复制
发表时间：
2019-12
期刊：
影响因子：
--
通讯作者：
M. Shibahara;Qiusheng Liu;K. Hata;K. Fukuda
中科院分区：
文献类型：
--
作者：
M. Shibahara;Qiusheng Liu;K. Hata;K. Fukuda
利用商用计算流体力学（CFD）软件PHOENICS . ver对小直径管内过冷水沸腾换热过程进行了数值模拟。2013. 模拟采用小直径管（d = 1.0-2.0 mm）。给出了管内壁面均匀的指数函数热流密度作为边界条件。内壁边界条件设为防滑。进气温度为302 ~ 312 K。d = 1.0 mm和d = 2.0 mm流速分别为9.29 m/s和2.34 m/s。在实验中，由于热流密度随时间增加，所以暂态分析是从非沸腾区进行的。在PHOENICS程序中采用有限体积法对包括能量方程在内的控制方程进行离散化。采用SIMPLE方法进行数值模拟。为模拟管内沸腾现象，采用腓尼基代码的相间滑移算法，采用欧拉-欧拉双流体模型。实验采用铂管作为实验管（d = 1.0-2.0 mm），通过直流电进行焦耳加热。蒸馏水和去离子水由增压器加压。用指数函数控制管的产热速率，得到非沸腾区瞬态传热特性。实验中，表面过热度随热流密度的增大而增大。数值模拟较好地预测了实验数据。当实验热流密度达到CHF点时，换热系数的预测值比实验值低3.5%左右。
Numerical simulation of boiling heat transfer for subcooled water flowing in a small-diameter tube was conducted using the commercial computational fluid dynamics (CFD) code, PHOENICS ver. 2013. A small-diameter tube (d = 1.0–2.0 mm) was modeled in the simulation. A uniform heat flux with an exponential function was given at the inner tube wall as the boundary conditions. The inner wall boundary condition was set to a non-slip. The inlet temperature ranged from 302 to 312 K. The flow velocities of d = 1.0 mm and d = 2.0 mm are 9.29 m/s and 2.34 m/s, respectively. The transient analysis was carried out from the non-boiling region since the heat flux increased with time in the author’s experiments. The governing equations including the energy equation were discretized using the finite volume method in the PHOENICS code. The SIMPLE method was applied for the numerical simulation. For modeling boiling phenomena in the tube, the Eulerian-Eulerian two-fluid model was adopted using the interphase slip algorithm of PHOENICS code. In the experiment, a platinum tube was used as the experimental tube (d = 1.0–2.0 mm) to conduct joule heating by direct current. The distilled and deionized water was pressured by the pressurizer. The heat generation rate of the tube was controlled with the exponential function to obtain the transient heat transfer characteristics from the non-boiling region. The surface superheat increased as the heat flux increased in the experiment. The numerical simulation predicted the experimental data well. When the heat flux of the experiment was reached to the CHF point, the predicted value of heat transfer coefficient was approximately 3.5 % lower than that of the experiment.