摘要
为提高涡轮级间篦齿的封严性能,设计了一种带分流通道的台阶型篦齿封严结构,研究了不同压比与转速下分流通道对篦齿封严泄漏特性和传热特性的影响规律,并与光滑篦齿进行对比。结果表明:在衬套上增设一个分流通道后,处于通道出口处的气流在离心力作用下被卷吸进通道,增强了篦齿的封严效果;分流通道的周向间隙率β越大,封严效果越佳,当周向间隙率β为3时,绝热壁面流量系数为13.47%,绝热壁面风阻温升为48.7 K;分流通道数量增加,篦齿的泄漏加剧;随着压比增加,密封泄漏量增大,分流通道可以引导核心射流回流,故泄漏量增速逐渐变缓而绝热壁面风阻温升增速加剧;当转速从3×1
非接触式环形篦齿封严通常用于多级旋转机械,如压气机和涡轮,目的是限制气体泄漏、调节各个腔室之间的压力平衡。这些密封件的封严性能可以直接影响发动机的效
对比诸多的密封技术,篦齿封严(Labyrinth seal, LS)具有低成本、维护简单和结构可靠等特点,经常被作为首选密封方
对于上述问题,Morgan
除了改变篦齿的流道来改善篦齿封严的泄漏性能,另一种方式则是设置流体屏障。流体屏障是一种以最小的几何变化来减少泄漏的流动控制方法,它主要通过高速射流阻断泄漏气流来实现密封性能的提
流体屏障方案需要从外部引入气体,为此,Ghaffari
本文在传统台阶型篦齿封严的基础上,提出了一种新的结构——带分流通道的台阶型篦齿封严,并通过对该结构进行数值仿真计算,研究分流通道周向间隙率、通道沿程数量、压比和转速等因素对篦齿封严性能的影响。
带分流通道的台阶型篦齿结构如

图1 篦齿模型示意图
Fig.1 Schematic diagram of the labyrinth model
Parameter | Value | Parameter | Value |
---|---|---|---|
c/mm | 0.3 | t/mm | 0.25 |
B/mm | 8 | H/mm | 3 |
b1/mm | 1 | b2/mm | 2.7 |
α1/(°) | 4.5 | α2/(°) | 1.5 |
基于商业软件Fluent Meshing进行多面体网格划分,在篦齿齿顶间隙、齿腔及分流通道处施加网格局部加密,采用Poly⁃Hexcore非结构体网格生成方法,使六面体网格与多面体网格共节点连接,提升整体网格中六面体网格的数量以节省网格量,同时提升计算效率和收敛速度。转子及静子壁面设置5层边界层网格,第一层棱柱网格高度为0.003 mm,增长率设置为1.2。经计算验证,转静盘壁面
台阶型篦齿腔内的气流流动可视为稳态的三维流动问题,篦齿转静盘壁面均设置为绝热。分流通道沿周向离散均匀分布,通道的周向旋转角度α1为4.5°,周向间隔角度α2为1.5°。为了保留在周向旋转特征下尽可能减少仿真计算网格量,模型周向旋转6°(周向仅包含一个分流通道)作为数值计算区域,计算域示意图见
针对光滑台阶型篦齿模型,选用近244.4万、301.4万、351.8万、398.1万、458.4万的网格数目,验证网格的无关性并绘制流量系数CD与风阻温升随网格数变化的关系曲线,如

