模型的转捩流动计算分析简

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

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

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

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

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

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

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

γ-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模型能正确地预测出三维流动中的横流不稳定性引起转捩的现象。     

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

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

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

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

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_θ模式应用于高速边界层转捩的研究

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

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

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

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

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

基于改进CST参数化方法和转捩模型的翼型优化设计

为提高翼型优化设计效率,增大设计空间,采用B样条基函数替代传统的形状类别函数(CST)方法中的Bezier多项式,增强了对翼型参数化表达的局部控制能力并提高了翼型局部表达精度。为了确保翼型在优化设计过程中的几何光顺特性和代理模型的准确性,采用小波分解技术提出了多分辨率翼型的局部光顺方法。采用基于本征正交分解(POD)的流场数值代理模型,并结合γ-Reθt转捩模型实现了快速准确的气动力与流动转捩预测。采用小波技术光顺的CST翼型参数化建模、POD流场数值计算代理模型以及γ-Reθt转捩模型,结合遗传算法建立了完整的翼型气动优化设计系统。针对低速层流翼型与超临界翼型进行优化设计,优化设计后的翼型升阻比分别提高了47.42%和45.85%,且对改进前后参数化建模方法的优化性能进行了对比,结果表明本文构建的翼型气动优化设计系统具备很高的优化效率。     

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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