第三章湍流模型_1.docx

《第三章湍流模型_1.docx》由会员分享，可在线阅读，更多相关《第三章湍流模型_1.docx（93页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三章湍流模型FLUENT讲义：第三章湍流模型第三章湍流模型第一节前言湍流流动模型很多，但大致能够归纳为下面三类：第一类是湍流输运系数模型，是Boussinesq于1877年针对二维流动提出的，将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即：u1u2u131tx2推广到三维问题，若用笛卡儿张量笛卡尔坐标系表示，即有：uiujuiuj232txikijxj3ij为DELT函数，一般i=j时为1，否则为0.模型的任务就是给出计算湍流粘性系数t的方法。根据建立模型所需要的微分方程的数目，能够分为零方程模型代数方程模型，单方程模型和双方程模型。模拟大空间建筑空气流动t=0.03874

2、vl模拟通风空调室内的空气流动比例系数由直接数值模拟的结果拟合而得，其中：v为当地时均速度，l为当地距壁面近期的距离。第二类是抛弃了湍流输运系数的概念，直接建立湍流应力和其它二阶关联量的输运方程。第三类是大涡模拟。前两类是以湍流的统计构造为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。实际求解中，选用什么模型要根据详细问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。参见：湍流模型的选择资料。FLUENT提供的湍流模型包括：单

3、方程Spalart-Allmaras模型、双方程模型标准-模型、重整化群-模型、可实现(Realizable)-模型及雷诺应力模型和大涡模拟。FLUENT讲义：第三章湍流模型Zero-EquationModelsOne-EquationModelsRANS-basedmodels包含更多物理机理Spalart-AllmarasTwo-EquationModels每次迭代Standardk-计算量增加RNGk-Realizablek-FLUENT提Reynolds-StressModel供的模型选择Large-EddySimulationDirectNumericalSimulation湍流模型

