涡轮叶片尾缘冷却特性数值研究

为了研究不同叶片尾缘结构对气膜冷却效果的影响规律,设计了三种尾缘结构,采用数值模拟的方法分别对其冷却效果进行了研究。研究结果表明:(1)三种尾缘结构的冷却效率沿壁面的分布有很大差异,针对模型Ⅰ,在小吹风比时,冷却效率沿尾缘壁面逐渐减小;而在大吹风比时,冷却效率呈现逐渐升高的趋势;(2)模型Ⅱ和模型Ⅲ的冷却效率沿壁面均呈现逐渐降低的趋势,但二者降低的规律又不相同;(3)在相同壁面位置,模型Ⅲ的冷却效率最高,模型Ⅰ的冷却效率最低,因此可以认为,模型Ⅲ所示的尾缘结构更有利于叶片尾缘的冷却。     

扰流柱对叶片尾缘对流换热特性的影响

为了研究扰流柱对涡轮叶片尾缘对流换热特性的影响规律,设计了4种不同结构的扰流柱,采用红外热像仪分别对不同形状和排列方式的扰流柱对叶片尾缘气膜冷却特性的影响规律进行了实验研究.研究结果表明,在吹风比较低时,水滴I型扰流柱叶片尾缘的对流换热特性较好;吹风比较大时,圆柱型扰流柱叶片尾缘的对流换热特性较好,水滴Ⅱ型扰流柱叶片尾缘的对流换热系数最小.      

流量分配对尾缘通道流动换热特性的影响

采用数值模拟方法,研究了不同进气比条件下双向进气叶片尾缘通道的流动与换热特性。结果表明,整个尾缘通道的全局平均换热强度随着进气比的增加而增大。小进气比情况下,尾缝出口流量沿径向分布较均匀,整个通道的对流换热都较弱,其中中间隔板附近对流换热最弱。大进气比情况下,尾缝出口流量分布沿径向变化较大,在通道顶部区域分布较多,通道顶部区域对流换热最强,中间隔板附近对流换热有所改善。     

下缘板孔对涡轮叶片尾缘内冷通道流动换热影响的数值研究

涡轮叶片下缘板出气孔对内冷通道的流动和换热性能有较大影响。通过数值模拟方法研究下缘板出 气孔对尾缘和下缘板双路出气涡轮叶片尾缘内冷通道内的流动和换热特性,对比分析孔径、孔形和孔位置对尾 缘溢流孔流量系数、尾缘出流比、尾缘通道内总压系数和尾缘内冷腔壁面换热特性的影响。结果表明:下缘板 出气孔孔径对流量系数分布的影响显著,孔径增大,尾缘溢流孔流量系数下降,尾缘出流比减小,尾缘内冷通道 内压力损失降低,内冷腔平均换热系数增大;孔形对上游内冷通道内流动和换热几乎没有影响;孔位置变化对 内冷通道壁面整体的换热系数影响很小,对局部影响较大。     

气膜-发散组合冷却结构换热特性的实验研究

为了研究气膜-发散组合冷却结构的冷却特征,保证相同的开孔率,设计了三种不同发散孔排布形式的组合冷却结构,采用实验的方法对气动参数和几何参数对绝热冷却效率和对流换热系数的影响规律开展了研究。结果表明:绝热冷却效率和对流换热系数沿主流方向先逐渐降低,达到最低点后沿流动方向二者基本保持不变;在研究参数范围内,主流雷诺数和吹风比对绝热冷却效率的影响不大,但对组合冷却结构的对流换热系数影响较大,随着主流雷诺数和吹风比的增加,对流换热系数均呈现逐渐增大的趋势;针对三种发散孔排布形式的绝热冷却效率和对流换热系数,流向间距大的气膜发散冷却结构最高,流向间距居中的气膜发散冷却结构次之,流向间距最小的气膜发散冷却结构最低。     

梯形突片气膜冷却特性的实验

为了研究突片对气膜冷却的影响规律,设计了3种不同堵塞比的梯形突片,采用红外摄像仪对不同形状突片冷却结构的壁面温度场进行了测试,并对气膜冷却特性的影响规律进行了实验研究.研究结果表明,与无突片的气膜冷却相比,突片的存在大大提高了气膜冷却效率和对流换热系数;冷却效率沿气流方向存在一个最佳吹风比值;对流换热系数在本实验范围内,堵塞比为0.214时气膜冷却对流换热系数最大.      

光滑楔形通道径向末端射流冲击冷却换热特性研究

