第三章 湍流模型
第一节 前言
湍流流动模型很多,但大致可以归纳为以下三类:
第一类是湍流输运系数模型,是Boussinesq于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即:
                                                    3-1
推广到三维问题,若用笛卡儿张量表示,即有:
                                    3-2
模型的任务就是给出计算湍流粘性系数的方法。根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。solowheel电动独轮车
实际求解中,选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT提供的湍流模型包括:单方程(Spalart-Allmaras)模型、双方程模型(标准κ-ε模型、重整化κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
Zero-Equation Models
One-Equation Models
    Spalart-Allmaras
Two-Equation Models
    Standard k-ε
    RNG k-ε上海二手车买卖
Realizable k-ε   
Reynolds-Stress Model
Large-Eddy Simulation
Direct Numerical Simulation
包含更多丰田威飒
物理机理
每次迭代
计算量增加
FLUENT提供的模型选择
RANS-based
models
湍流模型种类示意图

第二节 平均量输运方程
雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度,有:
                                                        3-3
其中,分别是平均速度和脉动速度(i=1,2,3)
类似地,对于压力等其它标量,我们也有:
                                                          3-4
其中,表示标量,如压力、能量、组分浓度等。
把上面的表达式代入瞬时的连续与动量方程,并取平均(去掉平均速度上的横线),我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式:
                                                  3-5
        3-6 
上面两个方程称为雷诺平均的Navier-Stokes(RANS)方程。他们和瞬时Navier-Stokes方程有相同的形式,只是速度或其它求解变量变成了时间平均量。额外多出来的项是雷诺应力,表示湍流的影响。如果要求解该方程,必须模拟该项以封闭方程。
如果密度是变化的流动过程如燃烧问题,我们可以用法夫雷(Favre)平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外,所有变量都用密度加权平均。变量的密度加权平均定义为:
                                                          3-7
符号~表示密度加权平均;对应于密度加权平均值的脉动值用表示,即有:。很显然,这种脉动值的简单平均值不为零,但它的密度加权平均值等于零,即:
Boussinesq近似与雷诺应力输运模型
为了封闭方程,必须对额外项雷诺应力进行模拟。一个通常的方法是应用Boussinesq假设,认为雷诺应力与平均速度梯度成正比,即:
                            3-8
Boussinesq雪佛兰乐风报价假设被用于Spalart-Allmaras单方程模型和双方程模型。Boussinesq近似的好处是与求解湍流粘性系数有关的计算时间比较少,例如在Spalart-Allmaras单方程模型中,只多求解一个表示湍流粘性的输运方程;在双方程模型中,只需多求解湍动能k和耗散率ε两个方程,湍流粘性系数用湍动能k和耗散率ε的函数Boussinesq假设的缺点是认为湍流粘性系数是各向同性标量,对一些复杂流动该条件并不是严格成立,所以具有其应用限制性。
另外的方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程,通常是耗散率ε方程。这就意味着对于二维湍流流动问题,需要多求解4个输运方程,而三维湍流问题需要多求解7个方程,需要比较多的计算时间,对计算机内存也有更高要求。
在许多问题中,Boussinesq近似方法可以得到比较好的结果,并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。但是,如果湍流场各向异性很明显,如强旋流动以及应力驱动的二次流等流动中,求解雷诺应力分量输运方程无疑可以得到更好的结果。
 
第三节 湍流模型
3.3.1 单方程(Spalart-Allmaras)模型
Spalart-Allmaras模型的求解变量是,表征出了近壁(粘性影响)区域以外的湍流运动粘性系数。的输运方程为:
              3-9
悦达起亚k5价格
其中,是湍流粘性产生项;是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少;是常数;ν是分子运动粘性系数。
湍流粘性系数用如下公式计算:
其中,是粘性阻尼函数,定义为:,并且
湍流粘性产生项,用如下公式模拟:
                                                      3-10
其中,,而。其中,k是常数,d是计算点到壁面的距离;S定义为:
                                                3-11
由于平均应变率对湍流产生也起到很大作用,FLUENT处理过程中,定义S为:
                                    3-12
其中,,平均应变率定义为:
什么是esp
                                                  3-13
在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这适合涡流靠近涡旋中心的
区域,那里只有“单纯”的旋转,湍流受到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变,估计的涡旋粘性系数产生项偏高。
湍流粘性系数减少项为:
                                                    3-14
其中,                                          3-15
                                            3-16
                                                    3-17
其中,是常数,。在上式中,包括了平均应变率对S的影响,因而也影响用计算出来的r
上面的模型常数在FLUENT中默认值为:
壁面条件
在壁面,湍流运动粘性设置为零。当计算网格足够细,可以计算层流底层时,壁面切应力用层流应力-应变关系求解,即:
                                                        3-18
如果网格粗错不能用来求解层流底层,则假设与壁面近邻的网格质心落在边界层的对数区,则根据壁面法则:
                                                  3-19