在转静干涉流动中的应用.pdfVIP

第 13 卷　第 2 期 航空动力学报 V o l13 N o 2 1998 年 4 月 Journa l of A ero space Power A p r. 　1998 B a ldw in - L om ax 模型 在转静干涉流动中的应用 北京航空航天大学　　吴先鸿 　陈矛章 【摘要】　本文以拉格朗日方式研究叶轮机内动、静叶间的相互作用, 通过直接跟踪流体涡团的时 ( ) 间历程, 来描述流体的非定常流动过程。发展了代数湍流模型 B a ldw in - L om ax 湍流模型 在尾迹 中的应用方法, 克服了其它数值方法中无法准确捕捉尾迹中心线的运动轨迹, 以及计算出的边界 层外的湍流涡粘性系数偏大的缺陷。用该方法对叶轮机内转子静子相互作用过程进行了数值模 拟, 与实验的对照表明, 该方法是可信的。数值结果还表明, 采用该模型能正确地模拟出叶轮机内 尾迹和各种旋涡等非定常因素的流动过程, 从而为研究叶轮机内非定常流动及掺混过程提供了一 个新途径。 　主题词: 　叶轮机械　湍流　数值模拟 　分类号: 　 2 11 V 1　引　言 　　已有的研究表明: 叶轮机内的非定常作用影响叶轮机各方面的性能, 如叶片载荷、级效率、 噪声特性和传热特性等等[ 1, 2 ] 。事实上, 叶轮机内部的非定常流动包含各种影响因素, 除湍流运 动固有的非定常性外, 还存在两种既有一定联系又有本质区别的非定常流动, 一种是由动、静 叶相互作用引起的; 另一种则是与流动失稳相联系的, 如旋转失速、喘振等。本文的内容主要涉 及前一种非定常运动。 本文中给出涡方法求解二维粘性问题的基本方程和所应用的湍流模型, 由于初始条件、边 界条件以及用涡方法求解动、静叶相互作用问题的计算步骤在文献[ 6 ] 中已有详细说明, 在此 不再赘述,B a ldw in - L om ax [3 ] 湍流模型在尾迹中的应用是本文的重点。 2　基本方程和湍流模型 2. 1　基本方程 　　由文献[ 6 ]二维扰动涡量的输运方程为: ′ d ′ ′ ′ ′ 2 ′ = - (u ) - ( U ) + (u ) + v ( 1) d t ( ) [3 ] 式 1 符号说明详见文献[ 6 ]。本文按B a ldw in - L om ax 模型计算湍流粘性系数。 1997 年 7 月收稿; 1997 年 10 月收到修改稿。 男　28 岁　博士　北京航空航天大学 404 教研室　100083 124 航 空 动 力 学 报 第　13　卷 2. 2　湍流模型 　　对湍流而言, 由于代数涡粘性模型非常简单, 计算量小, 而且无需判别边界层的边界和进 行积分运算, 因此在工程上应用很多。在叶轮机非定常的计算中, 也有很多工作采用这一模型, [4 ] [3 ] 如A rnon e an d P accian i 等的工作。但是B a ldw in - L om ax 湍流模型 在应用于动静叶相互 作用的非定常流动的数值计算中, 面临两个明显的困难: 首先由于尾迹的位置一般远离边界 层, 因而按照B a ldw in - L om ax 湍流模型计算出的y m ax 和 F m ax