摇摆对单相自然循环系统流动特性的影响分析

俞胜之，阎昌琪，王建军，郭雪晴

(哈尔滨工程大学 核安全与仿真技术国防重点学科实验室，黑龙江 哈尔滨 150001)

摇摆对单相自然循环系统流动特性的影响分析

俞胜之，阎昌琪，王建军，郭雪晴

(哈尔滨工程大学 核安全与仿真技术国防重点学科实验室，黑龙江 哈尔滨 150001)

针对摇摆对加热段为窄矩形通道的自然循环系统单相流动特性的影响，开展了实验研究。为了进一步分析摇摆对流动特性的影响机理，采用自主开发的系统分析程序，对实验工况下的系统行为进行了数值模拟。结果表明：在实验参数范围内，角加速度幅值是摇摆影响质量流量的本征量。当热工参数相同时，角加速度的幅值越大，摇摆运动下的时均流量与静止时的流量比值越小。基于静止条件下的阻力关系式开发的程序，可以较好地模拟摇摆条件下单相自然循环流动特性。摇摆引起的流量波动是摇摆对单相自然循环流动阻力特性影响的主要机制。

摇摆运动； 单相流动； 自然循环； 窄矩形通道； 流量； 数值模拟

由于受到海洋条件影响，核动力舰船、漂浮式电站等在海洋环境中运行时，会受到起伏、倾斜、摇摆等[1-2]的影响。一般而言，自然循环系统驱动压头小，自然循环能力受冷热源中心有效高度差影响，易受到系统装置运动的影响。研究表明，海洋条件对自然循环系统的热工水力特性有着较大的影响[3-5]。在20世纪90年代，我国已经开展海洋条件对自然循环系统的热工水力特性影响研究；受限于当时的条件，主要进行理论研究[1-2]。近十年来有关海洋条件对热工水力特性的影响研究，主要是在强迫循环条件下进行的[6-15]。因此，研究摇摆运动等海洋条件对自然循环系统热工水力特性的影响具有重要意义。

研究表明，摇摆会使单相自然循环的时均流量下降[2,16]，但是影响的机制还有待研究。摇摆条件下的自然循环系统的热工水力特性数值研究，多集中于计算程序的研发[1-2,4,17]，缺少与本体实验验证分析。因此，有必要开展摇摆对单相自然循环系统流动特性的影响的研究。

1 实验装置

如图1所示，实验装置主要由主循环回路，二次侧循环回路和采集系统组成。其中，主循环回路被安装在3.5 m×2.5 m矩形摇摆平台上，二次侧的主要供水管路通过软管与冷凝器二次侧的进出口进行连接。主循环回路主要由实验段，预热器，稳压器，启动泵，冷凝器等管路及阀门组成。实验段是尺寸为40 mm×2 mm的窄缝通道[13]，通过直流电进行加热。预热器的最大功率为45 kW。稳压器的安全压力为1.3 MPa。摇摆平台可以模拟船舶的简谐运动，其运动规律：

(1)

(2)

(3)

图1 实验装置示意图Fig.1 Scheme of experiment equipment

在摇摆条件下，自然循环的流量不仅与主循环回路自身几何尺寸相关，还与摇摆轴的距离及回路的三维布置密切相关[1]。表1列出了主循环回路与摇摆轴之间的距离。主循环回路的三维布置方式如图2所示。其中，冷凝器轴向与摇摆轴平行。

表1 管路几何尺寸和相对位置

在摇摆运动下的自然循环系统，整个环路会受到附加惯性力作用，流量会发生周期性波动[1-5,14]。研究表明，附加惯性力主要受β影响[1]。

因此，为了研究摇摆剧烈程度对自然循环系统的影响，可利用βmax来体现摇摆剧烈程度：

(4)

表2将实验中的摇摆工况，按照βmax的值进行了升序排列。

2 实验结果与分析

自然循环系统的驱动力较小，易受到摇摆产生的附加惯性力影响而产生流量波动[1-5,14]。在流量波动时，系统各个管路内的流体可能同时历经层流，过渡和湍流3个流态；也可能仅经历1或2个流态，这主要受初始流态(竖直条件)影响。在不同的流动状态下摩擦阻力特性有着较大区别，因此初始流动状态是一个关键影响因素。表3给出了实验工况的热工参数，其中，工况1和工况3初始流态处于过渡区，但工况3入口处于层流，工况2和工况4初始流态处于紊流区。