4、种类示意图大涡模拟启动需要用命令：(rpsetvarles-2d?#t)第二节平均量输运方程输运经过的粘滞系数、扩散系数和热传导率，故称为输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有：uiuiui33其中，ui和ui分别是平均速度和脉动速度i=1,2,3类似地，对于压力等其它标量，我们也有：34其中，表示标量，如压力、能量、组分浓度等。把上面的表达式代入瞬时的连续与动量方程，并取平均去掉平均速度ui上的横线，我们能够把连续与动量方程写成如下的笛卡儿坐标系下的张量形式：t(ui)35xiDuipuiuj2uluiuj36Dtxixjxj

5、xi3ijxlxj上面两个方程称为雷诺平均的Navier-StokesRANS方程。他们和瞬时Navier-Stokes方程有一样的形式，只是速度或其它求解变量变成了时间平均量。额外多出来的项uiuj是雷诺应力，表示湍流的影响。假如要求解该方程，必须模拟该项以封闭方程。假如密度是变化的流动经过如燃烧问题，我们能够用法夫雷Favre平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外，所有变量都用密度加权平均。变量的密度加权平均定义为：/37符号表示密度加权平均；对应于密度加权平均值的脉动值用表示，即有：。很显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零

6、，即：0，0Boussinesq近似与雷诺应力输运模型为了封闭方程，必须对额外项雷诺应力uiuj进行模拟。一个通常的方法是应用Boussinesq假设，以为雷诺应力与平均速度梯度成正比，即：uiuj2ui)ij38uiujtxi(ktxixj3Boussinesq假设被用于Spalart-Allmaras单方程模型和k双方程模型。Boussinesq近似的好处是与求解湍流粘性系数有关的计算时间比拟少，例如在Spalart-Allmaras单方程模型中，只多求解一个表示湍流粘性的输运方程；在k双方程模型中，只需多求解湍动能k和耗散率两个方程，湍流粘性系数用湍动能k和耗散率的函数。Boussine

7、sq假设的缺点是以为湍流粘性系数t是各向同性标量，对一些复杂流动该条件并不是严格成立，所以具有其应用限制性。另外的方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程，通湍流问题需要多求解7个方程，需要比拟多的计算时间，对计算机内存也有更高要求。在很多问题中，Boussinesq近似方法能够得到比拟好的结果，并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。但是，假如湍流场各向异性很明显，如强旋流动以及应力驱动的二次流等流动中，求解雷诺应力分量输运方程无疑能够得到更好的结果。粘性面板参数讲明：Cmu：(onlyforthestandardorRNG-modelortheRS

8、M)istheconstantthatisusedtocompute.C1-Epsilon：(onlyforthestandardorRNG-modelortheRSM)istheconstantusedinthetransportequationfor.C2-Epsilon：(onlyforthestandard,RNG,orrealizable-modelortheRSM)istheconstantusedinthetransportequationfor.C1-PS(onlyforRSM)istheconstantinEquation.C2-PS(onlyforRSM)isthecons

9、tantinEquationC1-PS(onlyforRSM)istheconstantinEquation.C2-PS(onlyforRSM)istheconstantinEquationPrandtlNumber(onlyfortheSpalart-Allmarasmodel)istheconstantinEquation10.3-1.TKEPrandtlNumber(onlyforthestandardorrealizable-model,thestandardorSST-model,ortheRSM)istheeffectivePrandtlnumberfortransportoftu

10、rbulencekineticenergy.ThiseffectivePrandtlnumberdefinestheratioofthemomentumdiffusivitytothediffusivityofturbulencekineticenergyviaturbulenttransport.TKE(Inner)Prandtl#(onlyfortheSST-model)istheeffectivePrandtlnumberforthetransportofturbulencekineticenergy,insidethenear-wallregion.SeeSection10.5.2fo

11、rdetails.TKE(Outer)Prandtl#(onlyfortheSST-model)istheeffectivePrandtlnumberforthetransportofturbulencekineticenergy,outsidethenear-wallregion.SeeSection10.5.2fordetails.TDRPrandtlNumberistheeffectivePrandtlnumberfortransportoftheturbulentdissipationrate,forthestandardorrealizable-modelortheRSM.Thise

12、ffectivePrandtlnumberdefinestheratioofthemomentumdiffusivitytothediffusivityofturbulencedissipationviaturbulenttransport.Forthestandard-model,theTDRPrandtlNumberistheeffectivePrandtlnumberforthetransportofthespecificdissipationrate,.SDR(Inner)Prandtl#(onlyfortheSST-model)istheeffectivePrandtlnumberf

13、orthetransportofthespecificdissipationrate,insidethenear-wallregion.SeeSection10.5.2fordetails.SDR(Outer)Prandtl#(onlyfortheSST-model)istheeffectivePrandtlnumberforthetransportofthespecificdissipationrate,outsidethenear-wallregion.SeeSection10.5.2fordetails.DispersionPrandtlNumber(onlyforthe-multiph

14、asemodels)istheeffectivePrandtlnumberforthedispersedphase,.SeeSection22.4.8fordetails.EnergyPrandtlNumber(foranyturbulencemodelexcepttheRNG-model)istheturbulentPrandtlnumberforenergy,Pr,inEquation10.6-20.(Thisitemwillnotappearforpremixedorpartiallypremixedcombustionmodels.)WallPrandtlNumber(foralltu

15、rbulencemodels)istheturbulentPrandtlnumberatthewall,PrinEquation10.8-5.(Thisitemwillnotappearforadiabaticpremixedcombustionorpartiallypremixedcombustionmodels.)Turb.SchmidtNumber(forturbulentspeciestransportcalculationsusinganyturbulencemodelexcepttheRNG-model)istheturbulentSchmidtnumber,Sc,inEquati

16、on13.1-3.PDFSchmidtNumber(fornon-premixedorpartiallypremixedcombustioncalculationsusinganyturbulencemodel)isthemodelconstantinEquation14.1-5.具体介绍请访问：:/jullio.pe.kr/fluent6.1/help/html/ug/node1141.htm第三节湍流模型3.3.1单方程Spalart-Allmaras模型，表征出了近壁粘性影响区域以外的湍流运Spalart-Allmaras模型的求解变量是动粘性系数。的输运方程为：D1()39DtGCb2

17、Yxjxjxj其中，G是湍流粘性产生项；Y是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少；和Cb2是常数；是分子运动粘性系数。湍流粘性系数用如下公式计算：tf1其中，f1是粘性阻尼函数，定义为：f1333，并且。C1湍流粘性产生项，G用如下公式模拟：GCb1310SSf2，而f21。其中，Cb1和k是常数，d是计算点其中，S2d2f1k1到壁面的距离；S2ijij。ij定义为：ij1ujui3112xixj由于平均应变率对湍流产生也起到很大作用，FLUENT处理经过中，定义S为：SijCprodmin(0,Sijij)312其中，Cprod2.0，ijijij，Sij2SijSij，平均应变率S

