多排叶轮机三维非定常流场的数值模拟方法研讨.pdfVIP

中国工程热物理学会 热机气动热力学 2006年学术会议 编号：062086 多排叶轮机三维非定常流场的 数值模拟方法研究 刘锡阳1，2赵晓路1 (1中国科学院工程热物理研究所，北京100080t 2中国科学院研究生院，北京100080) Tel：010E．mail：lxy922000@163．corn， 摘要：非定常性是叶轮机械内流动的本质属性．本文采用Jameson提出的双重时间步法．使用有 方程进行封闭，实现了多叶片排非定常流动的三维数值模拟．并使用该程序对某无导叶对转涡轮的流 场进行了计算，通过与商业软件Numeea计算结果的比较分析，初步校核验证了程序的准确性和可靠性， 为下～步工作提供了一个较好的平台． 关键词；双重时间步法；非定常流动；数值模拟 1 引言 由于叶片排相对转动带来的流场周期性变化是叶轮机非定常流动的主要特征之一。 在工程设计中的有着重要的意义．由于实验需要投入大量的财力，且准备的周期比较长， 人们一直希望通过数值计算的方法能够准确的预估它对叶轮机性能的影响。 近几十年来，众多的研究人员投入到叶轮机非定常数值模拟中，发展了很多描述叶 种方法，目前应用得最为广泛。 、 本文在原有定常计算程序的基础上13I，采用双重时间步法，实现了非定常流动时间 精确模拟。控制方程为Navier-Stokes方程，采用有限体积法进行空间离散，湍流模型采 用Baldwin．Lomax模型。本文还使用该程序模拟了某无导叶对转涡轮，并与商业软件 Numeca进行了比较，初步校核验证了程序的准确性和可靠性． 。 2数值模拟方法 2．1双重时间步法 。 ， 本文采用Jameson[41提出的双重时间步法，在求解过程中，将一个流动周期划分为几 十个物理时间步长。在每个物理时间步长内通过时间推进求解Navier-Stokes方程，时间 推进的步长称为虚拟时间步长l卵。 Navier-Stokes方程可以写成如下紧凑的形式： ‘． 豢=Jjc妙) (1) 研 其中u表示通量矩阵，等表示物理时间导数项，R(u)项包含了定常流动控制方 C玎 361 程中的非时间相关项，将物理时间导数项移到方程左边，并添加虚拟时间导数项，(1)式 可变为： ’(訇，=删一詈。 ．∞ 在某一物理时间层内，随着虚拟时间警进，当达到收敛标准时，有(期，哼。， 所以方程组(2)中添加虚拟时间导数项的做法并不违反物理本质。 2．2边界条件及计算网格 对转涡轮的叶型如图1所示，第一个叶片为静子，第二个叶片为第一级转子，第三 个叶片为第二级转子。涡轮叶片数均为36片。 52887。Sl流面上的网格分布如图2所示。 对转涡轮叶片示意图’ 、图2Sl流面上的网格分布 涡轮进口给定总温，总压，气流角；涡轮出口给定静压；固壁处采用无滑移壁面边 界条件，绝热；叶片前尾缘延长部分为周期性边界条件；叶排间网格交界处采用滑移界 面，上下游叶排间参数的传递采用插值得到。边界上的具体参数如表l所示。 表l边界条件 进口总压P·(Pa) 300000 进口总温p(K)