«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2019, Vol. 40 Issue (5): 938-943  DOI: 10.11990/jheu.201803099
0

引用本文  

田志涛, 郑群, 姜斌, 等. 氦压气机转子叶顶间隙结构研究[J]. 哈尔滨工程大学学报, 2019, 40(5): 938-943. DOI: 10.11990/jheu.201803099.
TIAN Zhitao, ZHENG Qun, JIANG Bin, et al. Investigation on the tip clearance structure of a helium compressor rotor[J]. Journal of Harbin Engineering University, 2019, 40(5): 938-943. DOI: 10.11990/jheu.201803099.

基金项目

国家自然科学基金项目(514776039,51676050);博士创新基本科研项目;中央高校基金项目(HEUCF180305,HEUCFP201847)

通信作者

郑群, E-mail:zhengqun@hrbeu.edu.cn

作者简介

田志涛, 男, 博士研究生;
郑群, 男, 教授, 博士生导师

文章历史

收稿日期:2018-03-28
网络出版日期:2018-12-28
氦压气机转子叶顶间隙结构研究
田志涛 , 郑群 , 姜斌 , 赵文峰 , 段昱     
哈尔滨工程大学 动力与能源工程学院, 黑龙江 哈尔滨 150001
摘要:为了解决高负荷设计方法导致高负荷氦气压气机叶顶泄漏问题,本文采用数值模拟的方法,对氦压气机转子叶顶间隙泄漏涡形成机理、转子上端壁间隙结构控制机理及其对叶顶间隙泄漏涡的控制作用进行了研究。结果表明:氦压气机转子叶顶前缘单侧槽对叶顶间隙泄漏损失的控制效果最为明显;槽的坡度为30°时,氦压气机转子的性能改善幅度最大;叶顶间隙越大,叶顶前缘单侧槽的有效性越强,当叶顶间隙为1 mm时,能够使压气机转子的绝热效率提高0.15%。该研究能有效控制叶顶间隙泄漏损失,提高压气机设计点性能。
关键词氦压气机    高负荷设计    数值模拟    单侧槽    叶顶泄漏    间隙结构    绝热效率    
Investigation on the tip clearance structure of a helium compressor rotor
TIAN Zhitao , ZHENG Qun , JIANG Bin , ZHAO Wenfeng , DUAN Yu     
College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China
Abstract: The highly loaded design method of helium compressors can effectively solve the compression problem arising from the use of helium in high-temperature gas-cooled reactors. However, it also leads to more leakage problems at the tip clearance of a highly loaded helium compressor. In this paper, the numerical simulation is used to study the formation mechanism of the tip clearance leakage vortex of a helium compressor rotor, the control mechanism of the rotor upper end wall gap structure, and its control function on the tip clearance leakage vortex. Results showed that the unilateral groove on the leading edge of the rotor tip of the helium compressor has the most obvious control effect on the tip clearance leakage loss, and the performance improvement of the helium compressor rotor is the greatest when the groove measures 30°. In addition, the greater the clearance, the stronger the effectiveness of the unilateral groove. Finally, when the tip clearance is 1 mm, the adiabatic efficiency of the compressor rotor increased by 0.15%. The single side slot on the tip clearance of the helium compressor rotor can effectively control tip clearance leakage loss and improve compressor performance at design point.
Keywords: helium compressor    high load design    numerical simulation    unilateral groove    tip clearance leakage    clearance structure    adiabatic efficiency    

