第3章,fluent湍流模型-1

24第三章，湍流模型第三章，湍流模型第一节，前言湍流流动模型很多，但大致可以归纳为以下三类： 第一类是湍流输运系数模型，是 Boussinesq 于 1877 年针对二维流动提出的，将速度脉 动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积即：3－121 21xuuut推广到三维问题，若用笛卡儿张量表示，即有：3－2ij ijji tjikxuxuuu32   模型的任务就是给出计算湍流粘性系数的方法根据建立模型所需要的微分方程的数t目，可以分为零方程模型（代数方程模型） ，单方程模型和双方程模型 第二类是抛弃了湍流输运系数的概念，直接建立湍流应力和其它二阶关联量的输运方程第三类是大涡模拟前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均大 涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的 Navier-Stokes 方程， 得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型 实际求解中，选用什么模型要根据具体问题的特点来决定选择的一般原则是精度要高， 应用简单，节省计算时间，同时也具有通用性 FLUENT 提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准 κ-ε 模型、重整化群 κ-ε 模型、可实现(Realizable)κ-ε 模型）及雷诺应力模型和大涡 模拟。

25Zero-Equation ModelsOne-Equation ModelsSpalart-Allmaras Two-Equation ModelsStandard k- RNG k-  Realizable k-   Reynolds-Stress ModelLarge-Eddy SimulationDirect Numerical Simulation包含更多包含更多 物理机理物理机理每次迭代每次迭代 计算量增加计算量增加FLUENT提 供的模型选 择RANS-based models湍流模型种类示意图 第二节，平均量输运方程雷诺平均就是把 Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分对于速 度，有：3－3iiiuuu其中，和分别是平均速度和脉动速度（i=1,2,3）iuiu类似地，对于压力等其它标量，我们也有：3－4其中，表示标量，如压力、能量、组分浓度等把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度上的横线） ，iu我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式：3－50)(i iuxt3－6 ji jll ij ijjijiiuuxxu xuxu xxp DtDu        32上面两个方程称为雷诺平均的 Navier-Stokes（RANS）方程。

他们和瞬时 Navier-Stokes 方程有相同的形式，只是速度或其它求解变量变成了时间平均量额外多出来的项是雷诺应力，表示湍流的影响如果要求解该方程，必须模拟该项以封闭方程jiuu如果密度是变化的流动过程如燃烧问题，我们可以用法夫雷（Favre）平均这样才可以 求解有密度变化的流动问题法夫雷平均就是出了压力和密度本身以外，所有变量都用密度 加权平均变量的密度加权平均定义为：3－7/~26符号～表示密度加权平均；对应于密度加权平均值的脉动值用表示，即有： 很显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零， ~即：， 0 0 Boussinesq 近似与雷诺应力输运模型近似与雷诺应力输运模型为了封闭方程，必须对额外项雷诺应力进行模拟一个通常的方法是应用jiuuBoussinesq 假设，认为雷诺应力与平均速度梯度成正比，即：3－8ij ii t ijji tjixukxuxuuu)(32    Boussinesq 假设被用于 Spalart-Allmaras 单方程模型和双方程模型。

Boussinesq 近似k 的好处是与求解湍流粘性系数有关的计算时间比较少，例如在 Spalart-Allmaras 单方程模型中， 只多求解一个表示湍流粘性的输运方程；在双方程模型中，只需多求解湍动能 k 和耗散k 率 ε 两个方程，湍流粘性系数用湍动能 k 和耗散率 ε 的函数Boussinesq 假设的缺点是认 为湍流粘性系数是各向同性标量，对一些复杂流动该条件并不是严格成立，所以具有其应t用限制性 另外的方法是求解雷诺应力各分量的输运方程这也需要额外再求解一个标量方程，通 常是耗散率 ε 方程这就意味着对于二维湍流流动问题，需要多求解 4 个输运方程，而三维 湍流问题需要多求解 7 个方程，需要比较多的计算时间，对计算机内存也有更高要求 在许多问题中，Boussinesq 近似方法可以得到比较好的结果，并不一定需要花费很多时 间来求解雷诺应力各分量的输运方程但是，如果湍流场各向异性很明显，如强旋流动以及 应力驱动的二次流等流动中，求解雷诺应力分量输运方程无疑可以得到更好的结果第三节，第三节，湍流模型湍流模型3.3.1 单方程（单方程（Spalart-Allmaras）模型）模型 Spalart-Allmaras 模型的求解变量是，表征出了近壁（粘性影响）区域以外的湍流运动~粘性系数。

的输运方程为：~3－9 YxCxxGDtDjb jj       ~~ )~(1~2 ~其中，是湍流粘性产生项；是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少；GY和是常数；ν 是分子运动粘性系数~2bC湍流粘性系数用如下公式计算：1~ ft27其中，是粘性阻尼函数，定义为：，并且1f3 1331  Cf~ 湍流粘性产生项，用如下公式模拟：G3－10~~ 1SCGb其中，，而其中，和 k 是常数，d 是计算点到222~~ fdkSS1211 ff1bC壁面的距离；S定义为：ijij2ij3－11   jiij ijxu xu21由于平均应变率对湍流产生也起到很大作用，FLUENT 处理过程中，定义 S 为：3－12), 0min(ijijprodijSCS其中，，，，平均应变率定义为：0 . 2prodCijijijijijijSSS2ijS3－13   jiij ijxu xuS21在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。

这适合涡流靠近涡旋中心 的区域，那里只有“单纯”的旋转，湍流受到抑止包含应变张量的影响更能体现旋转对湍 流的影响忽略了平均应变，估计的涡旋粘性系数产生项偏高 湍流粘性系数减少项为：Y3－1421~ dfCYww其中， 3－15 6/16 366 31  ww wCgCgf3－16)(6 2rrCrgw3－1722~~dkSr其中，，，是常数，在上式中，包括了平均应变率对1wC2wC3wC222~~ fdkSSS 的影响，因而也影响用计算出来的 rS~上面的模型常数在 FLUENT 中默认值为：，，，1335. 01bC622. 02bC3/2~，，，，1 . 71C~22 11/ )1 (/bbwCkCC3 . 02wC0 . 23wC41. 0k壁面条件在壁面，湍流运动粘性设置为零当计算网格足够细，可以计算层流底层时，壁面切~28应力用层流应力－应变关系求解，即：3－18yu uu如果网格粗错不能用来求解层流底层，则假设与壁面近邻的网格质心落在边界层的对数 区，则根据壁面法则：3－19   yuEkuuln1其中，k=0.419，E=9.793。

对流传热传质模型对流传热传质模型在 FLUENT 中，用雷诺相似湍流输运的概念来模拟热输运过程给出的能量方程为：3－20heffijj itpii iSuxT tckxpEuxEt       )(Pr)]([)(式中，E 是总能量，是偏应力张量，定义为：effij)(3－21ij ii eff jiij effeffijxu xu xu32)()(其中，表示粘性加热，耦合求解如果默认为分开求解，FLUENT 不求解处effij)(但是可以通过变化“粘性模型”面板上的湍流普朗特数（Prt） ，其默认值为 0.85effij)(湍流质量输运与热输运类似，默认的 Schmidt 数是 0.7，该值同样也可以在“粘性模型” 面板上调节 标量的壁面处理与动量壁面处理类似，分别选用合适的壁面法则 综上所述，Spalart-Allmaras 模型是相对简单的单方程模型，只需求解湍流粘性的输运方 程，并不需要求解当地剪切层厚度的长度尺度该模型对于求解有壁面影响流动及有逆压力 梯度的边界层问题有很好模拟效果，在透平机械湍流模拟方面也有较好结果。

Spalart-Allmaras 模型的初始形式属于对低雷诺数湍流模型，这必须很好解决边界层的粘 性影响区求解问题在 FLUENT 中，当网格不是很细时，采用壁面函数来解决这一问题当 网格比较粗糙时，网格不满足精确的湍流计算要求，用壁面函数也许是最好的解决方案另 外，该模型中的输运变量在近壁处的梯度要比中的小，这使得该模型对网格粗糙带来数k 值误差不太敏感 但是，Spalart-Allmaras 模型不能预测均匀各向同性湍流的耗散并且，单方程模型没有 考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合比如，平板射流 问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显3.3.2 标准标准模型模型k标准模型需要求解湍动能及其耗散率方程湍动能输运方程是通过精确的方程推导k29得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的该模型假设流动为 完全湍流，分子粘性的影响可以忽略因此，标准模型只适合完全湍流的流动过程模拟k标准模型的湍动能 k 和耗散率 ε 方程为如下形式：k3－22Mbk iktiYGGxk xDtDk       3－23kCGCGkCxxDtDbk ikti2231)(        在上述方程中，表示由于平均速度梯度引起的湍动能产生，是用于浮力影响引起kGbG的湍动能产生；可压速湍流脉动膨胀对总的耗散率的影响。

湍流粘性系数MY2kCt在 FLUENT 中，作为默认值常数，＝1.44，=1.92，，湍动能 k 与耗1C2C09. 0C散率 ε 的湍流普朗特数分别为＝1.0，＝1.3可以通过调节“粘性模型”面板来调节k这些常数值3.3.33.3.3 重整化群重整化群 κ-εκ-ε 模型模型重整化群 κ-ε 模型是对瞬时的 Navier-Stokes 方程用重整化群的数学方法推导出来的 模型模型中的常数与标准 κ-ε 模型不同，而且方程中也出现了新的函数或者项其湍动 能与耗散率方程与标准 κ-ε 模型有相似的形式：3－24Mbk ieffk iYG。