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

采用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软件的确认：DPWⅡ复杂组合体的数值模拟

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

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

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

大迎角翼—身组合体涡流特性计算

论述了用非线性离散涡法来模拟大迎角下翼－身组合体涡流绕流的计算模型及计算方法；给出了一个典型翼－身组合体的涡流流态计算结果及非线性气动特性和截面压强分布。     

湍流模型在民机跨声速绕流计算的应用

针对民机构型DLR-F6翼身组合体,研究了不同湍流模型对气动特性和流场特征的影响.控制方程采用全场三维可压、雷诺平均的Navier-Stokes方程;计算格式采用Jameson中心有限体积、多步Runge-Kutta时间步长格式进行求解;湍流模型分别采用EARSMk-ω模型、Spalart-Allmaras(S-A)模型、SST模型、Standardk-ω(STDk-ω)模型、Low-Reynolds Numberk-ω(LRNk-ω)模型.通过对巡航状态的计算研究,分析了五种湍流模型在翼身结合处气流分离的特点,为工程应用提供了一定的借鉴意义.     

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

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

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

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

半模试验翼身组合体垫块的数值模拟分析

采用数值模拟方法研究存在风洞侧壁边界层影响时垫块厚度和修型对翼身组合体气动特性的影响.本文采用Jameson中心有限体积、多步Runge-Kutta时间推进格式求解三维可压、雷诺平均Navier-Stokes方程,引入Baldwin-Lomax湍流模型,分析了垫块厚度和修型对DLR-F4翼身组合体半模型气动力和局部流场的影响.通过本文计算得出的结论对于跨声速半模型风洞试验垫块优化设计具有重要的指导意义,能更好地提高风洞试验数据的精准度和试验技术水平.     

Pω增强型k-ω湍流模型在三角翼旋涡流动的应用

在三角翼旋涡绕流数值模拟中,标准 Wilcox k-ω湍流模型生成项未考虑旋度的影响而导致预测的旋涡强度较弱。通过引入探测因子区分剪切层和涡核,在旋涡流动的高旋度区域增加ω方程生成项的方法,基于结构化网格上的 RANS 求解器,加入了 Pω增强型 k-ω湍流模型,对绕尖前缘三角翼亚声速和跨声速旋涡流场进行了数值模拟。计算结果与 NASA 的 NTF 风洞和 DLR 的 DNW-TWG 风洞试验数据进行了对比分析,结果表明：不论在亚声速还是跨声速自由来流条件下,Pω增强型 k-ω湍流模型计算的压力分布、涡破裂位置均与试验数据吻合良好,准确地预测出了三角翼上翼面的主涡、二次涡结构,特别是跨声速条件下激波干扰导致的涡破裂的临界迎角及涡破裂位置,表明 Pω增强型 k-ω湍流模型在绕三角翼旋涡流动数值模拟中具有良好的适用性。     

梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。     

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

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

飞行器大迎角非对称涡与气动特性的实验研究

针对飞行器大迎角绕流的非对称气动特性,采用变湍流度技术,通过对细长体和翼身组合体模型的实验研究,分析了湍流度和人工扰动对飞行器大迎角流动非对称性的影响规律。实验结果表明,湍流度和人工扰动对飞行器大迎角情况下的非对称流动特性影响显著。     

包含支撑装置和机翼变形的CRM-WB构型气动特性数值模拟

CRM(Common Research Model)是第四届~第六届AIAA阻力预测研讨会(DPW)选择的参考外形。第四届和第五届DPW会议气动特性数值模拟结果的统计分析表明试验与计算之间存在明显差异。采用CFD方法和流固耦合(FSC)方法数值模拟了带支撑装置的CRM翼身组合体(CRM-WBS)构型的气动特性,以评估模型支撑装置和静气动弹性变形对CRM翼身组合体(CRM-WB)构型气动特性数值模拟结果的影响。通过与不包含支撑装置的CRM-WB构型的数值模拟结果和NASA Langley NTF风洞试验结果的对比分析,CRM-WBS构型的CFD计算结果表明,支撑装置导致机翼上翼面激波位置前移,升力系数、阻力系数、低头力矩系数下降。CRM-WBS构型的流固耦合计算结果表明,静气动弹性变形主要影响机翼上表面激波位置并显著降低外侧机翼激波位置前的负压值,进一步导致升力系数、阻力系数、低头力矩系数下降。包含支撑装置和静气动弹性变形的CRM-WB构型气动特性数值模拟结果更加接近试验结果。     

