使用GAO-YONG模型对不可压圆管湍流的计算

基于有限体积法思想运用SIMPLER方法PLS格式求解轴对称柱坐标系下的GAO-YONG方程组,对管流边界层以及圆管湍流的自然转捩过程进行数值模拟。结果表明,方程不仅对于边界层流动的内层壁面律和外层亏损律能作出良好的预测,而且能够预报层流-湍流转捩过程中摩阻系数和速度型的变化。      

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

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

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

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

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

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

平板边界层转捩的数值模拟

对于自由湍流引起的压力梯度为零的平板边界层转捩过程的数值模拟,本文提出了一个单方程湍流模型 (VL M模型),对湍流粘性系数和湍流长度尺度都进行了修正以正确引入低雷诺数效应。数值计算表明,该模型在自由来流湍流度为 0.03～8%的大范围内能正确模拟边界层转捩全过程,包括转捩区起点位置和转捩区长度、边界层内的平均及脉动参数分布等。本文还提出了上述转捩过程的时均化五阶段模型。      

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

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

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

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

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

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

解决湍流圆/平射流异常现象的数值研究

圆射流/平面射流异常现象是湍流研究领域的一个著名的难题,是验证湍流微分方程模型的试金石。本文用GAO-YONG不可压湍流控制方程对圆射流进行了数值验证,成功地解决了该难题,进一步证明了该方程对真实湍流的适定性和通用性。      

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

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

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

Reflected shock/turbulent boundary layer interaction structure analysis based on large eddy simulation

In order to understand the physical phenomenon of the reflected shock/turbulent boundary layer interaction, the Large Eddy Simulation (LES) is conducted to investigate shock wave and turbulent boundary layer interaction in a $12°$ compression ramp with inlet high Mach number of 2.9. Rescaling/recycling method is used as inflow turbulence generation technique and validated on a supersonic flat plate turbulent boundary layer. The flow field of recycling plane in the plate computation domain is obtained to give the inlet boundary condition for the LES computation. This paper focuses on the reflected shock/turbulent boundary layer interaction region, where the fine flow structure and instantaneous flow field are analyzed in detail. It is found that the unsteady motion of the shock wave leads to the increase of wall pressure fluctuation.     

Recent experience with different turbulence models applied to the calculation of flow over aircraft components

An evaluation of ten turbulence models is made for compressible flows encountered in current aircraft applications. The Baldwin–Lomax and P.D. Thomas algebraic models, the Baldwin–Barth and Spalart–Allmaras one-equation models, five low-Reynolds-number k– models and the Menter SST blended k–/k–ω model are examined. A zonal, upwind, implicit, factored algorithm is used to solve both the mean flow equations and the turbulence model equations for three-dimensional, compressible turbulent flow. Calculations are presented for both internal and external flowfields including a two-stream mixing layer, a supersonic flat-plate boundary layer, a transonic supercritical airfoil, a shock wave/turbulent boundary layer interaction, an ejector nozzle, a highly offset diffuser, and a twin impinging jet flowfield. The influence of two modifications to the production of turbulent kinetic energy for the low-Reynolds-number k– models is evaluated, a vorticity-based strain rate and a production limiter. A compressibility correction for high speed shear layers is also examined. Comparisons of the results of the various turbulence models are made with experimental measurements. Significant differences are observed in the model predictions when applied to the same problem using the same computational mesh and mean flow solver. The algebraic models are unable to capture the physics of these complex flowfields, particularly for the internal flow calculations. The performance of each model is dependent on the application. No universal model is found for all flowfields examined. Each one-equation and two-equation model has specific strengths and weaknesses and the performance of each model is assessed.     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.