图2 主循环回路的三维布置图Fig.2 Three-dimensional layout map of main circulation loop

工况角度/(°)周期/sβmax/(rad·s-2)R00000R152000086128R2102000172257R351400175772R4152000258385R551000344514R6202000344514R7101400351545R8151400527317R9101000689028R105600956983R11151001033542

表3 实验工况热工参数

2.1 摇摆对时均流量的影响

如图3给出了不同条件下无量纲时均质量流量W+受不同摇摆参数的影响：

(5)

式中：W(t)表示任意时刻的质量流量，kg/s；W(t0)表示摇摆运动前的质量流量，kg/s。

图3 不同条件下的无量纲时均质量流量Fig.3 Time-average mass flow rate without dimensionless under different cases

由图3易知，在不同热工参数条件下，随着βmax增加W+减小。对比初始状态在湍流区的工况2和工况4与初始状态在过渡区的工况1，W+受βmax影响规律有着较为明显的区别。在工况2和工况4条件下，W+随βmax增加而减小的总体趋势中，存在一个变化减缓的趋势，在其他条件下没有出现相似的现象。这是因为在工况2和工况4条件下，除R10的工况以外，流量波动时所处的主要流态是在湍流与过渡区。这时流量变化引起的阻力系数变化较小，无量纲的波峰与波谷质量流量随βmax的变化率绝对值接近，如图4所示。

图4 不同条件下的无量纲波峰和波谷质量流量Fig.4 Peak and trough of mass flow rate without dimensionless under different cases

对比R2和R3，这两组βmax近似相等的摇摆工况下的W+，可以发现：在各个不同热工参数条件下，W+近似相等。再对比R5、R6和R7，这3组βmax近似相等的摇摆工况下的W+，可以发现在各个不同热工参数条件下，W+也近似相等。

2.2 摇摆对流量波峰和波谷的影响

图4给出了不同条件下的波峰和波谷的无量纲质量流量。此时，W+可计算得

(6)

图4显示，无量纲波峰质量流量随βmax增加而增加，无量纲波谷质量流量随βmax增加而减小。而且，βmax对流量波谷的影响比流量波峰大。引起这种现象的主要原因有两个：

1)流量的周期波动，使得ρu2/2也发生周期性波动。由于密度变化较小，流量与ρu2/2的关系可以近似为平方关系，具有较强的非线性关系。因此，在相等的流量变化条件下，与流量减小相比，流量增加引起ρu2/2的变化值更大。

2)在层流区流量的减小，雷诺数减小，而摩擦阻力系数却迅速增加。这也会使得在相等的流量变化条件下，流量增加引起的摩擦阻力系数增加量，大于流量减小引起的摩擦阻力系数减小量。

当然，在摇摆运动条件下，附加惯性力在垂直于流动方向的法相分力会使流体的速度分布发生改变，与静止状态有所不同[4,19]；同时，流量的波动使得流体的速度分布也会不同于稳态条件下的分布。

与无量纲时均质量流量相似，对比R2和R3，这两组βmax近似相等的摇摆工况下的W+，可以发现：在各个不同热工参数条件下，无量纲波峰与波谷质量流量近似相等。同时，再对比R5、R6和R7，在这3组βmax近似相等的摇摆工况下的W+，也可以发现：在各个不同热工参数条件下，无量纲波峰与波谷质量流量也是近似相等。

因此，摇摆条件下，βmax是影响流量变化的本质因素，无论是时均流量还是流量的波峰与波谷。在研究摇摆参数对流量的影响分析时，可以将摇摆幅度和摇摆周期的影响整合到一个变量中。

3 计算模型

为了进一步研究摇摆对自然循环系统流动特性产生影响的主要因素，基于竖直静止条件下的阻力和传热本构关系式，开发了可用于模拟摇摆条件下单相自然循环特性的系统分析程序。

3.1 简化与假设

1)实验段加热壁面径向划分5个节点，考虑入口效应，忽略轴向导热。

2)实验段出口至冷凝器入口之间绝热段的管路流体温度会发生波动，因此考虑该段管壁的蓄热变化。

3)冷凝器出口至实验段入口之间管路流体温度瞬态变化较小，忽略该段管壁的蓄热变化。

