非正交曲线坐标叶轮机械S2流面正问题计算机程序.pdfVIP

第 1卷第1期 工程热物理学报 Vol.1, No.1 1980年 2月 JOURNALOFENGINEERINGTHERMOPHYSICg Feb. 1980 非正交曲线坐标叶轮机械 S2流面正问题计算机程序术 丁云峰 邓子同 （中国科学院计算中心） 摘 要 对于非正交曲线坐标叶轮机械s：流面正问题 （亚声速情况）在计算机上实现了一种收敛 算法 ，对一 个高载 荷轴流 式两级风 扇 的验算说 明计算结果与实验数据相 比很一致 ， 主导方程采用新推得的非线性形式——方程(8),与传统沿用的线性形式——方程(13) 相 比 ， 当流 场 中高 理 声 速 区较 大 时 ，新 解 法 改 善 了数 值 求 解过 程 的 收 敛性 ． 一、 丹iJ 百.CN T 叶轮机械内的流动十分复杂．五十年代初，吴仲华教授n1提出的三元流动理论至今 仍是叶轮机械内部流动计算的基础．.IE P 文献[1】中推导了正交曲线坐标下沿S．流面和S2流面运动的主导方程，有非线性形 ET 式和线性形式．六十年代初，在[2-4]中又推导了非正交曲线坐标下沿Sl流面和＆流 J 面运动主导方程的线性形式．[1]认为，对于亚声速情况用线性形式的主导方程作数值 求解要方便些． 随着电子计算机的发展，到六十年代中期，三元流动理论开始了计算机程序实现阶 段．到现在十多年来，主导方程都传统沿用线性形式． 本文推导了非正交曲线坐标S2流面正问题主导方程的非线性形式，以此作为数值求 解时的出发方程．流函数掣的差分方程用超松弛迭代法求解，也用线性形式的主导方程 做了计算，两者相比说明新解法在高亚声速区域较大时改善了收敛性．同时，对于高参 数情况，采取了逐步提高参数和修正进气角的措施． 程序用TQ-16机-ALGOL-60语言编制，可用于分析叶轮机械气动性能 （亚声速情 况）．一次可计算两级叶片的全流场．计算结果与实验结局相比很一致． 二 、基本方程 在子午面上采用非正交曲线坐标(Xl，X2)．在这种坐标下，轮壳线，轮毂线及各列叶 片的前后缘在子午面上的投影均可取为z1线或X2线．因此不仅通用性强，而且便于处理 +本文曾于1978年11月在杭州举行的全国第二届工程热物理学会上宣读． 1期 丁云峰等：非正交曲线坐标叶轮机械s：流面正问题计算机程序 37 边界条件． 对于动叶，取随动叶以角速度扬旋转的动坐标(XI，x：，妒)．其余为静止坐标(xl，z：， p)．在 绝 热 、稳 定 的 假 设 下 ，S。流 面 上 流 动 所 遵 循 的方 程 [2-钉 为 ： 连续方程： { 言了 （加∥1~／瓦）+芸≯ （如∥2y／磊）一o (1) 运动方程： Xl向： W2 ~／丢嘉[（形lsin（pt一岛）+∥2)N／瓦]一岳[~／五 （矽1 a（Ver）亘￡ +∥2siri（p，-p。））]）+鼍予亘(3 axl 8xi +丁暑≥一一，．一暑}兰号手立 (2) X2向： 鲁 [（∥1+W2SiD (Ol-p：)）~／瓦]一吾 ≥[（∥1sin（口-一 日：） i a， 石r T.CN +驴）~／瓦]一每[一孚嘉 （矿。，）+丽一7砑一f,] (3) 能量方程： P.IE dI=0 (4)JET dt 即，沿流线 （在动叶区为相对流线）为常数．须注意，在动、静两区交界