引言
实验系统和方法
基于叶片弦长的叶栅进口雷诺数Re、出口马赫数Ma、质量流量比MFR定义如下
(1) |
(2) |
(3) |
式(1),(2),(3)中,ρ1,V1,μ1分别为叶栅入口的气流密度、速度和动力粘性系数;C为叶片弦长;p1t为叶栅前总压;p2为叶栅后静压;k为绝热等熵指数,取k=1.4;mc, mg分别是二次流质量流量和单个叶栅通道内的主流质量流量;mc1, mc2分别为叶片前腔和后腔的二次流质量流量。具体实验工况如表1所示。
1.4 数据处理和实验误差
本实验基于一维半无限大假设条件,采用瞬态测温技术测量叶片表面气膜冷却效率和换热系数,其定义如下
(4) |
(5) |
有气膜冷却时,叶片表面气膜冷却效率和换热系数的计算公式为
(6) |
式(4),(5),(6)中,qw为叶片表面热流量,Tgr为叶片表面的当地主流恢复温度,Tw为壁面温度,Taw为绝热壁温,Tc为二次流温度,qw可通过“脉冲响应法”求得。具体的测量原理和方法可参考文献[20]。
可以看出式(6)为一次线性方程,可写为以下形式
(7) |
式(7)中的斜率a对应表面换热系数h,与x轴的截距b对应气膜冷却效率η。其中,系数b通过系数a和系数c共同求解,均是进行线性拟合的待求量
(8) |
y的不确定度为
(9) |
系数a和c的不确定度为
(10) |
系数b的不确定度同样由a和c的不确定度求出,即
(11) |
结合式(9)和式(11)可以看出,实验数据点数n越多和越分散,拟合系数a的误差就越小,因此实验时采用高频(10kHz)进行数据采集,保证了在较短的实验时间内采集较多的数据。另外为了提高实验结果准确性,每个工况都分别在实验叶片常温和加热的情况下进行一次实验,保证了数据点的分散度。更详细的误差分析可参考文献[22],实验各参数的相对误差值见表2。
实验结果分析
2.1 光叶片换热结果
Fig. 4 Effect of Re on Nussult number
(12) |
式中h是叶片表面的换热系数,是进口空气的导热系数。
从图4中可以看出进口雷诺数的提高显著增强了叶片表面的换热。在压力面,三个雷诺数下Nu数都有一个从S/C=0处沿流向逐渐下降的过程,这是由于边界层逐渐变厚导致热阻变大;Re=6.4×105,9.0×105时,Nu数在S/C=-0.3附近位置明显上升,表明边界层发生了层流到湍流的转捩,在S/C<-0.4的区域湍流边界层逐渐发展导致换热沿流向逐渐增强。在吸力面,Re=3.0×105时的主流边界层在S/C=0.4附近区域发生了层流到湍流的转捩,Re=6.4×105,9.0×105时,边界层在前缘附近即发生了转捩。
本文的换热实验结果也以St数的形式与平板上层流边界层和湍流边界层St数的经验关联式进行了对比,以有助于更好地理解叶片表面边界层的状态。为便于表述,层流边界层的St数与当地Re数的关系用“层流经验式”表示,湍流边界层的用“湍流经验式”表示。
当地Re数和St数的定义式如下[23]
(13) |
(14) |
Bergman等[24]给出了当地Nu数的关系式
层流边界层
(15) |
湍流边界层
(16) |
式(13),(14),(15),(16)中,下标L代表当地参数,Pr是普朗特数。
图5是叶片表面St数和经验关联式的对比。对于吸力面而言,Re=3.0×105时的数据点最初保持层流经验式的下降斜率发展,与层流经验式不重合可能是由于压力梯度和主流湍流度造成的,在转捩点后数据点向湍流经验式移动。对于Re=6.4×105,9.0×105的工况,最开始的3个数据点下降速率和层流经验式一致,之后的数据点开始向湍流经验式移动,与Re=3.0×105不同的是St数没有一直增大,而是随着湍流经验式的趋势发展了一段长度,最后超过了湍流经验式,达到了旺盛湍流的状态,这和Nasir等[25]在涡轮导叶上得到的规律是明显不同的,表明对于高速来流工况,继续升高雷诺数对不带气膜孔叶片吸力面换热系数的发展具有显著的影响。
Fig. 5 St number compared to the flat plate correlations
2.2 质量流量比对气膜冷却效率的影响
Fig. 6 Effect of MFR on film cooling effectiveness
对于吸力面S/C>0.2的区域,MFR=5.5%时雷诺数对气膜冷却效率没有影响,这主要是由于MFR较小时,二次流对边界层的扰动较小导致二次流和主流的掺混耗散过程较弱,因此叶片表面的冷却效率取决于二次流流量和主流流量的比例,当MFR相同时气膜冷却效率也就基本相同。而对于中高MFR的工况来说,此时二次流由于动量较高对涡轮叶片表面的主流流动状态产生了较大的影响,根据前文分析可知高MFR时二次流发生了脱离壁面的现象,与主流掺混较为严重,而且由于本实验的主流湍流度较高,因此二次流很容易耗散到主流中去,对于相同MFR工况,雷诺数较高时的二次流绝对流量较多,抵抗耗散的能力较强,因此雷诺数越高气膜冷却效率越高。
而对于压力面S/C<-0.6的区域,MFR=5.5%和MFR=8.4%时,气膜冷却效率随着雷诺数的增大先降低后升高,MFR=12.5%时气膜冷却效率随着雷诺数增大一直降低。-0.6<S/C<-0.2区域内由于气膜孔较多导致其表面都处于湍流状态,因此雷诺数的影响规律更加复杂。与Fu等[27]的实验结果相比,本实验中MFR对气膜冷却效率的影响较小,这可能是由于本实验采用的涡轮叶片在压力面布置了多排孔。以上实验结果表明雷诺数对气膜冷却效率的影响规律是区域性的而且与MFR有关,整体来看在高MFR时雷诺数的变化对压力面的尾缘和吸力面具有较大的影响,因此涡轮进口雷诺数变化时应当适当调整二次流流量以保证气膜冷却效率不发生很大变化。
2.4 质量流量比对换热系数比的影响
图8给出了不同主流雷诺数工况下MFR对涡轮叶片表面换热系数比的影响,图中h为全气膜出流时叶片表面的换热系数,h0为相同主流条件下光叶片表面的换热系数。由于两个实验叶片的热电偶位置不完全对应,因此在数据处理时对光叶片表面的换热系数做插值处理。
从图中可以明显看出前缘(-0.1<S/C<0.1)的换热系数比整体较高,且和MFR有明显的正相关关系,最高值达到了4.5左右。根据光叶片的换热系数分布可知主流在前缘滞止后沿吸力面和压力面都经历了层流边界层的发展过程,边界层的存在会增大流体的传热热阻,而前缘区域布置了多排气膜孔,气膜孔密度很大,大量二次流出流打破了前缘附近的层流边界层状态并极大地扰动了前缘附近的流体运动,因此显著增强了换热。值得注意的是对于S/C=0的驻点区域,其换热系数比相对于其他前缘区域又相对较低,这是因为光叶片前缘驻点的换热本身就比较高,同时前缘气膜孔的二次流出流大都是向吸力面和压力面流动,对前缘驻点的扰动相对来说比较小。
在吸力面S/C>0.2的区域,MFR对换热系数比的影响规律和雷诺数有较大的关系,在Re=3.0×105和Re=6.4×105时,MFR对换热系数比的影响很小。而在Re=9.0×105时,MFR越高换热系数比越大。Re=3.0×105时,二次流出流导致S/C>0.2区域的换热系数比处于1.0和2.0之间且在气膜孔出口附近较高,造成这个现象的原因主要是由图4和图5可知小雷诺数时光叶片吸力面前半段S/C<0.4的区域为层流状态,二次流出流对边界层的扰动增强了壁面附近的湍流程度因此增强了换热。而换热系数比沿着流动方向逐渐降低,这有两方面的原因,一是二次流出流造成扰动逐渐被主流吸收,二是光叶片的换热沿主流流动方向是逐渐增强的。Re=6.4×105时,0.3<S/C<0.5的区域内换热系数比小于1,这主要是因为Re=6.4×105时,光叶片这个区域的边界层较薄,二次流出流在一定程度上加厚了边界层,导致换热系数降低,Eriksen等[3]在平板上,Fu等[27]和Newman等[28]在叶片吸力面上也发现了类似现象。
在雷诺数Re=9.0×105时,0.2<S/C<0.4区域内的换热系数比也小于1.0,这和Re=6.4×105是一致的,但是换热系数比受MFR影响明显,这应该是由于不同雷诺数时虽然MFR相同,但是气膜孔的实际出气量随着雷诺数的增大而增大,吸力面由于曲率大和具有由主流指向叶片的正压力梯度这两个特点,使得二次流会发生脱离壁面再贴附的现象[29],中低雷诺数时,不同MFR时的二次流更多地贴附到了壁面上,对主流扰动作用基本相同,而高雷诺数工况时,不同MFR时有较多的二次流在壁面附近与主流发生了掺混,较高的MFR对壁面附近流体的扰动较强,最终导致高雷诺数时吸力面的换热系数比整体上随着MFR增大而升高。值得注意的是,吸力面S/C=0.41的位置处换热系数比有极大值,这可能是由于吸力面二次流重新贴附壁面后对壁面附近流体产生较强的扰动引起的。
整体来看,压力面的换热系数比相比吸力面较高,每个气膜孔出口附近的换热系数比都比较高,这是由于其表面布置了多排气膜孔,二次流出流的叠加作用使得整个压力面附近的流动都处于湍流度比较高的状态,MFR越高对换热的增强作用就越明显。MFR相同时,叶片表面的换热系数比随着雷诺数升高逐渐降低,这主要是由于雷诺数升高使得光叶片的换热系数也明显升高,因此在对涡轮叶片压力面中间区域(-0.6<S/C<-0.2)进行气膜孔设计时,要重点考虑MFR对换热系数的影响来整体评估冷却效果,而对于压力面后半段(S/C<-0.6)以及吸力面,要重点考虑MFR对气膜冷却效率的影响。
结论
参考资料:
· end ·