4)采用静止条件的阻力与传热关系式。

3.2 守恒方程

1)质量守恒方程：

(7)

式中：W是质量流量，kg/s；l是沿流动方向的距离，m。

2)动量守恒方程：

(8)

式中：ρ是密度，kg/m3；u是流速，m/s；gpro是重力加速度沿流动方向的分量；apro是附加惯性加速度沿流动方向的分量。

3)能量守恒方程：

(9)

式中：i是比焓，J/kg；S是单位体积流体与加热壁面或热构件的热交换量，J/m3；

3.3 本构关系

1)附加惯性加速度

根据非惯性系下的流体受力分析，可以得到与摇摆引起的附加惯性加速度为

aroll=-ω×(ω×r)-β×r+2(ω×u)

(10)

apro=aroll·l

(11)

式中：r是回路某一位置在远动坐标系下的位置向量；aroll是摇摆运动引起的惯性加速度，l是流动方向在运动坐标系下的方向向量。

2)非惯性系下重力加速度

摇摆运动下，自然循环回路会产生倾斜。此时，坐标系建立在摇摆平台上，重力在该坐标系下，是一个变化的矢量，可表示为

(12)

(13)

式中：g是运动坐标系下的重力加速度。

(14)

式中：λ是摩擦阻力系数，Re是雷诺数。

4)传热计算

对于矩形通道的传热计算，采用文献[18]推荐的计算公式。而圆管层流区的传热计算采用等壁温边界条件处理。

4 计算结果与分析

4.1 与实验结果对比

1)流量。

图5～6给出了两组不同摇摆工况下的计算结果与实验结果对比。

图5 工况4和R9条件下的体积流量Fig.5 Volume of flow rate under case 4 and R9 conditions

图6 工况4和R1条件下的体积流量Fig.6 Volume of flow rate under case 4 and R1 conditions

结果显示：基于静止条件下的阻力计算关系开发的计算程序，能够较好地模拟摇摆条件下的单相自然循环系统的流量波动。

为了进一步验证计算结果的准确性，将矩形通道摩擦压降计算结果与实验结果进行对比。

2)实验段摩擦压降。

如图7～8所示，计算结果与实验结果符合较好。仅摩擦压降的平均值有些许差别，振幅和相位周期均能吻合较好。

4.2 摇摆对时均流量影响机理分析

在实验结果分析中，提及摇摆影响单相自然循环时均流量的因素有4个。这4个因素可以根据其影响机理分为两类：

1)摇摆引起的流量波动使管道流速发生改变，流速的变化会引起阻力的变化。因为流速与阻力的有较强的非线性相关性，会使摇摆运动条件下系统的平均阻力大于静止条件下的系统阻力。

2)摇摆引起的流量波动和附加惯性力会使流体流速分布发生变化，进而影响阻力特性。

在本文开发的程序中，仅考虑了第一类影响因素。因此，从与实验结果对比可以得出：第一类因素是摇摆对单相自然循环流动特性的主要影响机制。

图7 工况4和R9条件下的摩擦压降Fig.7 Friction pressure drop under case 4 and R9 conditions

图8 工况4和R1条件下的摩擦压降Fig.8 Friction pressure drop under case 4 and R1 conditions

为了更具有论证性，将工况4和R9工况下的系统驱动力，环路摩擦压降，环路局部阻力平均值与静止条件相对比，如图9～11所示。为了将摇摆前后的数据尽量在一个图中展示，同时保证摇摆稳定后的数据细节全部展示，将摇摆启动过渡到稳定周期波动的数据隐藏，在横坐标轴设置了个断点，将1 440～2 440 s期间数据隐藏。

图9 工况4和R9条件下的系统驱动力Fig.9 The driving head of system under case 4 and R9 conditions

图10 工况4和R9条件下的环路摩擦压降Fig.10 Friction pressure drop of system under case 4 and R9 conditions

图11 工况4和R9条件下的环路局部压降Fig.11 Local pressure drop of system under case 4 and R9 conditions

在工况4和R9条件下，流态主要是处于过渡区。在过渡区，流量降低，雷诺数减小，摩擦阻力系数下降，因此摩擦压降平均值是减小的。而局部压降平均值增加则是因为ρu2/2非线性变化引起的。

