跨音速轴流压气机级三维粘性流场全工况数值模拟.pdfVIP

中国工程热物理学会 热机气动热力学 第十届年会 编号：012028 跨音速轴流压气机级三维粘性 流场全工况数值模拟 张士杰z袁 新。叶大均： (1．清华大学鸯筑学院蹇筑技术科学系，北京100084：2．清华大学热能工程薯．北京ID0084) 文 摘：采用一种快速隶解三维粘性流场的计算方法求解跨音速轴流压气机缎 TVD迎风格 内部流场及全工况特性．该方法以LU-SGS-GE隐式格式和MUSCL 式为基础，结舍壁面函数方法和简单的混台长度湍流模型，对三维可压绾置诺平 均Navie-Stokes方程进行求解．叶列间参数的传递采用混台平面方法并应用了微 机网络并行计算技术．计算得到了NASA37号低展弦比、跨音速轴流压气机级 70％设计转速下的全工况性能曲线，并重点分析了其中一些典型工况下的内部 流场．计算与实验结果的对比表明此方法能快速得到三维粘性流场的流动特性 且计算精度较高，可用来模拟跨音速轴流压气机级内的全工况三维粘性流动． 关键词：全工况，跨音速压气机，f级j三堆流动，数值模拟 0引 言 叶轮机械多叶片排的三维粘性流计算方法已有许多报道．但这种计算往往需要大量 的时间。由于计算机硬件的限制，目前应用于工程实践还有一定的难度．发展一种既能 包含一些粘性效应又能快速求解的方法就成为当前研究的一个重要诹题． 粘性体积力方法【1】是Denton发展的一种粘性流场快速求解方法．但由于对方程 中的粘性力计算作丁大量的简化处理，实际对能量方程并不求解，在流动现象较为复杂 时。计算往往会出现大的偏差．基于Denton的思想并克服以上不足，本文使用有限差分 的隐式时间推进方法，发展了一种三维雷诺平均Navier-Stokes方程快速求牌方法用于模 拟多叶片排叶轮机械内部流场．叶列间的数据传递采用混合平面法．为突破计算机容量 和速度对叶轮机械多叶片排三维粘性流场计算的限制，利用现有微机资源，采用了租粒 度的微机网络并行计算技术． 由于计算方法采用了简单的壁面函数，固壁边界处网格较稀，对流场中存在强激设 37号 ／边界层干扰的模拟势必存在着不足．因此，本文只对具有详细实验结果的NASA 低展弦比轴流跨音速压气机级70％设计转速下的流场进行了全工况数值模拟，以验证计 算方法和计算程序对只存在徽跨音或高亚因流动的流场的数值模拟可靠性．在计算中通 过改变背压，得到了总压比、总温比及效率总特性曲线．对比实验数据，还重点分析了 其在一些典型工况下的内部流场． 1数理方程和数值方法 1．1控铡方程假设移动系统以常速度n围绕z-轴作定常旋转《对静叶n=0)，则 三维可压缩冒诺平均Navier—Stokes方程无量纲化后在相对参考坐标系中可以写成如下守 本文得到国寒重点基础研究专项经费(G1999022304)和国家自然科学基金(s0076019)资助 83 恒律形式 磐+篆+去口+s=o ㈣ 其中，q=I口l，啦…q5r，F。为对流通量矢量，口为粘性矢量，S为源项，矗e是雷 诺数． 假设靠近固体壁面网格上的第一个点在层流底层或对数率边界层内，则固壁边界处 J： 的剪切应力可由下面三个式子来求解【1 垂直距离，c，为摩擦系数，下标“2”表示离开壁面的第一个网格点．得到治壁面方向 的r之后，经过张量转换，即可得到计算所需的固壁上的各应力分量t，．由于这个显 式的壁面函数公式能较准确地估计固壁面处的剪切应力，壁面边界层内不再需要布置大 量的计算网格点数，因而计算工作量大大减少．而由于假设网格上的第一个点在边界层 内，在实际计算中固壁采用的是滑移边界条件． 由于本方法力求快速求解，我们采用简单的混合长