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

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

TRIP2．0软件湍流模块的开发与确认

基于TRIP2.0_SOLVER数值模拟软件,开发了湍流模型方法库,目前该CFD软件中集成的湍流模型主要包括B-L代数湍流模型、SA一方程模型和SST两方程模型.为了考核不同的湍流模型工程适用性,本文采用对接网格技术,通过有限体积法数值离散三维任意坐标系下的RANS方程组,应用LU分解、MUSCL差分格式和低雷诺数SA和SST两种湍流模型,数值模拟了二维NLR-7301两段翼型和三维DLR-F6翼身组合体的绕流流场,计算与试验比较的内容包括了表面的压力系数分布、典型站位的速度型和气动特性曲线等内容.通过计算与试验的比较,在本文的计算范围内,采用两种湍流模型均可以得到与试验结果相吻合的压力分布和升力曲线,但在边界层内的速度型、粘性阻力和力矩特性等方面,不同的湍流模型具有明显的差异.     

预旋进气的盘腔系统内流场的数值研究

用CFD软件(FLUENT6.0.12)对预旋进气的盘腔进行了数值模拟，计算中采用了多种湍流模型，包括Spalart—Allmaras单方程、K—epsilon和K-omega多种双方程模型，结果表明Spalart-Allmaras模型与结果符合得较好；通过计算研究了流场结构，流场主要由两个参数旋转比βp和湍流参数λT控制，λT主要决定源区的大小，计算结果验证了自由涡流的存在，并且盘腔压力分布的计算值与理论曲线符合得很好。     

不同湍流模型对增升装置气动特性计算结果的影响

为了研究湍流模型对增升装置气动特性计算的影响,分别采用了4种湍流模型对某民用飞机带增升装置翼身组合体构型进行了计算流体力学(CFD)数值模拟,并且将计算结果与风洞实验结果进行比较和分析。计算结果表明,不同湍流模型计算结果有所差异。在线性段,SST、Omega雷诺应力、BSL雷诺应力3种湍流模型计算结果较为接近,与实验结果符合的更好,计算阻力与实验值相比误差在6%以内,而Omega雷诺应力模型计算结果与实验结果相比差异较大。此外,这4种湍流模型对失速迎角、最大升力系数的预测也存在不小的差异。     

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

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

基于自主软件LSWS的典型地形大气CFD模式研究与应用

基于我国典型地形风资源评估、风电场微观选址及风电预测等的发展要求，根据风场实测数据，提出了两种风廓线模型：一种是利用实测数据修正的Gryning风廓线模型，另一种是结合分段函数及理论公式拟合的风廓线模型。考虑大气热效应、地转效应等大气效应，结合改进的雷诺时均k-ε两方程湍流模型，建立了能够模拟不同大气稳定度的典型地形大气CFD模式，开发了自主可控的风场CFD模拟软件LSWS(large scale wind simulation)。以锡林浩特平坦草原地形和山西丘陵地形的风场为例，比较了自主软件LSWS和商业软件Fluent的模拟精度。模拟结果显示：对于锡林浩特平坦草原地形风场，两个软件模拟结果与实测数据总体吻合良好；在稳定大气状态下，LSWS的计算精度较Fluent最大提升6.33%；在模拟山西丘陵地形时，LSWS的计算精度明显高于Fluent，特别是绕山脊的山后流场模拟更加合理，LSWS计算精度较Fluent最多提升40%以上。     

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

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

基于数值模拟的复杂地形风场风资源评估方法

与采用线性模型的WAsP软件相比,CFD方法具有可以逼真模拟复杂地形三维风场大气流动的优势,因此成为了复杂地形风场风资源评估的发展方向。为了把CFD软件包FINE/TURBO和带有壁面函数的k-ε湍流模型用于复杂地形地貌风电场的风资源评估和微观选址中去,本文以南澳岛风场作为研究对象,在中性大气条件假设下,对风场区域以30°风向为间隔进行数值模拟。研究采用不同入口边界条件对数值计算结果的影响。比对模拟结果发现,风场中的风加速与湍流强度基本上不随来流条件发生变化。基于这一规律,结合全年测风数据评估风场的风资源分布状况,建立了一套风资源评估的数值模拟及后处理方法,为工程实际提供了有力参考。     

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

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

某四级风扇特性计算的三维数值模拟

利用自主开发的程序对某四级风扇进行了特性计算,一方面是为了发动机整机的匹配,另一方面是为了校核程序对多级压气机的数值模拟精度.程序采用标准k-ε湍流模型计算湍流黏性系数,为了避免计算过程中由湍流模型方程带来的数值刚性问题,对湍流模型方程中的源项作了隐式处理.转静子之间的数据传递采用了掺混界面法.通过对一个单级风扇进行计算,并与常用的几种商用软件计算结果以及实验结果进行对比,初步验证了程序计算结果的可靠性.进一步对某四级风扇进行了计算分析,发现程序计算精度较好,可用于工程设计中.      

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

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

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

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

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

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

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

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

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

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

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

Recent experience with different methods of drag prediction

During the recent past the process of engineering design has changed entirely because of rapid developments in computational simulations. In aerodynamic design, computational fluid dynamics (CFD) methods are slowly superseding empirical methods and design engineers are spending more and more time applying CFD tools to analyze and predict the aerodynamic characteristics of vehicles. The purpose of this paper is to review recent experiences in CFD-based drag prediction with an emphasis on flow solutions governed by the Euler and Reynolds-averaged Navier–Stokes equations. For these types of flow solutions, various drag analysis methodologies are outlined and applied to determine the drag of components of as well as whole-body airplanes, helicopters, and ground-based vehicles at subsonic and transonic flow conditions. The review demonstrates that although significant progress has been made, CFD-based drag prediction still faces a number of hurdles that must be dealt with before it will become more widely accepted.