5 结论

1)在实验参数范围内，βmax是摇摆影响流量变化的本征因素；

2)实验参数范围内，热工参数相同时，βmax越大，时均流量与静止时的流量比值越小；

3)基于静止条件下的阻力计算关系式开发的系统分析程序，可以较好地模拟摇摆条件下单相自然循环流动特性。

4)摇摆运动下单相自然循环的流动阻力特性变化，主要是流量波动引起的。因此，在摇摆运动下的阻力计算，可以采用静止条件下的阻力关系式。

[1]高璞珍，庞凤阁，王兆祥. 核动力装置一回路冷却剂受海洋条件影响的数学模型[J]. 哈尔滨工程大学学报， 1997， 18(1)： 24-27. GAO Puzhen, PANG Fengge, WANG Zhaoxiang. Mathematical model of primary coolant in nuclear power plant influenced by ocean conditions[J]. Journal of Harbin Engineering University, 1997, 18(1): 24-27.

[2]苏光辉，张金玲，郭玉君，等. 海洋条件对船用核动力堆余热排出系统特性的影响[J]. 原子能科学技术， 1996， 30(6)： 487-491. SU Guanghui, ZHANG Jinling, GUO Yujun, et al. Effects of ocean conditions upon the passive residual heat removal system of ship reactor[J]. Atomic energy science and technology, 1996, 30(6): 487-491.

[3]宫厚军，杨星团，黄彦平，等.倾斜条件下一体化反应堆模拟回路单相自然循环实验与数值研究[J]. 核动力工程， 2014， 35(5)： 89-93. GONG Houjun, YANG Xingtuan, HUANG Yanping, et al. Experimental and numerical study on natural circulation of integrated reactor under inclination conditions[J]. Nuclear power engineering, 2014, 35(5): 89-93.

[4]李晓伟，杨明昆，宋宇，等. 摇摆对自然循环系统流量及载热能力影响分析[J]. 工程热物理学报， 2015， 36(2)： 414-418. LI Xiaowei, YANG Mingkun, SONG Yu, et al. Analysis of rolling effects on the flow and heat transport capacity of a natural circulation system[J]. Journal of engineering thermophysics, 2015, 36(2): 414-418.

[5]封贝贝，杨星团，姜胜耀. 起伏强度和周期对自然循环流动特性的影响[J]. 哈尔滨工程大学学报， 2015， 36(4)： 448-451. FENG Beibei, YANG Xingtuan, JIANG Shengyao. Effect of heaving strength and cycle on natural circulation flow characteristics[J]. Journal of Harbin Engineering University, 2015, 36(4): 448-451.

[6]TAN Sichao, WANG Zhanwei, WANG Chang, et al. Flow fluctuations and flow friction characteristics of vertical narrow rectangular channel under rolling motion conditions[J]. Experimental thermal and fluid science, 2013, 50: 69-78.

[7]YU Zhiting, TAN Sichao, YUAN Hongsheng, et al. Experimental investigation on flow instability of forced circulation in a mini-rectangular channel under rolling motion[J]. International journal of heat and mass transfer, 2016, 92: 732-743.

[8]WANG Chang, LI Xiaohui, WANG Hao, et al. Experimental study on single-phase heat transfer of natural circulation in circular pipe under rolling motion condition[J]. Nuclear engineering and design, 2014, 273: 497-504.

[9]WANG Chang, WANG Shaowu, WANG Hao, et al. Investigation of flow pulsation characteristic in single-phase forced circulation under rolling motion[J]. Annals of nuclear energy, 2014, 64: 50-56.

[10]XING Dianchuan, YAN Changqi, SUN Licheng. Effect of rolling motion on single-phase laminar flow resistance of forced circulation with different pump head[J]. Annals of nuclear energy, 2013, 54, 141-148.

[11]XING Dianchuan, YAN Changqi, SUN Licheng. Flow fluctuation behaviors of single-phase forced circulation under rolling conditions[J]. Ocean engineering, 2014, 82: 115-122.

[12]TIAN Daogui, YAN Changqi, SUN Licheng. Model of bubble velocity vector measurement in upward and downward bubbly two-phase flows using a four-sensor optical probe[J]. Progress in nuclear energy, 2015, 78: 110-120.