图2 网格无关性验证
Fig.2 Verification of grid-independence
研究结果表明,当网格数达到300万后,继续增加网格量,流量系数与风阻温升曲线均趋于稳定值,篦齿流量系数与风阻温升的变化均小于0.5%,可认为选取该网格数对计算结果的影响已经不大。故后续数值计算研究选取的网格数均在300万左右。
选用商业软件ANSYS Fluent进行数值计算求解,湍流模型选用Realizable k⁃ε模型。湍流方程中考虑黏性耗散项和曲率修正。计算过程中不考虑壁面传热。流体域工质为考虑可压缩性的理想气体。空气的物性参数考虑温度的影响,其中导热系数和比热容按Kinetic理论计算,黏性系数按Sutherland公式计算,自由度设置为5。动量方程、能量方程、湍流方程中的对流项采用二阶迎风格式,压力和速度耦合采用Simple⁃C算法。数值计算的收敛通过两个判据确定:(1)各方程残差小于1.0×1
进口截面为压力进口条件,给定总压、总温、轴向进气速度方向,出口截面为压力出口条件,给定静压。出口压力恒定为100 kPa,进口压力根据进出口压比给定,进口温度设定为600 K。计算利用冻结转子法处理衬套和转子的相对运动,旋转坐标系固定在转子上,转子相对于旋转坐标系静止,衬套绝对静止,流域在旋转坐标系下求解,旋转速度为15 000 r/min,篦齿转盘的最小旋转半径r=145 mm。
由于与本文计算模型对应的实验结果尚未得到,为验证计算方法的可靠性,暂选用文献[

图3 计算方法验证
Fig.3 Validation of computational methods
从
(1)篦齿雷诺数Re定义如下
(1) |
式中:m为篦齿泄漏空气的质量流量;μ为空气的动力黏度;R为篦齿的齿尖位置半径大小。
(2)篦齿压比(Pressure ratio, PR)定义为模型进口截面的总压与出口截面的静压之比,即
(2) |
式中:为进口总压;pout出口静压。
(3)篦齿风阻温升定义为
(3) |
式中:T为篦齿出口段截面的总温;T为入口段截面的总温。
(4)流量系数CD定义为实际泄漏流量m和理想流量mideal之比,即
(4) |
理想质量流量mideal定义为
(5) |
式中:Rg和k分别为气体常数和绝热指数;A一般取为级间篦齿齿顶的最小流通面积,常用来反映泄漏特性。
(5)周向间隙率β定义为
(6) |
式中:α1为分流通道的周向旋转角度;α2为分流通道的周向间隔角度。该参数反映了分流通道的周向分布情况。
本文主要研究分流通道对篦齿封严性能和风阻温升的影响,并对分流通道的周向间隙率、通道数量、旋转速度和进出口压比这4个参数进行研究。为此,所计算的工况为:旋转速度N=3×1

图4 篦齿封严流线图
Fig.4 Flow diagram of the labyrinth seal
考虑到分流通道沿周向设置为均匀离散分布,通道间隙率不同对其周向漩涡流场结构也会有不同的影响,故下面针对分流通道的周向间隙率进行研究。计算选取了5个不同周向间隙率的篦齿模型,其间隙率β分别为1/3、1/2、1、2、3,

图5 5种周向间隙率下通道出口处径向速度分布
Fig.5 Radial velocity distribution at the outlet of the channel under five circumferential clearance rates
气流流入篦齿进口段,受到转子壁面旋转拖拽从而获得了周向速度。在分流通道内部时,气流沿周向流动不断冲击分流通道的侧壁面1,造成通道出口处靠近侧壁面1区域出现低压,齿腔内的气流受到离心力和压差作用回流入通道内,而靠近侧壁面2处的气流则从通道内部泄漏到篦齿腔。
对比不同周向间隙率的径向速度分布,当β从1/3增加到1/2时,气流的回流区域有所减小,之后随着通道β的增大,出口截面的回流区域与泄漏区域均在增大,且回流区域明显大于泄漏区域。综上所述,分流通道在周向间隙率较大时能很好地起到抑制泄漏的作用。

图6 流量系数和风阻温升随通道周向旋转角度的变化
Fig.6 Variation of flow coefficient and wind resistance temperature rise with channel circumferential rotation angle

图7 流量系数和风阻温升随通道数量的变化
Fig.7 Variation of flow coefficient and wind resistance temperature rise with the number of channels
为分析流量系数在下游增设分流通道后大幅升高的原因,现以沿程通道数1和2时的第一级分流通道出口截面径向速度分布云图为对象进行分析。从

图8 第一级分流通道出口截面径向速度分布
Fig.8 Radial velocity distribution at the outlet section of the first stage diversion channel
通过上述对分流通道周向旋转角度以及分流通道沿程数量的影响研究,发现通道的旋转角度增大后篦齿流量系数逐渐下降,其密封效果增强,而在下游增加分流通道会降低篦齿封严性能。因此,下文压比与转速的研究均以分流通道周向旋转角度为4.5°(即周向间隙率β=3)、沿程分流通道数量为1的篦齿模型作为研究对象。

图9 不同转速下通道出口处径向速度分布(PR=1.5)
Fig.9 Radial velocity distribution at the outlet of the channel at different rotational speeds (PR=1.5)

图10 转速对流量系数和风阻温升的影响(PR=1.5)
Fig.10 Effect of rotational speed on flow coefficient and wind resistance temperature rise (PR=1.5)
对于篦齿风阻温升变化,随着篦齿的旋转速度增加,气流流经篦齿的风阻温升逐渐升高。在相同转速下,带分流通道的篦齿风阻温升要高于光滑篦齿。当篦齿的旋转速度为3×1

图11 压比对流量系数的影响(实线为带分流通道篦齿封严,虚线为光滑篦齿封严)
Fig.11 Influence of pressure ratio on flow coefficient(Solid lines represent the labyrinth seal with flow channels and dashed lines represent smooth labyrinth seal)

图12 压比对风阻温升的影响(实线为带分流通道篦齿封严,虚线为光滑篦齿封严)
Fig.12 Effect of pressure ratio on wind resistance temperature rise(Solid lines represent the labyrinth seal with flow channels and dashed lines represent smooth labyrinth seal)
针对分流通道对台阶式篦齿的影响,在同一转速下,当气流进出口的压比相同时,带分流通道的篦齿内部气流的风阻温升要远高于光滑篦齿,而且随着转速的增加,两个模型风阻温升的差异也更加显著。出现该现象主要是由于数值模拟设置的衬套和转子壁面均为绝热壁面,相较于光滑篦齿,带分流通道后,大量气流被卷吸进通道内,而通道内上下游气流互相碰撞、掺混产生的耗散热量又无法通过壁面传导出去,从而在通道处形成了局部高温区域,致使气流的风阻温升大幅升高。当转速升高时,分流通道内上下游气流之间的碰撞、掺混也随之增强,气流产生的耗散热量也随之增加,最终其与光滑篦齿风阻温升的差异也越大。
本文主要研究分流通道对篦齿封严性能和风阻温升的影响,并对分流通道的周向间隙率、通道数量、旋转速度和进出口压比这4个参数进行研究。为此,所计算的工况为:旋转速度N=3×1
(1)分流通道对篦齿封严的作用主要是通过离心作用使核心射流被卷吸进分流通道内,从而增大了气流流动阻力,降低了台阶型篦齿的流量系数。
(2)分流通道的周向间隙率β对篦齿的封严特性有所影响,当β从1/3增大到1/2,篦齿的流量系数有所升高,风阻温升有所降低;继续增大β,篦齿的流量系数逐渐降低、风阻温升逐渐升高。当β≥1时,增大周向间隙率β有助于提高篦齿的封严性能。
(3)随着分流通道数量的增加,篦齿的流量系数逐渐升高,分流通道的抑制泄漏能力逐渐恶化,其中当沿程分流通道从数量1增加到2时,流量系数变化非常剧烈,升高了4.7%。在下游增设分流通道会导致分流通道对台阶型篦齿的密封性能提升作用被抵消,不利于篦齿的封严。
(4)随着旋转速度从3×1
(5)在固定转速下,随着气流进出口压比从1.1增加到1.5,气流的风阻温升不断降低,且风阻温升随进出口压比的变化曲线近似线性;同时篦齿的流量系数随之升高,篦齿的封严性能不断恶化。
参考文献
LI Guoqing, ZHANG Qian, LEI Zhijun, et al. Leakage performance of labyrinth seal for oil sealing of aero-engine[J]. Propulsion and Power Research, 2019, 8(1): 13-22. [百度学术]
LIU Kuo, HAN Lingsheng, WANG Yongqing, et al. Theoretical and experimental research towards labyrinth sealing mechanism of liquid nitrogen in the cryogenic spindle[J]. Mechanical Systems and Signal Processing, 2022, 167: 108502. [百度学术]
RATHMANN U, OLMES S, SIMEON A. Sealing technology: Rub test rig for abrasive/abradable systems[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Montreal, Canada: ASME, 2007: 223-228. [百度学术]
UNTAROIU A, WOOD G H, MORGAN R N. Design of experiments to investigate geometric effects on fluid leakage rate in a balance drum seal[J]. Journal of Engineering for Gas Turbines and Power, 2016, 138(7): 072506. [百度学术]
CANGIOLI F, VANNINI G, PENNACCHI P, et al. Rotordynamic characterization of a staggered labyrinth seal: Experimental test data and comparison with predictions[J]. Journal of Engineering for Gas Turbines and Power, 2019, 141(1): 011009. [百度学术]
LI Zhigang, LI Jun, FENG Zhenping. Labyrinth seal rotordynamic characteristics part Ⅰ: Operational conditions effects[J]. Journal of Propulsion and Power, 2016, 32(5): 1199-1211. [百度学术]
LI Zhigang, LI Jun, FENG Zhenping. Labyrinth seal rotordynamic characteristics part Ⅱ: Geometrical parameter effects[J]. Journal of Propulsion and Power, 2016, 32(5): 1281-1291. [百度学术]
CHUN Y H, AHN J. Effects of geometric parameters of a staggered labyrinth seal on leakage flow[J]. Journal of Mechanical Science and Technology, 2023, 37: 2959-2968. [百度学术]
HILFER M, HOGG S, INGRAM G. Experimental validation of a curtain type fluidic jet seal on a turbine rotor shroud[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Montreal, Quebec, Canada: ASME, 2015: V05CT15A012. [百度学术]
MACCALMAN J, HOGG S, INGRAM G. Application of fluidic curtains to turbine rotor tip seal geometries[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Oslo, Norway: ASME, 2018: V008T29A017. [百度学术]
HOGG S I, GOMEZ RUIZ I. Fluidic jet barriers for sealing applications[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. [S.l.]: ASME, 2011: 779-788. [百度学术]
RAI A C, PRABHUDHARWADKAR D, MURTHY S, et al. Effect of air-curtains on labyrinth seal performance[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Seoul, South Korea: ASME, 2016: V05AT15A021. [百度学术]
唐萌, 赵晨旭, 翟世杰. 射流对不同封严篦齿特性的影响分析研究[J]. 机械制造与自动化, 2022, 51(3): 130-134. [百度学术]
TANG Meng, ZHAO Chenxu, ZHAI Shijie. Study on effect of jet on characteristics of different labyrinth seal[J]. Machine Building & Automation, 2022, 51(3): 130-134. [百度学术]
曹永华. 一种基于射流抑制泄漏的篦齿封严效果研究[D]. 南京:南京航空航天大学, 2017. [百度学术]
CAO Yonghua. A study based on using jet to suppress leakage flow in the labyrinth seal[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2017. [百度学术]
GHAFFARI P, WILLINGER R, BAUINGER S, et al. Impact of passive tip-injection on tip-leakage flow in axial low pressure turbine stage[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Montreal, Quebec, Canada: ASME, 2015: V02AT38A002. [百度学术]
GHAFFARI P, WILLINGER R. On the impact of passive tip-injection on the downstream flow field of a shrouded LP-turbine: CFD and experimental results[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Seoul, South Korea: ASME, 2016: V02BT38A009. [百度学术]
WASILCZUK F, FLASZYNSKI P, KACZYNSKI P, et al. Air curtain application for leakage reduction in gas turbine shroud sealing[J]. Aerospace Science and Technology, 2021, 112: 106636. [百度学术]
WASILCZUK F, FLASZYŃSKI P, KACZYŃSKI P, et al. Passive reduction of leakage by transverse jet[J]. Journal of Fluids Engineering, 2021, 143(1): 011204. [百度学术]
WASILCZUK F, FLASZYNSKI P, PYCLIK L, et al. Feasibility study of fluidic sealing in turbine shroud[J]. Materials, 2021, 14(13): 3477. [百度学术]
GU Qianlei, WU Zhengfeng, YANG Jiangang, et al. A flow control method for the leakage reduction of annular gas seals using the helical gas curtain generator[J]. Aerospace Science and Technology, 2023, 141: 108528. [百度学术]
ZHANG Li, ZHU Huiren, LIU Cunliang, et al. Experimental and numerical investigation on leakage characteristic of stepped labyrinth seal[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Seoul, South Korea: ASME, 2016: V05AT15A016. [百度学术]