典型航空分离流动的雷诺应力模型数值模拟

先进飞行器设计对CFD方法的边界层分离模拟能力提出了更高要求,传统雷诺平均纳维-斯托克斯(RANS)涡黏性模型因建模和构造层面的理论缺陷导致其分离流动预测可信度较低。雷诺应力模型由于未对雷诺应力及其分量关系进行建模,对湍流非平衡、旋转以及雷诺应力各向异性等流动现象具有天然的理论优势。为验证与确认雷诺应力模型对典型航空分离流动的预测能力,基于SSG/LRR-g模型,以NACA4412翼型大攻角分离、M6机翼跨声速分离以及F6翼身结合区分离流动为例,探讨了雷诺应力模型对逆压梯度、激波诱导分离、二次流动分离等典型航空分离流动预测的适应性。通过与k-ω剪应力输运(SST)模型模拟结果对比发现,雷诺应力模型对分离泡大小、速度型分布、雷诺应力分布和激波位置等关键特征的模拟精度较涡黏性模型显著提升,基本验证了雷诺应力模型可在翼身接合区角区流动和三维强激波诱导分离等问题中得到正确的流动特征,而SST模型在此类流动中基本失效,显示了雷诺应力模型在典型航空分离流动中较涡黏模型的优势。同时,发现k-ω SST模型所包含的Bradshaw假设在三维激波诱导分离较强时严重影响了模型预测的准确性,是预测结果偏离试...     

RAE2822翼型跨音速流动CFD计算的可信度分析

采用格点格式的有限体积法对RAE2822翼型的气动特性进行了粘流数值模拟,流动模型为雷诺平均NS（RANS）方程。在此基础上系统研究了不同网格类型以及不同湍流模型的计算精准度。通过与经典的试验数据进行确认研究,得出了适宜二维跨音速粘性流动计算的网格拓扑分布规律,并总结出了适宜的湍流模型。     

自由曲面修形技术在民机翼身组合体上的应用

翼身结合处的整流在民机气动中有很重要的作用,需要设计合适的外形以调整机翼上表面的压力分布及翼身结合处的气流,达到提高大型民用客机翼身组合体气动性能的目的.基于CATIA的自由曲面修形技术,采用自由曲线B样条的参数化方法,综合考虑工程约束(起落架布局),对民机构型翼身结合处进行了整流设计,并采用求解N-S方程的方法计算样...     

低速粘性分离流的数值模拟

本文介绍对一种采用涡控制技术的新概念兵器--环翼以及升力体航天飞机的低速气动力数值模拟.对这类外形复杂、流动状态也十分复杂的流场,需正确模拟出流场不同部分的流态及湍流运动参数.本文采用一种主要与流场平均运动涡量相联系的新的湍流模型,较好地解决了这一问题.数值计算给出了大量满足工程要求的气动参数,并解决了一些关键性气动问题.     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

前体涡诱导双极限环摇滚流动特性的实验研究

通过摇滚/PIV/压力同步测量实验,对翼身组合体前体涡诱导的双极限环摇滚过程中流动特性及演化规律进行了系统的研究,并分析了前体涡诱导翼-身组合体双极限环摇滚的流动机理。实验结果表明,前体涡与机翼翼面流动的相互作用使模型在正负滚转相位处分别出现极限环摇滚运动;正负滚转相位过渡是模型运动惯性与气动力共同作用的结果。