氦压气机是高温气冷堆动力转换单元的关键部件之一,其性能的优劣直接关系到高温气冷堆发电系统的功率输出与运行效率[1-2]。叶顶间隙泄漏对氦压气机性能有很大影响,如堵塞压气机主流流道、增加压气机气动损失以及降低压气机的运行稳定性。泄漏流与通道二次流相互作用产生泄漏涡,泄漏涡在后续发展过程中占据主流通道,降低了流道的流通能力,其强烈的分离效应是导致压气机失稳的主要因素之一[3-4]。因此,发展叶顶泄漏流控制技术对于提高压气机效率、扩展压气机稳定工作范围具有十分重要的理论价值和现实意义。目前控制压气机叶顶间隙泄漏损失的方法主要可以分为主动控制技术和被动控制技术2大类[5-6]。主动控制技术大多需要附属设备,能够通过附属设备调节控制量,调节准确性高,但是系统实现较为复杂。被动控制方法一般不需要附属设备,实现相对简单,被动控制技术一直发展至今。常见的被动控制技术主要通过机匣处理、减小叶顶间隙、改变叶尖结构等形式实现[7-8]。机匣处理作为一种主要的被动控制手段,主要的实现形式有周向槽、缝式处理槽、蜂窝和斜锥孔3种方式,能使喘振裕度有10%左右的提高,但是不同程度的降低了压气机的效率[9-10]。在高温气冷堆中的变工况主要通过改变闭式循环系统内部的工质密度实现,即质量流量改变而体积流量不变[11]。针对问题,本文提出了一种氦压气机转子叶顶前缘槽结构,采用流体力学计算软件ANSYS CFX 18.0作为数值计算工具,对此种前缘槽的控制机理进行了研究。

1 计算数值模型

本文的研究对象是某三级氦压气机第一级动叶,具体设计参数为进口总压2.442 MPa;进口总温308 K;转速5 540 rad/min;叶片数143;叶片高度h为53.8 mm;弦长c为33.6 mm;叶顶间隙τ为0.5、0.75、1 mm;稠度为1.8;轴向几何进气角及出气角为44.15°、14.08°。计算模型叶片的中弧线为抛物线,厚度分布为三次多项式,最大厚度的位置为弦长的35%处,中径处的最大厚度约为弦长的8%。

本文计算所采用的网格采用NUMACA软件的Autogrid 5模块对其进行结构化网格划分。其中,叶片通道采用HOH型网格拓扑结构,叶顶间隙采用蝶形网格。根据选取湍流模型的要求,对叶片表面边界层内网格进行加密,第1层网格厚度为3×10-6 m,网格如图 1所示。

Download:
图 1 计算网格 Fig. 1 Computational grid and local details

为了尽可能的保证计算的准确性,本文对网格的无关性进行了验证。在叶顶间隙为0.75 mm的情况下,选取了39万、45万、56万、73万和92万进行数值计算。计算结果如图 2所示,随着网格数量的增加叶顶间隙的相对泄漏量也有所增加,当网格的数量达到73万后,氦压缩机转子叶顶的泄漏量基本稳定,因此权衡计算的精度与计算的成本后选定73万网格为计算网格。

Download:
图 2 不同网格节点的计算结果 Fig. 2 Numerical results with different grid points

本文采用商用数值模拟软件ANSYS CFX 18.0进行计算。设计工况边界条件如下:进口总压设定为2.442 MPa,进口总温为308 K,进口湍流度为5%,出口平均静压值为2.368 MPa,壁面采用光滑、绝热、无滑移边界条件。计算工质为理想氦。通过比较不同湍流模型的计算结果,如图 3所示,SST湍流模型能够更为准确的贴合实验值[12],因此本文选择SST湍流模型封闭方程组。

Download:
图 3 不同湍流模型的计算值与试验值对比 Fig. 3 Predicted and experimental isentropic efficiency
2 压气机转子上端壁结构形式对叶顶泄漏的影响 2.1 转子上端壁槽的形式对叶顶泄漏的影响

为了研究氦压气机转子上端壁槽道形式对其叶顶泄漏的影响,本文选取了3种结构形式。如图 4所示,第1种为常规平顶结构,第2种为叶顶前缘凹槽结构,第3种为叶顶前缘与尾缘双侧凹槽结构。在进行此研究时保证叶顶结构间隙尺寸相同。

Download:
图 4 研究的叶顶示意 Fig. 4 Representation of the tip topologies under investigation

