DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

运输类飞机巡航阻力CFD计算分析

针对运输类飞机设计过程中的巡航阻力准确预测问题,采用基于多块结构网格求解三维雷诺平均Navier-Stokes方程的并行流动求解器“CCFD-MB”,对AIAA阻力预测小组提供的典型运输类飞机DLR—F6翼身构型进行了CFD计算分析,比较了小同规模网格及不同湍流模型对计算结果的影响;并通过与DLR-F4翼身构型计算结果...     

基于非结构平台的DLR－F6标模阻力预测

采用自研的非结构网格解算器UNSMB进行了AIAA第三届阻力会议提供的DLR－F6翼身组合体的阻力计算验证。重点分析了模型的网格收敛特性、升阻力曲线以及压力分布等，并把计算结果与阻力预测会议上各个软件的计算结果以及试验数据进行比较，在此基础上分析计算结果。分析结果显示，非结构混合网格解算器的计算结果与各个软件的计算结果以及风洞试验数据吻合度较好，一定程度上验证与确认了解算器的阻力预测精度。     

基于混合网格的DLR-F6模型跨声速阻力预测粘流计算

该文详细介绍了针对第二届AIAA阻力预测会议公布的标准模型及其网格进行的基于混合网格的纳维尔-斯托克斯粘流计算考题,我们将自主开发的WoF90程序的计算结果与DLR的TAU程序以及FOI的EDGE程序的计算结果进行了对比分析,发现三种程序的升力系数普遍比实验结果偏高,但偏高的程度大致相同,阻力计算的精度范围大体一致。结果确认了300万网格点的计算结果尚不足以得到网格收敛的效果,阻力预测需要更大规模的网格。     

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

NASA桘CRM阻力预测模型的计算验证

采用自研的非结构网格解算器UNSMB 进行了AIAA第四届阻力会议提供的NASA-CRM 翼身组合体（ WB）以及翼身组合加平尾（ WBT）两种构型的计算验证。重点分析了WBT模型的网格收敛特性、升阻力曲线、雷诺数效应、下洗效应以及压力分布等,并把计算结果与阻力预测会议上ONERA的计算结果进行了对比。分析结果显示,非结构混合网格解算器的计算结果与ONERA的计算结果吻合度较好,同时在一定程度上验证与确认了解算器的阻力预测精度。     

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

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

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

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

NASA共同研究模型翼身组合体数值模拟研究

采用结构网格求解器WiseManPlus软件,对DPW5提供的NASA共同研究模型CRM翼身组合体构型进行数值模拟研究.采用雷诺平均N-S方程,选择SA一方程湍流模型和SST两方程湍流模型,开展了网格收敛性分析及抖振分析计算研究,评估了该软件对此类民机构型的阻力预测能力.按照DPW5的要求,采用会议统一提供的基准结构网格进行计算分析,研究结果表明,计算的巡航设计点的阻力及大攻角下的气动力系数与试验值吻合较好,达到会议统计分析平均水平.计算结果表明,所采用的结构网格求解器WiseManPlus软件计算精度较高,适用于大型民机高速巡航及抖振现象的数值模拟研究.     

DLR-F4翼身组合体的阻力计算

为了考察自行研发的CFD软件的计算能力和阻力计算精度，本文采用LU—SGS方法、MUSCL差分格式和Baldwin—Lomax代数湍流模型，数值模拟了AIAA阻力计算工作室提供的DLR—F4翼身组合体的绕流流场，综合分析了easel和case2的气动力的计算结果，并与NASA Christopher L．Rumsey采用CFL3D6．0和AFRL/VAAC Don W．Kinsey采用Cobalt60提供的两组计算结果以及AGARD提供的两种不同风洞的测力试验结果作了比较。计算结果表明，本文计算精度与国外CFD软件相当。为了提高激波，边界层干扰的模拟精度，今后要重点加强湍流模型的应用研究。     

基于隐式嵌套重叠网格技术的阻力预测

 采用一种多层多块隐式嵌套重叠网格技术,对美国国家航空航天局通用化研究模型(NASA-CRM)翼身平尾(WBT)组合体进行了数值模拟与分析。多层多块隐式嵌套重叠网格技术是结合多层多块嵌套重叠网格处理策略和隐式切割方法,在建立重叠网格之间的流场信息传递关系时,基于网格单元切割准则选择"最优"重叠单元而无需人工设定插值边界。对美国AIAA委员会召开的第4届阻力预测研讨会(DPW-4)提供的CRM WBT组合体生成4种不同密度的结构化多层多块嵌套重叠网格,并采用计算流体力学(CFD)方法进行数值计算和阻力预测,计算结果与CFL3D和OVERFLOW的结果进行了对比。数值模拟结果表明:计算得到的压力分布和极曲线与CFL3D和OVERFLOW的结果几乎相同,说明了隐式嵌套重叠网格技术的有效性,同时也验证了流场求解方法与程序的可靠性。当迎角增大到3°左右时,在机身与机翼、尾翼连接处出现明显的分离涡,影响CRM WBT组合体的气动特性。在阻力预测方面,增加网格密度能够提高阻力预测的精度。采用不同的湍流模型会导致升、阻力系数的计算结果存在一定的差异,因此,湍流模型的选择也是阻力预测需要考虑的因素。     

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

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

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

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

栅格绕流数值模拟研究

栅格长宽比是影响栅格翼气动特性的重要参数,为研究它对栅格翼气动特性的影响规律,本文通过求解NS方程数值模拟研究单栅格长宽比对栅格气动特性的影响,并且对栅格边框厚度和形状对阻力的影响,在亚、跨、超声速多个马赫数下进行了数值模拟.结果显示,在所选择的长宽比范围内,长宽比对阻力的影响不大,对法向力有明显的影响,选择适当的栅格边框外形和厚度,可以大大减小栅格的阻力.     

基于RANS方法的民用飞机网格分析

基于RANS方程分别对代表民机巡航的DLR—F6模型和增升的NASATRAPWING模型进行多块结构化网格生成技术研究,通过网格拓扑、分布、质量参数等多项研究,得到影响民机外流场粘性计算需要的网格生成因子,形成一套适用高精度阻力预测的民机CFD网格生成指南。同时,为验证网格设计方法的可靠性及准确性,分析了不同网格生成准则对于CFD计算模拟精准度的影响,为民机计算策略研究及气动力设计提供了参考依据,具有重要的工程实用价值。