涡轮叶片前缘气膜冷却数值模拟.docxVIP

涡轮叶片前缘气膜冷却数值模拟 1 数值模拟方法 为了提高发动机的可靠性和输出，现代高性能航空发动机的电机输入温度增加。现代先进航空发动机涡轮进口温度普遍在2000K之上,这样的温度远远超过金属叶片的承受温度,因此必须采用高效的冷却方法对涡轮叶片进行冷却。通常燃气涡轮叶片前缘是热负荷最大的区域,因此这一区域的气膜冷却设计也是最困难的。 数值模拟是研究气膜冷却机理的重要手段,很多国内外研究者采用数值模拟的方法对涡轮叶片前缘气膜冷却问题进行了大量地研究。Garg和Rigby对前缘3排冷却孔的VKI叶片进行了数值模拟,计算域包含冷气腔,研究不同冷气流量下叶片表面绝热效率、冷却孔附近区域速度及温度的分布。Thakur等采用带有壁面函数的低Re数k-ε模型对带有气膜冷却的半圆形叶片前缘模型进行了数值模拟,分析了由主流和射流相互作用而引起的漩涡结构。Lin和Shih用一个用带有3排气膜冷却孔的半圆柱形前缘模型来研究气膜冷却,数值分析了前缘表面的气膜冷却效率分布、温度分布及速度矢量分布。颜培刚等采用具有TVD性质的三阶精度有限差分格式及多区网格技术数值模拟了在前缘滞止区域附近冷气射流的运动规律以及气膜冷却对涡轮气动参数的影响。安柏涛等采用标准k-ε模型价改进的壁面函数方法对某型燃气轮机透平静叶进行了详细的数值模拟,对前缘气膜冷却复杂三维流动以及整个叶片的气膜冷却特性进行分析。 本文基于全三维N-S方程隐式求解,以圆柱形前缘模型为研究对象,采用FNM(Full Non-Match,完全非匹配)形式分区重叠结构化网格来模拟冷却孔射流,研究前缘不同吹风比下交错孔排冷气射流对流场及气膜冷却效果的影响,并针对孔排进行优化设计。本文暂不考虑冷气腔及壁面热传导的计算。 2 物理模型和计算方法 2.1 数值计算模型 对于任意曲线坐标系ξi,可压缩Reynolds平均的三维Navier-Stokes方程和双方程湍流模型建立在旋转角速度为Ω的相对转动坐标系上,可写成 ??Q?t+??Fi?ξi+1Re?D+?S=0?Q??t+?F?i?ξi+1ReD?+S?=0(1) ?QQ?为未知变量的一般形式,?FiF?i为对流通量矢量,?D为粘性项,?S为源项,包含离心力和哥氏力。 ?Q=J[ρρu1ρu2ρu3ρErρ?1ρ?2]?Fi=J[ρUiρu1Ui+ξi,1pρu2Ui+ξi,2pρu3Ui+ξi,3pρΗrUiρ?1Uiρ?2Ui](2)?D=-??ξiJξi,j[0τj1τj2τj3τjkuk-qj?ωj(?1)?ωj(?2)]?S=-J[00ρΩ(Ωx2+2u3)ρΩ(Ωx3-2u2)0S?1S?2](3) 式中ρ为密度;u1,u2,u3为相对速度矢量的3个分量;Er,Hr分别为相对滞止内能和滞止转焓;J=?(x1,x2,x3)/?(ξ1,ξ2,ξ3)为坐标变换Jacobian矩阵;ξi,j=?ξi/?xj;Ui=ξi,juj为逆变速度。 由于q-ω模型的计算量相对较少,边界条件处理简单,又能适应粗糙的初始湍流流场,所以本文采用Coakley的q-ω双方程湍流模型。数值计算格式采用了LU-SGS-GE (Lower Upper Symmetric Gauss Seidel Gaussian Elimination method)隐式格式和改良型高分辨率MUSCL TVD(Monotone Upstream-centred Scheme for Conservation Laws, Total-Variation-Diminishing)格式以加快收敛速度以及提高解的精度和分辨率。 射流孔与主流区及冷气腔之间的数据交换采用FNM方法,文献和分别针对平面叶栅和平板气膜冷却进行了数值模拟,验证了该方法的可行性。采用FNM方法可以把主流场网格和射流孔网格分别独立生成,在射流孔网格改变的情况下不必重新生成主流场网格,为数值计算提供方便。 2.2 实验模型的建立 文本研究的物理模型是参照Cruse的实验模型,即一圆柱形前缘,在静止前缘设置了两排复合角度的圆柱形气膜冷却孔。前缘滞止线位置为x/d=0,射流孔每排9个,直径均为6.3mm,孔与叶高方向呈20°,第二排孔与第一排孔夹角为25°。孔与孔之间的距离为S/d=7.64。圆柱形前缘半径50.8mm,前缘壁厚25mm,孔长度L/d=12。前缘不考虑上下端壁的影响。前缘材料为聚苯乙烯,由于其导热系数相对较小,所以实验中假设为绝热壁面,具体实验装置参见文献。为简化计算模型,计算域如图1所示,滞止线附近的孔简化为4个半孔,第二排孔简化为3个孔,滞止线为对称边界。坐标轴x为原点位于滞止线上并沿曲面向下游方向,y为垂直壁面方向,z为叶高展向方向,z=0位于第二排第一个孔左端。 主流和射流进口均给定总温总压,出口给定背压,主流