3种不同叶顶结构形式的氦压缩机绝热效率计算结果分别为0.925 72、0.926 75及0.925 59。叶顶前缘单侧槽结构形式相比与平顶结构和前缘与尾缘双侧槽结构氦压缩机具有最高的效率,初步认为这是由槽结构减少叶顶泄漏损失与槽结构本身的附加损失2方面因素共同作用的结果。

图 5给出了在98%叶高处,叶片表面的静压系数分布,用以研究氦压气机转子叶顶结构形式对叶顶载荷分布的影响。叶片表面静压系数为:

$ {C_{\rm{P}}} = \left( {p - {p_{{\rm{in}}}}} \right)/\left( {{p_{{\rm{t}}{\rm{.in}}}} - {p_{{\rm{in}}}}} \right) $ (1)
Download:
图 5 98%叶高处的表面静压系数分布 Fig. 5 Static pressure coefficient distributions of rotor at 98% height

式中:CP为叶片表面静压系数;p为当地静压;pin为进口静压;pt.in为进口总压。如图 5所示,叶顶前缘凹槽使得氦压气机叶顶前缘处的静压发生跃升,同时相较于平顶结构前缘凹槽使得静压系数的整体水平有所提高。尾缘槽对氦压气机叶顶静压分布的影响不大,因此前缘槽是改善氦压气机叶顶性能的主要手段。

图 6为氦压气机转子叶顶前缘凹槽的作用机理图。如图 6所示,氦流经叶顶前缘凹槽后形成旋涡,此旋涡作用于转子叶片前缘使其降低了对冲角的敏感性,减小了叶片对流体的轴向切削作用,从而减小叶片对来流的剪切作用使得氦压气机转子叶顶泄漏涡减弱从而减少其叶顶泄漏损失。

Download:
图 6 前缘单侧槽作用机理 Fig. 6 Action mechanism of unilateral groove clearance

图 7为氦压气机转子叶顶端壁静压系数分布云图。如图 7所示,叶顶前缘与尾缘双侧槽的氦压气机的上端壁的低压区最小,单侧槽的次之,无槽的低压区最大。这说明叶顶双侧槽能够更为有效的减小叶顶泄漏涡,从而具有最小的叶顶泄漏损失。但是由于双侧槽形成了2个附加旋涡,使得附加损失增加较多,因此综合叶顶间隙泄漏损失以及叶顶槽所带来的附加损失的影响,叶顶前缘槽对氦压气机整体性能的提高最为有效。

Download:
图 7 上端壁表面静压分布 Fig. 7 Contour of shroud surface pressure distribution
2.2 叶顶前缘槽的角度对叶顶泄漏的影响

为了进一步研究叶顶前缘槽对叶顶泄漏的影响,在叶顶间隙为1 mm不变的前提下,本文选取了30°、45°以及60°三种坡度的槽作为研究对象,研究叶顶前缘槽的控制机理。槽的结构如图 8所示,槽的最深点位于氦压气机叶顶前缘正上方,槽的深度与叶顶间隙的大小一致,槽的坡长与槽的倾斜角度直接相关。

Download:
图 8 叶顶槽结构图示意 Fig. 8 Schematic diagram of groove structure

图 9为相同槽深不同坡度条件下的氦压气机转子叶顶叶片表面静压分布情况。如图 9所示,槽的坡度越大对叶顶叶片前缘的压力扰动越剧烈,这主要是因为坡度越大上端壁的后台阶作用越强从而形成更加强烈的旋涡使叶顶叶片前缘静压分布产生突变。在叶片压力面, 由于槽的深度相同所以压力跃升的最高值基本相同但是随着坡度的增加压力的回落值逐渐降低。而在叶片的吸力面,压力的跃升随着坡度的增加而减小。叶片从70%弦长到尾缘区间3种方案下叶片压力面静压分布差异明显,这主要是受人造旋涡的强度影响产生的,坡度越大产生的旋涡越强,对后续的流场扰动越大,在中间弦长位置载荷较大,因此表现出的差异较小。

