高温涡轮叶片表面流动与换热特性研究.docxVIP

高温涡轮叶片表面流动与换热特性研究 0 叶片表面边界层转合作 高温气和涡轮叶表面的换热系数是叶片冷却系统设计、温度场分析和热应力评估的基础。在满足叶栅气动设计的前提下,降低叶片外表面的换热系数,减少燃气向叶片的传热量,可以减少叶片内部冷却所需要的空气流量,同时减少燃气流动中的热损失,提高燃气轮机热效率。 叶片外表面换热系数的大小与表面流动的边界层发展过程有关。Dunn等在透平级叶片表面压力和表面传热的测量结果,证明透平叶片表面存在边界层转捩流动现象,包括三个阶段:层流、转捩和湍流。层流的对流换热系数低,湍流换热系数高,因此流动中的转捩位置决定了沿叶片表面换热系数的变化。假设叶片表面能保持层流,那么外表面换热系数应该是最低的。 Nealy等研究了C3X与Mark-II两种高负荷无气膜涡轮静叶表面的换热,分析了叶栅出口马赫数、出口雷诺数以及型线对叶片边界层发展及表面换热的影响。这两种叶栅具有几乎相同的进、出口流动条件,但是叶型完全不同。由于吸力面转捩特点的不同,结果两者表面换热系数分布存在很大的差异。Nicholon等在相同流动条件下,改变叶片的安装角,通过优化动叶片压力面型线,控制表面边界层发展,降低压力面热负荷。Ames等研究了来流湍流强度和湍流尺度对涡轮叶片边界层转捩和传热的影响,证明湍流度提高了传热强度,并使吸力面边界层转捩提前。Abuaf等、Hoffs和Bunker研究证明在粗糙表面上,流动转捩提前,叶片吸力面湍流区的传热增加。 本文工作的目的是对比两组几何设计的涡轮叶型,在出口马赫数和出气角相同的流动条件下,分析动静叶片不同叶型特点对边界层发展和表面换热的影响,评估表面换热与来流湍流度等相关性,为今后主动设计叶型降低热负荷提供理论基础。 1 数值方法 1.1 边界层流场模拟 湍流转捩模型是准确计算表面传热的关键。Medic认为k-ε湍流模型对涡轮叶片表面低雷诺数区域流场传热预测模拟能力不足,预测误差较大。全湍流假设的k-ε和k-ω模型只能求解湍流边界层流动,不能准确分辨层流及转捩状态,求解的壁面温度和换热系数误差都较大。本文选用SST湍流模型,并配以γ-Reθ转捩模型,能够很好地预测壁面边界层转捩位置,相对其它模型能较好地模拟附面层内的复杂流动过程和传热过程,对叶片表面热负荷分布预测更加准确。本文使用CFD软件Fluent计算流动传热,壁面条件为无滑移等温壁面。进口给定总温总压和湍流强度,出口给定静压。工质选用理想气体模型,黏性系数采用Sutherland公式表示。 表面换热系数定义为: 其中,Q为叶片表面热流密度;Tg为进口总温;Tw为叶片表面温度。 1.2 网格无关性验证 带有转捩模型的流动计算,对边界网格质量的要求远高于一般湍流计算。本文使用ICEM软件对计算域分区,划分四边形结构化网格。在叶片周围使用加密O型网格,划分边界层网格以提高计算精度。叶片型线外方向第一层网格高度满足y+1,同时以增长率1.005布置节点,保证边界层内布置足够网格层数。保持壁面处y+1,改变网格节点数,得到三种不同数量的网格,做网格无关性验证,网格数据列于表1。三种数量网格对压力预测没有影响,但是吸力面换热系数分布预测有所不同,中等网格与细网格预测换热系数分布结果差别很小,可以说明中等数量网格满足计算需要。本文后面计算全部选用中等数量的网格。 2 计算 2.1 表面等熵马赫数分布,3. MarkⅡ与C3X是高负荷涡轮导向叶栅常用的两个案例,具有丰富的实验测量结果。两个叶栅的主要几何参数如表2。如图2所示,两个叶型的主要差别在于叶栅前缘形状、相对节距和安装角,但是几何出气角相同。表3为叶栅流动的计算工况条件,其中出口雷诺数特征长度区叶片弦长。 图3和图4为两种叶栅的压力分布图。叶片表面压力分布预测值与实验结果一致,但是在吸力面喉口位置以后,预测压力低于实验值。两个叶型吸力面压力分布有明显差别,MarkⅡ叶型吸力面0.25位置后,有很大的逆压梯度,呈明显的V字形分布。 图5为两个叶型预测的表面等熵马赫数分布,在吸力面靠近前缘部分,两个叶栅的表面马赫数都快速增长,但MarkⅡ叶型叶栅表面马赫增长的更加剧烈,两个叶栅不同的表面马赫数分布特点由两种叶栅的吸力面型线特点所决定。MarkⅡ叶型叶栅喉口位置,在吸力面25%弧长处,表面马赫数都明显大于1,说明该处肯定有较强烈的激波产生。而C3X叶型叶栅吸力面表面马赫数都基本小于1。不同的流动特点导致表面换热分布也有明显的不同。 图6和图7为两种叶栅表面换热系数分布图,在该工况条件下,两个叶栅吸力面,在大约25%弧长的位置处,表面换热系数都有明显增高的现象,该处即为边界层转捩点。可以看出采用该湍流模型预测的换热系数,压力面和吸力面上的整体趋势都是正确的。 通过改变进口湍流条件来研究两个叶型表面换热的不同