基于SA和SSTγ-■模型对边界层转捩预测的比较

利用作者共同研发的In-house代码TRANS3D平台,在SA和SST两种常用湍流模型框架下构建了γ-■转捩模型,并以二维低速S809层流翼型和三维小后掠跨音速F5层流机翼为对象,比较了两种不同湍流模型构架下γ-■转捩模型预测的气动力特性和流场分布。结果表明:两种转捩模型均能预测航空转捩计算中常见的分离流转捩和自然转捩类型,明显改善了中低雷诺数流场下的预测精度,但由于选取基准湍流模型的不同,基于SA和SST的γ-■转捩模型在流动细节上依然存在着一些差异。     

γ-Reθ转捩模型在二维低速问题中的应用

介绍了基于当地关联的γ-Reθ转捩模型.通过Spalart建议的环境源项法对自由流湍流度衰减的有效控制,实现了该转捩模型的实用化.利用中国空气动力研究与发展中心的TRIP软件平台上实现的该转捩模型,对S809翼型、NLR7301两段翼型和30P30N三段翼型3个低速翼型绕流问题进行了计算,将不同迎角下的计算结果与试验结...     

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

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

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

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

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

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

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

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

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

采用γ-Re_(θt)模型的转捩流动计算分析

为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Reθt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序(HUNS3D)。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Reθt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明:γ-Reθt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

基于湍流模型的转捩预测方法研究

阻力的精确计算,一直是低雷诺数飞行器设计的关键.而与附面层分离和转捩相关的摩擦阻力的精确计算则是此类飞行器阻力精确计算的关键.因此,开展了基于湍流模型的转捩预测方法的研究,应用几种湍流模型对平板边界层进行计算分析,研究各湍流模型预测流动转捩的能力.由于各模型对扰动的过于敏感,引起转捩提前发生,使预测与实验差距较大,由于k-ωSST两方程湍流模型预测能力优于其它模型,在此模型基础上引入了间隙函数进行粘性系数修正,对平板和NACA 0012翼型的转捩重新进行了预测,计算结果表明改进后的模型对转捩位置具有较好预测能力.     

三维不可压缩层流、转捩、湍流边界层的数值计算

本文提出了一种计算三维不可压缩边界层的数值解法,其中包括对层流边界层的微分法,对湍流边界层的积分法以及利用Orr-Sommerfeld方程计算转捩的数值解,当扰动振幅达到e~(13.5)时发生转捩。 按照上述方法,本文计算了ONERA M6后掠角为30°的机翼,计算结果分别和Cebeci和Cousteix的结果作了对比,二者吻合较好。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

SST湍流模型改进研究综述

k-ωSST湍流模型是目前综合性能最佳的涡黏模型之一,近年来得到非常广泛的应用。然而,随着问题复杂性的增加和模拟精准度要求的提高,标准SST湍流模型在某些方面呈现出明显的局限性,大量学者对其进行了相应的改进研究。围绕旋转/曲率效应、可压缩效应、激波不稳定性效应、雷诺应力各向异性效应、应力-应变偏差效应和层流/湍流转捩效应等6个方面,对SST模型的改进研究进行了评述,同时对近年的基于数据驱动技术的模型改进也进行了简要介绍;梳理了各种改进研究的思路和发展趋势,阐述了各自的适用性和局限性,并分析了影响改进效果的原因和问题所在。最后,对未来工作提出了建议。     

[γ-Reθt]转捩模型在跨声速涡轮叶栅中的应用

