使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

基于湍动能输运的一方程转捩模型

传统转捩模型的构造十分依赖转捩机理的发现和理解,对机理的模化也存在误差。本文基于尺度自适应的KDO湍流模型,以输运变量r=μt/μ重新标定模型参数,使模型整体达到了流动结构自适应,可实现湍流/转捩一体化计算。根据长度尺度的不同,可细分为KDO-tran和CKDO-tran模型。CFD计算评估了经典的T3A,T3B平板边界层旁路转捩,T3A-平板边界层自然转捩,Aero-A翼型分离泡转捩,DLR-F5,6:1椭球横流转捩及可压缩平板边界层、尖锥转捩。该模型无法精确捕捉层流-湍流的过渡,但能准确捕捉转捩的起始位置(Transition onset),在高Re数、复杂流动的表现尤其突出。该模型捕捉转捩的机制在于流动结构的自适应和湍流输运特性的保存,不引入任何转捩机理却能捕捉多种类型的转捩现象,具有较好的适用性。     

电弧加热流场湍流度对尖锥边界层转捩影响的研究

电弧加热流场的热环境特性直接影响热防护系统的地面试验数据,由于电弧加热器高温气流和参数波动的原因,直接测量湍流度非常困难。为研究电弧加热流场湍流度对于边界层转捩的影响,采用红外热图热像仪,在电弧加热流场中进行了5°尖锥模型边界层转捩研究。结合数值计算,将试验结果与常规风洞的尖锥边界层转捩结果进行了比较。结果表明：马赫数影响的雷诺数转捩判别准则可以用于计算电弧加热流场的转捩雷诺数;电弧加热流场的尖锥边界层转捩雷诺数显著小于常规风洞的转捩雷诺数,表明在该试验条件下,电弧加热流场的湍流度显著大于常规风洞。     

变湍流度时翼型边界层及近场尾流的法向湍流特性初步试验 …

阐述了湍流度对边界层及近场尾流时均特性和边界层转捩位置的影响，讨论了流向和法向湍流特性分布规律，并着重分析了湍流度法法向湍流特性的影响，结果表明，法向湍流强度及雷诺正应力沿流向分段性更为明显，且法向与流向湍流参量之比的特征值有一定规律性。     

Breakdown过程长度及其前后边界层名义厚度的变化

采用数值模拟的方法,模拟了平板边界层的转捩过程及发展的湍流,在此基础上,研究了不同初始扰动对转捩持续时间和湍流边界层厚度影响的共性和差异.研究结果表明:在不同的波数情况下,计算得到不同的初始扰动引起的转捩的持续时间不同.在给定基本波数情况下,如果初始扰动在展向不是对称的,那么,它们引起的转捩后平均边界层厚度与转捩起始点处边界层厚度的比值δ_e/δ_b在相当长的时间内表现出相同的增长规律,这一增长规律与所加的初始扰动无关.然而,展向是非对称的表现出不同增长规律.      

一种基于湍流模式的转捩预测方法

通过对平板边界层进行计算分析,基于湍流模式来预测流动转捩是可行的,但预测出的转捩位置不准,比实际情况靠前。因此,本文从转捩流动的物理特征出发对湍流模式中的混合尺度作出了修正,提出了一种考虑间歇函数影响的k-ωSST两方程湍流模式。采用该模式来对平板湍流边界层和不同攻角下的NACA0012翼型绕流进行了转捩预测,计算结果表明,改进后的模式对转捩位置具有较好的预测能力。     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

较低湍流度范围湍流度对风洞实验结果的影响

变湍流变试验结果表明，湍流度对诸如洞壁边界层、平板边界层，各类翼型边界层和尾流的时均特性，湍流特性，转捩情况以及对翼型和旋成体压力分布等气动特性有显著影响。这种影响有一定规律，量值较大，不能忽视。     

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

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

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

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

使用GAO-YONG湍流方程组对翼型绕流的计算

将GAO-YONG湍流方程组应用到翼型复杂绕流场的计算中,研究了机械能方程中Du/Dt项和各向异性湍流黏性系数分量中的实度系数Cs的影响,给出了控制Du/Dt项作用区域和选择实度系数Cs值的判别函数,对NACA4412翼型最大升力迎角下的低速绕流算例、AS240翼型8°迎角下的转捩算例、AS240翼型19°迎角下的分离流算例进行了模拟,与实验结果比较,得到了比较好的结果.结果表明GAO-YONG湍流方程组能够准确预测翼型前缘层流到湍流的转捩问题,翼型尾迹流动松弛问题和翼型大迎角分离流.     

