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

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

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

腐蚀条件下基于Gerber模型的DFR推导及实测

为了研究等寿命曲线模型的选取对细节疲劳额定值计算结果的影响,针对六种典型航空材料对比了Gerber模型和Goodman模型对于高周疲劳数据的拟合精度;推导基于Gerber模型的DFR计算公式、腐蚀折算系数CC的表达式;针对2024-T3铝合金(表面阳极化)进行了预腐蚀0 h、6 h、12 h、24 h、36 h和72 h的疲劳实验并分析预腐蚀72 h的疲劳断口。结果表明:Gerber模型适用于LY12CZ等铝合金,并且在N_95/95>10~5次时,基于Gerber模型的DFR法才能发挥延性材料的潜能;随着预腐蚀时间增长,2024-T3铝合金DFR值下降,基于Gerber模型计算的DFR分别为84.251 MPa、84.721 MPa、79.683 MPa、80.745 MPa、77.026 MPa和74.996 MPa,腐蚀折算系数CC为1.006、0.946、0.958、0.914、0.890,拟合得到DFR随预腐蚀时长的变化曲线是DFR=84.251[lg(t+10)]^-0.15578;断口分析发现预腐蚀产生的蚀坑和材料中的夹杂物会加速疲劳裂纹的形成和扩展,导致结构的疲劳性能降低,但与裸材相比,阳极化过的试件的DFR在腐蚀环境中下降趋势减缓。     

TC4钛合金超声喷丸强化残余应力数值模拟分析

为了探究 TC4 钛合金超声喷丸过程中喷丸参数对残余应力的影响规律,基于 ABAQUS 建立了 TC4 钛合金超声喷丸强化的 3 维有限元模型,分别从超声喷丸模型中振动头振幅、弹丸数量和弹丸直径的变化对 TC4 钛合金表层及亚表层残余应力分布的影响进行了分析。结果表明:随着振动头振幅和弹丸直径的增加,试件表面及亚表面每层残余压应力分布范围及残余压应力值均增大,残余压应力层深度增加;随着弹丸数量的增加,仅能提高每层残余压应力的值,对残余压应力层深度影响较小。     

扰流片式推力矢量喷管气动特性数值模拟研究

为了满足工程应用中扰流片式推力矢量喷管控制率建模以及优化设计的需求,本文选取了三种几何外形的扰流片,并通过数值模拟手段,研究了扰流片几何形状对于轴对称喷管推力矢量气动特性的影响规律,提出了减小扰流片推力损失的设计方法。数值模拟结果表明,推力矢量角与推力损失系数都随着扰流片插入高度的增加而增加;对于矩形扰流片,可以通过增加扰流片宽度的方式减小推力损失,对于扇形扰流片,可以通过减小上圆弧圆心与扇形顶点距离的方式减小推力损失;在插入高度及面积一定时,对比不同形状的扰流片,弧顶矩形扰流片的推力矢量角及推力损失系数均为最大,圆形扰流片均为最小。     

第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。     

固定几何气动矢量喷管气动性能数值仿真

为掌握固定几何气动矢量喷管气动性能,通过CFD数值模拟的方法,研究了主流落压比、扩张段二次流落压比、扩张段二次流角度和引射对固定几何气动矢量喷管轴向推力系数的影响;主流落压比、扩张段二次流落压比和扩张段二次流角度对矢量角的影响;主流落压比、喉道二次流落压比和喉道二次流角度对喉道控制率的影响。结果表明:随主流落压比增大轴向推力系数增大,矢量角减小,喉道控制率减小;随扩张段二次流落压比增大推力系数减小,矢量角增大;随喉道二次流落压比增大,喉道控制率增大;随扩张段二次流角增大轴向推力系数减小,矢量角略有减小;随喉道二次流角增大喉道控制率增大;随引射方式增加喷管推力系数增大。     

复合固体推进剂细观结构建模及脱黏过程数值模拟

为研究复合固体推进剂损伤演化规律,基于分子动力学颗粒填充算法构建了HTPB(hydroxyl terminated polybutadiene)推进剂细观结构模型,通过在AP(ammonium perchlorate)颗粒/HTPB基体界面处引入黏接接触替代传统的黏接单元,并基于Hooke Jeeves的参数优化算法反演得到颗粒/基体界面处内聚力模型参数,利用双线性和自定义指数型损伤内聚力模型模拟了AP颗粒和HTPB基体黏接界面处损伤的萌生、发展、聚合直至宏观裂纹破坏的过程。通过数值仿真与实验结果对比发现,指数型损伤内聚力模型比双线性模型能更准确描述推进剂单轴拉伸过程中颗粒与HTPB基体界面间脱黏过程。最后对比了多阶段加载实验结果与仿真结果曲线,发现两者变化趋势基本一致,最大偏差仅为10%,验证了所建细观模型的可靠性及反演所得界面参数的准确性。      

涡桨飞机螺旋桨滑流气动干扰效应及流动机理

螺旋桨滑流对飞机各气动部件的干扰是涡桨飞机气动设计中面临的难点之一。以某双发涡桨支线客机为对象,采用数值模拟手段分析了滑流对机翼、平尾、垂尾的影响及其流动机理。计算采用基于动态面搭接网格技术的非定常方法,气动力与表面压力分布的计算结果均与实验值吻合良好。研究结果表明：在滑流影响下,全机升力和阻力有所增加,升阻比和纵向静稳定度有所降低,并在无侧滑状态下产生了滚转力矩和偏航力矩;在不同展向位置,滑流对机翼表面流动分离的影响存在显著的差异。在螺旋桨下行运动一侧,滑流的旋转减小了当地迎角,同时桨后气流速度较高,翼面流动分离被有效抑制。在螺旋桨上行运动一侧,滑流的旋转增大了当地迎角,而且桨后气流速度较低,因而翼面流动分离并未得到明显改善;在中小迎角下,滑流对背景飞机平尾当地的动压没有产生影响,但增加了下洗角变化率,进而导致平尾效率降低;垂尾的侧力与偏航力矩是由滑流的侧洗引起的,而滑流的侧洗又与左右两侧机翼升力分布的不对称性有关。     

非轴对称端壁设计的高负荷涡轮气热性能研究进展

非轴对称端壁设计因能够有效地减少涡轮叶栅的二次流损失和提高气动性能而在高负荷涡轮设计中得到应用。简要回顾了涡轮叶栅二次流模型和非轴对称端壁造型方法,重点综述了非轴对称端壁设计的高负荷涡轮气动性能研究进展和抑制端壁二次流的作用机制,介绍了非轴对称端壁设计的高负荷涡轮端壁气热耦合作用的冷却特性研究进展,总结了高负荷涡轮非轴对称端壁设计技术的应用成果,展望了非轴对称端壁设计在高负荷涡轮的高效气动和冷却布局应用方面需要深入研究的内容。