数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

发动机进排气效应对民机构型气动特性影响

基于非结构混合网格技术,通过数值求解Navier-Stokes方程,分析了航空发动机进排气效应对民机构型气动特性的影响。通过设置合适的进排气边界条件来模拟发动机进排气效应的影响,采用单独发动机短舱风洞试验模型,验证了计算方法的可靠性。在此基础上,分析了发动机进排气效应对翼吊式和尾吊式两种典型民机构型气动特性的影响,结果表明:翼吊式民机构型发动机进排气效应对升力的干扰主要由发动机尾喷流的引射效应对机翼的干扰引起;尾吊式民机构型随着发动机进气流量的增大,机翼上表面压力逐渐减小,激波位置逐渐后移;不同发动机进气流量会对飞机的阻力特性产生较大影响,在飞机详细设计阶段需充分考虑机体与动力装置之间的干扰影响。      

FlowStar:国家数值风洞(NNW)工程非结构通用CFD软件

计算流体力学（CFD）仿真软件是流体相关的数学物理知识和工程实践经验的数字化表达,是工业数字化转型的重要助推。然而,大型工业CFD软件研发难度极高,需要同时兼顾功能多样、系统稳定、性能优越、交互友好等特征。依托国家数值风洞（NNW）工程,研发出一款通用流场模拟软件NNW-FlowStar,并在航空、航天等工业部门大力推广使用。软件基于非结构有限体积求解方法和大规模并行计算技术开发,结合现代化软件工程思维设计,具备先进的数值方法、高效的计算效率和友好的用户操作界面,可满足各类复杂外形的高效气动模拟。独特的重叠网格技术配合六自由度运动模块,可帮助实现武器分离、舱门定轴转动等各类气动-运动协同仿真需求。多类标模案例和复杂工程应用表明,FlowStar软件算法鲁棒、精度可靠,是一款高精度、高效率、高可靠性的通用CFD仿真软件。通过对软件的架构设计和功能应用进行介绍,使相关从业人员能更好地了解FlowStar软件,最终促进国产自主CFD软件生态的良性发展。     

动量增升高升阻比飞行器横航向稳定性研究

考虑高超声速飞行器装填、防热及操稳等多约束条件下的实用化设计需求,提出一种新型下反式高升阻比滑翔飞行器气动布局。借鉴乘波体飞行器的高升力设计方法及动量增升原理,该新型飞行器采用下反式后掠翼构型,上表面为光滑倒圆“Λ”形设计,下表面为内凹的装填空间,为尽可能避免长时间超远距离高超声速飞行带来的防热负担,采用后缘体襟翼及体侧扩张方向舵面设计。采用数值计算方法对所提布局思路进行验证分析,计算结果表明：在飞行高度为40 km,Ma=10的条件下升阻比可以达到4.48左右,在一定迎角范围内均具备很高的气动效率,验证了所提布局的有效性。同时重点针对所提布局共性的横航向稳定性问题基于数值模拟方法探讨了3种不同优化改进方案的效果及可行性,并采用风洞试验对两侧翼梢V尾的控制方案进行横航向稳定性控制效果的试验验证,结果表明：采用两侧翼梢V尾的控制方案是实现横航向稳定性控制的较优方案。     

DLR-F6翼身组合体数值计算

使用自行研制的混合网格亚跨超声速流场解算器程序MFlow计算了AIAA第三届阻力会议提供的DLR-F6及其带整流装置FX2B情况下的阻力。重点分析了两者的网格收敛特性、阻力极曲线以及压力分布等,并与阻力会议的各个软件的计算结果进行比较。详细分析了DLR-F6翼身结合处后缘附近网格精度对分离气泡计算的影响。计算结果表明,MFlow程序的计算精度与国外软件相近,能够为阻力计算提供比较精确的结果。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

分离涡模拟方法在内埋弹舱机弹分离模拟中的应用

为了研究内埋弹舱非定常流场对武器投放分离特性的影响,基于分离涡模拟（DES）方法和非结构重叠网格技术,建立了内埋武器机弹分离高精度数值模拟方法。在此基础上,首先通过典型空腔流动和机弹分离算例验证了数值方法的有效性,然后采用平均雷诺应力模拟（RANS）和DES两种方法探究了机弹分离过程中弹舱腔内涡结构的演化历程,分析了两种方法在机弹分离流场结构模拟中的差异,研究了内埋弹舱流场结构对导弹分离特性的影响。计算结果表明：机弹分离过程中,高强度涡结构逐渐向弹舱腔体后缘聚集,同时弹舱剪切层流动结构被导弹破坏,导弹出弹舱阶段俯仰方向所受气动力矩急剧变化。相比RANS方法,DES方法捕捉到的高强度小尺度涡结构对导弹压力分布影响较大,使得分离后期导弹俯仰角、俯仰角速度产生明显差异,DES方法更适合内埋武器舱分离特性精细化研究。该研究结果可为内埋武器机弹分离研究的计算方法选择和安全分离策略提供一定的参考依据。     

多场耦合效应对高超声速进气道入口参数影响

为了研究多场耦合效应对高超声速进气道入口参数的影响,采用自主开发的热环境/热响应耦合计算分析平台FL-CAPTER,对吸气式高超声速进气道前体进行了数值仿真研究。介绍了采用的多物理场耦合分析策略及不同物理场求解方法,通过圆管和两级压缩楔外形,初步验证了多场耦合分析方法的可靠性。以此为基础,研究了进气道前体在长时间巡航飞行条件下的结构温升情况和宏观变形量,分析了进气道结构变形对入口参数的影响。结果表明:进气道前体迎风区域和背风区域不均匀的温度分布引起热应力变化,进气道前体压缩面在多场耦合效应作用下上翘约20mm,考虑变形影响后,进气道偏离设计状态,激波边界层干扰效应增强,喉道附近的分离区域有所增大,进气道入口的质量流量增加约4.2%,喉道平均马赫数降低,静压升高,总压恢复系数降低。