1、文章编号：1000-8055（2023）06-1340-10doi：10.13224/ki.jasp.20210652旋转条件下动叶前缘气膜孔排布局影响分析谷萌，谢刚，周志宇，孟龙（北京航空航天大学航空发动机气动热力国家级重点实验室，北京100191）摘要：研究了旋转动叶前缘两侧区域双排孔的孔排布局对气膜冷却特性的影响。采用数值模拟方法获得前缘气膜冷却效率分布以及流动特性。在前缘滞止线每一侧布置两排气膜孔，第 1 排气膜孔与滞止线夹角为10，通过调整第 2 排孔的径向位置和流向位置实现不同的单侧双排孔布局。结果表明：在前缘两侧区域，双排气膜孔布局对气膜冷却特性的影响规律并不相同。在前缘偏吸力面
2、侧区域，孔排布局对相邻孔排的射流掺混特征具有显著影响；而在前缘偏压力面侧区域，双排孔密布的布局形式可以获得更好的气膜冷却效果，这种效应在近孔区域更加明显。吸力面和压力面侧最佳布局的冷却效率提升分别为 0.07 和 0.02。关键词：前缘；气膜冷却效率；旋转叶片；孔排布局；涡轮中图分类号：V231.3文献标志码：AAnalysisonfilmcoolingholearrangementeffectforrotatingbladeleadingedgeGUMeng，XIEGang，ZHOUZhiyu，MENGLong（NationalKeyLaboratoryofScienceandTechnol
3、ogyonAero-EngineAero-thermodynamics，BeihangUniversity，Beijing100191，China）Abstract:The effects of hole arrangement of two rows of film holes on the film cooling per-formanceateachsideofrotatingbladeleadingedgewereinvestigated.Thefilmcoolingperformanceandflowfiledwereobtainedbyusingnumericalsimulatio
4、nmethod.Thefirstrowoffilmholeswasfixedat10downstreamstagnationlineandwasfollowedbythesecondrowofholes.Thespan-wisepositionandthestream-wiselocationofthesecondrowofholeswereadjustedtoachievedifferentholearrangement.Theresultsindicatedthattheeffectsofholearrangementoftworowsonfilmcoolingweredifferenta
5、ttwosidesoftheleadingedge.Attheleadingedgesuctionside，therelativepositionbetweentworowsofholessignificantlyaffectedtheinteractionbetweentheadjacentjetsfromtworows.Atthepressureside，closerdistancebetweentwoadjacentholesintworowsproducedbetterfilmcoolingperformance.Meanwhile，theeffectofrelativepositio
6、nbetweentworowsofholeswasmainlylocatedinthenearholeregion.Thefilmcoolingeffectivenessatthesuctionsideandpressuresideincreasedby0.07and0.02，respectively.Keywords:leadingedge；filmcoolingeffectiveness；rotatingblade；holearrangement；turbine现代燃气轮机为了获得更高的热效率，涡轮进口温度不断被提高，已经超过了涡轮部件的材料耐受极限。前缘区域的热负荷显著高于其余叶片表面1
7、。同时，前缘流场复杂，存在流动滞止、收稿日期：2021-11-15作者简介：谷萌（1983），男，助理研究员，博士生，主要从事航空发动机热端部件高效冷却技术等研究。通信作者：谢刚（1993），男，讲师，博士，主要从事高温旋转部件的流动与换热研究。E-mail：xiegang_引用格式：谷萌,谢刚,周志宇,等.旋转条件下动叶前缘气膜孔排布局影响分析J.航空动力学报,2023,38（6）：1340-1349.GUMeng,XIEGang,ZHOUZhiyu,etal.Analysisonfilmcoolingholearrangementeffectforrotatingbladeleadinge
8、dgeJ.JournalofAerospacePower,2023,38（6）：1340-1349.第38卷第6期航空动力学报Vol.38No.62023年6月JournalofAerospacePowerJune2023曲率大、压力梯度大等特点。因此，前缘采用了多种主动冷却方式。气膜冷却技术就是其中一种有效的方法，它利用从离散孔排出的冷气覆盖涡轮部件表面，通过隔绝高温燃气和降低换热温度的方式保护涡轮部件。目前的前缘气膜冷却研究大多基于简化模型，例如圆柱和半圆前缘结构。同时，这些研究的对象主要集中在气膜孔型、孔间距等结构参数以及主流马赫数、湍流强度等流动参数的影响。早期研究中，Luckey 等
9、2对带有一排斜孔的圆柱体进行换热测量。实验结果表明，孔倾角及其流向位置对气膜冷却有显著影响。Wadia 和 Nealy3利用相似的模型进一步研究了孔倾角和主流参数对气膜冷却的影响。随后，Mick 和 Mayle4测量了具有不对称三排气膜孔布局的半圆前缘气膜冷却效率和相应传热的详细分布，结果表明在大吹风比下，除近孔区域外整个前缘的表面热负荷都有所降低。Ou 等5和 Mehendale 等6详细研究了主流湍流强度对半圆柱前缘气膜冷却特性的影响。在小吹风比下主流湍流的影响更为显著。Cruse等7采用经典的半圆柱模型研究了主流湍流度以及滞止线位置对气膜冷却的影响，后者表现出显著的影响。为了获得更高的气
10、膜冷却性能，前缘气膜冷却一般采用多排孔。以往的研究表明，多排孔之间的射流相互作用对气膜冷却有显著影响。Salcudean 等8、Thakur 等9和 Kim 等10均通过简化前缘模型研究了两排气膜孔在叶高方向上相对位置变化对气膜冷却的影响。Salcudean 等8指出，孔排相对位置的变化对近孔区域有显著影响，顺排布置比交错排布冷却效果更好。Thakur 等9发现，孔排相对位置对第 1 排孔和第 2 排孔的气膜冷却效率有相反的影响。Kim 等10观察到当孔排展向相对位置越近时，孔排间的相互作用越强，冷却效率越高。最近，Anderson 等11在具有多排气膜冷却孔的导叶前缘模型上对叠加方法进行了评
11、估。实验结果表明，由于气膜相互作用并融合，叠加方法低估了气膜冷却效果，特别是在压力侧。由于旋转条件下数据测量难度较大，目前关于旋转条件下动叶的气膜冷却效率数据相对较少。Ahn 等12-13系统研究了主流进气攻角的影响，实验采用内窥镜拍摄和压敏漆（PSP）技术获得了动叶前缘表面的气膜冷却效率数据。结果表明，由于涡轮转速变化带来的主流进气攻角变化显著改变了前缘滞止线位置以及气膜尾迹偏转程度。最近，Li 等14和 Al-Zurfi 等15分别采用实验和数值方法对动叶压力面和吸力面的气膜冷却进行了研究，结果表明两侧区域不同的流场和曲率条件导致了不同的气膜冷却效果。上述研究表明两排气膜孔的相对位置对气膜
12、冷却性能有显著影响，同时这些研究多采用简化的前缘模型来研究这种影响。但是，前缘压力面侧和吸力面侧的流场和曲率明显不同，其气膜冷却特性也不同16-17。因此，本研究的目的是获得典型动叶前缘两侧区域的两排孔相对位置对气膜冷却效率影响的规律，这在公开文献中没有讨论过，因此将有助于前缘冷却设计。考虑到以往的前缘研究大多是在静止状态下进行的，本文的数值研究将在旋转状态下进行并讨论旋转对气膜冷却的影响。1数值计算方法1.1物理模型图 1 展示了用于数值模拟研究的叶片模型，模型通过真实叶片中截面径向拉伸获得，图中 为转速。动叶前缘布置有两排圆柱孔，每排布置有 10 个孔，两排孔均布置在前缘几何滞止线的一侧。
13、因此，两侧（压力面侧或吸力面侧）孔排布局研究的模型是独立的。图 2（a）和图 2（b）展示了前缘吸力面侧（SS）两排孔的布局示意图。压力面侧（PS）也采用了主流二次流(冷气)气膜孔标定点叶顶端壁图1带气膜孔的动叶模型Fig.1Testblademodelwithfilmcoolingconfiguration第6期谷萌等：旋转条件下动叶前缘气膜孔排布局影响分析1341类似的布置，这里不再展示。图中 s 为弧长距离，每个气膜孔的孔径 d 为 0.8mm，气膜孔轴线与流向和展向的夹角分别成 30和 90，其他符号含义见图 2 中标示。如图 2 所示，第 1 排气膜孔位于前缘滞止点下游 1=10处。
14、为了研究两排孔相对位置的影响，第 2 排孔位于前缘滞止点下游 2=25或40处。此外，两排孔在展向（径向）方向上的相对位置设置为 P2/P1=1/4，2/4 和 3/4。前缘的详细结构尺寸如表 1 所示。1.2边界条件图 3 展示了数值模拟的计算域，计算域模型根据课题组以往实验设置18，包括一级导叶流道和一级动叶流道。为了减少计算资源，数值模拟采用周期边界条件。考虑到导叶和转子叶片数量不同，使用了 frozenrotorframechangemodel 将导叶计算域和动叶计算域连接起来。其他边界条件根据之前的实验18进行设置，导叶主流进口采用速度条件，动叶的冷气进口指定固定的质量流量以保证面积
15、平均吹风比。计算采用理想空气为工质，比热比和导热系数采用温度拟合多项式计算，动力黏度采用 Sutherland 公式计算。为了模拟主流和冷气间的密度比 Rd，取 Rd=1.52，主流和冷气的温度分别设置为 460K 和 300K。本 研 究 的 动 叶 流 体 域 旋 转 速 度 设 定 为600r/min。因此，基于叶片旋转半径和转速的旋转雷诺数为 5.38105。相应的动叶计算域进口处的相对主流速度为 10.8m/s。在各个计算条件下，吹风比是基于气膜孔出口冷气平均参数和动叶入口主流的相对参数定义的。计算采用的工况和条件具体如表 2 所示。压力出口速度进口交界面绝热壁面绝热壁面导叶静子盘面
16、质量流量进口旋转周期边界外环静止面动叶图3计算流体域和边界条件Fig.3Computationaldomainsandboundaryconditions表2计算工况Table2Detailedcomputationconditions参数数值主流雷诺数 Reg0.8导叶入口湍流度 Tu/%5动叶进口相对速度 Vrel/（m/s）10.8转速/（r/min）10主流温度 Tg/K460冷气温度 Tc/K300密度比 Rd1.52吹风比 M0.5,1.0,2.0表1前缘结构参数Table1Configurationdimensionsofleadingedge参数数值孔径 d/mm0.8前缘直径
17、 D/d10前缘壁厚/d2.5第 1 排孔流向位置 1/（）10第 2 排孔流向位置 2/（）25,40展向倾角/（）30流向倾角/（）90单排孔间距 P1/d7.875两排孔相对径向间距 P2/P11/4,2/4,3/4叶片高度 H/mm80冷气腔压力面+吸力面s/d=01D2(a)前缘叶尖叶根s/d0dP1P2第1排孔几何滞止线第2排孔冷气(b)吸力面侧气膜孔排布图2前缘和吸力面侧气膜孔排布示意图Fig.2Schematicviewsofleadingedgeandholearrangementinsuctionside1342航空动力学报第38卷需要说明的是，本文的转速和雷诺数与真实发动
18、机的参数有一定差距，但是以往转速对于前缘气膜冷却效率分布的影响研究18表明，转速增加会显著提升气膜冷却效率数值，但是对于孔排间掺混规律影响较小。本文主要研究孔排相对间距对射流掺混的影响，因此相关研究结论对于动叶前缘气膜冷却设计依然具有一定指导意义。1.3数值方法本文数值模拟采用商业软件 ANSYSCFX16.0对稳态 RANS（ReynoldsaveragedNavier-Stoke）方程进行求解，其中湍流模型采用 SST（shearstresstransport）k-模型。求解过程中，所有守恒方程的对流项均采用具有 2 阶进度的高精度格式（highresolutionscheme）进行离散。
19、当计算残差小于 105并且监测的前缘面平均温度在 600 步迭代里保持稳定时，认为计算收敛。计算域采用 ICEM16.0 进行非结构网格划分，如图 4 所示。为了使用 SSTk-湍流模型精确求解近壁面区域的流动换热，所有壁面表面设置有20 层棱柱网格，其中第 1 层网格无量纲高度 y+控S制在 1 附近，边界层网格尺寸按照 1.2 的比例向主流方向增长。网格独立性测试在前缘两侧分别进行，测试采用 3 种密度的动叶计算域网格（3580万、2090 万和 1210 万）。测试选取前缘展向平均气膜冷却效率来评价网格独立性，如图 4 所示，其中横轴表示距离前缘滞止点的无量纲流向弧长 s/d。从图 4
20、中可以看出，在前缘的每一侧，网格数量较多的两个网格计算得到的冷却效率都比较接近，而最稀疏的网格计算获得的结果明显与另外两个有区别。考虑到计算精度和资源，本文数值模拟最终采用 2090 万单元的网格密度离散计算域。为了验证上述数值模拟方法，下面对利用上述方法获得的数值模拟结果与以往的实验数据进行比较18，如图 5 所示。需要说明的是，本文数值研究的叶片模型与参考文献 18 中的叶片模型完全相同，区别仅在气膜孔的布局。从图 5 中可以看出，数值模拟结果和实验数据有相似的值和趋势，特别是在吸力面侧的近孔区域。在压力面侧区域，数值模拟结果与实验数据的最大偏差约为 35%，但是依然能较好地捕捉雷诺数变化
21、带来的冷却效率变化规律。因此，本文数值模拟方法可以用于分析孔排布局对前缘气膜冷却的影响。0.50.40.3s-0.20.11050s/d510压力面侧吸力面侧实验数据,Reg=8.3104实验数据,Reg=16.6104数值模拟结果,Reg=8.3104数值模拟结果,Reg=16.6104图5数值模拟方法验证（M=0.5）Fig.5Validationofnumericalsimulationmethod(M=0.5)2结果分析本章主要讨论动叶前缘单侧两排气膜孔布局，包括流向和展向相对位置对前缘气膜冷却性能的影响。研究分为前缘吸力面侧和压力面侧进行并讨论。此外，讨论中采用数值模拟获得的流场结果
22、更好地理解孔排布局影响。2.1前缘吸力面侧气膜冷却结果图 6 展示了前缘吸力面侧在不同吹风比（M=0.50.40.3s-0.20.11086s/d4203580万网格2090万网格1210万网格第2排孔25第1排孔10(a)吸力面侧0.50.40.3s-0.20.11086s/d4203580万网格2090万网格1210万网格第1排孔10第2排孔25(b)压力面侧图4网格独立性测试结果（P2/P1=2/4）Fig.4Gridindependencetestresults(P2/P1=2/4)第6期谷萌等：旋转条件下动叶前缘气膜孔排布局影响分析13430.5，1.0 和 2.0）和孔排展向布局（
23、P2/P1=1/4，2/4 和3/4）下气膜冷却效率 的分布云图，其中横轴表示距离前缘滞止点的无量纲弧长 s/d，纵轴表示无量纲叶片高度 H/d。第 1 排孔位于前缘滞止点下游1=10处，第 2 排孔位于前缘滞止点下游 2=25处。气膜冷却效率定义如下式所示：=TgTwTgTcTgTcTw其中、和分别是主流温度、冷气温度和壁面温度。从图 6 中可以看出，由于射流的惯性以及主流/射流相互作用，气膜出流后向叶尖偏转。此外，在靠近叶尖区域，气膜尾迹的偏转较大，这是由于前缘滞止线偏转导致的。在叶尖区域，滞止点与气膜孔距离较近，主流速度越小，因此气膜尾迹的偏转量也随之增大。此外，随着吹风比的增大，气膜尾
24、迹的偏转也随之增大，这是因为大吹风比下射流的动量较大。图 6 还表明，不同孔排间射流的相互作用受两排孔展向相对位置（P2/P1）、吹风比（M）以及前缘滞止线位置影响显著。对于 P2/P1=2/4 的基础工况，在吹风比 M=0.5 下，从第 1 排孔出流的气膜与第 2 排孔的气膜尾迹形成完全交错分布，形成均匀的气膜覆盖效果。然而在靠近叶尖区域，两排孔的气膜尾迹发生融合现象，这是由于靠近叶尖区域的射流偏转较大引起的。随着吹风比继续增大，射流脱离壁面现象更加明显，气膜冷却效率显著降低。当吹风比进一步增大至 M=2.0 时，气膜覆盖效果出现明显的恶化，但是在叶高 45H/d50 区域。当第 2 排孔向
25、叶尖方向移动 1/4孔间距，即 P2/P1=3/4 时，从第 1 排出流的气膜尾迹与第 2 排气膜尾迹基本形成交错分布，只在近叶根区域，两排孔射流发生合并现象，这与 P2/P1=2/4 结构下的结果正好相反。为了更好地理解两排孔展向相对位置的影响，图 7 展示了前缘下游 s/d=3.5 处截面上的流向涡量 s分布，其中气膜孔从 1 到 4 进行编号。总的来说，不同排射流之间的相互作用随孔位置和气膜尾迹偏转变化而变化。图 7（a）的第 1 张图展示了吹风比 M=0.5 时，冷气出流后典型的不对称涡对结构。文献 19-20 指出具有涡量较强（逆时针方向）的主涡结构主要由射流-横流剪切形成，而这种涡
26、结构会将主流卷吸到壁面，并导致冷气覆盖效果变差。同时，孔下游也可以观察到涡量较小（顺时针方向）的涡结构，这是由冷却剂的孔内涡对形成的20。图 7 中气膜尾迹位置和上述较907050H/d301012 52s/dM=0.5907050301012 52s/dM=1.0907050301012 52s/dM=2.0907050H/d301012 52s/dM=0.5907050301012 52s/dM=1.0907050301012 52s/dM=2.0907050H/d301012 52s/dM=0.5907050301012 52s/dM=1.0907050301012 52s/dM=2.0
27、0.80.60.40.20.10.30.50.7(a)P2/P1=1/4(b)P2/P1=2/4(c)P2/P1=3/4 图6前缘吸力面侧不同孔排展向距离下气膜冷却效率分布云图Fig.6Contoursoffilmcoolingeffectivenessatvariousspanwiseholearrangementsinsuctionsideofleadingedge1344航空动力学报第38卷弱的涡的位置也是重合的。图 7（a）的第 1 张图表明在 P2/P1=1/4 时，两排孔的展向相对位置较近，不对称涡对发生合并现象，来自孔 4 的主涡将冷气从孔 3 卷吸到壁面，从而在下游区域形成更好
28、的气膜覆盖，获得较高的冷却效率，如图 6 所示。然而，在相同孔间距的 P2/P1=3/4 条件下，图 7（a）的第 3 张图中却没有出现射流合并，其原因有两方面：第 1 排射流偏转角比第 2 排孔大；相邻孔 2和孔 3 的相对位置。图 8 展示了前缘吸力面侧第 2 排孔流向位置改变时气膜冷却效率分布云图，其中第 2 排孔的流向位置从25变化至40。总的来说，当第 2排孔流向位置从25调整到40后，下游区域（s/d4）的气膜覆盖效果变好。然而，两排孔间出现大片无气膜覆盖区域，这对于前缘冷却是不利的。在较小的吹风比（M=0.5）下，当第 2 排孔从25移动至40后，第 1 排孔出流的气膜尾迹明显收
29、缩。随着吹风比增大，两种结构下的第 1 排孔出流气膜尾迹差异逐渐减小。但是，第 2 排孔位于25的结构在下游区域的气膜冷却效果明显好于第 2 排孔位于40的结构，这是由于吹风比变化导致射流偏转发生变化，从而产生类似图 7（b）的第 2 张图中的射流融合现象导致。a图 9 展示了不同孔排布局和吹风比下吸力面侧面积平均气膜冷却效率的结果，面积平均计算时去掉了气膜孔内数据。在第 2 排孔流向位置位于25时，从图中可以看到在第 2 排孔排展向距离 P2/P1=2/4 和 3/4 的结构的气膜冷却效率比较接近，且好于 P2/P1=1/4，特别是在大吹风比下，这与气膜冷却效率云图展示的规律一致，这主要是因
30、为 P2/P1=1/4 的结构发生了明显的射流融合现象，增强了射流非对称涡强度，降低了气膜覆盖效果。从图 9 中还可以看到，当第 2 排孔流向位置向下(a)M=0.5P2/P1=1/4P2/P1=2/4P2/P1=3/4213421342134213421342134(b)M=2.0P2/P1=1/4P2/P1=2/4P2/P1=3/40.80.70.60.50.40.30.20.110001000s/(1/s)截面0图7不同孔排布局下流向位置 s/d=3.5 处涡量分布Fig.7Streamwise-vorticitydistributionsonplanes/d=3.5withvaried
31、spanwiseholearrangements/dM=0.5s/dM=1.0(a)2=25907050H/d301012 52907050301012 52907050301012 5212 5212 5212 52s/dM=2.0(b)2=40s/dM=0.5s/dM=1.0101010907050H/d309070503090705030s/dM=2.00.80.60.40.20.10.30.50.7图8不同孔排流向距离下前缘吸力面侧气膜冷却效率分布云图Fig.8Contoursoffilmcoolingeffectivenessatvariousstreamwiseholearrang
32、ementsinsuctionsideofleadingedge第6期谷萌等：旋转条件下动叶前缘气膜孔排布局影响分析1345游移动到40时，其面积平均气膜冷却效率数值相比于流向位置为25的结构有明显下降。2.2前缘压力面侧气膜冷却结果图 10 展示了前缘压力面侧在不同吹风比（M=0.5，1.0 和 2.0）和孔排展向布局（P2/P1=1/4，2/4 和3/4）下气膜冷却效率的分布云图。其中横轴为无量纲流向距离 s/d，纵轴为无量纲叶片高度 H/d。第 1 排孔位于前缘滞止点下游 1=10处，第 2 排孔位于前缘滞止点下游 2=25处。从图 10 可以看出，压力面侧的气膜冷却效率要低于吸力面侧，
33、且前缘两侧的冷却效率分布特点有很大不同。在小吹风比时，气膜尾迹的偏转相比吸力面侧结果更大。这主要有两个原因：首先，如前所述，倾斜的前缘滞止线导致其在叶根处与气膜孔距离较近；其次，沿着叶高方向，叶尖附近泄漏流的影响越明显，从而增强了射流出流后的偏转。因此，整个压力面侧的气膜轨迹偏转较大。在大吹风比下，从云图中可以观察到展向覆盖较为均匀的缝状气膜分布特点。产生这一现象的原因如下：首先，凹面结构会导致流场不稳定，加剧主流-射流掺混程度；其次，射流出流偏转随着吹风比增加也会增大，相邻射流的绝对间距减小，冷气更容易连成一片；最后，由于大吹风比下冷气质量流量较大，结合前面的两点原因，形成了所谓的冷气积聚现
34、象16，最终呈现出缝状气膜分布特点。对于两排孔展向间距 P2/P1=2/4 工况，当吹风比 M=0.5 时，由于气膜尾迹偏转较强，第 1 排孔的气膜直接流向第 2 排孔处。当第 2 排孔向叶根移动 1/4 孔间距时，即 P2/P1=1/4，第 2 排孔会略微阻挡第 1 排孔气膜的路径，导致近孔区域的冷却效率较高。当第 2 排孔向叶尖方向移动 1/4 孔间距时，即 P2/P1=3/4，两排孔的气膜尾迹出现融合现象，孔下游区域的气膜冷却效果得到有效提升。随着吹风比进一步增加到 1.0，气膜尾迹偏转更加明显。对于 P2/P1=1/4 和 P2/P1=3/4 的结构，近孔区域的气膜覆盖效果更好。而对于
35、 P2/P1=2/4，则出现了较为明显的气膜脱离壁面现象。在吹风比为 2.0 时，三种孔排布局下均出现缝状气膜覆盖现象，其中 P2/P1=1/4 和 P2/P1=3/4 的结构下气膜覆盖效果相对较好。为了更好地理解孔排展向相对间距的影响，0.300.250.200.150.10a-0.050M=0.5M=1.0M=2.0P2/P1=1/42=25P2/P1=2/42=25P2/P1=3/42=25P2/P1=2/42=400.1850.1910.1910.1640.1780.1930.2020.1810.1140.1540.1440.0819图9前缘吸力面侧面积平均气膜冷却效率Fig.9Are
36、a-averagedfilmcoolingeffectivenessinsuctionsideofleadingedge907050H/d30102512s/dM=0.590705030102512s/dM=1.0(a)P2/P1=1/4(b)P2/P1=2/4(c)P2/P1=3/490705030102512s/dM=1.5907050H/d30102512s/dM=0.590705030102512s/dM=1.090705030102512s/dM=1.5907050H/d30102512s/dM=0.590705030102512s/dM=1.090705030102512s/dM=
37、1.50.80.60.40.20.10.30.50.7图10前缘压力面侧不同孔排展向距离下气膜冷却效率分布云图Fig.10Contoursoffilmcoolingeffectivenessatvariousspanwiseholearrangementsinpressuresideofleadingedge1346航空动力学报第38卷压力面侧孔排下游 s/d=3.5 和 s/d=7 处的流场截面被提取出来，如图 11 所示。图 11 中展示了截面上无量纲温度=（TgT）/（TgTc）的分布，其中 T 为温度。从图 11 中可以看出在小吹风比下，孔排间射流的掺混规律依然与吸力面类似，孔排径向间
38、距变化会产生不同的射流掺混特点，例如射流融合和交错分布等。随着吹风比增大到 1.5，如图 11（b）所示，冷气在近孔区域可以较好地贴附壁面。对比不同孔排结构还可以看出，P2/P1=1/4和 P2/P1=3/4 的结构射流贴壁效果更好。本文认为，这主要因为更近的孔排相对间距更有利于相邻射流形成冷气积聚现象，这可以将主流排除在外，获得更好的气膜展向覆盖效果。图 12 展示了前缘压力面侧第 2 排孔流向位置改变时气膜冷却效率分布云图，其中第 2 排孔的流向位置从 25变化至 40。在展向方向上，第1 排孔与第 2 排孔完全交错分布，即 P2/P1=2/4。总的来说，改变第 2 排孔流向位置对气膜冷却
39、效果的影响主要集中在近孔区域，两排孔流向间距越小，气膜冷却效果也越好。吹风比为 0.5的气膜冷却效率分布表明，当第 2 排孔流向位置从 25移动到 40时，两排孔出流的气膜尾迹出现明显的尾迹重叠现象。当吹风比增加到 1.0，气膜尾迹偏转角度增大，第 2 排孔移动到下游 40时两排孔的气膜尾迹出现交错排布的效果。因此，其孔排下游气膜冷却效率明显高于第 2 排位于下游 25的结构。当吹风比进一步增大到 2.0，第 2排孔移动到下游 40结构的冷气积聚现象并不明显。这可归因于当两排孔间距较大时，冷气难以形成积聚的效果。(a)M=0.75P2/P1=1/4P2/P1=2/4P2/P1=3/421212
40、1(b)M=1.5P2/P1=1/4P2/P1=2/4P2/P1=3/40.80.70.60.50.40.30.20.10.80.70.60.50.40.30.20.1截面图11不同孔排布局下流向位置 s/d=3.5 和 7 截面处无量纲温度分布Fig.11Non-dimensionaltemperaturedistributionsonplanes/d=3.5and7withvariedspanwiseholearrangement(b)2=40907050H/d30102512s/dM=0.590705030102512s/dM=1.090705030102512s/dM=2.090705
41、0H/d30102512s/dM=0.59070503010212s/dM=1.0(a)2=2590705030102512s/dM=2.00.80.60.40.20.10.30.50.75图12不同孔排流向距离下前缘压力面侧气膜冷却效率分布云图Fig.12Contoursoffilmcoolingeffectivenessatvariousstreamwiseholearrangementsinpressuresideofleadingedge第6期谷萌等：旋转条件下动叶前缘气膜孔排布局影响分析1347 a图 13 展示了不同孔排布局和吹风比下压力面侧面积平均气膜冷却效率的结果，面积平均计算
42、时去掉了气膜孔内数据。从图中可以看到，在大吹风比下，在第 2 排孔流向位置位于 25时，第 2 排孔排与第 1 排孔展向距离更近的结构（即P2/P1=1/4 和 P2/P1=3/4）的面积平均气膜冷却效率更高。此外，从图 13 中还可以看到，当第 2 排孔流向位置向下游移动到 40时，其面积平均气膜冷却效率数值相比于流向位置为 25的结构也有明显下降。0.300.250.200.150.10a-0.050M=0.5M=1.0M=2.0P2/P1=1/42=25P2/P1=2/42=25P2/P1=3/42=25P2/P1=2/42=400.1590.1940.2080.1590.1540.14
43、50.2040.1650.1520.2060.1660.181图13前缘压力面侧面积平均气膜冷却效率Fig.13Area-averagedfilmcoolingeffectivenessinpressuresideofleadingedge3结论本文研究了旋转叶片前缘两排气膜孔相对位置对气膜冷却性能的影响，其中第 1 排气膜孔固定在前缘滞止点下游10处。第 2 排孔与第 1 排孔在叶高方向交错布局，两排孔在叶高方向上的相对位置为 P2/P1=1/4，2/4 和 3/4，同时第 2 排孔在流向上的位置也不同，分别位于前缘滞止点下游25和40处。主要结论如下：1）由于泄漏流的存在和滞止线的位置，前
44、缘压力面侧的气膜尾迹偏转较大。同时，前缘吸力面侧的气膜冷却效率明显高于压力面侧。2）在吸力面侧，当两排孔完全交错布置时，气膜覆盖均匀均匀性最好。当第 2 排孔向叶尖方向移动时，两排孔间的射流相互作用发生了显著的变化，导致气膜孔下游的冷却效率较高，但是无气膜覆盖的区域较大。当第 2 排孔向下游移动时，孔排下游区域的冷却效率显著提升，但是两排间的气膜覆盖效果则有所下降。3）在压力面侧，更近的两排孔间距可以产生更好的冷却性能。在小吹风比下，气膜冷却效率受两排孔射流相互作用影响较为明显，类似于吸力面侧。在较大的吹风比下，冷气出流后会出现积聚现象，使气膜覆盖均匀，冷却效率提高。而两排孔间距越近，更有利于
45、冷气的积聚。4）以上研究表明动叶前缘两侧的气膜孔排布局需要采用不同的设计方法。在前缘压力面侧，通过合理调整孔排位置（如孔排相对距离和孔出口角等）和流动条件（吹风比和滞止线位置等），可以使得两排孔下游的气膜冷却效率分布更加均匀。对于前缘压力侧，在大吹风比下，两排孔之间的相对间距越近会带来较好的气膜覆盖效果，但是还需要进一步考虑密集布局带来的热应力影响，后者也是涡轮动叶设计中十分重要的参数。参考文献：HANJC，DUTTAS，EKKADSV.GasturbineheattransferandcoolingtechnologyM.BocaRaton，US:CRCPress，2012.1LUCKEYD
46、W，WINSTANLEYDK，HANUSGJ，etal.Stagna-tionregiongasfilmcoolingforturbinebladeleading-edgeapplica-tionsJ.AIAAJournalofAircraft，1977，14:494-501.2WADIAAR，NEALYDA.Experimentalsimulationofturbineairfoil leading edge film coolingJ.Journal of Turbomachinery，1988，110(2):226-232.3MICKWJ，MAYLERE.Stagnationfilmco
47、olingandheattransfer,includingitseffectwithintheholepatternJ.JournalofTurboma-chinery，1988，110(1):66-72.4OU S，MEHENDALE A B，HAN J C.Influence of high main-streamturbulenceonleadingedgefilmcoolingheattransfer:effectoffilmholerowlocationJ.JournalofTurbomachinery，1990，114(4):716-723.5MEHENDALEAB，HANJC.
48、Influenceofhighmainstreamtur-bulenceonleadingedgefilmcoolingheattransferJ.JournalofTurbomachinery，1992，114(4):707-715.6CRUSEMW，YUKIUM，BOGARDDG.Investigationofvari-ousparametricinfluencesonleadingedgefilmcoolingR.Orlan-do，US:ASMEInternationalGasTurbineandAeroengineCongr-essandExhibition，1997.7SALCUDE
49、AN M，GARTSHORE I，ZHANG K，et al.Leadingedgefilmcoolingofaturbineblademodelthroughsingleanddou-blerowinjection:effectsofcoolantdensityR.Hague，Netherlands:ASME International Gas Turbine and Aeroengine Congress andExposition，1994.8THAKURS，WRIGHTJ，SHYYW.Convectivefilmcoolingoverarepresentativeturbineblad
50、eleading-edgeJ.InternationalJour-nalofHeatandMassTransfer，1999，42(12):2269-2285.9KIMKS，KIMYJ，KIMSM.Enhancementoffilmcoolingper-formanceattheleadingedgeofturbinebladeR.Barcelona，Spain:ASMETurboExpo:PowerforLand，Sea，andAir，2006.10ANDERSONJB，WINKAJR，BOGARDDG，etal.Evaluationofsuperpositionpredictionsfor