Download:
图 9 不同坡度98%叶高处静压系数分布 Fig. 9 Static pressure coefficient distributions of rotor at 98% height for different slopes

图 10为不同槽坡度情况下叶顶上端壁的总压损失分布云图。总压损失系数为:

$ \omega = \frac{{{p_{{\rm{t}}{\rm{.in}}}} - {p_{\rm{t}}}}}{{{p_{{\rm{t}}{\rm{.in}}}} - {p_{{\rm{in}}}}}} $ (2)
Download:
图 10 压气机动叶上端壁总压损失系数分布 Fig. 10 Contour of shroud surface total pressure loss coefficient distribution

式中:ω为总压损失系数;pt为当地总压。如图 10所示,此氦压气机叶顶的损失主要由2方面组成:一方面为由凹槽产生的漩涡所造成的附加损失;另一方面为由叶顶泄漏所引起的叶顶泄漏涡的损失。凹槽所产生的漩涡主要集中在氦压气机叶片前缘,并且随着坡度的增加漩涡的强度增强。至于叶顶泄漏损失,槽的坡度为30°时所形成的漩涡对叶顶泄漏涡的抑制作用最好,不但推迟了叶顶泄漏涡的形成,而且降低了叶顶泄漏涡的强度,减小了叶顶泄漏涡所造成的损失。

为了更清晰的明确氦压气机叶顶槽坡度对其损失的影响,本文对总压损失系数沿轴向的变化规律进行了研究。如图 11所示,不同槽坡度的氦压气机转子总压损失系数的差异在凹槽处开始显现,随着坡度的增加总压损失系数也增加,并且其增加速度也随着坡度的增加而增加。这主要是叶顶泄漏损失与凹槽附加损失综合作用的结果。

Download:
图 11 压气机总压损失系数沿轴向的分布 Fig. 11 The distribution of the total pressure loss coefficient along the axial direction

图 12为不同坡度条件下,氦压气机的绝热压缩效率随流量的变化情况。如图 12所示,较小的坡度的压气机在所研究的流量范围内都具有最高的绝热压缩效率。在高于95%流量时。不同坡度槽的压气机的效率差异基本不变。而在流量低于95%设计流量时,30°坡度槽的压气机的效率优势更为明显。因此30坡度的氦压气机在较低流量的工况条件下对叶顶间隙泄漏损失的控制能力更强。

Download:
图 12 不同凹槽坡度情况下的绝热压缩效率 Fig. 12 Adiabatic compression efficiency at different slope
2.3 叶顶前缘槽对不同叶顶间隙的氦压球机的性能的影响

为了研究叶顶前缘槽在不同叶顶间隙情况下的作用效果,在凹槽坡度为30°的情况下,选取了0.5、0.75、1 mm 3种间隙尺度的氦压气机转子进行了研究。

图 13为设计点条件下,不同间隙尺度的有槽与无槽的压气机转子叶顶总压损失系数云图。如图 13所示,无论在什么间隙尺度条件下,有槽的叶顶上端壁高损失区都较无槽的小,并且在较大间隙尺度情况下这种变化反应的更为明显。至于凹槽所形成的附加损失由于凹槽的坡度相同所以在不同间隙尺度条件下差异较小,反应到上端壁的总压损失,间隙尺度大的转子叶顶对上端壁的影响较小。

Download:
图 13 压气机转子上端壁总压损失系数分布 Fig. 13 Contour of shroud surface total pressure loss coefficient distribution

图 14为设计工况条件下,不同间隙尺度的叶顶有槽与无槽的压气机转子的绝热压缩效率情况。在间隙尺度为0.5 mm时,有槽叶顶能够使压气机转子的效率提高0.05%;而当间隙的尺度为1 mm时,有槽的压气机转子能够使压气机转子的效率提高0.15%。这说明叶顶前缘凹槽结构在较大的间隙尺度下的控制氦压气机转子的叶顶泄漏的效果较好。

