转捩模型在翼型气动性能预测 中的应用研究

边界层转捩位置的准确预测对于提高飞行器气动性能的预测精度具有重要意义.选取与k-ω SST湍流模型相耦合的γ-Reθt模型,以零压力梯度平板为研究对象,通过求解基于有限体积法的雷诺平均N-S方程验证该模型自动捕捉流动转捩的准确性;将该模型应用于传统有压力梯度的NACA 0012翼型的流场特性和气动性能的研究中,并与原始k-ω SST模型的计算结果及全湍流试验数据进行比较.结果表明:远场边界距离对翼型阻力系数有较大的影响;与无转捩模型相比,γ-Reθt转捩模型对翼型阻力系数的预测精度有一定程度的提高;对于二维模型,γ-Reθt转捩模型难以捕捉翼型表面的三维效应和非定常分离特性.     

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

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

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

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

固定转捩和自由转捩流动预测方法验证

采用固定转捩位置和基于SST k-ω湍流模型的γ-Reθt转捩模型自由转捩模拟两种方法研究了RAE2822翼型流动转捩问题。先采用γ-Reθt转捩模型确定转捩位置,再采用所获得的转捩位置实施固定转捩位置计算,并与全湍流计算结果和其他软件计算结果进行了对比。表明两种方法均能有效捕捉超临界翼型流动转捩现象,固定转捩位置模拟结果基本上重现了自由转捩预测结果。     

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

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

Gamma-Theta经验转捩模型在DES中的应用

与RANS湍流模型相比,DES湍流模型的边界层外自由流部分存在LES区,该区域的涡粘性会有所降低。为了研究了这一特征对γ-Reθt转捩模型的影响,以平板实验为算例,选取基于k-w SST湍流模型的DES模型和DDES模型,检验γ-Reθt转捩模型在DES模型和DDES模型中的表现。结果表明：DES模型和DDES模型使转捩位置推迟,SST—DES模型由于对网格过度敏感而没有触发转捩,SST—DDES模型使转捩位置后移较少;DDES模型比DES模型更适合γ-Reθt转捩模型;转捩经验关系需要进行修正后才能将γ-Reθt应用于DDES模型。     

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

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

用转捩模型预测转捩及确定最佳粗糙带高度

采用SST k-ω二方程湍流模型,附加γ-Reθ转捩模型对不同迎角下的NACA0012翼型绕流进行了转捩预测,计算预测的转捩位置与实验结果符合较好,说明γ-Reθ转捩模型对转捩位置具有较好的预测能力。进一步用模型对带有粗糙带的翼型进行转捩预测,获得了不同雷诺数下刚好激发转捩又不引起过大局部压力扰动和附加阻力的最低粗糙带高度,对固定转捩风洞实验有指导意义。     

应用γ-Reθ湍流模型和PSE方法研究高超声速流动

采用考虑转捩的γ-Reθ湍流模型研究了高超声速复杂流动,结合PSE(parabolic stability equa-tion)稳定性分析方法和e-N方法,通过给出边界层最不稳定点的参数,用以控制γ-Reθ的转捩经验关联式,减少对经验参数的依赖.计算了锥体高超声速流场,得到的壁面换热量和边界层转捩位置与实验吻合较好.计算了平板激波/边界层干扰和高超声速拐角压缩流场,得到了准确的波系结构、压力分布和分离区的大小.结果表明,该方法与传统湍流模型相比具有较大的优越性.     

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

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

基于γ-Reθt转捩模型的翼型气动特性计算精度改进研究

将结合SSTk-ω湍流模型的γ-Reθt转捩模型与二维雷诺平均NS（RANS）方程数值求解相耦合,发展了考虑转捩影响的翼型气动特性计算方法。得到的计算结果与没有考虑转捩影响的计算结果的比较表明,与实验更加吻合,计算精度得到了明显提高。     

低压涡轮叶栅流动中转捩模型的校验及改进

为评估并提高现有转捩模型的预测精度,采用计算流体力学(CFD)软件FLUENT 12.1,选取层流模型、全湍流模型、剪切应力输运(SST)低雷诺数模型、k-kl-ω模型以及γ-Reθ模型对低压涡轮叶栅T106-EIZ进行数值模拟,通过与实验数据的对比校验了后3种模型对于转捩以及相关参数的模拟能力,并对结果以及模型的作用机理进行分析,校验结果表明所有模型都不能准确地预测分离流转捩以及尾迹诱导转捩.选取预测效果较好的γ-Reθ模型进行了修正,提出通过修改间歇因子输运方程中的参数Ca1和Ca2的方法来修正该模型,结果表明该方法可以提高模拟精度.     

γ-Re_θ模式应用于高速边界层转捩的研究

应用γ-Reθ转捩模式对超声速、高超声速边界层转捩进行数值模拟。γ-Reθ模式通过求解关于当地雷诺数和间歇因子两个输运方程给出转捩起始位置和转捩区长度等信息。本文对马赫数3.5至7范围内的四种算例进行计算,研究了来流雷诺数、攻角变化和头部钝化半径等关键参数的变化对γ-Reθ模式预测流动转捩性能的影响。给出了壁面摩擦阻力系数、热流值与实验值的对比以及壁面附近间歇因子等值线分布等计算结果。γ-Reθ模式能够正确预测攻角和钝化半径变化时转捩位置和流动参数的变化趋势;在较大雷诺数时计算结果与实验值吻合很好。γ-Reθ模式对于超声速或高超声速边界层转捩的模拟仍需修正和改进。     

基于Gama-Theta转捩模型的气冷涡轮气热耦合研究

基于非结构化混合网格,采用有限体积法求解了N-S方程,对流项格式为AUSM+,梯度的求解方法为加权最小二乘法,求解器采用LUSGS隐式求解方法,并进行了预处理和网格重排序。与导热方程耦合,搭建了气热耦合平台,采用直接耦合方法,交界面互相传递温度,并且采用能保证通量守恒的面积加权类的插值方式实现数据传递,将γ-Reθt两方程转捩模型应用在流场程序中。通过MARKII叶片4311和5411工况的实验结果进行了验证,结果表明采用转捩模型对压力分布的影响不大,与实验值吻合较好,通过间歇因子的分布可以看出,γ-Reθt模型成功地预测了从层流到湍流的转捩过程,在层流区采用γ-Reθt模型计算的涡粘系数与真实流动情况更加吻合,由于涡粘系数对传热的影响很大,从而有效地提高了边界层的层流区和转捩区的传热计算的精度,得到的温度和换热系数与实验值吻合更好,但是该模型是和SST模型耦合在一起的,由于SST模型的局限性,所以对激波和边界层干扰区域的模拟产生了7%左右的温度误差。     

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

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

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

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

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

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

压缩性修正对γ-Re_θ转捩模型的影响研究

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

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

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

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

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