18、ij定义为：Sij1ujui3132xixj在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这合适涡流靠近涡旋中心的区域，那里只要“单纯的旋转，湍流遭到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变，估计的涡旋粘性系数产生项偏高。湍流粘性系数减少项Y为：2YCw1fwd3141Cw61/6其中，fwg3315g6Cw63grCw2(r6r)316r317Sk2d22f2。在上式中，包括了平均应变率其中，Cw1，Cw2，Cw3是常数，SS2dk对S的影响，因此也影响用S计算出来的r。上面的模型常数在FLUENT中默认值为：Cb10.1335，Cb20.622，2/3，C1

19、7.1，Cw1Cb1/k2(1Cb2)/，Cw20.3，Cw32.0，k0.41。壁面条件在壁面，湍流运动粘性设置为零。当计算网格足够细，能够计算层流底层时，壁面切应力用层流应力应变关系求解，即：uuy318u假如网格粗错不能用来求解层流底层，则假设与壁面近邻的网格质心落在边界层的对数区，则根据壁面法则：u1lnEuy319uk其中，k=0.419，E=9.793。对流传热传质模型在FLUENT中，用雷诺类似湍流输运的概念来模拟热输运经过。给出的能量方程为： (E)ui(Ep)kcptTuj(ij)effSh320xiPrtxitxi式中，E是总能量，(ij)eff是偏应力张量，定义为：uju

20、2ui (ij)effeff(i)321xj3effijxixi其中，(ij)eff表示粘性加热，耦合求解。假如默以为分开求解，FLUENT不求解处(ij)eff。但是能够通过变化“粘性模型面板上的湍流普朗特数Prt，其默认值为0.85。Prt数：由流体物性参数组成的一个无因次数即无量纲参数群，表明温度边界层和流动边界层的关系，反映流体物理性质对对流传热经过的影响，它的表达式为：Pr=/=cp/k式中，为动力粘度；cp为等压比热容；k为热导率；为热扩散系数=/c单位：m2/s，v为运动粘度。湍流质量输运与热输运类似，默认的Schmidt数是0.7，该值同样可以以在“粘性模型面板上调节。Schm

21、idt数:表示动量和质量输运之间的关系：粘性系数与扩散系数的比值标量的壁面处理与动量壁面处理类似，分别选用适宜的壁面法则。综上所述，Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，并不需要求解当地剪切层厚度的长度尺度。该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果，在透平机械湍流模拟方面也有较好结果。Spalart-Allmaras模型的初始形式属于对低雷诺数湍流模型，这必须很好解决边界层的当网格比拟粗糙时，网格不知足准确的湍流计算要求，用壁面函数也许是最好的解决方案。另外，该模型中的输运变量在近壁处的梯度要比k中的小，这使得该模型

22、对网格粗糙带来数值误差不太敏感。但是，Spalart-Allmaras模型不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比拟大的流动问题不太合适。比方，平板射流问题，从有壁面影响流动忽然变化到自由剪切流，流场尺度变化明显。3.3.2标准k模型标准k模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过准确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟类似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响能够忽略。因而，标准k模型只合适完全湍流的流动经过模拟。标准k模型的湍动能k和耗散率方程为如下形式：dktkGkGbYM322dtxi

23、xikDt2C1k(GkC3Gb)C2323Dtxikxik在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数tk2C。在FLUENT中，作为默认值常数，C11.44，C2=1.92，C0.09，湍动能k与耗散率的湍流普朗特数分别为k1.0，1.3。能够通过调节“粘性模型面板来调节这些常数值。3.3.3重整化群-模型重整化的一般思想是：减少系统的自由度，并在这个缩减的空间中，通过特定的重整化技巧，在迭代经过中保持系统的自由度数不变，并使约化系统最终收敛到真正系统的低能态中。重整化群-模型是对瞬时的N

24、avier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准-模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准-模型有类似的形式：DkkeffkGkGbYMDtxixiD2effC1k(GkC3Gb)C2RDtxixik324325Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响，这些参数与标准-模型中一样。k和分别是湍动能k和耗散率的有效湍流普朗特数的倒数。湍流粘性系数计算公式为：2kd1.723d3261其中，Ceff/，C100对上面方程积分，能够准确得到有效雷诺数涡旋尺度对