为了研究在旋转状态下射流冲击冷却在楔形通道径向末端的冷却效果,在主流通道进口雷诺数为15000、旋转数为0或0.1的工况下,实验研究了不同射流孔位置、射流流量条件下的楔形通道沿程壁面换热规律。研究结果表明:射流孔附近的通道壁面换热系数被极大地增强,但冷却效果被限制在一定范围内;射流对通道内侧的影响范围要大于中部和外侧,但增大射流流量并不会增加内侧的换热系数;旋转工况下,射流影响区的换热对比静止基本上没有变化;在本实验工况下,射流不能强化沿通道径向射流孔后侧位置的换热;不同工况下的通道整体平均换热在10%以内,末端射流冲击对通道整体换热影响有限。     

基于磁共振测速的复合冷却涡轮叶片流动分析

针对传统光学测量技术难以完整地测量涡轮叶片复杂冷却结构内部三维流场的问题，采用磁共振测速（MRV）技术测量了一种复合冷却涡轮叶片内部的三维流场，并重点研究了叶片尾缘的流动特征。在通过MRV数据验证计算流体力学（CFD）准确性和可靠性的基础上，进一步通过CFD研究了全高、半高扰流柱对尾缘弦向出流的影响。结果表明，MRV成功地捕捉到了复合冷却涡轮叶片内部流动特征（如扰流柱和半劈缝处的漩涡结构等）。CFD与MRV数据在尾缘处有一致的出流趋势：气膜孔出流量沿展向逐渐减小，而半劈缝出流量沿展向逐渐增大。通过不同扰流柱设计的尾缘CFD对比发现，全高、半高扰流柱对尾缘流动的影响主要在于增大了流阻和出流的驱动压差，进而改变了出流情况，使得气膜孔整体出流量增大约2.8%，半劈缝整体出流量减小约2.8%。     

某涡轮导向叶片换热实验与计算

针对某涡轮导向叶片,实验测量了光滑叶片表面的压力系数和速比系数,并使用瞬态液晶测量技术获得了叶片全表面传热系数分布.分别使用shear stress transport（SST）,k-ω,k-ε和renormalization group（RNG）k-ε四种湍流模型模拟了相同结构尺寸的叶栅通道内的流动与换热,并与实验结果进行对比.结果表明:压力面压力系数沿弧长方向逐渐下降,吸力面上压力系数先快速下降达到最小值后缓慢上升（出现逆压梯度）.叶栅通道和叶片表面附近气流流动结构的复杂性导致叶片表面传热系数分布较为复杂.4种湍流模型对压力系数和速比系数的计算结果相互差别不大,计算数据也比较接近实验值.关于叶片表面传热系数,SST模型计算结果分布规律与实验接近,而其他3种湍流模型都没有能模拟出吸力面边界层分离对换热的影响.      

小间距梯形扰流柱通道内的流动换热数值计算

为了深入了解涡轮叶片尾缘扰流柱区域的流动结构和换热情况,通过数值计算研究了小间距梯形扰流柱通道内的流动和换热.数值计算结果表明:①数值计算得到的端壁静压分布、弦向出流量分布、总压系数与实验结果符合良好,端壁换热分布规律总体上与实验结果一致.②扰流柱区域的流动方向基本为弦向.沿弦向湍流度是逐渐升高的,靠近弦向出口扰流柱后的湍流度比较低.数值计算结果对涡轮叶片内部冷却计算具有重要的参考价值.      

尾缘冷却跨声速涡轮气动特性的数值模拟

为了明确跨声速涡轮中尾缘冷却措施的引入对其气动性能的影响,针对跨声速涡轮叶栅中压力面半劈缝尾缘冷却方式和尾缘全劈缝冷却方式条件下,叶栅性能和尾缘激波系结构的变化进行了数值研究。结果显示,两种尾缘冷却措施都降低了叶栅能量损失系数,压力面半劈缝尾缘冷却方式效果更好,最佳情况损失系数下降达到48%。冷却流量对作用效果有一定影响,存在最佳值。冷却措施的引入显著改变了尾缘激波系结构,尤其对PSEJ冷却方式,将尾缘激波系压力面分支由原本一道强激波分成三道弱激波。     

不同喷口修形的二元收敛喷管RCS数值模拟

用等效棱边电磁流(EEC)方法和物理光学迭代(IPO)法自主开发了计算程序,用来分析不同尾缘修形喷管的雷达散射特性;在计算中考虑喷管出口边缘的绕射场对雷达散射截面(RCS)的影响;提出了一种高效的判断面元间遮挡关系的算法和喷管出口轮廓线的辨别方法;在此基础上,研究了四种不同喷口修形的二元收敛喷管的RCS特性.结果表明:该开发程序具有较好的计算精度;适当的喷口修形可降低重点姿态角下喷管的RCS.      

波瓣混合器流场试验