使用GAO-YONG方程组对湍流边界层拟序结构特性的研究

与以往的湍流模型不同,GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。其级数形式的能量方程与非线性现象多尺度层次相对应,具备了描述湍流统计平均和拟序结构的双重功能。本文采用交错网格和SIMPLE方法求解GAO-YONG不可压湍流方程组,考察了GAO-YONG方程组中能量方程阶数以及计算网格对拟序结构的影响,证明了该方程在描述拟序结构方面具有一定的适定性。      

使用GAO-YONG方程组对可压湍流边界层的数值模拟

应用SIMPLE/SIMPLER方法及QUICK格式求解GAO-YONG可压湍流方程组,对二维零压力梯度可压平板湍流边界层进行了数值模拟。结果表明,不需要任何经验系数及壁面函数的GAO-YONG可压湍流控制方程组对较广马赫数范围(Ma=0.2～2.4)平板边界层各项参数(表面摩擦系数C_f,对数律,亏损律等)能做出良好的预测。      

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

GAO-YONG湍流模型对湍流传热的数值模拟

将GAO-YONG湍流模型应用于湍流传热的研究,分别计算了平板剪切湍流和二维平面冲击射流的湍流传热问题.边界层剪切湍流流动与换热的计算表明:与传统的湍流模型不同,GAO-YONG湍流模型不需要对近壁区域做任何特殊处理(比如壁面函数、低Reynolds数修正等)即可模拟出从壁面到主流区的全部流动与传热情况;另外,对于冲击射流Nusselt数的模拟也得到了与实验符合较好的计算结果,准确地捕捉到了2种冲击高度下流场换热的不同特征,表明了GAO-YONG湍流模型能够较高精度地计算湍流换热.      

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

运用GAO-YONG湍流模型对扩压器内跨声速流动的数值模拟

运用GAO-YONG可压缩湍流方程组,采用同位网格SIMPLE算法,对扩压器跨声速流动中的二维激波/湍流边界层干扰现象进行了数值模拟。将计算得到的流场的时均参数与实验值进行比较,数值模拟结果在激波强度、壁面压力分布以及分离点和再附点位置等方面,与实验值吻合较好,表明GAO-YONG可压缩湍流方程组能够比较准确的模拟较强激波/湍流边界层干扰流动,从而进一步为GAO-YONG湍流模型的正确性及其在可压缩流场模拟方面的适用性提供了佐证。     

湍流模型及施密特数对点火位置的影响研究

准确预测燃料和空气湍流混合水平对精确模拟发动机性能至关重要.通过对Burrows&Kurkov氢气顺喷经典算例的数值模拟,研究了湍流模型和湍流关键参数(施密特数)对氢气和空气的混合程度以及点火位置的影响.数值结果表明:湍流模型对点火的位置具有一定的影响,其中Menter SST k-w湍流模型计算结果与实验值总体吻合最好;湍流施密特数严重影响着氢气和空气的混合程度以及燃烧流场的点火位置,其中Sct=0.5能满足大部分超燃冲压发动机燃烧算例的数值模拟.     

锯齿尾缘翼型近场湍流试验研究

基于锯齿尾缘结构在航空发动机上的应用,对其降噪机理进行研究。通过3维热线风速仪测量2种尾缘结构的尾迹流场揭示锯齿降噪的流动本质,其结果显示出锯齿尾缘后流场的细微湍流结构变化规律,并在尾迹流场可见单个锯齿的齿峰和齿谷。结果表明:锯齿尾缘后尾迹中心线速度的衰减率比直尾缘的高;湍流峰值因为锯齿尾缘的存在出现在离翼型更远处,锯齿在近尾迹区产生了额外的马蹄涡。     

高超声复杂流动中湍流模式应用的评估

本文选择几个高超声速基准流动：二维可压缩拐角、锥柱裙组合体绕流、斜激波与湍流边界层干扰，采用几个常见的湍流模式，BL模式，CHκ-ε模式，κ-ε模式，SST模式，CMOTTκ-ε模式，SHIHκ-ε模式，通过将数值计算结果和实验结果进行比较，对有关的湍流模式地 评估，得到一些有意义的结论。