为了解涡轮叶栅在跨声速条件下的流动特性和准确预测涡轮叶栅外换热情况,对γ-Reθt转捩模型,依据零压力梯度时自由流湍流度衰减实验结果,给出了一种直接估算来流粘性比的方法,以保证叶片前缘附近具有正确的自由流湍流度分布,提高换热预测准确度;同时减少试算次数。对MARK II与VKI两种叶栅通道跨声速工况下的流动换热情况使用CFX软件,选取层流模型、SST k-ω模型以及缺省粘性比和设定合理粘性比的γ-Reθt转捩模型进行了数值模拟验证,计算结果与实验数据的对比表明:转捩模型优于其他模型;而采用本文方法给定进口粘性比,能准确预测转捩位置,同时显著改善γ-Reθt转捩模型对不同来流湍流度下涡轮叶栅表面换热的预测精度;当入口湍流度较高,相比采用缺省粘性比情况,压力面上换热系数的相对误差降低30%以上,控制在7%左右。     

基于S-A湍流模型和间歇因子输运方程的转捩流数值模拟

考虑到γ-Re_θt转捩模型中间隙因子输运方程并不依赖于具体的湍流模型,因而提出了耦合S-A湍流模型和间歇因子输运方程的转捩流模拟方法.其主要思想是构造与间歇因子相关的两个耦合函数,分别作用到S-A模型的生成项和耗散项用以控制湍流的产生与发展,从而实现层流到湍流的数值转捩.结果表明:该方法中模型方程具有完全当地性的特点,易于实施,可直接用以求解三维转捩流动.4个典型算例表明该方法模拟结果与γ-Re_θt模型准确度相当,但少求解了两个输运方程,计算耗时减少了15%.      

γ-Reθt转捩模型在微通道流动特性研究中的应用

应用ANSYS CFX对矩形微通道中的流动进行了数值研究。通过与实验数据的对比检验了层流模型、γ-Reθt转捩模型和剪切应力输运模型在微通道流动转捩模拟中的适用性,并在此基础上探究了微通道长径比以及截面宽高比对临界雷诺数的影响。结果表明:3种模型中,只有γ-Reθt转捩模型捕捉到了转捩的趋势,并较为准确地捕捉到了临界雷诺数。当微通道的长径比不小于200,进口段对整体流动特性产生的影响可以忽略。当微通道的当量直径和长径比一定,宽高比由1增加到5时,临界雷诺数呈逐渐增长趋势。      

基于改进SST模型的分离流动数值模拟

对SST湍流模型进行了改进,并通过对典型分离流动进行数值模拟,来检验和提高模型预测分离流动的能力.基于亚声速分离流动,提出了提高雷诺应力的模拟精度和分离流修正的改进方法,并进行了对比研究得出结论;在亚声速分离流动分析结论基础上,采用了可压缩性修正方法,开展了跨声速、超声速和高超声速激波/边界层干扰分离流动的数值模拟研究.结果表明:提高雷诺应力的模拟精度和采用分离流修正明显地提高了SST湍流模型对分离流动的模拟能力;适当地可压缩性修正对超声速和高超声速分离流动的计算精度有改善作用.     

转捩对压缩拐角激波/边界层干扰分离泡的影响

为了研究转捩对压缩拐角内分离泡结构的影响,进行了来流马赫数2.9,24°压缩拐角激波/转捩边界层干扰的直接数值模拟（DNS）。通过在拐角上游平板的不同流向位置处添加周期性吹吸扰动激发流动转捩,使得转捩不同阶段进入拐角入口,从而在拐角内产生激波/转捩边界层的相互干扰。计算得到的平均速度剖面、壁面压力分布以及分离泡大小与风洞试验及以往直接数值模拟结果吻合较好,验证了计算结果的可靠性。研究了转捩过程对角部干扰区内分离泡结构的影响规律,分析比较了不同转捩阶段下角部分离区内湍动能的生成、耗散和分配机制。研究结果表明：转捩初期的拟序涡结构对分离泡尺度及形状影响最大,发卡涡包在角部拐点附近发生展向融合,并在角部区域形成湍流斑,此时分离泡尺度最小,形状呈现中间高两边低的山峰型。随着转捩的发展,分离区内湍动能生成和近壁区的耗散逐步降低,此时输运项起到了主要的平衡作用。