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

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

非结构混合网格自适应并行技术

计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

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

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

烧蚀热响应计算中的不确定性传播分析方法研究

发展了基于代理模型的不确定度量化方法，以少量确定性模拟结果为样本，通过代理模型来构建目标变量和不确定性输入参数的响应关系，以此来分析目标变量的不确定性信息，包括均值、标准差、概率密度分布等统计信息以及各个输入参数对目标变量的敏感性大小。本文以烧蚀热响应计算为例，通过分析材料物性参数的不确定性对烧蚀热响应预测的影响，表明该方法具有准确度高、效率高的特点，为工程中的高维不确定度量化问题提供了很好的解决思路。     

基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。