为了获得加力燃烧室波瓣混合器内的三维流场,采用带热电偶的五孔探针对波瓣混合器内的流场进行了测量.通过试验,得到了波瓣混合器出口不同截面的速度场和温度场,以及波瓣混合器沿流向的热混合效率、总压恢复系数变化规律.试验结果表明:在波瓣混合器尾缘处会产生一对相互逆转的流向涡,流向涡的产生加强了内、外涵气流的对流混合;越靠近波瓣混合器尾缘,内、外涵掺混能力越强,热混合效率增加速度越快,总压恢复系数下降也越快.     

有间隙涡轮叶栅壁面流谱实验与数值模拟研究

对具有 3 .6 %相对叶顶间隙的涡轮叶栅进行了壁面流动显示和三维流场数值模拟 ,分析了大叶顶间隙涡轮叶栅的壁面流动特点 ,从数值模拟上特别对尾缘附近的壁面流动结构进行了详细讨论。实验与计算结果表明叶栅尾缘附近的流动是非常复杂的 ,同时由于叶顶间隙的存在 ,上下尾缘附近的壁面流谱明显不同     

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

前缘喷淋射流对导叶压力面冷却特性的影响

为了评估涡轮导叶的前缘喷淋射流对压力面多排气膜孔冷却特性的影响,在高亚声速风洞中进行了实验,获得了有无前缘喷淋射流时叶片表面的气膜冷却效率和传热系数。叶栅进口雷诺数(基于叶片弦长)范围为2.0×105~4.0×105,出口等熵马赫数为0.95,叶片前缘和压力面分别都包含6排圆形孔,质量流量比的范围分别为2.46%~4.57%和2.00%~3.71%。实验结果表明:在没有前缘喷淋射流时,压力面前半段的气膜冷却效率受质量流量比的影响较小,而后半段的气膜冷却效率随着质量流量比升高而增大。前缘喷淋射流使压力面多排气膜孔的冷却效率提高了20%~70%,并且使气膜冷却效率沿流向分布更均匀。不论是否有前缘喷淋射流,压力面的传热系数比都随质量流量比升高而增大,沿流向看,前缘喷淋射流提高了压力面前缘和尾缘区域的传热系数比而对压力面中间区域影响较小。     

带气膜出流的旋转叶片冲击冷却的实验研究

用实验方法研究了放大 4倍的涡轮叶片前缘模型在旋转状态下,带气膜出流时叶片内冲击面的换热特性。实验结果表明:雷诺数的增加显著地增强了换热,但叶片前缘和尾缘的增加是不同的,而浮升力对换热的影响是复杂的;在低雷诺数的情况下,旋转对换热影响不明显,而在高雷诺数时,会使换热降低。      

微气泡发射沸腾形成机理

为了探究具有超高换热性能的微气泡发射沸腾现象的形成机理,采用FLUENT软件对加热面上单个气膜周围的速度场进行数值模拟,并与实验结果进行对比.实验结果表明,对于水,微气泡发射沸腾现象发生时,加热壁面上会出现气膜破裂的过程,并且过冷度和壁面过热度的升高会加剧这一过程.对于酒精,微气泡发射沸腾现象很难发生.计算结果表明,在过冷条件下气膜周围存在marangoni对流,对于水而言,过冷度和壁面过热度的升高会增强气膜周围的marangoni对流过程,而在酒精气膜周围 marangoni对流相对较弱.因此由气膜周围强烈的marangoni对流过程引起的气液界面上的扰动可能造成气膜破裂,这可能是微气泡发射沸腾现象形成的原因之一.      

Experiments on unsteady vortex flowfield of typical rotor airfoils under dynamic stall conditions

A new experiment for airfoil dynamic stall is conducted by employing the advanced particle image velocimetry(PIV) technology in an open-return wind tunnel. The aim of this experimental investigation is to demonstrate the influences of different motion parameters on the convection velocity, position and strength of leading edge vortex(LEV) of airfoil under different dynamic stall conditions. Two different typical rotor airfoils, OA209 and SC1095, are measured at different free stream velocities, oscillation frequencies, and angles of attack. It is demonstrated by the measured data that the airfoil with larger leading edge radius could notably decrease the strength of LEV. The angle of attack(Ao A) of airfoil can obviously influence the dynamic stall characteristics of airfoil,and the LEV would be effectively inhibited by decreasing the mean pitch angle. In addition, the convection velocity of LEV is estimated in this measurement, and the results demonstrate that the influence of airfoil shape on convection velocity of LEV is limited, but the convection velocity of LEV would be increased by enlarging the oscillation frequency. Meanwhile, the convection velocity of LEV is a time variant value, and this value would increase as the LEV convects to the trailing edge of airfoil.