fluent教程湍流模拟解析课件

湍流模拟湍流模拟湍流：非定常，非周期性的三维速度脉动、

强化物质、动量和能量的输运.瞬时速度分解为平均速度和脉动速度:Ui(t)Ui+ui(t)压力、温度、组分浓度值具有类似的脉动什么是湍流?TimeUi(t)Uiui(t)湍流：非定常，非周期性的三维速度脉动、强化物质、动量和能平均量输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有：其中，和分别是平均速度和脉动速度（i=1,2,3）类似地，对于压力等其它标量，我们也有：把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式平均量输运方程雷诺平均就是把Navier-Stokes方程中平均量输运方程（续）上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程。

如果要求解该方程，必须模拟该项以封闭方程。

平均量输运方程（续）上面两个方程称为雷诺平均的Navier-湍流模拟的方法直接数值模拟（DNS）只适合低雷诺数流动。求解雷诺平均的

Navier-Stokes(RANS)方程:

其中(雷诺应力)

时间平均的湍流脉动用实验常数和流场平均速度信息来模拟.大涡模拟（LargeEddySimulation，LES）数值求解大涡，小涡用模型.(定常,不可压缩流动有/无

体积力)湍流模拟的方法直接数值模拟（DNS）只适合低雷诺数流动。(定湍流输运模型涡旋粘性系数模型涡旋粘性系数模型是根据湍流应力和分子运动引起的粘性应力相似而提出来的雷诺应力模型

雷诺应力模型则是从基本方程出发，直接推出雷诺应力的方程，但在雷诺应力的方程中包含有更高阶的相关项，对这些更高阶的相关项再建立相应的模型。湍流输运模型涡旋粘性系数模型湍流应力湍流应力在三维空间中，下标和分别可取为1，2和3，所以湍流应力有9个分量组成，是一个二阶张量，若用矩阵形式表示可写为：主对角线上的三个分量，和称为湍流正应力，其余的六个分量称为湍流切应力。三个湍流正应力之和是湍流脉动动能的两倍湍流应力湍流应力在三维空间中，下标和分别可取湍流应力若以主对角线上的三个分量作为对称轴，则对称的两个切应力分量是相等的，很显然，这是一个对称的二阶张量。很容易可以证明，在各向同性湍流中，湍流正应力的三个分量相等，即湍流应力若以主对角线上的三个分量作为对称轴，则对称的两个切应判断湍流外流内流自然对流沿表面绕流其中其中其它因素，如自由流湍流度,表面条件,扰动可能导致流动从层流向湍流转捩L=x,D,Dh,etc.判断湍流外流内流自然对流沿表面绕流其中其中其它因素，如自由流湍流特点额外应变率流向曲率测向分离加速或减速有旋回流(或分离)二次流3D振荡流动Transpiration(吹风/吸气)自由湍流剪切层相互作用湍流特点额外应变率需要作出选择湍流模型

和近壁处理物理流体精度要求计算资源计算时间要求计算网格需要作出选择湍流模型物理流体精度要求计算资源计算时间要求Zero-EquationModelsOne-EquationModels

Spalart-AllmarasTwo-EquationModels

Standardk-eRNGk-e

Realizablek-e

Reynolds-StressModelLarge-EddySimulationDirectNumericalSimulation湍流模型IncludeMorePhysics每次迭代增加计算时间FLUENT5的湍流模型基于雷诺平均（RANS）的模型Fluent6Zero-EquationModels湍流模型Includ涡旋粘性系数模型

涡旋粘性系数模型是Boussinesq于1877年最早提出来的。他假定在近似平行的剪切流中，湍流应力张量中的切应力分量和平均速度在横向方向的梯度成正比，其比例系数称为涡旋粘性系数，以后把这个假定再推广到三维的流动。即湍流应力张量和平均流场应变率之间有关系：仿照分子运动引起的粘性系数，我们可以假定涡旋粘性系数和湍流中含能涡旋的特征长度和速度尺度成正比，即：而根据确定速度尺度V和长度尺度L方法的不同，又可以细分为各种不同的模型。而其中最简单，V和L的确定方法又是一致的模型应该是双方程模型。涡旋粘性系数模型涡旋粘性系数模型是Boussinesq于1涡旋粘性系数模型在双方程模型中，假定特征速度V和湍流动能k的平方根成正比，特征长度由湍流动能k和另外一个辅助的量确定。比如在k-e模型中，辅助的量选为湍流动能的耗散率，根据量纲分析得长度尺度，模型涡旋粘性系数模型在双方程模型中，假定特征速度V和湍流动能k的涡旋粘性系数模型辅助的量是湍流频率，根据量纲分析得长度尺度模型涡旋粘性系数模型辅助的量是湍流频率，模型推导我们都假定流体不可压，即Favre平均和雷诺平均方程完全相同。瞬时的动量方程减去平均的动量方程得脉动速度的方程

乘以并求平均，利用湍流动能的定义

推导我们都假定流体不可压，即Favre平均和雷诺平均方程完全上式右端第二项可以重新整理成：这里湍流动能的方程，方程右端各项依次为输运项，产生项和耗散项上式右端第二项可以重新整理成：这里湍流动能的方程，方程右端各各项对求导数，乘以，并求平均，可得耗散率的方程

各项对求导数，乘以，并求平均，可得耗散率的方程目前采用的标准k－e模型方程为如果用湍流频率代替湍流动能耗散率，频率的模型方程为：目前采用的标准k－e模型方程为如果用湍流频率代替湍流动能耗RANS方程需要对雷诺应力进行封闭.对于单方程模型，从修正的粘性系数输运方程简介求解湍流粘性系数。对于双方程模型,湍流粘性系数根据湍动能(TKE)和耗散率TKE来确定.TransportequationsforturbulentkineticenergyanddissipationratearesolvedsothatturbulentviscositycanbecomputedforRANSequations.ReynoldsStressTermsinRANS-basedModelsTurbulentKineticEnergy:DissipationRateofTurbulentKineticEnergy:湍流粘性系数:BoussinesqHypothesis:(isotropicstresses)RANS方程需要对雷诺应力进行封闭.ReynoldsSt湍流模型0方程模型单方程模型双方程模型雷诺应力模型大涡模拟湍流模型0方程模型湍流粘性系数:

求解的输运方程:Theadditionalvariablesarefunctionsofthemodifiedturbulentviscosityandvelocitygradients.单方程模型:Spalart-AllmarasGenerationDiffusionDestruction湍流粘性系数:单方程模型:Spalart-Allmaras单方程模型:Spalart-Allmaras用于计算航空中考虑边界作用的高速流动问题可以考虑边界层逆压梯度叶轮机械流动壁面网格可粗可细网格细了，可以考虑低雷诺数流动，考虑边界层影响.网格粗时候，可以得到相对较好的流动结果.不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合。比如，平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显。单方程模型:Spalart-Allmaras用于计算航空中应用举例无粘流结果单方程湍流模型结果应用举例无粘流结果无粘流结果单方程湍流模型结果无粘流结果双方程模型:标准k-

模型湍动能方程耗散率输运方程实验常数(equationswrittenforsteady,incompressiblefloww/obodyforces)ConvectionGenerationDiffusionDestructionDestructionConvectionGenerationDiffusion双方程模型:标准k-模型湍动能方程耗散率输运方程实验双方程模型:标准k-

模型“Baselinemodel”(Two-equation)最广泛用于工程流动与换热计算结果比较透彻，优点、缺点明显半经验性k方程经过严格推导得到

方程通过物理推理得到只对完全湍流适用(大雷诺数流动)对于多数的湍流流动问题有较为合理的结果一般工程流动问题换热问题双方程模型:标准k-模型“Baselinemode双方程模型:Realizablek-与标准k-模型的区别:湍流粘性系数

是变量，其中(A0,As,和U*是速度梯度的函数)雷诺正应力为正;Schwarz不等式成立新的耗散率输运方程:GenerationDiffusionDestructionBuoyancy双方程模型:Realizablek-与标准k-模型与标准的k-模型具有相同的湍动能输运方程在如下流动中，比标准的双方程模型更具有优越性:平板射流、圆射流边界层有强的逆压力梯度，分离有旋（rotation）,回流（recirculation）强曲率影响双方程模型:Realizablek-与标准的k-模型具有相同的湍动能输运方程双方程模型:R双方程模型:RNGk-湍动能方程耗散率输运方程ConvectionDiffusionDissipationGeneration其中用RNG理论推导得到(equationswrittenforsteady,incompressiblefloww/obodyforces)Additionaltermrelatedtomeanstrain&turbulencequantitiesConvectionGenerationDiffusionDestruction双方程模型:RNGk-湍动能方程耗散率输运方程Conv双方程模型:RNGk-k-两个方程都从瞬时N－S方程出发，采用重整化群方法(RenormalizationGroupMethod)严格推导得到.与标准的k-方程有类似的形式、但:在方程中多出一项用于分析快速应变的流动问题考虑了旋流对湍流的贡献采用了分析Prandtl数的分析公式有效粘性系数公式不同对如下流动有较好模拟结果:高流向曲率与应变率流动过渡态流动壁面传热传质问题双方程模型:RNGk-k-两个方程都从瞬时N－S方雷诺应力模型ReynoldsStressModel产生项压力应变项耗散项TurbulentDiffusion(modeled)(relatedtoe)(modeled)(computed)(equationswrittenforsteady,incompressiblefloww/obodyforces)雷诺应力输运方程.Pressure/velocityfluctuationsTurbulenttransport雷诺应力模型ReynoldsStressModel产生项雷诺应力模型（RSM，ReynoldsStressModel）RSM封闭了雷诺平均的Navier-Stokes方程，求解雷诺应力输运方程。方程严格推导得到。（瞬时量方程－雷诺平均方程）×脉动量需要求解耗散率输运方程，用以封闭方程组无需涡粘性各向同性假设雷诺应力输运方程有需要模拟量（项）.RSM有高精度模拟复杂湍流流动的潜质.考虑了曲率,旋流,旋转和高应变率影响气旋流动,有旋燃烧室内流动问题旋转流动通道,二次流雷诺应力模型（RSM，ReynoldsStressMod大涡模拟（LES，LargeEddySimulation）大涡:决定了动量、能量和其它标量场的特性、对平均流动场起主导作用.各向异性,可能是湍流结构历史影响，也许是流动场影响，大涡受到流畅、边界调节和流动参数影响。.小涡:趋于各向同性，受流动影响较小比大涡更易于模拟.LES直接结算(求解)大涡，小涡用模型来模拟(Subgrid-ScaleModeling).计算量较大网格点,NLES

非定常计算大涡模拟（LES，LargeEddySimulationRANS湍流模型比较RANS湍流模型比较近壁处理

靠近壁面区域、多数情况下

k-

和RSM湍流模型的计算结果都不够精确.需要近壁处理.标准壁面函数（Standardwallfunctions）非平衡壁面函数（Nonequilibriumwallfunctions）双层区模型（Two-layerzonalmodelV5）增强的壁面处理（Enhancedwalltreatment，V6）Boundarylayerstructure近壁处理靠近壁面区域、多数情况下k-和RSM湍流模标准壁面函数（StandardWallFunctions）平均速度温度其中whereandPisafunctionofthefluidandturbulentPrandtlnumbers.热边界层厚度标准壁面函数（StandardWallFunctions非平衡壁面函数－NonequilibriumWallFunctionsLog-lawissensitizedtopressuregradientforbetterpredictionofadversepressuregradientflowsandseparation.RelaxedlocalequilibriumassumptionsforTKEinwall-neighboringcells.Thermallaw-of-wallunchangedwhere非平衡壁面函数－NonequilibriumWallFuTwo-LayerZonalModel应用于低雷诺数流动或复杂流动中的近壁流动现象模拟。用点到壁面举例y为特征长度求当地雷诺数，区分区域（层）对于高雷诺数流动，湍流核心区域用k-e

模型来模拟.双层区模型中，在粘性影响的区域只求解k

的输运方程.

通过计算得到（与特征长度相关）.分区是动态过程，并根据计算结果调节分区区域属性.Two-LayerZonalModel应用于低雷诺数流动近壁处理方法的比较近壁处理方法的比较计算网格划分注意事项WallFunctionApproachTwo-LayerZonalModelApproachFirstgridpointinlog-lawregionAtleasttenpointsintheBL.Bettertousestretchedquad/hexcellsforeconomy.Firstgridpointaty+

1.

Atleasttengridpointswithinbuffer&sublayers.Bettertousestretchedquad/hexcellsforeconomy.计算网格划分注意事项WallFunctionApproa估计第一个网格点的位置根据经验与实验公式，估计表面摩擦系数（skinfrictioncoefficient）:平板-管流-计算摩擦速度:反过来计算第一网格点到壁面的距离:壁面函数 •

Two-layermodel用post-processing检查近壁网格的空间分辨率y1=50/uy1=/u估计第一个网格点的位置根据经验与实验公式，估计表面摩擦系数（边界条件设定给出湍流inlets&outlets(potentialbackflow)参数k-

模型需要k和Reynoldsstress模型需要Rij

和

根据不同输入选择，给出自己熟悉或方便给出的量Turbulenceintensity和lengthscale（湍流强度与长度尺度）长度尺度与包含主要能量的大涡旋结构尺度相关.对于边界层流动: l0.4d99流体流过栅板/孔板: l

openingsizeTurbulenceintensityandhydraulicdiameter（湍流强度与水力学直径）适合通道与管内流动Turbulenceintensityandturbulentviscosityratio（湍流强度与粘性比）外流:直接给出k

和

(用户自定义函数，可以给出分布).边界条件设定给出湍流inlets&outlets(p湍流模型用户界面Define

Models

Viscous...湍流模型选项，TurbulenceModeloptions近壁处理，NearWallTreatments无粘，层流，湍流Inviscid,Laminar,orTurbulent其余湍流选项，AdditionalTurbulenceoptions湍流模型用户界面DefineModelsViscExample:ChannelFlowwithConjugateHeatTransferadiabaticwallcoldairV=50fpmT=0°FconstanttemperaturewallT=100°Finsulation1ft1ft10ftPPredictthetemperatureatpointPinthesolidinsulationExample:ChannelFlowwithCon湍流模拟检查是否是湍流

ReDh=5,980低雷诺数下的发展湍流与壁面边界层会导致压力梯度，因此选择：RNGk-

模型，用nonequilibriumwallfunctions处理近壁流动.网格划分简单几何形状

quadrilateralcells估计在水平壁面的垂直方向上梯度比较大，近壁第一个网格要落在对数区（log-lawregion.）流向网格渐大，用于捕捉边界层增长过程.根据求解结果，调节网格，用于提高求解温度边界层（温度梯度）湍流模拟检查是否是湍流ReDh=5,980速度等值线温度等值线上下边界层加速了核心区域的流体流动动量&热边界层预测结果温度边界层清晰可见P速度等值线温度等值线上下边界层加速了核心区域的流体流动动量举例:圆柱绕流wallwall1ft2ft2ft

airV=4fpsComputedragcoefficientofthecylinder5ft14.5ft举例:圆柱绕流wallwall1ft2ft2ft检查是否湍流

ReD=24,600钝体绕流尾迹中，可以看到非定常的涡旋脱落现象。侧下游的分离点计算和靠近壁面的近壁影响，所以需要考虑用,RNGk-

湍流模型，2-layerzonalmodel处理边界.网格划分简单几何形状&BLsquadrilateralcells.边界层网格，2-layermodelfinemeshnearsurface&firstcellaty+

=1.湍流模拟步骤检查是否湍流ReD=24,600湍流模拟步骤圆柱绕流的网格圆柱绕流的网格湍流模拟，涡旋脱落（VortexShedding）有效粘性系数等值线图eff=+tCD

=0.53 StrouhalNumber=0.297其中湍流模拟，涡旋脱落（VortexShedding）有效粘性不同边界结果处理结果比较标准壁面函数非平衡壁面函数双层区壁面模型不同边界结果处理结果比较标准壁面函数非平衡壁面函数双层区壁面小结:湍流模拟成功的湍流模拟需要对实际过程作出如下的正确判断:流体物理特性计算机资源项目要求精度时间可用的湍流模型与近壁处理方法可以先用k-

得到初步结果，不满意改用RNG或Realizablek-.对于高旋流数流动，选择RSM模型.除非是求解低雷诺数流动，或者关心复杂流动的近壁流体特性，否则壁面处理选择壁面函数.小结:湍流模拟成功的湍流模拟需要对实际过程作出如下的正确判湍流模拟湍流模拟湍流：非定常，非周期性的三维速度脉动、

强化物质、动量和能量的输运.瞬时速度分解为平均速度和脉动速度:Ui(t)Ui+ui(t)压力、温度、组分浓度值具有类似的脉动什么是湍流?TimeUi(t)Uiui(t)湍流：非定常，非周期性的三维速度脉动、强化物质、动量和能平均量输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有：其中，和分别是平均速度和脉动速度（i=1,2,3）类似地，对于压力等其它标量，我们也有：把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式平均量输运方程雷诺平均就是把Navier-Stokes方程中平均量输运方程（续）上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程。

如果要求解该方程，必须模拟该项以封闭方程。

平均量输运方程（续）上面两个方程称为雷诺平均的Navier-湍流模拟的方法直接数值模拟（DNS）只适合低雷诺数流动。求解雷诺平均的

Navier-Stokes(RANS)方程:

其中(雷诺应力)

时间平均的湍流脉动用实验常数和流场平均速度信息来模拟.大涡模拟（LargeEddySimulation，LES）数值求解大涡，小涡用模型.(定常,不可压缩流动有/无

体积力)湍流模拟的方法直接数值模拟（DNS）只适合低雷诺数流动。(定湍流输运模型涡旋粘性系数模型涡旋粘性系数模型是根据湍流应力和分子运动引起的粘性应力相似而提出来的雷诺应力模型

雷诺应力模型则是从基本方程出发，直接推出雷诺应力的方程，但在雷诺应力的方程中包含有更高阶的相关项，对这些更高阶的相关项再建立相应的模型。湍流输运模型涡旋粘性系数模型湍流应力湍流应力在三维空间中，下标和分别可取为1，2和3，所以湍流应力有9个分量组成，是一个二阶张量，若用矩阵形式表示可写为：主对角线上的三个分量，和称为湍流正应力，其余的六个分量称为湍流切应力。三个湍流正应力之和是湍流脉动动能的两倍湍流应力湍流应力在三维空间中，下标和分别可取湍流应力若以主对角线上的三个分量作为对称轴，则对称的两个切应力分量是相等的，很显然，这是一个对称的二阶张量。很容易可以证明，在各向同性湍流中，湍流正应力的三个分量相等，即湍流应力若以主对角线上的三个分量作为对称轴，则对称的两个切应判断湍流外流内流自然对流沿表面绕流其中其中其它因素，如自由流湍流度,表面条件,扰动可能导致流动从层流向湍流转捩L=x,D,Dh,etc.判断湍流外流内流自然对流沿表面绕流其中其中其它因素，如自由流湍流特点额外应变率流向曲率测向分离加速或减速有旋回流(或分离)二次流3D振荡流动Transpiration(吹风/吸气)自由湍流剪切层相互作用湍流特点额外应变率需要作出选择湍流模型

和近壁处理物理流体精度要求计算资源计算时间要求计算网格需要作出选择湍流模型物理流体精度要求计算资源计算时间要求Zero-EquationModelsOne-EquationModels

Spalart-AllmarasTwo-EquationModels

Standardk-eRNGk-e

Realizablek-e

Reynolds-StressModelLarge-EddySimulationDirectNumericalSimulation湍流模型IncludeMorePhysics每次迭代增加计算时间FLUENT5的湍流模型基于雷诺平均（RANS）的模型Fluent6Zero-EquationModels湍流模型Includ涡旋粘性系数模型

涡旋粘性系数模型是Boussinesq于1877年最早提出来的。他假定在近似平行的剪切流中，湍流应力张量中的切应力分量和平均速度在横向方向的梯度成正比，其比例系数称为涡旋粘性系数，以后把这个假定再推广到三维的流动。即湍流应力张量和平均流场应变率之间有关系：仿照分子运动引起的粘性系数，我们可以假定涡旋粘性系数和湍流中含能涡旋的特征长度和速度尺度成正比，即：而根据确定速度尺度V和长度尺度L方法的不同，又可以细分为各种不同的模型。而其中最简单，V和L的确定方法又是一致的模型应该是双方程模型。涡旋粘性系数模型涡旋粘性系数模型是Boussinesq于1涡旋粘性系数模型在双方程模型中，假定特征速度V和湍流动能k的平方根成正比，特征长度由湍流动能k和另外一个辅助的量确定。比如在k-e模型中，辅助的量选为湍流动能的耗散率，根据量纲分析得长度尺度，模型涡旋粘性系数模型在双方程模型中，假定特征速度V和湍流动能k的涡旋粘性系数模型辅助的量是湍流频率，根据量纲分析得长度尺度模型涡旋粘性系数模型辅助的量是湍流频率，模型推导我们都假定流体不可压，即Favre平均和雷诺平均方程完全相同。瞬时的动量方程减去平均的动量方程得脉动速度的方程

乘以并求平均，利用湍流动能的定义

推导我们都假定流体不可压，即Favre平均和雷诺平均方程完全上式右端第二项可以重新整理成：这里湍流动能的方程，方程右端各项依次为输运项，产生项和耗散项上式右端第二项可以重新整理成：这里湍流动能的方程，方程右端各各项对求导数，乘以，并求平均，可得耗散率的方程

各项对求导数，乘以，并求平均，可得耗散率的方程目前采用的标准k－e模型方程为如果用湍流频率代替湍流动能耗散率，频率的模型方程为：目前采用的标准k－e模型方程为如果用湍流频率代替湍流动能耗RANS方程需要对雷诺应力进行封闭.对于单方程模型，从修正的粘性系数输运方程简介求解湍流粘性系数。对于双方程模型,湍流粘性系数根据湍动能(TKE)和耗散率TKE来确定.TransportequationsforturbulentkineticenergyanddissipationratearesolvedsothatturbulentviscositycanbecomputedforRANSequations.ReynoldsStressTermsinRANS-basedModelsTurbulentKineticEnergy:DissipationRateofTurbulentKineticEnergy:湍流粘性系数:BoussinesqHypothesis:(isotropicstresses)RANS方程需要对雷诺应力进行封闭.ReynoldsSt湍流模型0方程模型单方程模型双方程模型雷诺应力模型大涡模拟湍流模型0方程模型湍流粘性系数:

求解的输运方程:Theadditionalvariablesarefunctionsofthemodifiedturbulentviscosityandvelocitygradients.单方程模型:Spalart-AllmarasGenerationDiffusionDestruction湍流粘性系数:单方程模型:Spalart-Allmaras单方程模型:Spalart-Allmaras用于计算航空中考虑边界作用的高速流动问题可以考虑边界层逆压梯度叶轮机械流动壁面网格可粗可细网格细了，可以考虑低雷诺数流动，考虑边界层影响.网格粗时候，可以得到相对较好的流动结果.不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合。比如，平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显。单方程模型:Spalart-Allmaras用于计算航空中应用举例无粘流结果单方程湍流模型结果应用举例无粘流结果无粘流结果单方程湍流模型结果无粘流结果双方程模型:标准k-

模型湍动能方程耗散率输运方程实验常数(equationswrittenforsteady,incompressiblefloww/obodyforces)ConvectionGenerationDiffusionDestructionDestructionConvectionGenerationDiffusion双方程模型:标准k-模型湍动能方程耗散率输运方程实验双方程模型:标准k-

模型“Baselinemodel”(Two-equation)最广泛用于工程流动与换热计算结果比较透彻，优点、缺点明显半经验性k方程经过严格推导得到

方程通过物理推理得到只对完全湍流适用(大雷诺数流动)对于多数的湍流流动问题有较为合理的结果一般工程流动问题换热问题双方程模型:标准k-模型“Baselinemode双方程模型:Realizablek-与标准k-模型的区别:湍流粘性系数

是变量，其中(A0,As,和U*是速度梯度的函数)雷诺正应力为正;Schwarz不等式成立新的耗散率输运方程:GenerationDiffusionDestructionBuoyancy双方程模型:Realizablek-与标准k-模型与标准的k-模型具有相同的湍动能输运方程在如下流动中，比标准的双方程模型更具有优越性:平板射流、圆射流边界层有强的逆压力梯度，分离有旋（rotation）,回流（recirculation）强曲率影响双方程模型:Realizablek-与标准的k-模型具有相同的湍动能输运方程双方程模型:R双方程模型:RNGk-湍动能方程耗散率输运方程ConvectionDiffusionDissipationGeneration其中用RNG理论推导得到(equationswrittenforsteady,incompressiblefloww/obodyforces)Additionaltermrelatedtomeanstrain&turbulencequantitiesConvectionGenerationDiffusionDestruction双方程模型:RNGk-湍动能方程耗散率输运方程Conv双方程模型:RNGk-k-两个方程都从瞬时N－S方程出发，采用重整化群方法(RenormalizationGroupMethod)严格推导得到.与标准的k-方程有类似的形式、但:在方程中多出一项用于分析快速应变的流动问题考虑了旋流对湍流的贡献采用了分析Prandtl数的分析公式有效粘性系数公式不同对如下流动有较好模拟结果:高流向曲率与应变率流动过渡态流动壁面传热传质问题双方程模型:RNGk-k-两个方程都从瞬时N－S方雷诺应力模型ReynoldsStressModel产生项压力应变项耗散项TurbulentDiffusion(modeled)(relatedtoe)(modeled)(computed)(equationswrittenforsteady,incompressiblefloww/obodyforces)雷诺应力输运方程.Pressure/velocityfluctuationsTurbulenttransport雷诺应力模型ReynoldsStressModel产生项雷诺应力模型（RSM，ReynoldsStressModel）RSM封闭了雷诺平均的Navier-Stokes方程，求解雷诺应力输运方程。方程严格推导得到。（瞬时量方程－雷诺平均方程）×脉动量需要求解耗散率输运方程，用以封闭方程组无需涡粘性各向同性假设雷诺应力输运方程有需要模拟量（项）.RSM有高精度模拟复杂湍流流动的潜质.考虑了曲率,旋流,旋转和高应变率影响气旋流动,有旋燃烧室内流动问题旋转流动通道,二次流雷诺应力模型（RSM，ReynoldsStressMod大涡模拟（LES，LargeEddySimulation）大涡:决定了动量、能量和其它标量场的特性、对平均流动场起主导作用.各向异性,可能是湍流结构历史影响，也许是流动场影响，大涡受到流畅、边界调节和流动参数影响。.小涡:趋于各向同性，受流动影响较小比大涡更易于模拟.LES直接结算(求解)大涡，小涡用模型来模拟(Subgrid-ScaleModeling).计算量较大网格点,NLES

非定常计算大涡模拟（LES，LargeEddySimulationRANS湍流模型比较RANS湍流模型比较近壁处理

靠近壁面区域、多数情况下

k-

和RSM湍流模型的计算结果都不够精确.需要近壁处理.标准壁面函数（Standardwallfunctions）非平衡壁面函数（Nonequilibriumwallfunctions）双层区模型（Two-layerzonalmodelV5）增强的壁面处理（Enhancedwalltreatment，V6）Boundarylayerstructure近壁处理靠近壁面区域、多数情况下k-和RSM湍流模标准壁面函数（StandardWallFunctions）平均速度温度其中whereandPisafunctionofthefluidandturbulentPrandtlnumbers.热边界层厚度标准壁面函数（StandardWallFunctions非平衡壁面函数－NonequilibriumWallFunctionsLog-lawissensitizedtopressuregradientforbetterpredictionofadversepressuregradientflowsandseparation.RelaxedlocalequilibriumassumptionsforTKEinwall-neighboringcells.Thermallaw-of-wallunchangedwhere非平衡壁面函数－NonequilibriumWallFuTwo-LayerZonalModel应用于低雷诺数流动或复杂流动中的近壁流动现象模拟。用点到壁面举例y为特征长度求当地雷诺数，区分区域（层）对于高雷诺数流动，湍流核心区域用k-e

模型来模拟.双层区模型中，在粘性影响的区域只求解k

的输运方程.

通过计算得到（与特征长度相关）.分区是动态过程，并根据计算结果调节分区区域属性.Two-LayerZonalModel应用于低雷诺数流动近壁处理方法的比较近壁处理方法的比较计算网格划分注意事项WallFunctionApproachTwo-LayerZonalModelApproachFirstgridpointinlog-lawregionAtleasttenpointsintheBL.Bettertousestretchedquad/hexcellsforeconomy.Firstgridpointaty+

1.

Atleasttengridpointswithinbuffer&sublayers.Bettertousestretchedquad/hexcellsforeconomy.计算网格划分注意事项WallFunctionApproa估计第一个网格点的位置根据经验与实验公式，估计表面摩擦系数（skinfrictioncoefficient）:平板-管流-计算摩擦速度:反过来计算第一网格点到壁面的距离:壁面函数 •

Two-layermodel用post-processing检查近壁网格的空间分辨率y1=50/uy1=/u估计第一个网格点的位置根据经验与实验公式，估计表面摩擦系数（边界条件设定给出湍流inlets&outlets(potentialbackflow)参数k-

模型需要k和Reynoldsstress模型需要Rij

和

根据不同输入选择，给出自己熟悉或方便给出的量Turbulenceintensity和lengthscale（湍流强度与长度尺度）长度尺度与包含主要能量的大涡