基于改进牛顿法的航空发动机起动建模技术

为进行涡扇发动机起动过程性能仿真研究,在发动机慢车以上部件级模型的基础上,通过指数外推法获得低转速部件特性,在建立起动附件模型的同时,修正起动过程各部件总压恢复系数,从而建立起包含地面起动过程的全状态性能模型。在求解发动机各部件共同工作非线性方程组时,采用基于传统牛顿拉夫森迭代求解方法的改进方法,即在第1次完成中心差分得到Jacobian矩阵后,使用修正项迭代更新上一步的Jacobian矩阵。结果表明:该建模方法可以有效地建立起动过程性能模型,并大幅改善模型的执行效率、提高模型仿真求解速度。     

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

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

热冲击后气膜冷却叶片结构强度流固热耦合仿真

研究了冷气流量对气孔周围热应力的影响,为气膜冷却叶片可靠性设计提供参考。改变气孔的孔径,并建立有限元模型,结合有限元/边界元理论,通过流固热三场耦合技术获得热冲击后的叶片最大温度、温度不均衡程度及最大热应力。研究表明:增加冷气量有利于改善叶片冷却效率降低叶片温度,但也会使叶片温度不均衡程度增加,加剧尾缘气孔内的热应力载荷;增加前缘气孔直径可提升66%的平均冷却效率,有利于减缓气孔内的热应力,增加尾缘气孔的直径对冷却效率及热应力的影响均较小。此外,数值计算结果与试验及解析解较为吻合,对于气膜冷却叶片结构设计具有参考价值。      

凹腔对一体化支板火焰稳定器燃烧性能的影响

在来流温度为780~850℃、来流马赫数为0.16及油气比为0.002~0.006的条件下,试验研究了凹腔对喷油/稳定一体化支板火焰稳定器燃烧效率及熄火性能的影响,并结合数值模拟进行辅助分析。结果表明:在不同油气比条件下,带凹腔的一体化支板火焰稳定器均能实现稳定高效燃烧;不带凹腔的一体化支板火焰稳定器燃烧效率始终低于带凹腔的一体化支板火焰稳定器,随着油气比的增加,两者燃烧效率差距逐渐缩小;带凹腔的一体化支板火焰稳定器较不带凹腔的一体化支板火焰稳定器有更好的熄火性能;凹腔结构促进了燃油雾化与蒸发,从而提高一体化支板火焰稳定器的燃烧性能。      

三维并联式TBCC发动机排气系统设计与实验

为实现三维并联式TBCC排气系统的设计,利用基于二维最大推力理论的密切法设计了圆转矩变截面非对称冲压发动机喷管,并在冲压发动机喷管上壁面为三维曲面的条件下,利用转动加滑动共用面的调节方式实现了涡喷发动机喷管的内型面设计及其喉道面积的调节,从而完成了排气系统设计。随后对设计结果进行了冷流实验验证。结果表明,这种排气系统设计方法具有可行性。在冲压单独工作状态,壁面压力分布的数值计算与实验结果吻合较好,相对方均根偏差低于10.8%;在涡喷单独工作状态及共同工作状态,上壁面压力在激波附面层干扰区域内偏差较大,最大相对偏差可达30.6%。但壁面压力分布的数值计算与实验结果的整体变化趋势一致。      

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

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

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

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

结冰风洞中液滴相变效应数值模拟

为明晰结冰风洞中液滴相变效应,发展了基于Euler法的气液两相传质传热耦合流动计算方法,探索了3 m×2 m结冰风洞主试验段构型相变效应对液滴传热过程的影响,开展了参数影响研究,评估了试验段内液滴过冷状态。结果表明：构型内液滴经历了先蒸发后凝结两个阶段,蒸发效应促进了准一维传热阶段液滴温度的下降趋势,使液滴温度趋于湿球温度,而凝结效应则抑制了三维收缩阶段液滴温度的下降趋势,进而增大了试验段液气温差,影响液滴过冷状态;增大初始相对湿度和试验段气流速度,会导致蒸发效应减弱而凝结效应增强,进而增强了相对湿度和气流速度对液滴过冷状态的影响程度,与此相反,初始液滴尺寸的增加,则会导致蒸发效应和凝结效应均减弱,进而减弱了液滴尺寸对液滴过冷状态的影响程度;在典型工况下,小尺寸液滴（液滴直径范围为40~100 μm）在高风速时（试验段气流速度大于164 m/s）将偏离过冷状态（液气温差超过3℃）。     

感性耦合夹层等离子体隐身天线罩电磁散射分析

设计了一种石英夹层感性耦合等离子体（ICP）隐身天线罩模型,采用有限元与Z变换时域有限差分（ZT-FDTD）联合仿真的方法,建立了ICP放电的流体模型,得到不同气压及功率下与电磁散射相关的电子密度空间分布,在此基础上建立Z变换时域有限差分法模型,对石英夹层等离子体隐身天线罩的宽频段后向散射进行计算,同时利用微波干涉法及XFDTD软件对算法及程序进行验证。结果表明：感性耦合等离子体可产生较高密度等离子体,能有效实现雷达散射截面（RCS）的减缩,在气压为2 Pa时,碰撞衰减较弱,等离子体密度分布较均匀,衰减带宽集中在等离子体振荡频率附近,功率增加会使衰减峰值向高频方向移动,气压为20 Pa时,碰撞衰减增强,且等离子体密度分布有较大梯度,衰减带宽有效增加,同时RCS曲线的波动特性加强。