基于IDDES方法和γ-Re_θ转捩模型的粗糙颗粒诱导高速边界层强制转捩模拟

采用γ-Re_θ转捩模型与IDDES (Improved Delayed Detached Eddy Simulation)相结合的方法对BAM6QT (Boeing/AFOSR Mach-6 Quiet Tunnel)中的马赫6来流条件下粗糙颗粒诱导转捩情况进行了数值模拟研究,通过与试验测量的压力脉动均方根值、脉动主频和边界层内的皮托压分布的定量对比及与文献中DNS (Direct Numerical Simulation)流场结构的定性对比,表明该方法可以捕捉粗糙单元诱导出的流向涡,能够模拟颗粒前缘分离激波和弓形激波之间的震荡现象,能够模拟流向涡向下游的发展失稳过程及其脉动发展过程。计算结果表明流向涡结构在粗糙颗粒下游40倍颗粒直径位置开始破碎,非定常扰动能也在该点附近增长达到最大值。粗糙单元诱导出了明显的条带涡结构,其中低速条带处于边界层上层,在向下游发展的过程中逐渐扩散至边界层外缘并耗散掉。高速条带位于边界层底部,在向下游发展的同时往展向两侧拓展,最终展向上多条高速条带接触并互相耦合,导致最后条带涡结构的破碎和尾迹区边界层的转捩。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

C-γ-Reθ高超声速三维边界层转捩预测模型

针对国家数值风洞（NNW）工程高超声速三维边界层转捩预测需求,开展了高超声速边界层横流转捩判据及模型研究。采用线性稳定性分析e^N方法对高超声速转捩数据进行扩展,结合横流强度与表面粗糙度构造当地化的高超声速横流转捩判据,对基于Chant 2.0计算平台的高超声速修正γ-Re_θ转捩模型进行了横流模式拓展,建立了适用于高超声速三维边界层横流转捩预测的C-γ-Re_θ转捩模型。采用构建的转捩模型对多状态下的高超声速尖锥进行横流转捩预测,取得了与试验结果符合较好的预测效果。     

γ-Re_θ转捩模型在高超声速下的应用及分析

为了初步研究γ-Reθ转捩模型在高超声速领域的适用性,在自行开发的CFD程序中添加了该模型,利用T3系列平板算例对该模型的实现进行了验证;针对高超声速流动,选取压缩面绕流和双圆锥扰流等算例,分别利用该转捩模型、全层流模拟、以及传统的湍流模型进行了数值模拟计算,并与实验结果进行对比。计算结果表明,尽管该转捩模型是在低速流动的基础上发展的,它仍能准确预测高超声速流动下的转捩现象。     

基于新的经验关联公式的γ-Re_θ转捩模型在高超声速流动中的应用

为提高γ-Reθ转捩模型在高超声速流动下的转捩预测精度,参考部分高超声速风洞的实验数据对该模型进行改进,在原始经验关联公式的基础上提出了转捩动量厚度雷诺数基于马赫数、湍流度、压力梯度因子的经验关联公式;选取高超声速平板绕流、高超声速圆锥绕流等算例对提出的经验关联公式进行了验证,并将原始经验关联公式和提出的经验关联公式计算得到的重要参数进行对比,分析了提出的公式的作用机理。结果表明:基于低速不可压流动经验关联公式的γ-Reθ转捩模型在高超声速流动中计算得到的流场中Reθt分布偏低(一般只有真实值的40%~60%),进而导致边界层内Reθt偏低而过早达到转捩触发准则,与实验值相比所预测的位置可能提前10%(以参考长度为基准)甚至更多。新的经验关联公式有效地改善了这一问题,计算结果与实验值符合得更好。     

基于IDDES框架的γ-Reθ转捩模型

针对高超声速流动的γ-Re_θ转捩模型在模拟强制转捩时存在捕获边界层内扰动不足的缺点,将RANS（Reynolds Averaged Navier-Stokes）框架改为IDDES（Improved Delayed Detached Eddy Simulation）框架,既能像基于RANS框架的转捩模型一样模拟复杂构型的自然转捩,又能发挥IDDES能够捕获更多脉动信息的优点,较为准确地模拟粗糙颗粒诱导强制转捩。通过对一系列简单构型的自然转捩及来流马赫数为6条件下平板上单个粗糙颗粒诱导强制转捩的模拟表明,模型既能体现γ-Re_θ转捩模型的优点,在自然转捩模拟中具有较强的鲁棒性,能够反映雷诺数等因素对转捩位置的影响规律;也能体现IDDES方法的优点,能够捕捉粗糙颗粒诱导的扰动及涡结构,从而较为准确地刻画出强制转捩的整个流程。     

