1. 引言
蒸汽调节阀是船舶动力系统中的流量调节元件,其产生的噪声是整个系统噪声中不可忽视的部分 [1]。调节阀内部的气动噪声分为两类:一是阀内节流结构在节流过程中产生的气动噪声;另一类是阀门内不稳定汽流与管道相互用作激发的噪音。本文关注的是前者,蒸汽阀内部结构复杂,阀套是其核心部件之一,用于控制气体的压力及流动状态。在高压差条件下,高速蒸汽流经调节阀内阀门套筒时,节流必然伴随着压力损失,部分流体机械能转换成声能,产生噪声 [2] [3]。
针对高速射流的噪声控制问题,已有学者进行了大量的研究。Zeng [4] 对汽轮机控制阀进行了实验与模拟,发现振动和噪声明显与流型密切相关,环形流和芯流的变化会引起声突变,与芯流相比,环形流具有更低的噪声和更稳定的流型。Peng Wang [5] 将调节阀主阀加装滤网,发现倾斜的涡流对在具有过滤器的构造中被矫直而变得对称,圆形过滤器有效地控制了主阀中复杂的蒸汽流,从而大大减小了主阀主轴上的横截面力。王干新 [6] 对不同内构件调节阀进行了实验研究,研究发现内构件的降噪关键在于如何将压力突变转化为压力渐变,盘片内构件为调节阀的主要降压区间,大大缓解了喉部的压力,从而降低喉部处产生的噪声。马大猷 [7] 结合实验测量和理论分析研究了小孔与小孔群对喷流噪声的影响。研究表明小孔可抑制射流冲击噪声,把射流噪声频谱推向高频范围,在相同压差条件下,单位面积小孔射流所辐射的噪声功率随孔径减小迅速下降,这正是小孔群控制射流噪声的基础。本文利用该原理将调节阀套筒窗型泄流孔改型为小孔群结构,并采用计算流体力学和计算气动声学方法计算调节阀内非定常湍流场和流动诱导的噪声,对比分析小孔群型阀套对蒸汽调节阀流场与声场分布的影响,为低噪声蒸汽调节阀的设计提供依据。
2. 数值计算
本文先对阀门内部流动进行定常模拟计算,收敛后,以定常场为初场进行非定常计算,获得阀内非定常流动特征及壁面非定常压强分布。最后,以壁面压强非定常压强为激励源,采用声学边界元方法计算阀内流噪声。
2.1. 流动控制方程及湍流模型
高温高压的过热蒸汽在调节阀内做三维可压缩的非定常湍流运动,可通过求解非定常雷诺时均N-S方程获得阀内流场数值解。在直角坐标系下,控制方程组的通用式 [8] 为:
(1)
式中:ρ为流体密度;t表示时间;U是速度矢量;
为通用变量,可以代表u、v、w、T、k和ε等待求变量;
表示广义扩散系数;
是广义源项,通过将式(1)对不同变量进行展开,可以得到相应的连续性方程、动量方程以及能量方程。
考虑到调节阀内复杂的剪切流与涡流流动特征,本文选择RNG k-ε两方程湍流模型对方程(1)进行封闭。李新一 [9] 的工作表明RNG k-ε模型包含低雷诺数流动效应和旋流修正子模型,对复杂剪切流动、旋流和分离流具有较好的模拟精度。获得调节阀内定常湍流解后,本文采用大涡模拟(Large Eddy Simulation, LES)方法计算阀内非定常湍流场。速度场与压力场的耦合计算则基于SIMPLEC算法。本文由于多采用四面体网格,故梯度法使用Green gauss node based,流动控制方程对流项采用具有二阶精度的有界中心差分格式进行离散。
2.2. 计算模型
图1是本文研究的蒸汽调节阀剖视示意图。可见,蒸汽进入阀门后,经阀套上的泄流孔绕阀芯折转向下,最后经阀芯与阀座间的间隙喷射进入出口管道。在此过程中,蒸汽速度大小、方向发生剧烈变化,伴随较为强烈的非定常流动现象。图1中A为沿调节阀出口段中心轴线布置的流速及噪声观测点,M1、M2、M3及M4则是四个压力脉动监测点,Surface 1为压损计算面,plane a、b、c则为湍动能观测面。
图2给出了调节阀阀套及泄流孔的示意图。原型窗型孔阀套的泄流孔由四组直径分别为25 mm和10 mm的圆孔组成,如图2。为控制蒸汽流经阀套时的噪声,在保证泄流孔面积基本不变的条件下,将窗型孔改型为由320个小孔组成的孔群型泄流孔,如图2所示。图中Z型泄流孔是孔径为3 mm的等截面圆柱孔,S型泄流孔是进口直径为3 mm、面积比为4/9的突缩孔,两者的气流流向均是由阀套外壁指向内壁面,K型泄流孔则与S型相反,气流由小孔经大孔流出。
![](//html.hanspub.org/file/42-2570495x11_hanspub.png?20220606104303529)
Figure 1. Schematic diagram of regulating valve structure and monitoring position
图1. 调节阀结构及监测位置示意图
![](//html.hanspub.org/file/42-2570495x12_hanspub.png?20220606104303529)
Figure 2. Flow hole pattern of valve sleeve
图2. 阀套泄流孔孔型图
2.3. 网格划分及边界条件
图3为蒸汽阀流体域示意图,其中进出口分别延长2倍和3倍,流体域开度取50%。为提高计算网格质量,本文将计算域分为进口域、核心域和出口域三个子域,分别单独生成网格。由于蒸汽阀核心域几何结构复杂,存在众多小孔和曲面结构,对其进行四面体网格划分,为保证网格质量,在核心域局部加密网格,同时对近壁区网格进行加密处理,近壁面第一层网格高度为0.008 mm,设置12层,增长率为1.05。数值模拟计算结果表明,调节阀壁面y+小于1,满足LES湍流模型计算要求。
![](//html.hanspub.org/file/42-2570495x13_hanspub.png?20220606104303529)
Figure 3. Computational domain and grid schematics
图3. 计算域及网格示意图
本文涉及的物理边界条件包括计算域进口、出口以及固体壁面边界。进口给定来流蒸汽总温552 K和总压6.3 MPa,出口给定蒸汽静压5.18 MPa,固壁均采用绝热、无滑移边界条件。考虑到人耳所能听到的频率上限为20,000 Hz,根据采样频率公式设置时间步长为2.5 × e−5。为消除初场影响,只取后2000个时间步作为可用数据。
2.4. 网格无关性验证
为验证网格数量对数值模拟计算准确性的影响,以窗型孔套筒模型为例,在阀门开度为10%的条件下,对五套不同的网格离散方案进行了数值计算。表1给出了五套计算网格获得的蒸汽质量流量以及观测点A处的速度大小,网格数量在300万到1000万之间。从表1可见,随着计算网格总量的不断增加,点位流速逐渐上升,误差趋于稳定。针对阀体内流动,进出口流量间差值均小于5%,可以认为数值解结果可靠。MSH-D与MSH-E对应的进出口质量流量和流速相差较小,考虑到MSH-E比MSH-D的网格数量多了近150万,为合理利用计算资源,在保证近壁面网格尺度与质量的前提下采用850多万的计算网格。
![](Images/Table_Tmp.jpg)
Table 1. Each curve represents the number of mesh
表1. 各曲线代表网格数
3. 流场计算结果分析
3.1. 流量特性及压力损失系数
图4给出了进出口压差Δp为1.12 MPa时,窗型孔阀套、Z型孔群阀套、S型孔群阀套和K型孔群阀套调节阀模型的流量系数随开度的变化。流量特性曲线 [10] 是在阀门两端压差保持恒定的条件下,介质流经调节阀的流量系数与它的开度之间的关系曲线。流量系数为衡量阀门流通能力的关键参数,常用
表示,计算公式见式
(2)
式中:Q——体积流量(m3/h);ρ——流体密度(kg/m3);
为阀门进出口的总压损失(102 Pa)。
从图中可见,三种阀的流量系数
随着阀门开度的增大均有显著提高。这种变化趋势与阀门内流动的基本规律吻合。当开度从10%到50%时,流量变化较为剧烈,而开度在50%到90%之间时变化相对平缓。从流通能力方面看,窗型孔最强,S型比Z型略高,K型孔最差。
图5是进出口压差为1.12 MPa时,压力损失系数ξ随调节阀开度的变化曲线,压力损失系数计算公式为:
(3)
式中:
——阀门进口与Surface 1的总压损失(Pa);ρ——流体密度(kg/m3);v——进口速度(m/s)。
由图5可知,随着阀门开度的提升,阀门压力损失系数逐渐下降,尤其在10%到50%段降低幅度较大,近83%,在50%到90%区间内变化幅度较小,这也与图4中的变化有所对应。而窗型孔损失最小,K型最大,Z型与S型相差较小。
![](//html.hanspub.org/file/42-2570495x20_hanspub.png?20220606104303529)
Figure 4. Variation diagram of flow coefficient of regulating valve with hole type and opening of valve sleeve
图4. 调节阀流量系数随阀套孔型与开度的变化图
![](//html.hanspub.org/file/42-2570495x21_hanspub.png?20220606104303529)
Figure 5. Influence of valve sleeve hole type on pressure loss coefficient of regulating valve (Δp = 1.12 MPa)
图5. 阀套孔型对调节阀压力损失系数的影响(Δp = 1.12 MPa)
3.2. 阀套孔型对阀内瞬态流场的影响
图6为Δp为1.12 MPa,阀门50%开度时ZX截面上某瞬时的速度流线及速度云图。由图6可知,蒸汽气流经窗型泄流孔通道后,速度增加迅速,时速超过240 m/s,在阀座区域形成了高速对流区域并在阀腔右侧出现大回流区和左侧的小涡。这种非对称分离流使得中心流域变窄。由于Z型孔、S型孔和K型孔是多孔孔群结构阀套,使得进入下游腔室的流动空间减少,极大地降低了喉口段及下游的蒸汽流速,但Z型孔没有改变局部的涡系结构的范围,而S型孔型与K型孔型的流线形更趋向于附壁流 [11],这可能是对降低噪声的有利因素。S型孔内发生了明显的节流现象,蒸汽气流进入小孔后加速通过,使气流相对有序平稳地进入阀芯下方,此处的射流区最大速度与高速冲击范围均比Z型孔有一定的减小,且S型孔右侧的涡旋区相较Z型孔又进一步减小了。
![](//html.hanspub.org/file/42-2570495x23_hanspub.png?20220606104303529)
Figure 6. The velocity distribution on the ZX plane at some instant
图6. 某瞬时ZX面的速度分布
由于流体在掺混过程中会发生能量的损失,为了更好地分析上述流场规律,图7给出了三处截面流道内的平均湍流动能分布,分别是plane a、b、c (见图1)。从plane a可以看出,窗型在左侧进口方向窗型套筒壁面内外两侧均出现较大的湍流动能梯度分布,且在阀套内圈流域产生低速的漩涡区。更换成孔群结构阀套后,射流入口段的高动能尾迹长度明显缩短,而K型孔更削减了尾迹动能的厚度。在K、S型的变截面小孔内有比大孔更强的湍动能。在plane b阀座段,窗型中部区域动能明显较大,Z型孔则减小了中心流域的湍流动能,并分成了三处小涡,K型与S型进一步降低了阀内流体的不稳定性 [12],平面内涡与高湍流区向四周分散。从截面a、b、c的纵向比较而言,与流场情况较接近,离冲击核心段越远湍动能越小。
阀门管道因90度弯角形成分离流向涡,涡束的发展导致下游壁面的压力波动,是管道的宽频噪声的主要激励源。阀门内腔的压力脉动是阀门重要的非定常特性指标,在阀门壁面处设置监测点(图1),监测其压力脉动情形。图8分别是阀座段左右两侧监测点M1、M2与扩压段M3、M4压强脉动的频谱分布。由图8可见,阀套泄流孔从窗型孔改为小孔群后,监测点的压力脉动幅值明显下降。M1点和M2处,窗型孔模型分别在980 Hz和1480 Hz处存在压强脉动峰值,其余三型模型频谱曲线中该峰值已完全被抑制,其中K型与S型孔压力脉动均方根差值相较Z型降低67.6%和73%。因为监测点M3、M4远离脉动核心区,相比于M1和M2点,整体脉动幅值扩压段比阀座段低。且由于扩压段流场波动趋于稳定,四种孔型间的脉动幅值差距较小。
![](//html.hanspub.org/file/42-2570495x24_hanspub.png?20220606104303529)
Figure 7. Multi-section turbulent kinetic energy diagram
图7. 多截面湍动能图
(a)
(b)
(c)
(d)
Figure 8. Pressure pulsation diagram of measuring point
图8. 测点压力脉动图
图9为阀芯0.1 s内所受非定常流体振荡力时域及频域图。图(a)表现了在核心域流体的紊乱流动下,蒸汽作用在阀芯的力随时间脉动,其中以窗型的波动最为剧烈。窗型、Z型、S型和K型的阀芯Z方向载荷脉动均方根值分别为6765 N、6681 N、6644 N和6636 N,窗型孔载荷在6721~6885 N内波动,Z型在6635~6716 N内波动,震荡力波动幅值减少约16.5%,S型在6619~6669N内波动,幅值减少约48.5%,K型在6608~6654 N内波动,幅值减少约52.6%。图(b)为对时域信号进行傅里叶变换后的非定常力频域信号,512 Hz之前载荷力先随着频率升高而增大,1024 Hz之后力随频率迅速衰减。窗型阀芯Z方向在1280 Hz存在峰值为8.3 N,在此频率下Z型孔幅值为4.2 N,K型孔为2.2 N,S型孔仅为0.98 N。由此可见小孔孔群阀套减轻了蒸汽阀内部的力学载荷强度 [13],提高了阀门内部的震荡稳定性,可能使阀门的疲劳损伤使用寿命延长。
4. 声场计算结果分析
本文以作用在调节阀内壁的声源作为激励源,采用间接边界元方法计算阀门声场。如图10所示,在距离阀门核心区中心点1 m处的环形每隔15˚布置一个噪声接收点布置点位共24个接收点,通过接收点监测阀门不同方向处声压级,分别计算不同阀套情况的声场变化,定量探究其中变化规律。
(a)
(b)
Figure 9. Time domain and spectrum diagram of unsteady fluid oscillation force in Z direction of valve core
图9. 阀芯Z方向非定常流体振荡力时域及频谱图
![](//html.hanspub.org/file/42-2570495x31_hanspub.png?20220606104303529)
Figure 10. Directional fixed measuring point position diagram
图10. 指向性固定测点位置示意图
声压级指向性分布如图11所示,可以看出阀门声场基本没有明显的极性特征,S型与K型降噪效果较好,但数值十分接近,S型孔在15到300度范围内声压级均略低于K型孔,仅在315到360度高于K型孔。图12为阀座段中心A监测点(见图1)的声压频谱图可见,阀门声压级呈现随频率的增大而减小的趋势,阀门的噪声基本为宽频噪声,在640 Hz之后的声谱频带三孔群阀套较原型有明显的下降,在此监测点原型声压级为154.1 dB,Z型声压级为147 dB,S型声压级为144.8 dB,K型声压级为141.8 dB。由此可见K型孔阀套降噪效果最佳。这验证了多孔结构能极佳地吸收声波以减小噪声的干扰。S型孔的二级降压结构增加了声阻 [14],Z型孔的抗性消声设计可进一步削弱声波能量达到吸收声波的效果。
![](//html.hanspub.org/file/42-2570495x32_hanspub.png?20220606104303529)
Figure 11. Sound pressure level directivity diagram
图11. 声压级指向性图
为更好地了解阀门流体域内部的声压分布,图13展示了100 Hz时的四孔型的声压云图。可以从中看出距离阀门节流区位置越远,气动声源声辐射的传播变得愈加衰弱,说明阀门内部声源来自于流体对阀套与阀芯的剧烈冲击。四类孔型均呈现出强声辐射区由核心噪声区通过管道传递至下游壁面附近的规律,但孔群孔型阀套对下游声传播有一定的遏制作用,其中S型孔与K型孔有更佳的效果。从下游管外的声辐射观察,窗型喉口与扩压段右侧声压辐射区声压较大,Z型通过孔群的节流大大缓解了喉部的压力,降低了喉部产生的声辐射,K型与S型又在此处进一步减小了噪声的量级,其中K型的扩张孔降噪效果最为突出。
5. 结论
本文采用流体可压缩模型和声学边界元方法研究了小孔孔群的流动及孔结构特征对阀门的噪声影响,得到以下结论:
![](//html.hanspub.org/file/42-2570495x34_hanspub.png?20220606104303529)
Figure 13. ZX plane sound pressure distribution diagram
图13. ZX平面声压分布图
1) 阀门流量随开度的增大而增大且增涨幅度也越大。孔群阀套相较于窗型孔阀套的流通能力有所减小,压力损失增大。小孔孔群阀套减轻了蒸汽阀内部的力学载荷强度,S型与K型节流孔将来自上游流体的剧烈冲击大大减缓,同时使气流进入下游的腔室更为有序,同时提高了阀的稳定性。
2) 窗型孔喉部及扩压段蒸汽气流速度较大,通过使用孔群阀套结构,不仅使得局部速度大大减小,速度分布更为均匀,同时控制了阀门壁面的分离涡,让气流更趋于附壁流。在阀内使用S型孔群阀套相比Z型孔可降低RMS压力脉动能级67.6%,K型孔阀套可降低73%。
3) 孔群阀套结构可有效降低了阀内部的气动噪声,其中由抗性消声器原理设计的K型孔降噪效果最为突出,将阀门管道中的噪声能量吸收,总声压级比原型降低13 dB。S型孔则利用二次降压与单次降压的Z型孔相比降低2.2 dB。综合考虑阀门性能与噪声之间的平衡将尤为关键。
基金项目
国家自然科学基金资助项目(51106099, 50976072)。