1. 引言
核反应堆中板状燃料元件由于具有换热面积大、结构紧凑等优势而被广泛关注,研究人员针对矩形通道中气液两相流的两相流型、压降及空泡份额分布开展了相关实验和数值研究,其中实验大多集中于非绝热工况 [1] 。由于矩形通道窄缝结构,气泡间的相互碰撞空间受到限制,碰撞和聚合以沿板方向的横向作用为主 [2] 。
采用高速摄影机能直观地观察到流动沸腾时气泡的产生、生长、脱离、混合以及两相流流型的变化,结合流动沸腾的传热特性,能更好地理解其传热机理。张永豪等 [3] 对矩形通道的速度边界层进行了空气–水可视化实验研究,分析了矩形通道边界层内速度分布和雷诺切应力等流场特性,探究了通道间隙对边界层的影响规律。Liu等 [4] 根据水平方向窄缝矩形通道空气–水可视化实验分析了一维气泡流动的特性界面输运特性,评估了输运方程的适用性。Lu等 [5] 根据可视化的矩形通道气泡分布特性分析了气泡流动与核速率和聚结速率的影响特性。居一伟等 [6] 基于蒸汽加热垂直矩形窄通道实验分析了其传热特性,并以Kandlikar公式为基础改进拟合了新的经验关系式。国内外已有许多两相流可视化实验研究,受限于高压可视化实验本体的设计,工况大都在常压或者低压范围。为了解矩形窄缝通道加热工况下的过冷沸腾气泡分布特性,本文建立了蒸汽加热垂直矩形窄通道沸腾传热特性实验台,分析入口温度、压力、流量和壁面热流密度下的气泡形态。
2. 实验装置
实验回路(图1)包含一回路和二回路两个回路系统、仪表系统、电气系统。一回路主要由主泵、稳压器、预热段、实验段、过滤器、换热器等设备及相应管道、阀门构成。回路设计压力2.5 MPa,主泵流量3.5 t/h。补水箱中的去离子水由柱塞泵驱动,经预热器预热到设定温度后,自下而上进入实验本体,再经换热器冷却后回到补水箱,形成循环。二回路由冷却水箱、冷却水泵、换热器及相连管道、阀门组成。换热器为套管式逆流换热器,实验运行时使用氮气稳压器调节实验回路压力,实验本体通过低电压大电流直流发电机直接进行加热。
2.1. 矩形通道可视化试验段
矩形通道试验段用于开展0~1 MPa压力条件下水工质的流动沸腾核化现象、沸腾工况下相界面形态及行为观测、空泡份额测算等可视化实验。实验段采用一面加热三面可视窗设计,包含2个方向的观测可视窗,以实现高速相机从流道正面和侧面进行拍摄,正面窗口用于观测壁面核化行为,侧面为左右通透的透视可视窗,除可观测近壁核化行为外,还可对沸腾过程中相界面沿径向、轴向演化过程进行观测研究。可视化实验本体总体设计图如图2(a)所示,测点布置图如图2(b)所示,试验段的几何参数见表1。距加热起点240 mm的可视窗的总高度为400 mm,宽度比同位置处的流道宽2 mm,能够保证全流场宽度的图形观测和图形数据记录。实验段测点包括进口/出口流体温度测点、加热壁面外壁面温度测点(沿轴向分别距加热段起点283 mm、355 mm、427 mm、499 mm、571 mm布置5个外壁温测点),用于实验过程中采集实验段的热工参数。
![](//html.hanspub.org/file/5-3150271x8_hanspub.png?20230724104856847)
Figure 1. Schematic of rectangular channel experiment loop
图1. 矩形通道实验回路示意图
![](Images/Table_Tmp.jpg)
Table 1. The geometry parameter of the test channel
表1. 实验段几何参数
(a) 结构组成
(b) 测点布置
Figure 2. Schematic of rectangle test channel
图2. 矩形通道实验段结构示意图
2.2. 实验工况
为探究矩形通道内在不同参数范围下过冷沸腾的气泡生长特性,流场可视化测量工况选取覆盖了泡状流区域到弹状流区域工况,如表2所示。
![](Images/Table_Tmp.jpg)
Table 2. Experimental parameter range
表2. 实验工况参数范围
实验段进出口温度分别为Tin1和Tout,流量为F,给定功率为P,根据热量公式,
,得到全段有效功率
为公式(1)。由于实验段的散热,热效率(
)在30%~50%,温度升高,热效率升高。热效率测量结果见表3。
(1)
![](Images/Table_Tmp.jpg)
Table 3. Experimental section thermal efficiency
表3. 实验段热效率
实验回路通过调节稳压器和背压阀调节系统压力,通过回流阀和实验支路调节阀调节实验段进口流量,通过预热器调节实验段进口温度。当实验段进口压力、流量、温度达到预设工况后,开展沸腾传热实验,逐步增加实验段功率,分别通过热电偶和高速相机监控壁温和流型演变规律,直至达到指定工况,记录实验数据,降低功率,切换工况。
高速摄像机在5 min内拍摄1500张图片进行空泡份额的测定,根据可视化图像后处理程序测得空泡份额为时均空泡份额,根据测点位置,生成某一位置存在空泡的时间占比,将各个位置的时均空泡份额进行成像,得到实验段内空泡份额分布。
3. 实验结果与分析
矩形窄缝通道由于壁面的限制作用,区别于常规圆管通道,带来了气泡形变、局部湍流变化等问题,会对气泡的聚合、破碎、形变等计算造成严重影响。为了解气泡在矩形窄缝通道内的生长情况,图3给出了不同工况下,气泡在不同压力、入口温度和不同功率下的可视化图片。
![](//html.hanspub.org/file/5-3150271x18_hanspub.png?20230724104856847)
![](//html.hanspub.org/file/5-3150271x17_hanspub.png?20230724104856847)
![](//html.hanspub.org/file/5-3150271x16_hanspub.png?20230724104856847)
![](//html.hanspub.org/file/5-3150271x15_hanspub.png?20230724104856847)
2.34 kW2.66 kW 3.07 kW 3.47 kW
工况一:压力0.506 MPa,进口温度120.3℃,流量91.4 kg/(m2∙s),功率逐渐增加
3.1. 数据处理
根据可视化图像后处理程序测得空泡份额为时均空泡份额,方法是基于测点位置,测量5 min时间内某一位置存在空泡的时间占比,将各个位置的时均空泡份额进行成像,得到实验段内空泡份额分布。当实验工况为压力0.767 MPa,进口温度138.3℃,87.8 kg/m2s,功率2.848 kW,高速摄像仪5 min长时间拍摄模式得到1500张照片,经过可视化图像后处理程序处理后,得到可视化实验段时均空泡份额(图4),以及加热板不同距离空泡份额(图5)。
![](//html.hanspub.org/file/5-3150271x25_hanspub.png?20230724104856847)
Figure 4. Time-average avoid fraction in experimental section
图4. 实验段时均空泡份额
![](//html.hanspub.org/file/5-3150271x26_hanspub.png?20230724104856847)
Figure 5. Average avoid fraction of heating plate at different distances in experimental section
图5. 实验段加热板不同距离平均空泡份额
3.2. 参数分析
流体入口段处于过冷状态,通道右侧为加热板,功率较低时实验段内没有气泡产生。随着功率的增加,加热板附近开始出现气泡,并向上流动。功率持续增加,实验段内气泡变大,数量变多呈现出泡状流,气泡位置贴近加热板。随着功率继续增大,实验段内开始出现弹状流。从空泡份额的云图可以看出(图4),气泡从壁面产生,随后向主流区扩散,气泡产生后,在同一轴向位置上,由于壁面润滑力的作用,空泡份额最高的地方不是壁面,而是离壁面一定距离的地方。气泡在离开成核位置后,总是沿着受热面滑动,并且在滑动过程中,滑动气泡的直径增大(图3)。
比较不同热流密度可以发现,随着热流密度的增大,壁温升高的速度变快,过冷沸腾起始点随着热流密度的增加而提前,且气泡的脱离直径随加热功率的增大而增大。质量流速越大,壁温越低,较不容易发生过冷沸腾,需要更高的流体温度才会发生过冷沸腾。
4. 结论
本文通过可视化实验研究了压力、壁面过热度、质量流速和流体过冷度对矩形窄缝通道气泡分布特性的影响,得到如下结论:
1) 三面可视窗加热实验的工况范围为系统压力0.5~1.0 MPa,入口流量50~200 kg/(m2∙s),过冷度小于40 K,沸腾段壁面热流100~500 kW/m2,观测流型范围为泡状流区域到弹状流区域。
2) 基于可视化图像后处理程序可以得到试验段二维分布的时均空泡份额。
3) 实验观测到在受热面上成核的气泡离开成核点后,会沿着垂直加热壁面向上滑移,最高空泡份额点出现在壁面附近区域,且在滑动过程中气泡的直径增大。
4) 过冷沸腾起始点随着加热功率的增大而提前,气泡的脱离直径随加热功率的增大而增大。质量流速越大,壁温越低,较不容易发生过冷沸腾。
本实验观测的气泡形态可为矩形通道气泡流行判断和聚合破碎模型提供数据支持。
基金项目
中国核动力研究设计院核反应堆系统设计技术重点实验室资助项目(HT-KFKT-02-2020002);广东省基础与应用基础研究青年基金资助项目(2020A1515110512)。
NOTES
*通讯作者。