压缩性修正对γ-Re_θ转捩模型的影响研究北大核心CSCD

将基于流场当地变量的γ-Reθ转捩模型加入可压缩RANS方程计算程序,对平均流场控制方程、湍流和转捩控制方程进行基于LU-SGS的隐式紧耦合求解。结合三种压缩性修正方法对高超声速平板、双楔和圆锥流动进行了数值模拟,给出了壁面斯坦顿数、热流值和实验值的比较,以及湍动能和间歇因子分布云图。数值计算结果表明,相同来流湍流度下不同压缩性修正方法得到的转捩位置差别较大。相比于未修正的模型,基于当地马赫数压缩性修正后的转捩位置最靠后,其对带有分离的双楔流动预测较为准确;基于来流马赫数的压缩性修正仅使得转捩位置稍有延迟;对湍动能源项的压缩性修正可提高模型在转捩后湍流段的热流预测精度。     

γ-Re_θ转捩模型的标定研究

介绍了γ-Reθ转捩模型,补充了相应的关联函数。基于平板的试验数据在中国空气动力研究与发展中心(CARDC)的trip软件平台上对该关联函数进行了标定,并研究了关联函数的变化对转捩结果的影响规律。通过修正与SST湍流模型的结合,改进了转捩位置后壁面摩擦阻力系数偏低的状况。     

高超声速边界层的转捩问题

取40km高空处的气体参数,对来流马赫数为8、10和12等三种不同情况的平板边界层和楔角为20°的楔型体边界层,做了转捩发生位置的预测。所用方法为改进的eN方法。结果发现,壁面条件对是第一还是第二模态波决定转捩位置有很大关系。如对于传统的eN方法,如果取N为10作为转捩判断标准,对马赫数为12的楔体边界层,等温壁条件下由第二模态波决定的转捩位置距前缘13m左右,而平板边界层中两模态波所决定的转捩位置距前缘均超过50m。而采用改进的eN方法,平板边界层中两模态波所决定的转捩位置距前缘均超过10m,而对马赫数为12的楔体边界层,两模态波所决定的转捩位置距前缘均出现小于5m的情况。     

改进k-ω-γ转捩模式对不同雷诺数下H IA D的转捩预测

高超声速充气式柔性减速器（HIAD）在气动力作用下会变形为波纹状,从而促进流动转捩为湍流,准确预测其转捩位置和壁面热流对热防护系统的设计至关重要。拓展了分离诱导转捩预测性能的改进k-ω-γ模式,同时具备对第1模态、第2模态、横流模态以及流动分离失稳的预测能力。本文将其应用于不同雷诺数下壁面波纹变形的HIAD边界层转捩预测,并与原始k-ω-γ模式的预测结果进行了对比,以评估和验证其对复杂转捩现象的预测性能。在此基础上,细致剖析了改进k-ω-γ模式的转捩预测机制。结果表明,改进k-ω-γ转捩模式可准确预测不同来流雷诺数下HIAD的转捩起始位置、转捩阵面形态和壁面热流分布。波纹壁面波峰处的转捩预测主要由构造的分离间歇因子猝发。而在波谷位置,第1模态、横流模态以及流动分离的贡献都很重要。以上研究显示了改进k-ω-γ模式在复杂外形中的应用潜力,可为多重不稳定耦合作用下的转捩预测方法发展提供参考。     

高超声速圆锥边界层转捩数值研究

以类似于处理湍流脉动的方式来模型化非湍流脉动,考虑湍流的间歇性特征及引起超声速流动失稳机制,研究包括转捩起始位置及转捩区范围的高超声速流动转捩问题,发展可反映扰动模态和可压缩性影响的两方程SST转捩/湍流模式.针对尖锥模型,在马赫数6、攻角O°～2°时,采用第二模态(Second-mode)机制,预测并比较完全层流态、完全湍流态,以及自然转捩发展过程中的表面热流、表面摩阻等气动特性.     

