动静叶相互干扰非定常流动特性的数值研究.pdfVIP

动静叶相互干扰非定常流动特性 的数值研究 王正明，王嘉炜 (中国科学院工程热物理研究所北京100080) Tel：(010 Email：wzm(苞,etpServer．etP．ac,cn 摘要：本文对动静叶栅相互干扰的非定常流动特性进行数值研究。研究内容主要涉及两个特定的方面，即 一方面对叶型优化前后在动静叶干扰的非定常流动条件下的性能进行比较，另一方面对非设计工况有大分 离流动的动静叶干扰的特殊性进行研究。 关键词：动静叶干扰：非定常流动；叶型优化；分离流动 1 引言 叶轮机械中的流动实质上是十分典型的非定常流动，灏年来，髓着计算机与计算技术的 不断进步，为了更精确地预测叶轮机械的气动性能，非定常流动的数值模拟已越来越被重视 (如文献[1卜[5])。叶轮机械中的非定常流动现象是多方面的，造成非定常的因素也是多方 面的。对单排叶栅而言，在非设计工况下有大分离的流动，有尾缘涡脱落现象的流动，都具 有十分明显的非定常流动特征”’”，对于多排叶栅，由动静叶交替排列而产生的非定常特性 则是叶轮机械中最根本的非定常流动特性，对于任何一排动叶来说，任何造成其进出口流场 不均匀的因素都是产生非定常流动的根源。 容易理解，对于一名设计者来说，如何设计得更好和更准是其考虑最多的问题。多年来， 反问题方法作为叶片优化设计的重要手段之一已经被广泛用于新型叶轮机械的气动设计过 程中，然而人们不难发现，迄今为止所有反问题优化设计方法均是在定常流动的假设下构思 形成的，而且可以预见，人们还不可能在短时间内找到在非定常流动祭件下求解反问题的思 路。于是退而求其次，研究在定常条件下优化得出的叶型在非定常流动条件下特别是在动静 叶干扰条件下的性能，应该是一个有现实意义的课题。在非设计工况下，如果存在大攻角下 的分离流动，则动静叶之间的干扰将具有双重非定常流动特性，即-方面是由大分离流动本 身引起的，另一方面是由动静叶相互作用引起的。本文基于任意曲线坐标下的N—S方程，采 用分区计算的数值方法，将对关于动静叶干扰的以上两个特定方面的问题进行研究。 2基本方程与数值方法 本文求解动静叶干扰非定常粘性流动的计算采用非正交曲线坐标与非正交速度分量下 张量形式的守恒型N—S方程“1： 日一 日．‘ 连续方程：鼍(√g力+-暑●·；呈p矿。)=0 (1) 8 哆qgu 动量方程的分方程： 鲁晤一∽吲唔彘∥∥删狒，， +砉(而巾+专(一瓜掣v一=。 本工作得到国家重点基础研究项目(G1999022309)的资助 争参肿∽％，+≥彘加Ⅵ砩川 +善(函内+砉(一瓜。^矾zⅢ 能量方程： 刍√虱∥一p)】+妥(粤P矿‘』) “ 蕊恼 (4) +嘉c一瓜峻砉)+考(一庙O’zjk睾=。 两步显式时间推进的数值离散格式与B叱湍流模型被用于数值求解的过程中。 本文采用分区计算的方法计算多通道动静叶相互作用的非定常流场，分区后每个区的计 算十分类似于单排叶栅单通道的计算。按本文的分区方法，所有边界条件将如下处理：在前 排叶片的进口边界给定总压，总温和气流角：在后排叶片的出口边界给定静压；在叶片表面 均同时采用无滑移边界条件和绝热壁面条件；在叶片前后延伸段的周向边界处，相邻两区的 共同边界上采用流场连续性条件，而多通道的外边界按动静叶干扰的物理特性应满足空间周 期性条件；前排叶片出口与后排叶片进口的交界面上应满足流场连续性条件，具体处理方法 是将每个时间步后排叶片计算区进口的静压作为前排叶片计算区出口的背压，将前排叶片区 出口的温度和两个速度分量作为后排叶片计算区的进口值，计算中由于动叶与静叶在不同瞬 时的相对位置不同，故交界面上的数据传递需通过插值计算来完成。 3．舁1刘及分机 本文采用两个不同算例来分别探讨前面提到的关于动静叶干扰的两个特定方面的问题。 算例一是涡轮叶栅算例，用一排静叶栅和一排动叶栅来进行计算。由于静叶与动叶叶片