叶栅内定常及非定常粘性流动数值模拟.pdfVIP

申田工程蕾■曩拳台 犏号l20幛2D16 燕枫气动燕力擎学术畚说 吐接内定常及非定常糙丝麴邀焦垄苣塾 刘建军 中国科学院工程热物理研究所 北京市2706信箱，100080 Email：jj-@eq蝽er嘲je单．∞．血 ● 摘要 奉文蛤出了一种模拟叶栅内准三堆定常和非定常粘性流动的数值方法．对于定常 ● 潼动．采用1、rDk-W哪dro廿格式和代数措湔薛型求解霄诺平均N州盯．stokcs方程，使 ● 用当地时伺步长和多同格技术使计算加速收敛到定常状志；对于非定常流动．使用双时间 步长和垒童式离散．采用与球解定常流动相似的多罔格方法求解隐式离教方程．文中描述 ● 了具体方法井给出了VⅪ透平叶栅内的定常流结果和1．5擐透平叶橱内的非定常数值结 集． 1、弓I言 计算流体力学在叶轮机械气动设计及其内部流场分析中发挥着越来越重要的作 用．基于定常流动假设的二维、准三维和全三维数谴方法，在叶型／叶片设计、叶栅变 工况损失特性分析以及多级叶片摔之间的气动匹配等方面得以广泛应用．由于转动和 静止部件的变错捧列．叶轮机械的内部流动实际上是非定常的。得益于计算机技术和 数僮计算方法的持续发展，现在人{|’1可以甩非定常方法，研究上游叶片捧尾迹对下游 叶片捧的毒响l}I及耜邻叶片撵之间的势流相互作用：研究非定常流动条件下才有的静 叶片捧之间的相对罔向位置和叶片捧之间的轴向位置变化对叶轮机效率的影响；也可 以研究时胸乎蚜漉蓊与定常流场韵差异，以及在定常瀛场数值模拟中考虑非定常因素 从而藏少这种差异的数学模型和方法。 除了稿流模型带来的不确定性，采用求解雷诺平均№vi稍toI。∞方程数值模拟叶 樯内定常及非定常粘性流动的方法是现实可行的。作者在近十年前研究开发了用于模 拟叶搪内准三维定常无粘流场的数值方法和软件刑，本文是把该方法向模拟叶橱内准三 ● 维定常和非定常粘性流场的推广．下面先给出使用多重弼格T、，D格式数值模拟定常粘 性流场的方法．然后介绍使用酶式离散和双时间步长求解菲定常流动的方法。最后给 ● 出v】Ⅱ蟠59遗平_叶栅的定常流场敷值结果及与实验结果的对比．和1．5级透平叶栅内 的非定常流场数值结果。 2、定常流动教幢解法 本文使用的方法求解任意回转面上叶栅内的流动，可以考虑流管厚度和回转半径 变化的影响．这里以二维控制方程为例，简要说明求解定常粘性流动的方法。二维无 量纲形式的雷诺平均N¨kr．S幻kes方程可以写为如下形式： (1) 塑+塑+盟：上f堡+堑1 … af 融 砂 ＆。I盘 砂J 一…一——————F_—可瓤虿 式(1)中u为守恒变量，F和G为无粘通量，R和q为粘性通量，黜．为基于无量纲特 征变量的雷诺数。使用ddwin．Lom“代数湍流模型估计湍流粘性系数，并假定流动 为全湍流，工质为完全气体。控制方程中各项的具体表达式不赘述。 采用显式有限体积T、巾hx．wendm圩格式和多重网格加速收敛技术数值求解方程 (1)，使用时间推进方法得到定常解的过程可表示为： ，，“+I