叶轮机通流计算的时间推进方法.pdfVIP

航空动力学报JOURNAL OF AEROSPACE POWER1999年　第14卷　第1期　VOL 航空动力学报 JOURNAL OF AEROSPACE POWER 1999年　第14卷　第1期　VOL.14　No.1　1999 叶轮机通流计算的时间推进方法 　　季路成　孟庆国　周　盛 摘要：本文首次将时间推进思想应用于通流计算，采用四级龙格—库塔方法求解，根 本上解决了目前出现的轴向流动跨音高负荷叶轮机设计困难，初步工作取得了较好结 果。作者们预期，时间推进通流计算方法还将弥补现有许多通流计算方法的不足，将 在工程设计中发挥更大作用。 主题词：跨音速　高载荷　叶轮机械　气动力计算 分类号：V231.3 TIME-MARCHING METHOD FOR THROUGH-FLOW COMPUTATION OF TURBOMACHINERY Ji Lucheng，Meng Qingguo，Zhou Sheng 1　引　言 ［1］ 　　五十年代，吴仲华教授的开创性工作 使叶轮机设计方法得以实现本质提高， 设计者开始可以将叶轮机内全三维非定常粘性流动近似为两种流面内准三维定常流动 （即叶片截面(S1)和子午面(S2)），并在两种流面解间迭代，从而初步解决了叶轮机设 计和分析问题，完成了叶轮机设计体系由一维方法向准三维方法的转变。随着计算技 术和计算机技术的飞速发展，近些年来，全三维N-S解已发展成熟，叶轮机设计体系正 在实现由准三维设计体系向基于分析三维设计体系的过度。虽然如此，准三维设计体 系建立起来的通流计算仍扮演着不可忽视的角色，通流设计仍是确定设计目标范围和 设计基本参数的有力工具。因而，即便在三维设计体系中，通流计算仍是不可缺少的 部分。这样，完善和发展通流设计方法仍有一定的现实意义。 　　从六十年代中后期至今，已有流线曲率、流函数法、有限元方法成功用于通流设 计，它们在近年来叶轮机设计中起到举足轻重的作用。然而，随着人们对叶轮机关键 file:///E|/qk/hkdlxb/hkdl99/hkdl9901/990106.htm（第 1／5 页）2010-3-22 22:44:04 航空动力学报JOURNAL OF AEROSPACE POWER1999年　第14卷　第1期　VOL 性能参数要求的日益提高，某些流动的轴向马赫数局部已突破音速，流动轴向分量场 同时有超音和亚音区存在，通流控制方程不再属椭圆型，而是混合型的，流线曲率、 流函数法、有限元方法已不再有效。轴向速度跨音问题最早出现在高负荷涡轮中，为 ［2］ ，这种求解方法人为认知因 解决问题，研制者常采用人为限制等方法，如Denton 素很大，得到的叶排内超音膨胀不正确，因而错误给出给定压 比下出口气流参数。另外，近年来在国内超音通流风扇设计提到日程，在设计工况因 为流动在整个域内是轴向超音的，只要对原有流线曲率方法通流程序加以考虑流动特 ［3］ 征的修改即可，但计算速度却由于极低的松弛而大大降低 。对于超音通流风扇某 些非设计工况，由于高背压，激波可能出现于风扇转子或静子，原有通流方法用于分 析非设计工况时便遇到了根本困难。另外，在轴向超音来流激波转子或激波静子的组 合风扇研究探索中，由于轴向流动跨音，原有通流计算方法即便在设计工况也是不可 能的。 　　基于时间推进思想解决跨音流动问题的有效性，本文保留了非定常通流控制方程 中时间导数项，使轴向流动跨音通流控制方程统一为双曲型，从而用时间推进法完成 通流计算。 2　通流计算气动方程组 　　描述叶轮机内流动的模型方程是N-S方程不同层次上简化的结果，本文采用如下 假设：(1)周向平均假设；(2)流动无粘；(3)绝热流动；(4)理想