Download:
图 14 不同叶顶间隙时的绝热压缩效率 Fig. 14 Adiabatic compression efficiency at different tip clearances
4 结论

1) 氦压气机转子叶顶单侧槽能够通过形成旋涡减小转子叶片与来流的对冲作用有效地控制泄漏涡的形成,减小叶顶泄漏损失。

2) 30°的叶顶单侧凹槽坡度由于较好的权衡了凹槽的附加损失与其对叶顶泄漏的控制作用使其对压气机叶顶间隙泄漏损失的综合控制效果最好。

3) 叶顶前缘凹槽结构在较大的间隙尺度的条件下的作用效果更为明显。当间隙的尺度为1 mm时,有槽的压气机转子能够使压气机转子的效率提高0.15%。

参考文献
[1]
TIAN Zhitao, ZHENG Qun, JIANG Bin, et al. Research on the design method of highly loaded helium compressor based on the physical properties[J]. Journal of nuclear science and technology, 2017, 54(8): 837-849. (0)
[2]
TIAN Zhitao, ZHENG Qun, MALIK A, et al. Numerical investigation of the effect of highly loaded design on the tip leakage in helium compressor rotors[J]. Progress in nuclear energy, 2018, 105: 263-270. (0)
[3]
SCHLEER M, SONG S J, ABHARI R S. Clearance effects on the onset of instability in a centrifugal compressor[J]. Journal of turbomachinery, 2008, 130(3): 031002. DOI:10.1115/1.2776956 (0)
[4]
黎琨, 单晓明, 吴伟力, 等. 某型压气机失速先兆研究[J]. 测控技术, 2015, 34(7): 46-49.
LI Kun, SHAN Xiaoming, WU Weili, et al. Research on the precursors of compressor stall[J]. Measurement & control technology, 2015, 34(7): 46-49. DOI:10.3969/j.issn.1000-8829.2015.07.013 (0)
[5]
BOTROS K K, HENDERSON J F. Developments in centrifugal compressor surge control-a technology assessment[J]. Journal of turbomachinery, 1994, 116(2): 240-249. DOI:10.1115/1.2928358 (0)
[6]
BOOUGH R M. Advanced controls for airbreathing engines. NASA CR-189205[R]. Washington: NASA, 1993. (0)
[7]
洪林.机匣处理对压气机叶尖间隙流动影响的数值研究[D].沈阳: 沈阳航空航天大学, 2017.
HONG Lin. Numerical study of the effect of casing treatment on compressor tip clearance flow[D]. Shenyang: Shenyang Aerospace University, 2017. http://cdmd.cnki.com.cn/Article/CDMD-10143-1017046585.htm (0)
[8]
HARTMANN M J, BENSER W A, HAUSER C H. Fan and compressor technology. NASA SP-259[R]. Washington: NASA, 1970. (0)
[9]
BAILY E E. Effects of grooved casing treatment on the flow range capability of a single stage axial-flow compressor. NASA TM X-2459[R]. Washington: NASA, 1972. (0)
[10]
刘志伟, 张长生, 时静珣, 等. 关于周向槽机匣处理的若干观测[J]. 西北工业大学学报, 1985, 3(2): 207-216.
LIU Zhiwei, ZHANG Changsheng, SHI Jingxuan, et al. Some observations on circumferentially grooved casing treatment[J]. Journal of Northwestern Polytechnical University, 1985, 3(2): 207-216. (0)
[11]
SCHLEER M, SONG S J, ABHARI R S. Clearance effects on the onset of instability in a centrifugal compressor[J]. Journal of turbomachinery, 2008, 130(3): 031002. (0)
[12]
TIAN Zhitao, ZHENG Qun, JIANG Bin. Effect of Reynolds number on supercritical helium axial compressor rotors performance in closed Brayton cycle[J]. Energy, 2018, 145: 217-227. (0)