《一种快速模拟振荡叶栅非定常流的数值方法.》.pdfVIP

第 32卷 第 4期 空 气 动 力 学 学 报 Vol_32．No．4 2014年 O8月 ACTA AERoDYNAM ICA SINICA Aug．，2014 文章编号：0258—1825(2014)04～0481—07 一 种快速模拟振荡叶栅非定常流的数值方法 施永强，杨青真，黄秀全，郭 霄 (西北工业大学 动力与能源学院，陕西 西安 710072) 摘 要 ：为了提高叶轮机 内周期性非定常流数值计算效率 ，发展 了一种计算振荡叶栅非线性非定常流的谐波平衡 方法 ，将周期性非定常流动参数采用傅里叶级数表示 ，从而将传统的定常流动数值计算方法应用于非定常流 问题 求解 中，大大提高了非定常流计算效率 。计算结果表明，发展 的谐波平衡方法与线性方法、非线性时域法计算结果 吻合较好，且其计算效率比时域法高 1个数量级。另外 ，系统研究了谐波阶次对计算效率与计算精度的影响，认为 在气流未出现严重分离时，应采用低阶谐波平衡方法 ，以减少计算耗时。此外，从压气机振荡叶栅非定常流场计算 结果可以看 出，叶片吸力面诱导涡的生成频率与振荡频率并不一致 ，涡生成频率要高于振荡频率 。 关键词 ：振荡叶栅 ；谐波平衡方法 ；非定常流 ；双时间法 ；形状修正法 中图分类号 ：V221．3 文献标识码 ：A doi：10．7638／kqd1xxb一2012．0146 常流动计算效率 。与时间推进方法相 比，这种方法计 O 引 言 算速度是传统方法 的7～9倍 ]，因此被广泛应用于 叶轮机 内流场具有高度复杂的非定常性 。传统 非定常流动数值模拟。 的二维轴对称气动设计理论和全三维定常流场计算 虽然双时间方法能够加速非定常流动数值计算， 方法 ，都无法反映这一真实的流动特性 。用这种定常 提高计算效率 ，但对于复杂的叶轮机械非定常流动来 流体系处理具有强烈非定常特征 的叶轮机设计问题 说 ，其计算工作量还是相当巨大 的。为此 ，国外学者 时，从一般意义上讲 ，会导致对非定常流 “潜力”的某 在这方面做了大量的研究工作，发展了高效的频域计 种忽略而使性能有所损失 。叶轮机械中的叶片是弹 算方法求解非定常流动 。Verdon ，WhiteheadL4和 性体 ，在气动力、结构等因素影响下，会发生振动。强 Hall[5等人分别提 出了模拟二维叶栅非定常流动的 的振动现象会导致 叶片在短 时间内断裂 ，造成严重的 时间线化方法 ；Hall和Lorenee发展 了三维线化非定 后果 。为 了深入研究叶栅非定常流动，探索叶轮机械 常流动计算方法 ；2005年 ，Ekici和 Hall采用时间 强迫振动 、颤振的发作机理，就必须对振荡 叶栅非定 线化方法分析了多级叶轮机械颤振问题l7]。时间线 常流场以及叶片所产生 的周期性非定常气动力 的性 化方法在求解叶轮机械非定常流动问题时，仅需要模 质给予准确、高效地模拟与预测 。 拟单个通道的流动状况，且其计算时间与定常流动计 目前 ，利用计算流体力学方法对振荡叶栅非定常 算时间相当，因此 ，被广泛应用于叶轮机非定常流动 绕流及叶片所受气动力进行数值模拟 ，进而分析振荡 数值模拟中。然而 ，时间线化方法也存在着不可忽视 叶片与流场间的相互作用 ，已成为研究叶片强迫振动 的缺点——忽略了非定常流动的非线性现象。为此 ， 与颤振 的重要手段 。在振荡叶栅非定常流动数值模 Hall提出了谐波平衡方法用于模拟 叶栅中非线性非 拟研究方面 ，传统 的分析方法是 时间推进