叶轮机械全环非定常大规模并行模拟程序设计

为了更加精确地模拟叶轮机械内部复杂的三维非定常流动,开展了全环网格模拟方法探索,提出了针对滑移面的守恒变量插值和高效的大规模并行计算解决方案。基于千万亿次高性能计算平台,自主开发了非定常数值模拟程序,采用单级压气机、1.5级涡轮和3.5级高负荷压气机三个算例进行了程序的验证与确认,以及并行效率测试工作。模拟采用的最大网格数为2亿,计算核心数为4096个。结果表明:程序计算得到的整体气动性能与试验数据吻合良好;程序所采用的全环非定常计算模型能保证交界面处物理量的连续,削弱非物理熵增,得到更加可信的流场解;所采用的METIS的网格分区方法以及MPI并行策略使程序具有良好的负载均衡和并行效率;同时程序对超大规模的复杂问题具备良好的可扩展性和适应性,能够满足工程实际问题的流动分析和叶轮机械的精细流动研究要求。     

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

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

武器分离及舱门开启过程数值模拟研究

基于非结构重叠网格技术,发展了一套能较好模拟空腔流动和武器分离问题的数值计算方法.通过模拟典型空腔流动和WPFS模型分离轨迹,并与试验结果比较,说明本文方法可较好模拟空腔复杂流动和武器分离轨迹.采用发展的数值方法对某内埋弹舱舱门开启过程和内埋武器分离轨迹进行了数值模拟,研究了舱门运动的非定常效应对舱门气动特性的影响.研究表明:一定条件下,舱门运动产生的非定常效应对舱门气动特性影响较大,舱门和舱门控制系统设计时必须考虑；在本文给定条件下,弹舱内的内、外侧导弹均能安全分离.通过研究,为内埋弹舱的舱门设计和内埋武器分离安全评估提供了技术手段和依据.     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放。该方法在重叠网格"挖洞"时,以"离壁面距离"为挖洞判据;在建立重叠区域插值关系时,对Neighbor to Neighbor(N2N)和Alternating Digital Tree(ADT)两种方法的效率进行了比较。采用机翼/挂架/带舵外挂物模型开展了验证计算。在此基础上,完成了某飞行器外挂投放的初步计算。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

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

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

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

基于混合平面法的跨声速轴流压气机三维数值模拟研究北大核心CSCD

利用自主开发的三维流场求解器MFlow,基于混合平面法技术对航空发动机压气机流场进行了数值模拟。在转/静子交界面处的不同的展向位置,将流动周向平均后进行传递,使得非定常流动问题简化成为"准定常"问题,从而将计算过程简化为只需要对一个叶片流道进行计算,大大减小了计算量。通过计算得到了压气机的总体工作特性和几个关键气动参数的展向分布曲线,并将结果与试验数据进行了对比分析。计算结果和试验数据基本吻合,验证了计算程序的适用性和可靠性。