DPWⅢ机翼和翼身组合体构型数值模拟

采用亚跨超CFD软件平台(TRIP)数值模拟了DPW Ⅲ提供的DPW_W1/W2两种机翼构型和DPW-F6/F6_FX2B两种翼身组合体构型,主要目的是通过两种机翼构型和两种翼身组合体构型的数值模拟,研究网格密度对运输机构型气动特性计算结果的影响。数值模拟采用的多块对接网格来自AIAA CFD Drag Prediction Workshop Ⅲ(DPW Ⅲ),采用National Transonic Facility(NTF)的试验结果和CFL3D的计算结果作对比。详细研究了网格密度对两种机翼构型和翼身组合体的总体气动特性和压力分布的影响。采用SST两方程模型计算两种构型均得到了网格收敛结果,网格密度主要影响压差阻力而对摩擦阻力影响较小,计算结果较好地预测了机翼和翼身组合体外形优化前后总体气动特性的变化量。     

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

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

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翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

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

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

TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟

 采用“亚跨超声速计算流体力学软件平台”（TRIP2.0）数值模拟了阻力预测小组（AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ）翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果。重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致。采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响。     

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

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

运输机组合体绕流非结构混合网格N-S方程数值模拟

本文利用通用网格生成软件生成运输类飞机翼身组合体（DLR-F6）非结构混合网格，在边界附面层区域内，生成三棱柱形网格，在其他流动区域采用四面体网格进行填充。粘性流场求解器采用基于N-S方程的CFX5计算流体软件求解，通过计算结果与风洞试验所得到的数据对比分析，验证了网格生成方法及网格数据的正确性和实用性。本文研究结果也为运输类飞机高速风洞试验提供了对比数据。     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

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

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

结冰对翼型及翼身组合体气动特性影响研究

飞机结冰是影响飞机飞行安全的重要因素之一,研究结冰对飞机气动性能的影响具有重要意义。应用多块结构化网格生成方法,通过求解雷诺平均N-S方程分别对翼型及DLR-F4翼身组合体干净外形和三种不同带冰外形的流场进行计算,分析不同冰型对翼型及翼身组合体绕流流场及气动特性的影响。     

基于非线性气动力的飞机飞行载荷计算

为了研究气动力非线性对飞机静气动弹性特性及飞行载荷的影响,发展了一种基于非线性气动力的飞机静气弹及飞行载荷计算方法。通过引入CFD计算结果,实现对线性方法中的气动力影响系数矩阵的非线性修正,采用模态法求解静气弹配平方程得到飞机非线性的气弹变形及飞行载荷。以某翼身组合体构型为计算算例,分别针对刚性模型和柔性模型,基于非线性气动力完成了飞机静气动弹性分析及载荷计算,并与线性结果进行了对比分析。结果表明,非线性方法可对大迎角范围内飞机静气弹变形及飞行载荷做出较合理的预测。     

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

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

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

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

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

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

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

某运输机全机构型雷诺数效应阻力修正研究

采用CFD方法分析了雷诺数效应对某运输机全机构型阻力的影响。基于对某运输机的翼/身/架/舱组合体绕流进行的计算分析,归纳总结出了一种适应于复杂构型的低雷诺数到高雷诺数的阻力修正方法。通过验证可知,该修正方法修正效果良好,为工程中的阻力预测提供了一种简洁有效的方法。     

民机翼身组合体参数化研究

飞机方案设计阶段,迫切需要一个能够快速评估、计算不同参数组合的翼身组合体基本气动力特性的分析平台。为满足这一实际工程设计需要,本文提出运输类飞机参数化概念和具体研究方法,开发了快速分析计算翼身组合体气动力特性分析平台,结合具体算例验证了方法的可行性和工程实用价值。     

翼身融合布局与传统布局飞机结冰特性对比分析

翼身融合（BWB）布局是提升未来民机综合性能的重要布局方式，其结冰飞行安全问题逐渐受到人们重视。针对翼身融合布局飞机与常规布局飞机的结冰问题，本文通过数值模拟的方法来开展结冰特性研究。本文提出基于Navier-Stokes方程对空气流场进行求解，并应用欧拉法计算水滴撞击特性，之后采用Shallow-Water结冰热力学模型的结冰计算方法。首先，通过将翼身融合布局飞机和传统布局飞机的计算结果与风洞试验数据进行对比，验证空气流场计算的正确性，并将两者进行对比分析；其次，数值预测飞机表面冰形的特征，将两种布局飞机结冰特性进行对比，结果表明，结冰对两种布局飞机气动外形的破坏程度从后掠翼翼根至翼尖逐渐变大，但传统布局飞机结冰只发生在机翼前缘和机头处，而翼身融合布局飞机前部几乎都发生了结冰，可为相关的结冰特性研究及防除冰设计提供技术参考。     

考虑机翼尾流影响的运输类飞机后体气动外形优化设计（英文）

在某典型运输机翼身组合体的构型上,进行了考虑机翼尾流影响的机身后体气动外形优化设计研究。基于翼身组合体构型建立了考虑尾流影响和部分工程约束的优化设计系统,并对后体构型在巡航状态下进行了优化设计。以比较适合描述后体变形的NURBS样条基函数为空间控制体属性引入FFD自由变形技术,通过在FFD控制框架对该运输机后体进行了空间属性构建。采用无限差值动网格技术提高空间网格的更新效率并保证网格质量。利用改进Kriging代理模型、量子粒子群优化算法提高优化效率和全局寻优能力。为减小巡航阻力,对某运输机后体的上翘角、截面形状等参数进行了优化设计,优化结果显示,设计后的机身气动特性明显提高。     

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

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