γ-Re_θ转捩模型的标定与应用

简单介绍了γ-Reθ转捩模型,并在mbns3d内部软件上实现。应用零压力梯度平板试验数据对关联函数进行了标定,得到了一组关联函数。对A-aerofoil迎角13.1°及13.3°进行了验证计算,结果表明该模型可以比较准确模拟转捩效应,提高了升力、阻力计算精度。随后研究了NLR7301两段翼型,通过应用环量修正及转捩模型提高了阻力计算精度,阻力最大误差控制到4%以内。     

高超声速飞行器前体边界层强制转捩数值模拟

边界层强制转捩是保证超燃冲压发动机进气道正常启动的关键技术之一,k-ω-γ转捩模式是适用于高超声速边界层转捩预测的方法。为研究该方法对边界层强制转捩的预测性能及不同转捩带对边界层强制转捩的影响特征,对原始的k-ω-γ转捩模式进行了壁面温度影响修正,采用修正的k-ω-γ转捩模式对高超声速飞行器进气道前体边界层强制转捩进行数值分析,计算了光滑外形、钻石型转捩带外形和斜坡型转捩带外形在马赫数Ma=6,7条件下的边界层转捩,并与试验结果进行了对比。研究结果表明,修正的k-ω-γ转捩模式对边界层强制转捩具有较好的预测能力,计算得到的转捩起始位置与试验结果基本吻合。两种转捩带强制转捩效果明显,其中:钻石型转捩带产生的扰动强于斜坡型转捩带,且转捩区长度较斜坡型转捩带短;斜坡型转捩带在控制边界层流动分离、减小流动横向溢出效果上优于钻石型转捩带。     

基于改进的k-ω-γ转捩模式预测高超声速飞行器气动特性

为验证边界层转捩对高超声速飞行器气动特性的影响,采用改进的k-ω-γ转捩模式对类X-51A高超声速飞行器进行了全机边界层转捩预测.不仅系统分析了飞行攻角和雷诺数对边界层转捩的影响规律,同时研究了边界层转捩对飞行器气动力和进气道性能的影响.发现边界层转捩对飞行器升力系数和俯仰力矩系数影响较小,对阻力系数影响较大.针对本文...     

从总体设计角度透视高超声速飞行器边界层转捩问题

边界层转捩对飞行器总体性能影响很大,按照弹道式再入、再入机动、高超声速滑翔、高超声速巡航四类飞行器进行分类,与边界层转捩相关的典型问题可归结为气动力、气动热及掺混效率三类,相应带来飞行稳定性与配平能力、落点散布、防热风险、飞行器减重、推进系统优化等问题.针对边界层转捩控制问题,从总体设计方面可采取弹道设计、翼载荷控制、...     

高超声速三维边界层横流转捩的数值研究

横流失稳是高超声速三维边界层转捩的主要机制。然而到目前为止,由横流失稳导致高超声速三维边界层转捩的数值研究还不多见。本文采用直接数值模拟方法计算了马赫数6、后掠角45°钝板边界层中定常横流涡的演化,在此基础上引入高频二次失稳模态,计算了二次失稳模态的非线性演化,直至湍流发生。计算结果表明:横流定常涡的非线性作用引起平均流修正,可使壁面摩擦系数曲线有一定程度的抬升。而高频二次失稳波的增长导致低频及定常扰动快速增长,促使壁面摩擦系数急剧抬升,同时饱和横流涡结构破碎,最终触发转捩发生。     

再入飞行器迎风面边界层转捩的数值预估

基于求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,首先数值预估了高超声速三维流动的边界层特性;采用已有的高超声速三维边界层转捩准则,即在转捩位置Ｒｅθ／Ｍｅ＝常数,以及通过修正由ＯＮＥＲＡ提出的代数转捩模型,建立了预估三维边界层转捩的半经验数值技术并应用于三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的求解。利用ＨＡＬＩＳ模型及相关的实验数据,首次进行了在Ｍ∞＝６时不同雷诺数及迎角的验证试验。预估结果表明：在不考虑