化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

细长体飞行器自由边界热模态试验与仿真

随着耐热承载一体化材料在新型高超声速飞行器上的应用,承力结构的工作温度不断提高,各类热模态特性逐渐得到关注。针对非平面形状的细长体飞行器自由边界条件下的热模态特性开展了研究。通过研究模拟气动加热条件的圆筒形加热笼、模拟自由边界的耐高温柔性支撑、非接触激光多普勒测振、耐高温激振杆激励等试验方法,获得了细长体结构自由边界条件下随温度变化的前3阶模态变化情况。结果表明：对此类薄壁长圆筒类结构,温度升高对模态频率影响可以超过6 Hz。开展有限元仿真,并与试验取得的热模态结果的变化规律进行对比。结果表明：建立考虑温度对结构弹性模量、热应力影响的壳单元模型,能够较好地预测出前3阶模态频率在全部受热时间范围内的最大下降量,可为高超声速飞行器控制系统设计时的拉偏范围提供参考。     

轮缘封严气流与转子干涉损失机理的数值研究

为了研究涡轮转静盘腔中轮缘封严气流与下游转子的干涉机理和损失机制,在不同封严流量下,针对轮缘封严气流对下游转子气动性能和流场分布的影响进行了数值模拟。结果表明,随着封严流量增加,涡轮级效率不断降低。在轮毂到10%叶高和30%～60%叶高,流量系数增加。封严气流主要通过在转子进口封严出流不同的分布和在通道输运过程中与转子通道二次流的交互作用对转子的气动性能产生影响。同时,封严气流引起下游转子气动性能的降低,存在三个损失机制:一是封严腔体摩擦泵效应减弱,增加了封严出流与主流之间周向速度差所引起的粘性剪切损失;二是不同分布的封严出流造成进气负攻角和堵塞效应,降低了转子的做功能力;三是不同分布的封严出流在输运过程中与转子通道二次流的交互作用,导致了二次流损失加剧。     

对转压气机外伸激波对叶顶泄漏流影响

摘要:为了揭示对转压气机下游转子外伸激波对上游转子泄漏流的影响规律,针对上游转子叶顶间隙分别为0.2、0.5、0.8 mm的对转压气机开展了非定常数值模拟研究。研究发现:受下游转子外伸激波掠扫影响,上游转子尾缘附近压力面会形成弱压缩波,且随上游转子泄漏流增强而逐渐减弱;而该外伸激波在上游转子尾缘附近吸力面,会形成与型线切向相垂直的较强压缩波,且其位置基本不受叶顶间隙大小影响;外伸激波使上游转子尾缘附近吸、压力面压差增大,叶顶泄漏流增强,进而导致其损失增大;随着叶顶间隙增大,上游转子叶尖区弦长前半段压力波动的频率,由通道激波转为叶顶泄漏流主导,且呈现减小的趋势,而弦长后半段压力波动的频率主要由外伸激波主导,且基本不变。      

三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

弹性空腔流致噪声/结构振动特性试验

高速空腔中经常存在高强度且多频率分量的流致噪声,空腔噪声与结构振动之间耦合效应严重,甚至可能发生结构共振。为此,在0.6m×0.6m高速风洞中,通过调整空腔底板厚度,改变其结构固有频率,模拟空腔流致噪声/振动相互作用。利用脉动压力和振动加速度测试技术,获取亚跨声速条件下,弹性空腔流致噪声特性及其结构振动响应特性。马赫数变化范围为0.6~1.2。结果表明,当振动强度较弱时,结构振动对空腔噪声影响较小,而空腔噪声对结构振动影响较大,在噪声载荷主频位置,振动谱出现峰值并且噪声/振动相关性达到最强;此外,空腔结构振动还与其固有频率特性密切相关,振动主要以低阶模态为主。     

一种涡轴发动机控制系统自适应扭振抑制方法研究

为了解决直升机机动飞行时涡轴发动机控制系统扭振不稳定性问题,根据扭矩动力传递链特点,设计了一种基于最小均方差(LMS)的自适应滤波器,以抑制自由涡轮转速中的扭振动态,并与常规的陷波滤波器效果进行了对比。仿真结果表明,在定旋翼转速下,陷波滤波器和基于LMS的自适应滤波器均能抑制1.90Hz对应的固定扭振分量,但后者能使扭振幅值降低至5%以下,滤波效果明显优于前者;同时,变旋翼转速飞行时,发动机端的扭振基频率会随着转速变化在1.30~2.20Hz变动,此时常规陷波滤波器不再适用,而基于LMS的自适应滤波器仍可显著抑制扭振动态。     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

锥形间隙孔型阻尼密封动力特性分析及抑振机理

应用非定常动网格技术建立了锥形间隙孔型阻尼密封动力特性多频椭圆涡动求解模型,研究了锥形度和涡动频率对孔型阻尼密封泄漏特性与动力特性的影响,分析了锥形间隙孔型阻尼密封的抑振机理。结果表明:锥形间隙增强了孔型阻尼密封的泄漏直通效应,增加了密封泄漏量。在平均间隙相等的情况下,收敛间隙孔型阻尼密封泄漏量小于发散间隙孔型阻尼密封;在同一涡动频率下,刚度系数随着锥形度的增加而增大,阻尼系数随着锥形度的增加而减小;收敛间隙孔型阻尼密封具有较大的正直接刚度,增加了转子系统的固有频率,发散间隙孔型阻尼密封产生较大的负直接刚度;收敛间隙孔型阻尼密封具有较大的有效刚度和穿越频率,较小的有效阻尼,发散间隙孔型阻尼密封具有较大的有效阻尼,较小的有效刚度和穿越频率;在锥形间隙孔型阻尼密封结构设计时,需要同时考虑转子系统固有频率和同心度的影响。