2. 上海电力股份有限公司闵行发电厂,上海 201100
2. Minhang Power Plant, Shanghai Electric Power Co. Ltd., Shanghai 201100, China
汽封由于可减少动静部件之间的间隙,减少漏气并提高机组运行效率,被广泛应用于汽轮机、航空发动机及燃气轮机等透平机械中[1-3]。在传统迷宫汽封[4]基础上,国内外学者先后提出了蜂窝[5]、刷式[6]、新型结构[7]等汽封型式,并通过大量实验与仿真研究密封泄漏特性,但如何有效提高汽封泄漏性能仍是一个难点。
汽封型式对汽封的泄漏特性起关键作用。Childs等[8]通过实验测试了螺旋齿汽封不同螺旋角度、转速、压比及预旋对密封泄漏量的影响,发现随着螺旋角度增大泄漏量逐渐增加。Zhou等[9]提出了一种新型交错螺旋齿结构迷宫汽封,通过数值模拟分析了不同密封间隙、压降、齿数与齿宽等对新型交错螺旋齿结构汽封泄漏量的影响规律,结果表明密封泄漏量随着压降和密封间隙增大而增加,随着齿数增加而减少,当齿数与密封长度一定时齿宽变大会增加泄漏量,而腔宽增加会减少泄漏量。张万福等[10]设计了一种新型径向环形汽封并通过Fluent软件分析了径向环形汽封与梳齿汽封在不同间隙与压比下的泄漏系数,结果表明在同一间隙下径向环形汽封的泄漏系数比梳齿汽封的降低约20%,且系统稳定性得到了提高。Jia等[11]研究了一种T型迷宫密封并建立了整圆计算模型,计算结果表明:在轴向长度与密封间隙相同时T型迷宫密封泄漏量比直通式迷宫密封的减少约23.6% ~ 25.3%,比交错式迷宫密封的增加7.4% ~ 8.5%。
结构参数同样对迷宫汽封泄漏量有着重要影响。Yücel等[12]通过实验与理论计算方法对阶梯式迷宫汽封泄漏量进行了研究,发现密封泄漏量随着压差增大而上升,且给出了阶梯迷宫汽封中可压缩性流体流动的连续性方程与周向动量方程。康乐嘉[13]通过建立三维数值模型分析了不同齿形结构对迷宫汽封泄漏量的影响,结果表明:在齿数足够且密封间隙一定时高低齿汽封比平齿汽封的泄漏量降低约10%,且侧齿汽封相较于平齿汽封的密封泄漏量降低约5%。宁哲等[14]建立了迷宫汽封三维数值模型并数值分析了压比、密封间隙、齿宽与台阶高度对阶梯式迷宫汽封的影响,研究发现汽封泄漏量随压比、密封间隙增大而增加,且上升速率逐渐降低,齿宽越小泄漏量越低,而台阶高度变化对汽封泄漏量影响较小。
本文基于此建立了梳齿汽封三维数值模型,通过实验与计算流体力学(CFD)对比的方法分析了在不同背压阀门开度、进口压力和转速下梳齿汽封泄漏量的变化规律。
1 实验装置与方法 1.1 实验装置图1为梳齿汽封泄漏密封实验装置。该装置主要由梳齿汽封实验台、背压装置及供气装置等组成。梳齿汽封实验台主要由滑动轴承、平衡板、气缸、转子、联轴器、增速器及无极变频电机构成,其中,空载下电机转速最高可达
|
图 1 梳齿汽封泄漏密封实验装置 Fig.1 Experimental device of labyrinth gas seal for its leakage |
图2为梳齿汽封二维几何结构,其中:d为转子直径;Cr为密封间隙;H为齿高;t为齿尖宽;w1为齿根宽;w2为腔宽;N为齿数。经压缩的高压空气由气缸中间进气,向两端排出,依次通过梳齿密封、背压腔、背压阀门和背压梳齿密封。其中,梳齿密封段长度为70.5 mm,背压梳齿密封段长度为20 mm,具体参数如表1所示。
|
图 2 梳齿汽封几何结构 Fig.2 Geometry of labyrinth gas seal |
|
|
表 1 梳齿密封件几何尺寸 Table 1 Geometric dimensions of seal |
梳齿汽封泄漏特性实验主要分为实验预处理、测试及数据后处理三个阶段。开机前,检查实验装置台是否有异物、油路系统是否正常以及冷却电机。汽封泄漏测试过程中,测试梳齿汽封在不同转速、进口压力与背压阀门开度下密封泄漏量的变化。其中,背压阀门开度为0代表阀门全关,背压阀门开度为1.0代表阀门全开。
表2为汽封泄漏特性的实验工况。梳齿汽封泄漏特性测试工质为空气,通过调节进气口阀门控制进口压力Pin,其范围为0.25 ~ 0.6 MPa。由于梳齿密封段出口为背压腔室,通过调节背压阀门控制背压P,泄漏量一部分通过背压阀门流向大气,一部分流经背压梳齿密封件。电机转速ω通过变频器调节,范围为0 ~ 1 500 r·min−1。
|
|
表 2 实验工况 Table 2 Experimental conditions |
使用前处理器GAMBIT建立梳齿汽封三维模型,为提高精度,采用结构化网格。为节约计算资源,对梳齿汽封件进行简化,模拟单侧流体流动,出口段为背压,计算域如图3所示。图4给出了梳齿汽封件的网格分布示意图。
|
图 3 计算域 Fig.3 Computational domain |
|
图 4 网格分布示意图 Fig.4 Schematic diagram of grid |
表3为数值计算时采用的求解方法。基于k−ε两方程湍流模型和三维RANS方程求解方法,采用CFD方法分析不同进出口压力、转速对梳齿汽封泄漏特性的影响,并与相同工况下的实验结果进行对比。计算时工作介质为空气;固壁面视为绝热、光滑壁面,同时近壁面处通过壁面函数法进行处理;采用高精度离散模型。当动量方程与连续性方程的残差数量级均小于10−6,且进、出口流量差小于0.1%时可认为计算收敛。
|
|
表 3 汽封泄漏特性数值计算方法 Table 3 Numerical method of leakage characteristics of the labyrinth gas seal |
在转速为500 r·min−1、背压阀门开度为1.0、进口压力为0.25 MPa工况下,选择4种网格数,分别为150万、215万、285万和365万,通过比较梳齿汽封的泄漏量进行网格无关性验证,计算结果如图5所示。
|
图 5 网格无关性验证 Fig.5 Grid-independent verification |
当网格数由150万增加至365万时,密封泄漏量逐渐变大且增长速率骤降。同时,当网格数增加而密封泄漏量偏差小于0.5%并趋于稳定时,则视为网格无关。
2.2 数值计算方法有效性验证为了验证实验方法的可行性,测试了在转速为500 r·min−1、背压阀门开度为1.0工况下不同进口压力对密封泄漏量的影响,并将实验值与CFD模拟值进行对比,结果如图6所示。
|
图 6 不同进口压力对汽封泄漏量的影响 Fig.6 Effect of inlet pressure on the leakage of labyrinth gas seal |
由图6可知,梳齿汽封泄漏量的实验值与模拟值均随进口压力变大而增加,且整体变化趋势近似,两者相对误差均小于5%,具有较高的吻合度,从而验证了该方法的可行性。
3 结果分析 3.1 阀门开度对背压的影响为了解背压腔室内背压的变化规律,通过调节背压阀门开度,测试了不同转速、背压阀门开度下进口压力对背压的影响,结果如图7所示。
|
图 7 不同转速与阀门开度下进口压力对背压的影响 Fig.7 Effect of inlet pressure on back pressure at different rotation speed and valve opening |
由图7(a)可知,进口压力对背压起关键作用,背压随着进口压力增加而变大,且上升速率逐渐放缓。背压阀门开度对背压有显著影响,在同一进口压力与转速下,背压阀门开度越大背压越小。在低进口压力下,背压幅值变化随着阀门开度变化甚微;在高进口压力下,背压幅值变化随着阀门开度变化剧烈。
由图7(a) ~ (d)可知,转速对背压幅值变化影响较小。在进口压力与阀门开度一定时,随着转速由0 r·min−1上升至1 500 r·min−1,背压幅值缓慢增大且增加速率越来越小。
3.2 进口压力对汽封流场特性的影响为了解梳齿汽封腔室内部流场变化情况,通过CFD方法分析了在转速为500 r·min−1、背压阀门开度为0.5工况下不同进口压力对汽封腔室压力场的影响,结果如图8所示。
|
图 8 不同进口压力下汽封腔室压力场 Fig.8 Effect of inlet pressure on pressure field in labyrinth gas seal chamber |
由图8可知,在不同进口压力下汽封各级腔室内的压力场变化均不明显,而在汽封入口段压降变化剧烈,在腔室四周形成小漩涡(中间压力低而四周压力高),这是由于在节流作用下气体速度骤升而压力急剧下降;压力在入口段与出口段变化显著而中间变化较小,同时腔室内压力沿轴向逐级递减;随着进口压力增大汽封形成的小漩涡在第一级腔室内逐渐消失,且缓慢向下一级腔室转移。
3.3 进口压力对汽封泄漏量的影响为了解梳齿汽封泄漏特性变化规律,通过实验测试了不同转速、阀门开度及进口压力下的梳齿汽封泄漏量,并将实验值与CFD模拟值进行对比,结果如图9所示。
|
图 9 不同进口压力对汽封泄漏量的影响 Fig.9 Effect of inlet pressure on the leakage of labyrinth gas seal |
由图9中的实验值可知,梳齿汽封泄漏量随着进口压力增大而上升,且近似呈线性关系;在同一进口压力下,转速与阀门开度对汽封泄漏量影响较小,随着转速升高、背压阀门开度减小,密封泄漏量缓慢减少。由于实验过程中存在误差,实验值略高于模拟值。
由图9中的模拟值可知,在低进口压力下,背压阀门开度对汽封泄漏量的影响甚微;在高进口压力下,汽封泄漏量随着背压阀门开度增大而显著增加。相较于进口压力的影响,转速对汽封泄漏量影响较小,随着转速提高,泄漏量缓慢下降且近似保持恒定。模拟值与实验值整体变化趋势相同且吻合度较高。
4 结论通过实验测试了不同背压阀门开度、转速及进口压力对梳齿汽封泄漏量的影响,并与CFD模拟值进行对比,主要结论为:
(1)梳齿汽封泄漏量的实验值与模拟值整体变化趋势相同且吻合度较高,两者相对误差较小。
(2)同一背压阀门开度下,梳齿汽封泄漏量随着进口压力增大而增加,近似呈线性关系;而随着转速升高,泄漏量缓慢下降且近似保持恒定。
(3)在进口压力和转速一定时,汽封泄漏量随着背压开度的变大而增加,且在背压阀门开度为1.0时,汽封泄漏量达到峰值。
| [1] |
CHUPP R E, HENDRICKS R C, LATTIME S B, et al. Sealing in turbomachinery[J]. Journal of Propulsion and Power, 2006, 22(2): 313-349. DOI:10.2514/1.17778 |
| [2] |
史进渊, 李军, 刘霞, 等. 我国大型汽轮机技术研究进展与展望[J]. 动力工程学报, 2022, 42(6): 498-506. |
| [3] |
王玉明, 杨惠霞, 姜南. 流体密封技术[J]. 液压气动与密封, 2004(3): 1-5. |
| [4] |
李志刚, 郎骥, 李军, 等. 迷宫密封泄漏特性的试验研究[J]. 西安交通大学学报, 2011, 45(3): 48-52,64. |
| [5] |
盛伟, 刘旭东, 卫运钢. 蜂窝汽封流场及泄漏特性的数值模拟研究[J]. 汽轮机技术, 2009, 51(3): 201-203,240. |
| [6] |
李尊平, 丁常富, 刘志博. 刷式汽封和梳齿汽封泄漏流动特性的数值模拟研究[J]. 汽轮机技术, 2016, 58(5): 348-350,369. |
| [7] |
吴强. 迷宫式汽封和薄叶式汽封的数值模拟与对比分析[D]. 哈尔滨: 哈尔滨工业大学, 2008.
|
| [8] |
CHILDS D W, GANSLE A J. Experimental leakage and rotordynamic results for helically grooved annular gas seals[J]. Journal of Engineering for Gas Turbines and Power, 1996, 118(2): 389-393. DOI:10.1115/1.2816601 |
| [9] |
ZHOU W J, ZHAO Z B, WANG Y F, et al. Research on leakage performance and dynamic characteristics of a novel labyrinth seal with staggered helical teeth structure[J]. Alexandria Engineering Journal, 2021, 60(3): 3177-3187. DOI:10.1016/j.aej.2020.12.059 |
| [10] |
张万福, 杨建刚, 曹浩, 等. 径向环形汽封流动特性研究[J]. 中国电机工程学报, 2012, 32(20): 115-119,148. |
| [11] |
JIA X Y, ZHENG Q, JIANG Y T, et al. Leakage and rotordynamic performance of T type labyrinth seal[J]. Aerospace Science and Technology, 2019, 88: 22-31. DOI:10.1016/j.ast.2019.02.043 |
| [12] |
YÜCEL U. Calculation of leakage and dynamic coefficients of stepped labyrinth gas seals[J]. Applied Mathematics and Computation, 2004, 152(2): 521-533. DOI:10.1016/S0096-3003(03)00574-5 |
| [13] |
康乐嘉. 迷宫式汽封内部流动的数值模拟与结构优化[D]. 北京: 华北电力大学, 2006.
|
| [14] |
宁哲, 杨荣祖. 迷宫汽封泄漏特性数值模拟[J]. 热力发电, 2017, 46(12): 87-92. |
2024, Vol. 40

