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

为了提高叶轮机内周期性非定常流数值计算效率,发展了一种计算振荡叶栅非线性非定常流的谐波平衡方法,将周期性非定常流动参数采用傅里叶级数表示,从而将传统的定常流动数值计算方法应用于非定常流问题求解中,大大提高了非定常流计算效率。计算结果表明,发展的谐波平衡方法与线性方法、非线性时域法计算结果吻合较好,且其计算效率比时域法高1个数量级。另外,系统研究了谐波阶次对计算效率与计算精度的影响,认为在气流未出现严重分离时,应采用低阶谐波平衡方法,以减少计算耗时。此外,从压气机振荡叶栅非定常流场计算结果可以看出,叶片吸力面诱导涡的生成频率与振荡频率并不一致,涡生成频率要高于振荡频率。     

两级对转风扇非定常特性分析

为研究两级对转风扇的非定常特性,对该风扇进行了设计转速下的定常和非定常数值模拟。结果表明:非定常得到的压气机裕度(15.3%)比定常(14.5%)的更高,主要原因在于近失速时,上游尾迹和泄漏流能改善叶尖吸力面分离;非失速工况时,上游尾迹和泄漏流与下游泄漏流、激波之间复杂的非线性作用会引起熵增,从而导致叶尖流动损失的增加;此外,所有工况下,上游尾迹对叶片根部角区分离都起到了抑制作用,降低了损失。     

局部进气ATR涡轮性能影响因素数值模拟研究

为研究局部进气条件下ATR涡轮内的流动特性,对局部进气涡轮进行定常和非定常数值模拟.定常计算采用冻结转子方法,改变动叶与静叶的周向相对位置研究涡轮内流场特征;改变进气扇区布局、进口参数和局部进气度,研究其对涡轮性能的影响.研究结果表明,对于给定局部进气度,进气扇区整体布局时,涡轮性能的优越性最大;在局部进气ATR涡轮的设计中,需要综合考虑流量、效率和功率3个因素来决定局部进气度.不同位置定常流动结果的平均值与非定常流动的时均值存在差异,非定常计算更好地展示了静叶尾迹和流动积累效应的影响,因而更加可靠.     

跨音速压气机轮毂型线优化设计

针对跨音速转子Rotor37,采用Numeca软件中的Design3D模块对其叶片区域的轮毂型线进行了优化设计计算。模块结合了CFD、人工神经网络和遗传算法,优化效果显著。轮毂型线优化后的Rotor37相比原始设计,在整个设计转速范围内的等熵效率得到了明显提高,而裕度与优化前相比略微提高。通过流场分析,其性能的改善主要源于叶片根部区域的激波前移和尾缘处附面层分离的减弱。     

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

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

振荡叶栅三维非定常流动数值分析

以全三维、非线性、非定常欧拉方程计算作用于振荡叶栅的非定常气动载荷。欧拉方程的数值离散采用了具有时、空二阶精度的MacCormack显式分裂的两步格式,用时间推进方法得到其数值解。应用动网格技术捕获叶片的运动,保证叶片振动过程中网格边界与叶片表面保持一致。以相移周期条件处理具有不同相位的叶片振动问题,同时发展了伪边界形状修正方法以简化周期边界条件处理。     

激波对跨声速风扇叶片颤振特性的影响

为了研究激波对叶片气动弹性稳定性的影响,基于谐波平衡计算方法,结合结构动力学分析,采用单向耦合方式下的能量法对跨声速风扇的颤振特性进行了数值研究。计算结果表明:设计点工况下,二阶弯曲模态时叶片的最小气动阻尼系数小于0,数值为-0.338%,意味着叶片存在失稳风险;此外,100%等转速线不同工况下叶片的颤振特性不尽相同,激波位置变化对气流做功的影响是导致叶片气动弹性稳定性差异的主要诱因。     

针对叶轮机流场分析的不同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湍流模型不仅能极大地提高计算的压比和流量范围,而且能提高计算的稳定性和失速工况下的计算准确性,但计算的堵塞流量和效率较小。     

振动叶栅中尾迹/激波相干的数值研究

研究了尾迹/激波相干对振动叶栅非定常流动的影响.数值方法利用了先进的谐波方法.该方法基于傅立叶变换,把叶栅非定常流变量分解为一时间平均量和若干谐波分量,继而,时域非定常N-S方程被转化为频域定常傅立叶系数方程.定常傅立叶系数方程的求解采用了伪时间推进方法.为加速数值解收敛,采用了本地时间步长和多重网格技术.通过NASA67转子算例结果表明:上游尾迹对振动叶栅非定常流动具有显著影响.当存在尾迹相干时,叶片通道内的激波在轴向方向上表现出较大的漂移运动,叶片所受的非定常力幅值变大,同时失去了时间周期性.      

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

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