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模型。     

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

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

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

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

γ-Re_(θt)-CF转捩模型在Spalart-Allmaras湍流模型中的推广及验证

横流(CF)不稳定性是三维流动中诱发转捩的一项非常重要的因素,考虑到γ-Reθt-CF转捩模型对流向Tollmien-Schlichting波和横流波不稳定性引起转捩的判定均是完全基于当地变量,且Spalart-Allmaras(SA)湍流模型计算效率高,因而将γ-Reθt-CF转捩模型与SA湍流模型相结合,并将其引入开源Standford University Unstructured(SU2)计算流体力学分析平台。为了考察和验证模型的预测精度,分别使用原始γ-Reθt模型和γ-Reθt-CF-SA模型,对NLF(2)-0415后掠翼型和标准6∶1椭球模型进行了转捩预测数值模拟。算例结果表明,γ-Reθt-CF-SA模型的计算结果与试验数据吻合程度远远优于原始γ-Reθt模型,γ-Reθt-CF-SA模型能正确地预测出三维流动中的横流不稳定性引起转捩的现象。     

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

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

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

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

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

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

尾迹对涡轮叶栅边界层转捩的影响

为研究周期性尾迹对下游涡轮叶栅边界层转捩的影响机理,并验证γ-Reθ转捩模型在非定常计算中的适用性,在自行开发的CFD程序上实现了该模型,对带有尾迹发生器的T106D-EIZ涡轮叶栅进行数值模拟;对比分析了实验数据与数值模拟的结果,并对流动机理进行阐述。主要结论如下:γ-Reθ转捩模型对周期性尾迹影响下的边界层转捩预测得较为准确,转捩位置随时间的变化规律及幅值均与实验值符合得很好;尾迹的时均作用是使吸力面转捩位置向尾缘移动,相比于无尾迹情况其移动距离约为5%轴向弦长;尾迹使吸力面阻力系数时均值较无尾迹情况增加约40%,同时减弱了分离泡的强度,但对分离泡的位置影响不大。     

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

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

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

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

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

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

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

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

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

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

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

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

压气机静叶栅层流分离泡转捩与角区分离数值模拟与实验

采用γ-Reθ转捩模型对某可控扩散叶型(CDA)平面叶栅全攻角范围进行了三维数值计算，通过对比数值计算结果与叶栅实验吻合较好。在此基础上，分析了进口来流湍流度和雷诺数变化对叶栅表面层流分离、转捩以及角区分离的影响。结果表明:进口湍流度低于5%时，吸力面存在层流分离，当进口湍流度大于5%后，层流分离移除，但转捩会一直存在；随着进口湍流度或雷诺数增加，吸力面和压力面转捩位置均会前移；随着进口湍流度增加，吸力面角区分离会有所减小，雷诺数增加对角区分离的影响不大。      

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

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

涡轮叶栅流动和传热耦合计算

通过求解三维雷诺平均N-S(Navier-Stokes)方程和带转捩模型的二方程SST(Shear stresstransport)湍流模型,完成了带简单冷却的MARKⅡ涡轮导向叶栅流动和传热耦合计算.计算结果与试验数据的对比表明,发展的带转捩模型的数值方法明显地提高了叶栅温度和外换热系数的计算精度.同时研究发现,湍流模型对叶栅表面压力分布影响很小,而对叶栅表面温度和外换热系数影响很大;在其它边界条件相同的情况下进口湍流度对壁面压力几乎无影响,但对叶栅温度和外换热系数影响较大.      

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

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

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

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

跨声速涡轮导叶吸力面换热特性数值研究

为了了解涡轮导叶吸力面在跨声速条件下的换热特性,采用数值模拟的方法,分析了出口马赫数对平面叶栅内流场与换热特性的影响,以及气膜出流对吸力面气膜冷却特性的影响。结果表明,跨声速条件下,斜激波导致的逆压梯度导致了吸力面层流边界层分离和转捩;亚声速条件下,吹风比从0.5增大至1.5时,转捩位置前移了约0.1倍弦长;跨声速条件的转捩位置随吹风比增大未发生变化,但是边界层分离现象被抑制,分离泡的尺寸明显变小。在吸力面小吹风比更容易获得更高的冷却效率;边界层的分离导致冷效率分布不同于亚声速条件,在分离区冷却效率迅速降低,在吹风比0.75时降低约50%。