高负荷轴流压气机叶栅二次流诱发损失的数值模拟.docxVIP

高负荷轴流压气机叶栅二次流诱发损失的数值模拟 1 端壁流动诱发损失的数学模型 随着高压气机的开发，叶片流量损失控制技术已成为国内外许多科学家的研究重点，尤其是在流量损失机理的基础研究方面。 在压气机通道中, 除壁面摩擦作用外, 损失的产生往往与二次流现象相关。近年来的不少研究[1~4]均显示, 受到横向压差作用而生成的通道横向二次流能使低速气体堆积于吸力面角区, 诱发流动分离并带来损失;此外, 大量有关端壁流动控制技术的研究, 如非轴对称端壁造型、涡发生器、端壁导流片等, 均通过消弱横向流动或通道涡流以便对损失进行有效控制。然而, 详细对比上述文献就会发现:虽然各人结论中, 关于端壁横向流动诱发损失的定性认识是明确的, 但大多数文献止步于对通道总损失参数的探讨;在分析损失产生的机理时, 也仅仅是定性地将其与观测到的流动现象进行联系, 缺乏对各类机理主导性的判断与更加精确的定量关联。 事实上, 压气机近端壁的流动现象是极为复杂的, 即使是同一种二次流现象, 也会由于不同的机理诱发多种损失。当外部环境发生变化时, 各类二次流出于不同机理而诱发的损失的能力也会不同程度的改变。通道损失的总体变化, 正是多种损失源此消彼长, 综合作用的结果。因此, 弄清各类二次流诱发损失机理的主导因素, 建立二次流与损失之间的定量关联, 对于理解压气机流场损失的产生机理及对端壁流动控制方案的合理利用具有重要的意义。 本文通过对一高负荷轴流压气机平面叶栅开展数值研究, 试图达到以下两点研究目标:其一是要对通道内二次流诱发损失的机理进行定性判别, 并进一步推导各种机理与损失的定量关系;其二是分析在改变攻角情况下, 各类损失的变化规律。 2 对象与数值法的研究 2.1 高压压气机出口静子级去预旋静子叶片的设计 本文研究对象为一高负荷轴流压气机叶栅。叶栅为双圆弧叶型, 按照高压压气机出口静子级的流动特征设计为去预旋静子叶片。表1给出了该叶栅的部分几何及气动设计参数, 有关该叶栅的设计方法与对应的实验风洞的详细介绍详见文献。 2.2 计算值的方法 2.2.1 叶栅定常流场计算 数值模拟采用的网格由Auto Grid生成, S1流面贴体网格拓扑为H-O4H-H, 结构如图1所示。进口边界位置同实验中叶栅上游附面层测点, 即叶栅上游1.2倍轴向弦长处 (为与轴流压气机中表述一致, 本文定义S1流面内垂直叶栅额线方向为轴向) , 出口边界选取尾缘下游2倍弦长位置。另外, 为节约计算空间, 引入周期性与轴对称假设, 因此实际求解的计算域仅为单个叶栅通道半叶高流动。近壁面第一层网格无量纲尺度y+1, 计算域网格数合计约76万。 计算采用商用软件NUMECA的Euranus求解器, 通过时间推进RANS方法获得叶栅定常流场。湍流模型选用Wilcox k-ω模型。时间项推进采用四级四阶Runge-Kutta格式, 空间项则采用Jameson的有限体积中心格式进行离散。计算中另使用多重网格技术与局部时间步方法加速收敛。 边界条件按照实验给定, 进口边界给定均匀静温、及速度展向分布 (如图2所示) 。出口给定均匀背压边界。叶表与端壁均给定绝热无滑移边界条件, 在50%Span对称面上设置镜面边界条件。 2.2.2 数值结果对比 图3, 4给出了在-1°攻角下, 实验和数值计算结果的对比。各图中的无量纲损失Loss, 依照实验定义为 式中pt, p分别表示总压与静压, 显然Loss是以叶展中部的总压与动压为参照进行构造的。实验采用7孔总压探针及不对向测量法测量流场中各点总压与气流方向。经实测及换算, 七孔探针所测角度误差约±0.5°, 总压损失系数误差约±0.3%。 由图3 (a) 可知, 数值计算结果对高叶展区域的损失略有过估计, 损失核心区的最大值比实验偏高, 但损失区域的边沿区及总体分布与实验呈现较好的一致性;图3 (b) 、 (c) 关于平均参数的计算结果说明:数值模拟所得损失系数在展向范围内分布趋势基本准确, 但定量来看, 在实验损失峰值处过估计约0.034;至于气流角, 数值模拟结果在全叶高范围内高于实验1°~2°, 于10%Span左右的欠偏转现象定位准确, 除在定量方面略显不足, 其展向分布的总体趋势与实验吻合。 图4中给出了叶表流动结构的实验与计算对比图。由图4 (a) 标注, 端壁面内的主要临界点预测准确, 吸力面分离区比试验略有偏小, 流动拓扑与实验基本相符。图4 (b) 中吸力面叶根前缘区鞍点 (S2) 的轴向位置预测准确, 角区分离除径向扩展范围略大于实验外, 其拓扑与流动趋势均与实验符合较好。 综合以上关于出口参数与叶表流动结构的对比, 认为数值模拟结果基本反映了实验中的真实流动现象。在以二次流与损失关联性探讨为目标的研究中, 本文所用计算方法的精度是