蒸汽发生器致畸变入流对核主泵流动性能的影响
Influence of distorted inflow caused by steam generator on flow properties of reactor coolant pump
通讯作者:
收稿日期: 2018-09-11
Received: 2018-09-11
作者简介 About authors
王悦荟(1992—),女,博士生,从事微流体及流体机械研究.orcid.org/0000-0003-0956-8551.E-mail:
为了探究由蒸汽发生器(SG)引起的畸变入流对核主泵(RCP)流动性能的影响,采用计算流体力学方法,建立核主泵与下封头的联合计算模型;利用多参考坐标系模型和滑移交界面技术,开展运行工况下的全三维稳态和瞬态数值模拟;采用入流畸变度和平均偏流角公式定量表征畸变入流的流场速度分布,通过正则化螺旋度法捕捉核主泵叶轮的涡流流动特征. 研究结果表明:与均匀入流相比,畸变入流复杂化了核主泵进口流场,导致核主泵的湍动能和湍耗散增大,降低了主泵的水力效率;在泵的进口流场形成了明显的非对称分布的涡核区域,引起叶轮各流道流态产生差异,造成叶轮内部压力和速度分布不均;导致各流道流量分配不均,加剧叶轮受载波动;降低其运行的稳定性和安全性.
关键词:
A united model for the reactor coolant pump (RCP) and the channel head was generated by using the method of computational fluid dynamic, in order to analyze the influence of the distorted inflow caused by the steam generator (SG) on the flow property of RCP. Three-dimensional steady and transient numerical simulations in normal operating conditions were performed by using multiple frames of reference and sliding interface technology. The velocity distribution of distorted inflow was quantitatively characterized by inflow distortion degree and mean drift angle, and the regularization helicity method was used to capture the flow characteristics of vortex core for the impeller of RCP. Results show that with the comparison of uniform inflow, the inflow field of the RCP is complicated by the distorted inflow, and the turbulent kinetic energy and turbulent eddy dissipation are increased while the hydraulic efficiency of the RCP is decreased. Obvious asymmetric vortex core regions form in the RCP inlet, causing differences between the flow patterns of the impeller flow channels, and then causing the non-uniform distribution of the pressure and velocity within the impeller. Furthermore, due to the distorted inflow, the flow distribution of each flow channel is uneven, increasing the load fluctuation. The distorted inflow has an adverse effect on the RCP, reducing the stability and safety of the operation of RCP.
Keywords:
本文引用格式
王悦荟, 刘聪, 王鹏飞, 许忠斌, 阮晓东, 付新.
WANG Yue-hui, LIU Cong, WANG Peng-fei, XU Zhong-bin, RUAN Xiao-dong, FU Xin.
通常在泵的设计中均假定入流均匀,但是在工程实际中,受泵的安装方式和结构限制[3-4],入流畸变是泵领域普遍存在的现象. 目前广泛研究的关于泵的入流畸变主要是由进水弯管引起的,其引起的畸变与管路系统相耦合,继而加剧了振动、噪声等问题[5]. 与弯管相比,第3代先进非能动压水堆(如AP1000)的核主泵的入流结构更复杂,核主泵的入口管被取消,泵直接与上游的蒸汽发生器相连. 由于缺乏过渡管道,从蒸汽发生器流出的高温高压流体未经整流而直接倒灌进核主泵,核主泵的入流不再是设计时假定的均匀入流,而是非均匀入流,即畸变入流. 这种由蒸汽发生器下封头(channel head of steam generator,SGCH)引起的畸变入流使得入口流场情况更加复杂[6],半圆形下封头到短管道的结构一方面导致流道转变产生离心力,另一方面导致出流截面突缩引起脱流,两者叠加导致畸变程度加强.
综上所述,由于蒸汽发生器出流的影响,第3代先进非能动压水堆的核主泵入流存在明显的流场畸变,且在长期使役周期内始终存在. 为了避免核主泵在长期使役过程中发生由畸变入流导致的叶片疲劳等问题,对核主泵在由蒸汽发生器引起畸变入流的情况下的泵内流场压力、速度分布情况、湍动能、涡演变情况等采用计算流体力学(computational fluid dynamic,CFD)仿真分析进行研究,为核主泵的设计提供参考.
1. 数值计算方法及验证
1.1. 计算模型
图 1
图 1 核主泵与蒸汽发生器下封头的位置关系图
Fig.1 Position diagram of reactor coolant pump and channel head of steam generator
图 2
1.2. 网格划分与边界条件
采用ICEM软件对各部件网格进行独立划分,其中对叶轮水体和导叶水体采用结构网格划分,对进口段、环形压水室水体和出水段采用非结构四面体网格划分,网格划分如图3所示.
图 3
图 3 下封头与核主泵联合模型网格划分图
Fig.3 Computational domain grids of united model for RCP and SGCH
在建模完成后,将核主泵的扬程作为检验的主要指标,对网格进行无关性验证,如图4所示. 图中,H为核主泵扬程,N为网格数. 当网格数达到535万后,所计算的核主泵扬程趋于稳定,网格达到要求. 同时,对模型网格进行局部细化和全局光顺化处理,使得网格质量满足计算要求,所采用的网格质量高于0.2,网格的最小角大于18°,以保证核主泵模拟的准确性.
图 4
图 4 下封头与核主泵联合模型的网格无关性验证
Fig.4 Grid independence verification of united model for RCP and SGCH
标准k-ε模型为两方程模型,其中湍动能输运方程由理论推导求得,耗散率方程由物理推理和数学模拟相似原型方程求得,表达式分别为
式中:ρ为密度;xi为坐标位置;ui为坐标位置对应方向的速度;p为压力;
1.3. 试验验证
采用课题组自主设计的主泵水力模型进行仿真模型,其主要设计参数如表1所示. 试验用泵为设计泵的缩尺模型(1∶2.5),通过相似换算得到模型泵的测试值[19].采用相似理论对原型泵的参数进行变换,然后加工模型泵开展试验[20-21]. 试验利用第三方单位的水力试验台进行测量,试验台主要包括控制台、水箱、缩尺模型泵、电机以及管道仪表等. 在缩尺泵进出口设置压力变送传感器,检测压差即可得到泵的扬程;在缩尺泵的出口设置流量计,用于测量泵的工作体积流量. 泵的转速由扭矩仪测量得到,试验回路如图5所示. 对比体积流量qV和扬程的曲线,如图6所示. 图中,Hd为标准体积流量下的扬程,
表 1 核主泵叶轮几何参数
Tab.1
设计参数 | 数值 |
流量/(m3∙h−1) | 17 886 |
扬程/m | >111.3 |
同步转速/(r·min−1) | 1 800 |
设计压力/MPa | 17.2 |
设计温度/K | 616 |
比转速 | 428 |
图 5
图 6
图 6 核主泵体积流量-扬程的模拟和实验数据对比图
Fig.6 Comparison diagram of volume flow-head of RCP between numerical and experimental results
2. 泵进口流场的畸变入流
式中:e为入流畸变度;ω为平均偏流角;
图 7
图 7 核主泵入口处不同剖面的位置示意图
Fig.7 Diagram of location of different sections of RCP inlet
图 8
图 8 入流畸变度和平均偏流角的变化规律
Fig.8 Variation of inflow distortion degree and mean drift angle
3. 畸变入流对主泵流场的影响
模拟标准工况下,在均匀入流与畸变入流2种工况下,核主泵性能的对比与差异.
3.1. 核主泵的叶轮进口涡结构辨识
正则化螺旋度的表达式为
式中:Hn为正则化螺旋度,v为速度矢量,w为涡量矢量.
如图9、10所示,分别为在均匀和畸变入流下,核主泵进口流场一个周期T内的正则化螺旋度分布图,图中截取的平面为如图7所示的Lb1、Lb2、Lb3面. 可以看出,在主泵进口处,由蒸汽发生器引起的畸变入流,在进口管内形成明显的3块旋向不同的涡核区域,涡核的存在影响了叶轮内的流场,在叶轮内产生了分布不均的涡核区域. 这说明畸变入流对主泵流场产生了不利影响,叶轮流道内流场的对称性被打破,从而降低了叶轮的稳定性. 正则化螺旋度法的优点是算法强健,兼顾计算精度和计算资源,可以较准确地捕捉到旋转叶片流体机械的涡流流动特征. 然而,二维结构的涡流特征只能反映各垂直轴平面上的涡的种类,无法体现各涡流流动对流体机械内部流动的影响.
图 9
图 9 均匀入流下核主泵叶轮进口流场的正则化螺旋度分布图
Fig.9 Regularization helicity distribution of impeller inlet flow field of RCP in uniform inflow condition
图 10
图 10 畸变入流下核主泵叶轮进口流场正则化螺旋度分布图
Fig.10 Regularization helicity distribution of impeller inlet flow field of RCP in distorted inflow condition
3.2. 核主泵的叶轮水力特性分析
图 11
图 12
图 13
图 14
在运行工况下,叶轮旋转对流体做功,流体在叶轮流道内的速度逐渐增大,在叶轮出口边侧速度达到最大. 在均匀入流工况下,由于进口流场的对称性,核主泵叶轮的各流道速度分布相对均匀. 在畸变入流工况下,来流以不同的速度流入叶轮各流道,使得各流道的入流速度产生差异. 如图15所示,在主泵叶轮内流道速度流场的差异体现在叶轮各流道的流体体积流量分配中. 图中,B1~B7为核主泵叶轮的7个叶片.
图 15
图 15 核主泵叶轮各流道体积流量分配图
Fig.15 Volume flow distribution diagram for each flow passage of impeller of RCP
与均匀入流工况不同,在畸变入流工况下,各流道内叶轮体积流量分配明显不均衡,不符合原本设计时假定的叶轮受力的对称性,极易造成运行过程中泵的受载不均以及振动、噪声增强等问题.
3.3. 核主泵的湍动能分析
图 16
图 17
4. 结 论
(1)由蒸汽发生器引起的畸变入流使得核泵的入口流场工况复杂化,造成叶轮内部各流道内压力和速度的差异,从而导致各流道体积流量分配不均,最终引起叶轮受载不均,降低主泵叶轮受载的稳定性和安全性.
(2)采用正则化螺旋度法对主泵进口段和叶轮的流场进行分析,可以看出,与均匀入流相比,由蒸汽发生器引起的畸变入流在核主泵进口流场形成了明显的涡核区域,导致叶轮各流道产生不对称的、甚至旋向相反的涡核区域. 这说明畸变入流会对主泵流场产生不利的影响,打破叶轮流道内流场的对称性,降低叶轮的稳定性. 对于涡内部的演化过程,还须进行更高精度的计算以得到三维涡结构,从而进一步分析畸变入流产生的涡在主泵内部流场的演化规律.
(3)由蒸汽发生器引起的畸变入流,会引起局部湍动能、湍流耗散高值,增大核主泵湍动能和湍耗散,造成核主泵的扬程下降. 若要深入了解畸变入流下的主泵性能变化,须进行畸变入流下的模型主泵试验,以揭示由蒸汽发生器所引起的核主泵畸变入流的形成机制.
参考文献
浅谈压水堆核电站主泵
[J].
Discussion on the main pump of PWR nuclear power plant
[J].
AP1000核主泵叶轮典型质量问题及监造预防措施
[J].
Typical quality problems and preventive measures for AP1000 nuclear main pump impeller
[J].
大型立式轴流泵叶片进口流场及其对水泵影响研究
[J].DOI:10.3321/j.issn:0577-6686.2005.04.006 [本文引用: 1]
Study on flow field in blade inlet of large vertical axial-flow pump and its influence on pump
[J].DOI:10.3321/j.issn:0577-6686.2005.04.006 [本文引用: 1]
进口弯管及前置扰流子对离心泵性能影响
[J].DOI:10.3969/j.issn.1002-6819.2013.20.011 [本文引用: 1]
Effect of import bend and forward turbolator on performance of centrifugal pump
[J].DOI:10.3969/j.issn.1002-6819.2013.20.011 [本文引用: 1]
Study on the characteristic of self-excited cavitation oscillation in pump-piping system
[J].DOI:10.4028/www.scientific.net/AMR.268-270.148 [本文引用: 1]
AP1000核岛主泵流场数值模拟
[J].
Simulation of AP1000 reactor coolant pump flow field CFD
[J].
Research on the non-uniform inflow characteristics of the canned nuclear coolant pump
[J].DOI:10.1016/j.anucene.2018.02.007
蒸汽发生器下封头对核主泵入口流场影响
[J].DOI:10.3969/j.issn.1674-8530.15.0207 [本文引用: 2]
Influence of steam generator channel head on reactor coolant pump inflow field
[J].DOI:10.3969/j.issn.1674-8530.15.0207 [本文引用: 2]
基于额定参数的核主泵惰转工况计算模型
[J].DOI:10.7538/yzk.2014.48.08.1435 [本文引用: 1]
Coast-down model based on rated parameters of reactor coolant pump
[J].DOI:10.7538/yzk.2014.48.08.1435 [本文引用: 1]
蒸汽发生器换热管流量分配及其对核主泵入口流场的影响
[J].
Flow distribution of heat exchanger tubes in a steam generator and its effect on flow field at entrance of reactor coolant pump
[J].
Westinghouse AP1000 advanced passive plant
[J].
湍流模型对混流泵性能预测的影响
[J].DOI:10.6041/j.issn.1000-1298.2013.11.008 [本文引用: 1]
Influence of turbulence model on performance prediction of mixed-flow pump
[J].DOI:10.6041/j.issn.1000-1298.2013.11.008 [本文引用: 1]
不同湍流模型在混流泵数值模拟中的应用与评价
[J].DOI:10.11662/j.issn.1673-5196.2014.05.012 [本文引用: 1]
Application of different turbulence models in numerical simulation of mixed-flow pump and its evaluation
[J].DOI:10.11662/j.issn.1673-5196.2014.05.012 [本文引用: 1]
Numerical investigation on the " S” characteristics of a reduced pump turbine model
[J].DOI:10.1007/s11431-011-4295-2 [本文引用: 1]
导流件和叶轮强作用涡轮流量计的CFD仿真方法
[J].DOI:10.3969/j.issn.1000-3932.2013.10.019 [本文引用: 1]
CFD simulation method for turbine flowmeter with strong action from guide vane and impellor
[J].DOI:10.3969/j.issn.1000-3932.2013.10.019 [本文引用: 1]
Unsteady flow structure and global variables in a centrifugal pump
[J].DOI:10.1115/1.2234782 [本文引用: 1]
核主泵启动过程压力脉动和径向力研究
[J].
Study on pressure fluctuation and radial force during startup of reactor coolant pump
[J].
核主泵升速过程最高转速设定方法研究
[J].
Study on maximum speed setting standard during speed up process of reactor coolant pump
[J].
大型泵站进水流道技术改造优选设计
[J].DOI:10.3969/j.issn.1003-1243.2006.02.011 [本文引用: 1]
Numerical simulation and experimental study on bell-shaped suction box under non-symmetric inflow conditions
[J].DOI:10.3969/j.issn.1003-1243.2006.02.011 [本文引用: 1]
非对称入流工况下钟形进水流道数值模拟试验研究
[J].
Numerical simulation and experimental study on bell-shaped suction box under non-symmetric inflow conditions
[J].
Graphical visualization of vortical flows by means of helicity
[J].
/
〈 |
|
〉 |