25、湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。对于高雷诺数，上面方程能够给出：tCk2，C0.0845。这个结果非常有意思，和标准-模型的半经历推导给出的常数C0.09非常近似。在FLUENT中，假如是默认设置，用重整化群-模型时候是针对的高雷诺数流动问题。假如对低雷诺数问题进行数值模拟，必须进行相应的设置。重整化群-模型有旋修正通常，平均运动有旋时候对湍流有重要影响。FLUENT中重整化群-模型通过修正湍流粘性系数来考虑了这类影响。湍流粘性的修正形式为：tt0f(s,k)327其中，t0是不考虑有旋计算出来的湍流粘性系数；是FLUENT计算出来的特征旋流数；s是旋流常数，不同值表

26、示有旋流动的强度不同。流动能够是强旋或者中等旋度的。FLUENT默认设置s0.05,针对中等旋度的流动问题，对于强旋流动，能够选择较大的值。湍动能及其耗散率的有效湍流普朗特数倒数的计算公式为：0.63210.36791.39292.3929mol3281.39292.3929eff式中，01，在高雷诺数流动问题中，mol/eff1，k1.393。湍流耗散率方程右边的R为：C3(1/0)2R13k329其中，Sk/，04.38，0.012。为了更清楚体现R对耗散率的影响，我们把耗散率输运方程重写为：D2C2*2xieffxiC1(GkC3Gb)C2k330Dtkk则：C2*C3(1/0)C213

27、331在0的区域，R的奉献为正；C2*大于C2。以对数区为例，3，C2*2.0，这和标准-模型中给出的C21.92接近。因而，对于弱旋和中等旋度的流动问题，重整化群-模型给出的结果比标准-模型的结果要大。重整化群模型中，C11.42，C21.68。3.3.4可实现-模型可实现-模型的湍动能及其耗散率输运方程为：DktkGkGbYM332DtxjxjkDt2C1SC2C1kC3Gb333Dtxjtxjk其中，C1max0.43,，Sk/5在上述方程中，Gk表示由于平均速度梯度引起的湍动能产生，Gb是用于浮力影响引起的湍动能产生；YM可压速湍流脉动膨胀对总的耗散率的影响。C2和C1是常数；k，分别

28、是湍动能及其耗散率的湍流普朗特数。在FLUENT中，作为默认值常数，C11.44，C2=1.9，k1.0，1.2。可实现-模型的湍动能的输运方程与标准-模型和重整化群-模型有一样的形式，只是模型参数不同。但耗散率方程有较大不同。首先耗散率产生项方程右边第二项不包含湍动能产生项Gk，如今的形式更能体现能量在谱空间的传输。另外的特色在于耗散率减少项中，不具有奇异性。并不象标准-模型模型那样把K放在分母上。该模型合适的流动类型比拟广泛，包括有旋均匀剪切流，自由流射流和混合层，腔道流动和边界层流动。对以上流动经过模拟结果都比标准-模型的结果好，十分是可实现-模型对圆口射流和平板射流模拟中，能给出较好的