[13]LI Shaodan, TAN Sichao, XU Chao, et al. Visualization study of bubble behavior in a subcooled flow boiling channel under rolling motion[J]. Annals of nuclear energy, 2015, 76: 390-400.

[14]曹夏昕，阎昌琪，孙中宁. 摇摆状态下竖直管内弹状流单元物理模型解析 [J]. 核动力工程， 2007， 28(3)： 47-50． CAO Xiaxin， YAN Changqi， SUN Zhongning. Analysis of slug flow unit model in rolling vertical tube[J]. Nuclear power engineering, 2007, 28(3): 47-50．

[15]阎昌琪, 金光远, 孙立成,等. 摇摆状态下矩形通道内弹状流压力模型[J]. 原子能科学技术, 2015， 49(6): 1062-1068. YAN Changqi, JIN Guangyuan, SUN Licheng, et al. Pressure model of slug flow in rectangular channel in rolling motion[J]. Atomic energy science and technology, 2015, 49(6): 1062-1068.

[16]谭思超, 高璞珍, 苏光辉. 摇摆运动条件下自然循环流动的实验和理论研究[J]. 哈尔滨工程大学学报， 2007， 28(11)： 1213-1217. TAN Sichao, GAO Puzhen, SU Guanghui. Experimental and theoretical study on natural circulation flow under rolling motion condition[J]. Journal of Harbin Engineering University, 2007, 28(11): 1213-1217.

[17]YAN Binghuo, YU Lei, LI Hongquan. Research on operational characteristics of passive residual heat removal system under rolling motion[J]. Nuclear engineering and design, 2009, 239: 2302-2310.

[18]马建，黄彦平，刘晓钟. 窄缝矩形通道单相流动及传热实验研究[J]. 核动力工程， 2012， 33(1)： 39-45. MA Jian, HUANG Yanping, LIU Xiaozhong. Experimental study on single-phase flow and heat transfer in a narrow rectangular duct[J]. Nuclear power engineering, 2012, 33(1): 39-45.

[19]XING Dianchuan, YAN Changqi, WANG Chang, et al. A theoretical analysis about the effect of aspect ratio on single-phase laminar flow in rectangular ducts[J]. Progress in nuclear energy, 2013, 65: 1-7.

本文引用格式：

俞胜之，阎昌琪，王建军，等. 摇摆对单相自然循环系统流动特性的影响分析[J]. 哈尔滨工程大学学报， 2017， 38(7): 1065 -1071.

YU Shengzhi, YAN Changqi, WANG Jianjun, et al. Analysis of the effect of rolling motion on flow characteristics of a single-phase natural circulation system[J]. Journal of Harbin Engineering University， 2017， 38(7): 1065 -1071.

Analysis of the effect of rolling motion on flow characteristics of a single-phase natural circulation system

YU Shengzhi, YAN Changqi, WANG Jianjun, GUO Xueqing

(Fundamental Science on Nuclear Safety and Simulation Technology Laboratory, Harbin Engineering University, Harbin 150001, China)

Experiments were conducted to investigate the effect of rolling motion on single-phase flow characteristics in a natural circulation system with heated narrow rectangular channel. Using the self-developed code, the system′s behaviors were numerically simulated in order to further understand the mechanism of effect of rolling motion on flow characteristics. Results indicate that, within the range of experimental parameters, the amplitude of angular acceleration is the eigenvalue of the influence on mass flow rate by rolling motion. For the given thermal parameter combination, the ratio between the mass flow rate under static condition and the average mass flow rate under rolling motion decreases with the increase of the angular acceleration amplitude. According to the calculation results, under static condition, the single-phase natural circulation flow characteristics can be preferably simulated by the code, which has been developed based on the resistance relationship. Furthermore, fluctuant flux has been identified as the main mechanism of the effect of rolling motion on single-phase natural circulation flow resistance characteristics.

rolling motion; single phase flow; natural circulation; narrow rectangular channel; flow rate; numerical simulation

2016-05-17.

日期：2017-04-26.

国家自然科学基金项目(11675045).

俞胜之(1990-), 男, 博士研究生； 阎昌琪(1955-), 男, 教授, 博士生导师.

阎昌琪,E-mail：changqi_yan@163.com.

10.11990/jheu.201605055

TL334

A

1006-7043(2017)07-1065-07

网络出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20170426.1802.080.html