湿蒸汽流场下介入式探针振动数值模拟
Numerical simulation for vibration of interferometric probein wet steam flow field
通讯作者:
收稿日期: 2018-04-10
Received: 2018-04-10
作者简介 About authors
胡展豪(1995—),男,硕士生,从事两相流光学测量技术研究.orcid.org/0000-0002-8812-4635.E-mail:
介入式探针在汽轮机末级与湿蒸汽流场相互影响,采用双向流固耦合技术对该三元流动问题进行数值分析,研究探针在湿蒸汽流场下的受力情况和振动特性. 利用Workbench平台搭建流固耦合模型,采用基于湍流剪应力输运的k-ω模型以提高计算精度. 计算结果表明:探针在X、Y、Z方向的振动周期略有波动,最大振幅分别为0.28 mm、4.75和5.50 μm. 由于光学探针结构特殊,探针前端测量区域受到压差作用力很小,前端的形变量增长率有所降低. 在振动过程中,探针背流方向约2/3探针长度处应力值最大,该值远低于探针材料的许用应力,且在计算工况下不会引发共振效应.
关键词:
The bidirectional fluid solid coupling technique was adopted to analyze the three-element flow problem that the interferometric probe interacts with the wet steam flow field at the end stage of the steam turbine. The vibration characteristics and stress of the probe under wet steam flow field were studied. The fluid solid coupling model was built by using Workbench platform. The k-ω model based on turbulent shear stress transport was used to improve calculation accuracy. Results show that the vibration period of the three directions X, Y and Z of probe fluctuate slightly. The oscillation amplitudes in directions X, Y and Z were 0.28 mm, 4.75 μm and 5.50 μm, respectively. Due to the special structure of the optical probe, the pressure differential force acting on the tip of the probe is very small, leading to a decrease in the growth rate of the shape variable at the front end. During probe vibration, the maximum stress position of the probe is about 2/3 of the point in the direction of back flow of the probe, the maximum stress is far below the allowable stress of the probe material, and the resonance effect is not induced in the calculation condition.
Keywords:
本文引用格式
胡展豪, 冯俊涛, 盛德仁, 陈坚红, 李蔚.
HU Zhan-hao, FENG Jun-tao, SHENG De-ren, CHEN Jian-hong, LI Wei.
核电站大型汽轮机组的全部级和大型凝汽式汽轮机的末几级在湿蒸汽下工作,湿蒸汽对机组的安全性和经济性有重要影响. 一方面,湿蒸汽的存在会降低汽轮机的热效率;另一方面,湿蒸汽流过叶片时形成的二次水滴易对汽轮机造成严重水蚀破坏,甚至引发叶片断裂事故. 为了解湿蒸汽的形成机理、性质特点、危害影响,国内外学者进行了不少理论和实验研究,主要集中在湿蒸汽两相流模拟和探针实验两方面.
White等[1-3]发展了二维叶栅非平衡凝结流动的计算理论. 韩中合等[4]提出了一种基于传热传质平衡耦合求解方法与水滴生长的修正模型. 余兴刚等[5]研究了1 000 MW火电汽轮机叶片表面粗糙度对叶栅通道内蒸汽自发凝结过程和级的性能的影响. 在实验方面,基于光学原理的湿蒸汽测量方法发展迅速,主要有角散射法、消光法、图像法、光脉动法等[6-11]. 光学法测量需要借助探针系统将光敏元件布置到汽轮机流道内,国内外近年来研制出不少新型探针并投入到实验中[12-14]. 汽轮机低压流道内湿蒸汽进行着复杂的三元凝结流动,细长的探针相对于狭小的动静叶间隙来说具有很大的几何尺度,如1 000 MW汽轮机末级叶片高度约1.3 m. 探针经过低压缸的外缸与内缸进入汽轮机末级叶片后的测量区域,长达2.2 m,以悬臂梁的形式一端固定在外缸上. 杨荣菲等[15-17]评估测试了探针支杆对压气机气动性能的影响,得出支杆直径越大,压气机特性线偏移量越大;但其主要着眼于气动探针对压气机流场及转子性能的影响,未考虑流场对探针的冲击以及两者的双向耦合作用,并且汽轮机末级测量探针长度远大于压气机测量探针(0.1 m),其在流场中的振荡与耦合作用远比气动探针激烈.
测量探针在汽轮机末级和湿蒸汽流场相互影响. 一方面,探针受到高速汽流的冲击,在压差、水滴、汽流激波等因素的作用下发生力与热的交换产生形变和振动,探针的强度、刚度与振动频率等对测量的安全性有重要影响;另一方面,探针的介入会引发两相高速流体的钝体绕流现象,进而产生边界层分离、漩涡脱落、相变等现象,导致压力场、温度场、湿度场等变化. 探针振动产生的位移又会对流场产生进一步的作用,探针介入对流场测量精度的影响也有待研究.
本文通过建立湿蒸汽探针的流固耦合计算模型,探究在汽轮机末级湿蒸汽流场中探针与流场双向流固耦合作用下探针的振动规律和受力情况.
1. 数值模拟方法
1.1. 研究方案
通过基于双向流固耦合的数值模拟方法,研究湿蒸汽流场对探针的影响,数值计算流程如图1所示. 根据实际中的探针外形结构和尺寸数据,进行合理的简化,选取探针附近一定区域为模拟流场区域,在Inventor平台上建立固体和流场的物理模型. 将物理模型导入到ANSYS Workbench中拆分为固体部分和流体域,分别在该软件自带的Meshing模块和专业流体网格划分软件ICEM CFD中进行网格划分. 在CFX进行动网格设置和求解器设置,固体域与流体域之间反复迭代,直到结构和流体域迭代计算收敛为止,然后再进行下一时间的耦合计算,保证流体和固体域能实现数据交换,达到流固耦合分析求解的目的.
图 1
图 1 双向流固耦合数值计算流程图
Fig.1 Flow chart of bidirectional flow solid couplingnumerical calculation
1.2. 计算模型选择
选用Euler/Euler两相流动模型作为湿蒸汽两相流动的数值计算基础,建立并求解欧拉坐标系中各变量的控制方程以充分实现汽液两相耦合. 流体控制方程由矢量形式的可压缩Navier-Stokes方程给定,离散格式采用二阶精度以确保计算精度.
凝结相质量份额β的输运方程如下:
式中:ρ为湿蒸汽密度,v为湿蒸汽在z方向的分速度,Γ为凝结核蒸发作用产生的质量生成率.
每单位容积液滴密度变化数量δ的输运方程如下:
式中:I为凝结成核率.
在推导湿蒸汽相变模型时,根据湿蒸汽的相变特性和液滴经典成核理论,考虑凝结过程中存在潜热形式的能量转移,加上等温修正,得到凝结过程的成核率公式:
式中:qc为凝结系数,θ为非等温影响修正系数,ρν为蒸汽密度,ρl为液相密度,σ为表面张力,m为单分子的质量,r为开尔文−赫姆霍兹临界液滴半径,kB为波尔兹曼常量,Tν为蒸汽热力学温度.
湍流模型的选择综合考虑了计算要求和计算资源的可靠性,对比LES模型、DES模型等模型后,选择基于湍流剪应力输运(SST)的k-ω模型:
式中:ηt为流体动力黏度,k为湍流动能并用湍流脉动速度的变化来定义,ω为湍流脉动频率,νt为流体运动黏度,α1为常数,S为应变率的常量,F为将限幅器限定在壁面边界层上的混合方程.
1.3. 物理模型及网格
图 2
图 3
图 4
图 4 局部网格加密处理的流场结构化网格
Fig.4 Flow field structured meshes with local meshes densification
1.4. 边界与初始条件
流场域由7个边界面包围,将上、下壁面设置为无滑移型边界条件,流体域两侧设置为对称性周期边界条件,探针表面设置为流固耦合交界面,进口边界面给定汽轮机末级额定工况,蒸汽进口总温为308.15 K,干度为0.9,平均流速为200 m/s,出口边界面给定平均绝对静压为5 kPa.
2. 计算结果分析
2.1. 流场稳态数值计算
图 5
图 5 稳态计算探针表面的应力场
Fig.5 Stress field on surface of probe on steady statecalculation
图 6
图 6 稳态计算探针附近流体的速度场
Fig.6 Velocity field of fluid near probe on steady state calculation
2.2. 网格无关性验证
在Block分块的基础上,采用4种网格方案,分别进行稳态数值计算,网格节点数分别为4×105,8×105,1.3×106和2.1×106,得到的探针表面最大应力如图7所示,其中,σm为探针表面最大应力,nw为网格节点数. 由图可得,计算的网格达到1.3×106时,已基本保证了网格的无关性,因此在瞬态数值计算过程中,采用网格节点数为1 338 965的网格方案.
图 7
图 7 探针表面最大压力与网格节点数的关系
Fig.7 Relationship between maximum pressure on probesurface and number of mesh nodes
2.3. 理论验证
由于探针外形独特,且模拟流场环境位于汽轮机末级内部,目前没有与本文一致的相关实验结果. 但垂直探针结构上除了前端消光法测量区域外,可近似视为圆柱体,此部分流固耦合与经典流体力学中流体绕流圆柱模型类似,可通过其漩涡脱落理论验证计算结果.
图 8
图 8 雷诺数约为30 000时的湍流尾迹
Fig.8 Turbulent wake form with Reynolds number ofabout 30 000
图 9
图 9 稳态计算得到的探针周围流场流线图
Fig.9 Flow field diagram around probe of steady statecalculation
2.4. 探针振动分析
2.4.1. 振动形态分析
以稳态计算结果为初场,设置探针在流场中的所有表面为流固耦合面,建立流场与固体间的信息交换,设置动网格模型,取时间步长为0.000 1 s,进行双向流固耦合计算. 数值计算过程中在探针前端设置监测点监测固体振动情况,监测点位置如图10所示,其中X方向为正对来流方向.
图 10
图 10 用于监测探针振动情况的监测点位置
Fig.10 Location of monitoring point used for monitor probe vibration
在整个计算过程中,监测点的位移-时间曲线如图11所示,其中ΔL为监测点的位移,t为监测时间. 可以看出,探针3个方向的振动周期略有波动,约为0.018 s. 计算结果显示X、Y、Z三个方向的最大振幅分别为0.28 mm、4.75和5.50 μm. 其中X方向由于正对来流,承受大部分流体冲击,其振幅远远大于Y和Z方向的振幅,且振幅波动较大. X方向由于受到流体冲击其振动基线偏离0点,Y和Z方向振动基线比较接近0点.
图 11
图 11 探针在3个方向上的位移-时间曲线
Fig.11 Displacement-time curve of probe in three directions
图 12
图 13
图 13 探针中部截面周围流场压力云图
Fig.13 Pressure cloud picture of flow field around centralsection of probe
图 14
图 14 探针测量部分周围流场压力云图
Fig.14 Pressure cloud picture of flow field around measurement section of probe
2.4.2. 探针振动安全性分析
图 15
图 15 探针应力最大时在XZ平面的应力云图
Fig.15 Probe’s stress cloud picture on XZ plane at maximum stress
图 16
图 16 探针应力最大位置在一个振动周期内的应力
Fig.16 Stress of probe’s maximum stress position within one vibration period
对研究的探针而言,应力最大位置在振动过程中基本保持不变,实际测量时应注意此处是否为结构连接点、是否存在刚度不足、长期测量时此处的疲劳损耗等问题. 直杆探针应力最大处大多处于迎流方向,而本研究探针的应力最大位置则出现在背流方向,根本上是由图像法测量探针的独特外形造成的独特振动形态导致的.
计算得探针最大应力为16.5 MPa,远低于探针材料的许用应力137.0 MPa. 为验证是否会引发共振效应,计算得到探针1阶振动固有频率约为6.8 Hz,2阶振动固有频率约为42.4 Hz,3阶振动固有频率约为117.8 Hz. 计算出的探针振动频率为55.6 Hz,与探针振动固有频率偏离较大,不会引发共振.
在其他参数不变的情况下,悬臂梁的振动固有频率与其长度L的平方成反比. 在实际测量过程中,需要通过调整探针支点的位置来满足不同深度的测量需求,悬臂梁长度随之改变。计算得该探针悬臂梁长度L=1.135 m时,探针的二阶振动固有频率与振动频率接近,极有可能引发共振现象。测量过程中必须避开该支点位置,从而避免发生安全事故.
3. 结 论
(1)探针在X、Y、Z方向上均存在周期性振动,振动频率约为55.6 Hz,X方向正面受到流体冲击,振幅较大且有波动,最大振幅为0.28 mm. X方向上的振动基线偏离0点,Y、Z方向上的振动基线接近0点.
(2)由于所用探针尾段有光学测量部分,流场在此部分压差不大,削减了探针前端部分流场前、后压差平均值,使探针形成特殊的振动形态.
(3)探针应力最大位置为探针背流方向约2/3探针长度处,最大应力远低于探针材料的许用应力.
(4)在额定工况下,探针振动频率与探针振动固有频率偏离较大,不会引发共振现象. 当探针悬臂梁长度为1.135 m时,探针振动频率与探针振动固有频率接近,极有可能引发共振现象,为保证测量安全性,探针应避免在悬臂梁长度接近1.135 m的状态下工作.
参考文献
Experimental validation of condensing flow theory for a stationary cascade of steam turbine blades
[J].DOI:10.1098/rsta.1996.0003 [本文引用: 1]
An analytical solution for the Wilson point in homogeneously nucleating flows
[J].
On the performance of a cascade of turbine rotor tip section blading in nucleating steam: I: surface pressure distributions
[J].
湿蒸汽两相凝结流动中水滴生长模型研究
[J].
Droplets growth model in wet steam two-phase condensation flow
[J].
汽轮机叶片表面粗糙度对自发凝结级性能影响的数值研究
[J].
Numerical investigation of the effect of blade roughness on the performance of steam turbine nucleating stage
[J].
基于激光的后向角散射法测量蒸汽湿度
[J].DOI:10.3969/j.issn.1001-5078.2016.01.004 [本文引用: 1]
Research on steam humidity measurement by the backward angle scattering
[J].DOI:10.3969/j.issn.1001-5078.2016.01.004 [本文引用: 1]
Probes for industrial particle measurement: application for water droplet sizing in wet steam turbines
[J].DOI:10.1016/0021-8502(89)90866-5
Studying the characteristics of wet steam in turbine flow sections
[J].
试验汽轮机内湿蒸汽测量研究
[J].DOI:10.3969/j.issn.1672-5549.2012.02.003
Research on the measurement of wet steam in a test steam turbine
[J].DOI:10.3969/j.issn.1672-5549.2012.02.003
布拉格光纤光栅测量湿蒸汽两相流温/湿度的理论数学模型
[J].DOI:10.3321/j.issn:0258-8013.2006.07.008
Theoretic mathematics model based on fiber Bragg grating sensor measurement in two-phase wet steam flow temperature/moisture
[J].DOI:10.3321/j.issn:0258-8013.2006.07.008
Low-cost relative humidity sensor based on thermoplastic polyimide-coated fiber bragg grating
[J].DOI:10.1016/j.snb.2007.05.007 [本文引用: 1]
湿蒸汽联合探针的风洞标定和数值模拟的比较
[J].
Comparison between calibration and numerical simulation of an integrated probe for wet steam measurement
[J].
Unsteady flow field and coarse droplet measurements in the last stage of a low-pressure steam turbine with supersonic airfoils near the blade Tip
[J].
330MW低压汽轮机湿蒸汽二次水滴测量实验研究
[J].
Experimental research on the coarse water droplets measurement in a 330 MW LP steam turbine
[J].
进口探针支杆对1.5级压气机气动性能的影响
[J].
Effects of inlet probe support on aerodynamic performance of 1.5-stage compressor
[J].
探针支杆对压气机转子性能及流场影响的数值模拟研究
[J].
Numerical investigation of effects of probe support on performance and flow field of compressor rotor
[J].
叶型探针对跨声压气机性能影响的数值模拟
[J].
Numerical simulation of the effect of a leaf probe on the performance of a transonic compressor
[J].
/
〈 |
|
〉 |
