基于高效特征值分析方法的旋转失速先兆预测

旋转失速限制了压气机稳定工作的范围,对其进行深入理解并实现准确预测是控制失速、提高裕度的关键。现有模型大多基于一定程度的几何和流动简化,不考虑三维效应及流动复杂性,直接应用于三维压气机失速预测时仍面临巨大挑战。同时,尽管试验测量和模拟仿真水平不断提高,试验和数值模拟多为唯象研究,缺乏对压气机流动失稳根本原因的揭示。此外,由于三维复杂流动精细化测量和高保真模拟的复杂性,大多数失速研究针对某一压气机若干孤立工况开展,缺乏系统的参数化研究,难以提炼出旋转失速关键影响因素。为弥补试验测试空间分辨率低和非定常流动模拟成本高的缺陷,提出了一种基于三维流动方程高效特征值求解的全局稳定性分析方法。一方面可以获得试验测量难以达到的空间分辨率,另一方面能够以比非定常模拟小2～3个量级的成本获得丰富的三维流场小扰动发展过程。针对某典型跨声速压气机环形叶栅,所发展的分析方法计算成本仅为定常特性线计算的28%,相比于非定常计算实现了约155倍的加速,为压气机旋转失速准确快速预测和机理研究提供了重要的研究工具。     

自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

叶轮机周期性非定常流高效降阶数值计算方法研究进展

简要介绍了叶轮机内非定常流对叶轮机性能的重要影响,总结了叶轮机内动静相干和叶片振动诱发的非定常流的时空周期性特点,并指出这些时空周期性特点可以加以利用, 发展分析这些非定常流的高效降阶数值计算方法。本文根据不同数值降阶方法的特点对其进行了归类,详细分析了各自的优缺点及在气动弹性力学领域的应用,并进一步展望了降阶数值方法在未来的几个关键技术及发展趋势。     

叶轮机气弹耦合方程求解的时间推进方法评估

为了提高叶轮机颤振预测时的计算效率,同时兼顾求解的稳定性,有必要选取最适合的时间推进方法。基于降阶结构 动力学方程的双向流固耦合方法常用于压气机气弹性能分析,针对气动弹性控制方程中结构动力学方程考查了多种时间推进方 法对计算结果的影响。通过单自由度弹簧系统的数值求解探讨了不同时间推进方法的特点;选取实际压气机算例NASA Rotor 67 进行颤振流固耦合分析,进行大量数值试验找到了各时间推进方法需要的最大时间步长,以此为依据对比出不同时间推进方法的 计算效率。结果表明：采用单自由度弹簧系统和压气机颤振耦合求解得出的结论基本一致,说明时间推进方法的性能更多地与方 法本身的数学性质有关,但是流固耦合求解时流场和结构的信息交互会带来一定影响。以实际颤振算例叶片瞬态响应振幅的对 数衰减率作为衡量精度的标准,经典龙格-库塔方法是颤振耦合计算达到相同精度耗时最少的方法,其次是4阶隐式Adams方法 和Newmark方法。     

针对叶轮机流场分析的不同Spalart-Allmaras湍流模型评价

Spalart-Allmaras（SA）湍流模型在Reynolds-Average Navier-Stokes（RANS）方程求解中得到了广泛的应用。针对叶轮机数值模拟，国内外学者提出了很多修正的SA模型，但它们之间的对比少有研究。为此，综述了原始的SA（SA-standard），SA-neg，SA-Helicity，SA-Ning及SA-R五种不同SA湍流模型。针对NASA Rotor 67和NASA Rotor 37，考察这五种湍流模型对数值计算稳定性的影响及对叶轮机内部流动细节捕捉的能力。研究表明：SA-standard和SA-neg湍流模型计算结果几乎相同；相比于SA-standard湍流模型，SA-Ning湍流模型计算的压比、效率和堵塞流量均较大；SA-R湍流模型模拟的激波/边界层干涉更严重，分离区域更大，计算的压比小，模型修正效果不理想；SA-Helicity湍流模型不仅能极大地提高计算的压比和流量范围，而且能提高计算的稳定性和失速工况下的计算准确性，但计算的堵塞流量和效率较小。