基于广义Richardson外插方法阻力预测精度分析.PDFVIP

航 空 学 报 Jan. 25 2014 Vol.35 No.X: XX-XX Acta Aeronautica et Astronautica Sinica ISSN 1000-6893 CN 11-1929/V hkxb@ DOI: 10.7527/S1000-6893.2015.0005 基于广义 Richardson 外插方法的阻力预测精度 分析 1 1, *，张宇飞 1 李钊 ，陈海昕 1. 清华大学 航天航空学院，北京 100084 摘 要：利用CFD方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二 维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA0012翼型无粘绕流和平板湍流 边界层两个算例中验证了NSAWET程序和广义Richardson外插方法，分别得到数值算法预测压差阻力和摩擦阻力能达到 的名义精度。进而模拟三维CRM翼身组合体绕流，得到的阻力名义精确值在DPW5 的统计误差带范围之内；综合DPW5 的计算结果来看，不同CFD解算器的结果之间存在一定差别，阻力预测精度总体上不符合二阶。可见，标准Richardson 方法采用的二阶精度假设难以普遍适用，有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度， 采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性，可以进一步降低网格对阻力 预测的影响。 关键词：广义Richardson外插方法；阻力预测；网格收敛性；验证与确认；计算流体力学；CRM翼身组合体 中图分类号：V211.3 文献标识码：A 文章编号：1000-6893 （2014 ）XX-XXXX-XX 近年来，计算流体力学（Computational Flu- 发展和应用的基本问题之一。国内也非常重视 id Dynamics ，CFD ）作为一种气动性能预测的有 CFD 的验证和确认工作[5-7] ，并且主要结合前四 力工具，在飞行器气动设计中得到了广泛的应 次的 DPW 的工作开展了阻力预测的相关研究[8- 用。而对 CFD 方法的可信度研究在其中发挥了 10] 。本文则针对第五次阻力预测会议（DPW5 ） 关键和基础的作用，也是目前工业界和学术界共 的三维翼身组合体进行阻力预测的精度分析。 同关注的研究热点之一[1] 。近年来，美国美国航 DPW5 延续了前一届阻力预测会议采用的通 空航天学会（American Institute of Aeronautics and 用 研 究 模 型 （ Common Research Model ， Astronautics ，AIAA ）发起举办了五次阻力预测 CRM ），其关键的不同之处在于统一了计算网格 会议[2] （Drag Prediction Workshop ，DPW ）、两 [11] 。主办方针对 CRM 翼身组合体的简单构型提 次增升装置预测会议[3] （High Lift Prediction 供了包括多块结构网格、重叠网格和非结构网格 Workshop ，HiLiftPW ）和两次推进气动会议[4] 在内的共同网格系列，对固定升力巡航状态的阻 （Propulsion Aerodynamic Workshop ，PAW ）等 力预测能力进行检验。藉此减小不同网格引入的 一系列旨在验证和确认 CFD 方法的工作会议。 不确定性，阻力预测的分布范围和标准偏差也 特别地，CFD 方法对飞机阻力的预测能力，已经 由 DPW4 的 40.9 和 8.1 个阻力单位（Count，阻 成为认识和衡量 CFD 发展水平并推动其进一步 力系数 0.0001 ）分别减小到 27.