郝鑫,连文磊,尚登阅
(南京航空航天大学,能源与动力学院,江苏 南京 210016)
预旋系统作为航空发动机气流冷却系统的重要组成部分,极大地影响了涡轮叶片冷却气流的冷却效果。如何优化预旋系统的结构设计,进一步对冷却气流进行优化利用,更加有效地满足涡轮叶片的冷却要求是涡轮叶片冷却技术后续发展的一个重要研究方向。
KILIC M[1]利用Navier-Stokes方程对盘腔内结构下的流动和传热特性进行了数值模拟。SNOWSILL G D[2]分析了转静交界界面的处理方式对数值计算结果的影响,证明Frozen Rotor处理方法耗时更少,并可以取得很好的结果。肖竞雄等[3]研究了径向进气装置进气孔数量、尺寸的变化,总压恢复系数与畸变指数的变化规律。王锁芳等针对随着盘腔进气径向位置[4]、预旋喷嘴的径向角度[5]对气流旋流比、径向预旋系统温降系数及总压损失系数变化影响规律进行了研究。
以上及其他相关研究大多集中在盘腔传热的基础研究上,对涡轮盘腔内部结构对流动的影响涉及得较少。本文将针对涡轮盘腔内部对冷却气流温降增压效果的影响展开研究,并分析影响温降增压效果的主要因素、相应的评价标准以及各影响因素的影响规律。
本文根据实际航空发动机涡轮盘腔的实际设计尺寸,进行了涡轮盘腔的几何模型构建如图1所示,并根据涡轮引气系统实际内部流场,对模型进行了适当的简化,构建了如图2所示的涡轮盘腔内部的流体域计算模型。由于该模型具有周期性,选取了模型的1/7作为计算域。
图1 涡轮盘腔的几何模型
图2 涡轮盘腔内部流体域
本文研究的内容是稳态流动,各物理量可以取稳态下的时均值。计算模型的边界条件由实际航空发动机的实际工况来确定,预旋喷嘴的进口设置为压力进口边界,给定进口总压、总温和出口压力,出口设置为流量出口。转静交界面采用Frozen Rotor法处理,采用RNG 湍流模型作为计算模型,并选用CFX软件进行计算。壁面的y plus在30~300之间;入口处的湍流强度为5%,涡流黏度比为10,密度、动量、湍动能、湍动能耗散率以及能量的离散格式均为1阶迎风,压力-速度耦合采用SIMPLE算法。收敛标准为各项残差均<10-4,出口相对总温的残差变化趋于稳定。在相关文献中也利用试验结果对该模型进行了验证[6],证明了计算方法的可行性(图3)[7]。
图3 温降计算值与试验值的对比
由于涡轮叶片冷却气流主要服务于高速旋转的涡轮叶片的冷却工作,旋转系下的相对总温对其实际工作效果更具有参考意义[8],故选定无量纲温降θ为温降效果的参考量,同时气流压比也影响着冷却气流的冷却效果。上述性能参数的定义如下:
1)无量纲温降
无量纲温降θ是有关进口总温和旋转坐标系下出口相对总温的无量纲参量,定义如下:
(1)
(2)
2)压比
在绝热的旋转盘腔中,静压增加,相对总压保持不变,总压变化等于动压变化和静压的变化。
π=Pout/Pnozzle
(3)
3)旋流比
为研究关键性能参数的影响变化规律,现选取旋流比为设计过程中的无量纲设计参量,旋流比的定义如下:
(4)
式中:Vφ为周向速度,m/s;ω为轴向角速度,rad/s。将上式代入相对总温的定义式(2),则有
(5)
对预旋系统来说其静温是固定不变的,所以当旋流比β接近于1的时候,相对总温最低。
针对涡轮盘腔内叶片设计参数对冷却气流的影响规律,本节针对叶型设计、叶片长度以及叶片构型等主要参数展开分析,同时选择旋流比为设计过程中的主要无量纲参考量。
根据叶片的流动规律和实际加工制造的可行性,初步选定叶片进口位置为距涡轮盘腔中心半径100mm处,出口安装位置为半径130mm处。并基于直叶片、弯曲叶片和弯扭叶片设计了3种基本叶型,各叶型下的盘腔内部气流流场如图4、图5所示。
图4 不同叶型下的盘腔内部的流场(旋转系下)
图5 不同叶型典型工况下旋流比分布图
由流场图可以看出,相对于空盘,添加叶片的各盘腔出口处的周向相对速度均显著降低。这是由于导流叶片使得盘腔内部气流的周向速度上升,与旋转盘腔的相对速度差逐步减小,从而实现了相对速度的降低,进一步使相对总温下降。其中弯扭叶片的流动分离区域最小,所造成的压力损失相对较小,结合叶片所带来的增压效果较好。
对比不同叶型的旋流比分布图可以发现,布置导流增压叶片对整个盘腔内部旋流比分布均会产生影响。空盘内部旋流比分布呈层状分布,受半径位置旋流比逐步下降的影响。叶片的存在改变了旋流比变化规律,明显提升了高半径出口位置的旋流比,并使盘腔内部大部分区域的旋流比接近于1,从而提升静压和总温降,提升了导流系统总体性能。
图6、图7为典型工况下3种基础叶型的涡轮盘腔与空盘的无量纲温降及压比变化曲线图。3种基础叶型的无量纲温降、压比相差不大,但均远好于空盘,弯扭叶型的温降效果最为优秀,计算结果与流场及旋流比分析的预估结果一致。
图6 典型工况下各叶型无量纲温降
图7 典型工况下各叶型压比
为避免因长度变化带来的叶片叶型变化所造成的影响,对比不同叶片数时均选取直叶片进行数值计算,选择典型工况作为计算工况,盘腔内叶片数均选为7片,选择叶片的出口半径为130mm,叶片的进口半径分别取为80mm、90mm、100mm、110mm。图8、图9为不同叶片长度与不同叶片旋流比对流场及分布的影响。
图8 不同叶片长度下的盘腔内部流场(旋转系下)
图9 不同叶片长度的旋流比分布图
图8中叶片长度越短,出口处的流动分离区影响范围越大。增加叶片长度主要影响了低半径处的气流流动方向,使得气流提前偏转,强化了对气流的增压效果,同时随着叶片长度缩短,流动分离尾迹区逐步增大,造成出口位置的流动不均匀,会产生较大的压力损失。
从图9中可以看出随着叶片的长度的增加,高旋流比分离区逐步减小,盘腔内部大部分区域旋流比均更为接近于1。总体来看,盘腔内部旋流比变化较大的区域主要集中在叶片附近,其他区域旋流比变化不大,加长叶片主要对底半径处叶片附近区域的旋流比产生了影响。
图10、图11中的数值计算结果与流场及旋流比分析结果一致。虽然不同叶片长度长度间的无量纲温降相差较小,但是随着叶片长度增加,无量纲总温降依旧有着升高的趋势。因此,可以认为叶片在低半径处可以提高系统的温降效果,但是总体来说提升力有限,可以认为加长叶片长度能够在一定范围内提升温降增压效果。
图10 不同叶片长度无量纲温降变化
图11 不同叶片长度压比变化
本文通过对空气预旋系统涡轮盘腔进行数值模拟,对其内部结构采用添加叶片的方式进行优化设计,并进一步研究了叶型设计、叶片长度、叶片数对冷却气流的温降增压效果的影响,得出了以下结论:
1)涡轮盘腔的作用机理是通过导流叶片的导流增压效果使得盘腔内部气流的周向速度升高,从而使得盘腔旋转系下气流的相对速度下降,达到降低相对总温的效果,同时可提高气流的压比。优化效果的主要影响因素是气流旋流比,全局流场旋流比均接近1时,优化效果最为明显。
2)叶型设计对涡轮盘腔优化效果有着显著影响。在不同工况下,布置叶片均能对温降增压效果产生收益,其中弯扭叶型带来的收益最高。
3)叶片长度对温降增压效果影响较大,加长叶片长度能有效增强对冷气流的导流增压效果,但随叶片加长而带来的收益会逐渐降低。