29、射流扩张角。湍流粘性系数公式为tCk2，这和标准-模型一样。不同的是，在可实现-模型中，C不再是个常数，而是通过如下公式计算：C1334A0AsU*K其中，U*SijSijijij，ijij2ijkk，ijijijkk，ij是isthemeanrate-ofrotationtensorviewedinarotatingreferenceframewiththeangularvelocityk。模型常数A04.04，As6cos，而：1SijSjkSkjSijSij1ujui)arccos(6W)，式中W，S，Sij(xj3S2xi我们能够发现，C是平均应变率与旋度的函数。在平衡边界层惯性底层，

30、能够得到C0.09,与标准-模型中采用底常数一样。双方程模型中，无论是标准-模型、重整化群-模型还是可实现-模型，三个模型有类似的形式，即都有和的输运方程，它们的区别在于：1，计算湍流粘性的方法不同；2，控制湍流扩散的湍流Prandtl数不同；3，方程中的产生项和Gk关系不同。但都包含了一样的表示由于平均速度梯度引起的湍动能产生Gk，用于浮力影响引起的湍动能产生Gb；可压缩湍流脉动膨胀对总的耗散率的影响YM。湍动能产生项Gkuj335uiujxiGbgitT336Prtxi式中，Prt是能量的湍流普特朗数，对于可实现-模型，默认设置值为0.85。对于重整化群-模型，Prt1/，1/Prk/Cp

31、。热膨胀系数1，对Tp于理想气体，浮力引起的湍动能产生项变为：Gbgit337Prtxi在FLUENT程序中，假如有重力作用，并且流场里有密度或者温度的梯度，浮力对湍动能的影响都是存在的。浮力对耗散率的影响不是很清楚，因而，默认设置中，耗散率方程中的浮力影响不被考虑。假如要考虑浮力对耗散率的影响，用“粘性模型面板来控制。浮力对耗散率影响是用C3来体现。但C3并不是常数，而是如下的函数形式：C3tanhv338uv是平行于重力方向的速度分量；u是垂直于重力方向的速度分量。假如流动速度与重力方向一样的剪切流动，C31，对于流动方向与重力方向垂直的剪切流，C30。对于高马赫数的流动问题，可压速性对湍

32、流影响在YM中体现。YM2Mt2其中，Mt是马赫数，定义为：MtkRT是声速。a2a默认设置中，只要选择可压速理想气体，可压速效应都是考虑的。在上述的双方程模型中，对流传热传质模型都是通过雷诺类似湍流动量输运方程得到的。能量方程形式为：(E)uiEpkeffTuj(ij)effSh339xixitxi式中，E是总的能量，keff是有效导热系数；(ij)eff是偏应力张量，定义为： (ij)effeffujui2effuiij340xixj3xi (ij)eff表示的是粘性加热，耦合求解时总是计算。假如不是耦合求解时候，作为默认设置，并不求解该量。假如有需要，需在“粘性模型面板中设置。对于重整化

33、群-模型，有效导热系数为：keffcpeff341用328计算，式中，01/Prk/Cp。事实上，随着mol/eff的变化而变化，这是重整化群-模型的一个优点，由于实验中证实，湍流普朗特数随分子普朗特数及湍流而变化。湍流质量输运处理经过与能量输运经过类似。对于标准-模型和可实现的-模型，默认的Schmidt数是0.7，重整化群模型中，是通过方程328来计算的，其中，1/Sc，Sc是分子Schimidt数。Schmidt数:表示动量和质量输运之间的关系：粘性系数与扩散系数的比值3.3.5雷诺应力模型RSM雷诺应力：湍流动量输送的切向应力雷诺应力模型是求解雷诺应力张量的各个分量的输运方程。详细形式

34、为： (uiuj)(Ukuiuj)txk对流项Cijuiujukp(kjuiikuj)xkuiujxkxk湍流扩撒项DijT分子扩散DijLuiukUjujukUigiujgjuixkxk应力产生项Pij浮力产生项目Gijuiuj2uiujpxixkxkxj压力应变项ij耗散项ij2kujumikmuiumjkm342系统旋转产生项Fij上面方程中，Cij，DijL，Pij，Fij不需要模拟，而DijT，Gij，ij，ij需要模拟以封闭方程。下面简单对几个需要模拟项的模拟。DijT能够用DelayandHarlowL38的梯度扩散模型来模拟，即：TCskukuluiujDij343xkxl但这

35、个模型会导致数值不稳定，因而FLUENT程序中采用标量湍流扩散模型：DijTtuiujxkkxk式中，湍流粘性系数用t344Ck2来计算，根据LienandLeschzinerL98，k0.82，这和标准-模型中选取1.0有所不同。根据GibsonandLaunderL58,FuL55,LaunderL88,L89,压力应变项ij能够分解为三项，即：w345ijij,1ij,2ijij,1，ij,2和wij分别是慢速项，快速项和壁面反射项。C1uiuj2，常数C11.8。ij,1ijkk3ij,2C2PijFijGijCij2ijPGC，C20.60，P1Pkk，32G1Gkk，C1Ckk。2

36、2壁面反射项用于重新分布近壁的雷诺正应力分布，主要是减少垂直于壁面的雷诺正应力，增加平行于壁面的雷诺正应力。该项模拟为：wijC1ukumnknmij3uiuknjnk3ujukninkk3/2k22CldC2km,2nknmij3ik,2njnk3k3/234622jk,2ninkdCl式中，C10.5，C20.3，nk是xk在垂直于壁面方向上的单位分量，d是到壁面的距离；ClC3/4/k，C0.09，k=0.41。默认设置时候，FLUENT不计算wij。假如需要计算时候，在“粘性模型面板中设置。线性压力应变模型对于小雷诺数流动，十分是用双层模型求解近壁流动问题时，FLUENT中通过改良模型

37、常数C1，C2，C1和C2来改良压力应变项LaunderL91。这一经过只要在选择双层流模型时候，在“粘性模型面板上调节。C112.58AA21e(0.0067Ret)2C20.75AC12C11.6732C21C2max36,0C2其中，Retk2/()，参数A和张量不变量A2,A3定义为：A9(A2A3)18A2ikkiA3ikkjji式中，ij是雷诺应力张量各向异性部分，定义为：uiuj2kijij3347k二阶压力应变模型二阶压力应变模型由SpezibleL157等人提出。ij(C1C1*P)bijC2(bikbkj1bmnbmnij)(C3C3*bijbij)kSij3中国科学技术大

38、学FLUENT讲义：第三章湍流模型C4k(bikSjkbjkSik2bmnSmnij)C5k(bikjkbjkik)3483式中，bij是雷诺应力各向异性张量，定义为：uiuj2kijbij33492k平均应变率Sij定义为1ujui；1ujui；模型常数:Sijxixjijxixj22C13.4，C1*1.8，C24.2，C30.8，C3*1.3，C41.25，C50.4。二阶压力应变模型不需要考虑壁面反射影响去模拟对数区湍流边界层经过。浮力对湍流的影响浮力引起的产生项模拟为：GijtgiTT350xjgjPrtxi其中，Prt是能量的湍流普朗特数，默认设置值为0.85。对于理想气体，把热膨

39、胀系数的定义代入上式，得：Gijtgigj351Prtxjxi耗散项ij的模拟耗散张量ij模拟为：2(YM)352ij3ij式中，YM2Mt2，Mt是马赫数；标量耗散率用标准k-模型中的采用的耗散率输运方程求解。雷诺应力模型的边界条件在流场进口，雷诺应力模型需要各个雷诺应力分量和湍动能耗散率的值。这些值能够直接输入，可以以湍流强度和特征长度来计算。在壁面，雷诺应力模型通过壁面函数，给出各个雷诺应力分量和耗散率的值。雷诺应力模型的能量与质量输运方程在雷诺应力模型中，对流传热传质模型都是通过雷诺类似湍流动量输运方程得到的。能量方程形式为：FLUENT讲义：第三章湍流模型(E)uiEp(kcpt)T

40、Sh353xixiPrtuj(ij)efftxi式中，E是总的能量；(ij)eff是偏应力张量，定义为： (ij)effeffujui2effuiij354xixj3xi (ij)eff表示的是粘性加热，耦合求解时总是计算。假如不是耦合求解时候，作为默认设置，并不求解该量，并且Prt=0.85。假如有需要，需在“粘性模型面板中设置。3.3.6大涡模拟LES大涡模拟,英文简称LES(Largeeddysimulation),是近几十年才发展起来的一个流体力学中重要的数值模拟研究方法。它区别于直接数值模拟(DNS)和雷诺平均(RANS)方法。其基本思想是通过准确求解某个尺度以上所有湍流尺度的运动，

41、进而能够捕捉到RANS方法所无能为力的很多非稳态，非平衡经过中出现的大尺度效应和拟序构造，同时又克制了直接数值模拟由于需要求解所有湍流尺度而带来的宏大计算开销的问题，因此被以为是最具有潜力的湍流数值模拟发展方向。由于计算消耗仍然很大，目前大涡模拟还无法在工程上广泛应用，但是大涡模拟技术对于研究很多流动机理问题提供了更为可靠的手段，可为流动控制提供理论基础，并可为工程上广泛应用的RANS方法改良提供指导。湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。LES的基本假设：1，动量、能量、质量及其它标量主要由大涡输运；2，流动的几何和边界条件决定了大涡的特性，而流动特性主要在大涡中体现；3，小尺度涡旋受几何和边界条件影响较小，并且各向同性；大涡模拟经过中，直接求解大涡，小尺度涡旋模拟，进而使得网格要求比DNS低。3.3.6.1大涡模拟的控制方程LES的控制方程是对Navier-Stokes方程在波数空间或者物理空间进行过滤得到的。过滤的经过是去掉比过滤宽度或者给定物理宽度小的涡旋，进而得到大涡旋的控制方程。过滤变量上横线定义为： (x)D(x)G(x,x)dx355其中，D表示流