第三章湍流模型
第三章 湍流模型
第一节 前言
湍流流动模型很多,但大致可以归纳为以下三类:
第一类是湍流输运系数模型,是Boussinesq于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即:
u2tu1
u1
3-1 x2
2
kij 3-2 3
推广到三维问题,若用笛卡儿张量表示,即有:
uuj
uiujti
x
jxi
模型的任务就是给出计算湍流粘性系数t的方法。根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT提供的湍流模型包括:单方程(Spalart-Allmaras)模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
包含更多 物理机理
RANS-based models
每次迭代 计算量增加
提的模型选
Direct Numerical Simulation
湍流模型种类示意图
第二节 平均量输运方程
雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度,有:
uiuiui 3-3
其中,ui和ui分别是平均速度和脉动速度(i=1,2,3)
类似地,对于压力等其它标量,我们也有:
3-4 其中,表示标量,如压力、能量、组分浓度等。
把上面的表达式代入瞬时的连续与动量方程,并取平均(去掉平均速度ui上的横线),我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式:
(ui)0 3-5 txi
DupiDtxixj
uiuj2ul
ijuiuj 3-6 xjxi3xlxj
上面两个方程称为雷诺平均的Navier-Stokes(RANS)方程。他们和瞬时Navier-Stokes方程有相同的形式,只是速度或其它求解变量变成了时间平均量。额外多出来的项uiuj是雷诺应力,表示湍流的影响。如果要求解该方程,必须模拟该项以封闭方程。
如果密度是变化的流动过程如燃烧问题,我们可以用法夫雷(Favre)平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外,所有变量都用密度加权平均。变量的密度加权平均定义为:
~
/ 3-7
~
。很显然,这种脉动值的简单平均值不为零,但它的密度加权平均值等于零,
即:
符号~表示密度加权平均;对应于密度加权平均值的脉动值用表示,即有:
0, 0
Boussinesq近似与雷诺应力输运模型
为了封闭方程,必须对额外项雷诺应力uiuj进行模拟。一个通常的方法是应用Boussinesq假设,认为雷诺应力与平均速度梯度成正比,即:
uiuj
uiujt
x
jxi
2
(ktui)ij 3-8 3xi
Boussinesq假设被用于Spalart-Allmaras单方程模型和k双方程模型。Boussinesq近
例如在Spalart-Allmaras单方程模k双方程模型中,只需多求解湍动能k和耗散率ε两个方程,湍流粘性系数用湍动能k和耗散率ε的函数。Boussinesq假设的t是各向同性标量,对一些复杂流动该条件并不是严格成立,所以具有其应用限制性。
另外的方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程,通常是耗散率ε方程。这就意味着对于二维湍流流动问题,需要多求解4个输运方程,而三维湍流问题需要多求解7个方程,需要比较多的计算时间,对计算机内存也有更高要求。
在许多问题中,Boussinesq近似方法可以得到比较好的结果,并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。但是,如果湍流场各向异性很明显,如强旋流动以及应力驱动的二次流等流动中,求解雷诺应力分量输运方程无疑可以得到更好的结果。
第三节 湍流模型
3.3.1 单方程(Spalart-Allmaras)模型
~,表征出了近壁(粘性影响)区域以外的湍流运Spalart-Allmaras模型的求解变量是
~的输运方程为: 动粘性系数。
~~~D1~Y 3-9 G()Cb2Dt~xjxjxj
其中,G是湍流粘性产生项;Y是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少;
~和Cb2是常数;ν是分子运动粘性系数。
湍流粘性系数用如下公式计算:
~f t1
~其中,f1是粘性阻尼函数,定义为:f13,并且。 3
C1
湍流粘性产生项,G用如下公式模拟:
~~
3-10 GCb1S~~
其中,SS22f2,而f21。其中,Cb1和k是常数,d是计算点
kd1f1
3
到壁面的距离;S
2ijij。ij定义为:
ujui1ij
2xixj
3-11
由于平均应变率对湍流产生也起到很大作用,FLUENT处理过程中,定义S为:
SijCprodmin(0,Sijij) 3-12
其中,Cprod2.0,ij
ijij,Sij2SijSij,平均应变率Sij定义为:
Sij
u1jui2xixj
3-13
在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这适合涡流靠近涡旋中心的区域,那里只有“单纯”的旋转,湍流受到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变,估计的涡旋粘性系数产生项偏高。
湍流粘性系数减少项Y为:
~
YCw1fw 3-14
d
2
6
1Cw3
其中,fwg6 3-15 6
gCw3
grCw2(r6r) 3-16
~
r22 3-17
Skd
1/6
~
其中,Cw1,Cw2,Cw3是常数,SS22f2。在上式中,包括了平均应变率
kd
~
对S的影响,因而也影响用S计算出来的r。
上面的模型常数在FLUENT中默认值为:Cb10.1335,Cb20.622,~2/3,
C17.1,Cw1Cb1/k2(1Cb2)/~,Cw20.3,Cw32.0,k0.41。
壁面条件
~设置为零。当计算网格足够细,可以计算层流底层时,壁面切在壁面,湍流运动粘性
应力用层流应力-应变关系求解,即:
uyu
3-18 u
如果网格粗错不能用来求解层流底层,则假设与壁面近邻的网格质心落在边界层的对数区,则根据壁面法则:
uyu1
lnE 3-19 uk
其中,k=0.419,E=9.793。
对流传热传质模型
在FLUENT中,用雷诺相似湍流输运的概念来模拟热输运过程。给出的能量方程为:
cptT
(E)[ui(Ep)]kxuj(ij)effSh 3-20 txixiPrti
式中,E是总能量,(ij)eff是偏应力张量,定义为:
(ij)effeff(
ujxi
uiu2
)effiij 3-21 xj3xi
其中,(ij)eff表示粘性加热,耦合求解。如果默认为分开求解,FLUENT不求解处,其默认值为0.85。 (ij)eff。但是可以通过变化“粘性模型”面板上的湍流普朗特数(Prt)
湍流质量输运与热输运类似,默认的Schmidt数是0.7,该值同样也可以在“粘性模型”
面板上调节。
标量的壁面处理与动量壁面处理类似,分别选用合适的壁面法则。
Spalart-Allmaras模型是相对简单的单方程模型,只需求解湍流粘性的输运方程,并不需要求解当地剪切层厚度的长度尺度。该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果,在透平机械湍流模拟方面也有较好结果。
Spalart-Allmaras模型的初始形式属于对低雷诺数湍流模型,这必须很好解决边界层的粘性影响区求解问题。在FLUENT中,当网格不是很细时,采用壁面函数来解决这一问题。当网格比较粗糙时,网格不满足精确的湍流计算要求,用壁面函数也许是最好的解决方案。另外,该模型中的输运变量在近壁处的梯度要比k中的小,这使得该模型对网格粗糙带来数值误差不太敏感。
但是,Spalart-Allmaras模型不能预测均匀各向同性湍流的耗散。并且,单方程模型没有考虑长度尺度的变化,这对一些流动尺度变换比较大的流动问题不太适合。比如,平板射流问题,从有壁面影响流动突然变化到自由剪切流,流场尺度变化明显。
3.3.2 标准k模型
标准k模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到,但耗散率方程是通过物理推理,数学上模拟相似原形方程得到的。该模型假设流动为完全湍流,分子粘性的影响可以忽略。因此,标准k模型只适合完全湍流的流动过程模拟。
标准k模型的湍动能k和耗散率ε方程为如下形式:
DkDtxiDDtxi
t
k
kxGkGbYM 3-22 i
t2xC1k(GkC3Gb)C2k 3-23
ki
在上述方程中,Gk表示由于平均速度梯度引起的湍动能产生,Gb是用于浮力影响引
起的湍动能产生;YM可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数
tC
k2
。
在FLUENT中,作为默认值常数,C1=1.44,C2=1.92,C0.09,湍动能k与耗散率ε的湍流普朗特数分别为k=1.0,=1.3。可以通过调节“粘性模型”面板来调节这些常数值。
3.3.3 重整化群κ-ε模型
重整化群κ-ε模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准κ-ε模型不同,而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准κ-ε模型有相似的形式:
Dkk
keffGkGbYM 3-24
DtxixiD2R 3-25 effC1(GkC3Gb)C2Dtxixikk
Gk表示由于平均速度梯度引起的湍动能产生,Gb是用于浮力影响引起的湍动能产生;YM可压速湍流脉动膨胀对总的耗散率的影响,这些参数与标准κ-ε模型中相同。k和分
别是湍动能k和耗散率的有效湍流普朗特数的倒数。 湍流粘性系数计算公式为:
~2k~ 3-26 1.72dd
31C
~/,C100 其中,
eff
对上面方程积分,可以精确得到有效雷诺数(涡旋尺度)对湍流输运的影响,这有助
于处理低雷诺数和近壁流动问题的模拟。
,C0.0845。这个结果非常有
意思,和标准κ-ε模型的半经验推导给出的常数C0.09非常近似。
在FLUENT中,如果是默认设置,用重整化群κ-ε模型时候是针对的高雷诺数流动
对于高雷诺数,上面方程可以给出:tC
k2
问题。如果对低雷诺数问题进行数值模拟,必须进行相应的设置。
重整化群κ-ε模型有旋修正
通常,平均运动有旋时候对湍流有重要影响。FLUENT中重整化群κ-ε模型通过修正湍流粘性系数来考虑了这类影响。
湍流粘性的修正形式为:
tt0f(s,,) 3-27
其中,t0是不考虑有旋计算出来的湍流粘性系数;Ω是FLUENT计算出来的特征旋流数;s
是旋流常数,不同值表示有旋流动的强度不同。流动可以是强旋或者中等旋度的。FLUENT默认设置s=0.05,针对中等旋度的流动问题,对于强旋流动,可以选择较大的值。
湍动能及其耗散率的有效湍流普朗特数倒数的计算公式为:
k
1.392901.0.6321
2.392902.0.3679
mol
3-28 eff
式中,0=1,在高雷诺数流动问题中,mol/eff1,k1.393。
湍流耗散率方程右边的R为:
k1
其中,Sk/,04.38,0.012。
R
C3(1/0)2
3
3-29
为了更清楚体现R对耗散率的影响,我们把耗散率输运方程重写为:
D22*
3-30 C2effC1(GkC3Gb)C2Dtxixikkk
C3(1/0)*
则:C2C2 3-31 3
1
**
在0的区域,R的贡献为正;C2大于C2。以对数区为例,3,C22.0,
这和标准κ-ε模型中给出的C2=1.92接近。因此,对于弱旋和中等旋度的流动问题,重整化群κ-ε模型给出的结果比标准κ-ε模型的结果要大。 重整化群模型中,C11.42,C21.68。
3.3.4 可实现κ-ε模型
可实现κ-ε模型的湍动能及其耗散率输运方程为:
DkDtxjDDtxj
t
ktt
k
xGkGbYM 3-32 j2CSCCC3Gb 3-33 121xkkj
其中,C1max0.43,
,Sk/ 5
FLUENT讲义:第三章 湍流模型
起的湍动能产生;YM可压速湍流脉动膨胀对总的耗散率的影响。C2和C1是常数;k,
在上述方程中,Gk表示由于平均速度梯度引起的湍动能产生,Gb是用于浮力影响引分别是湍动能及其耗散率的湍流普朗特数。在FLUENT中,作为默认值常数,C1=1.44,
C2=1.9,k=1.0,=1.2。
可实现κ-ε模型的湍动能的输运方程与标准κ-ε模型和重整化群κ-ε模型有相同的形式,只是模型参数不同。但耗散率方程有较大不同。首先耗散率产生项(方程右边第二项)不包含湍动能产生项Gk,现在的形式更能体现能量在谱空间的传输。另外的特色在于耗散率减少项中,不具有奇异性。并不象标准κ-ε模型模型那样把K放在分母上。
该模型适合的流动类型比较广泛,包括有旋均匀剪切流,自由流(射流和混合层),腔道流动和边界层流动。对以上流动过程模拟结果都比标准κ-ε模型的结果好,特别是可实现κ-ε模型对圆口射流和平板射流模拟中,能给出较好的射流扩张角。
湍流粘性系数公式为tC
k2
,这和标准κ-ε模型相同。不同的是,在可实现
κ-ε模型中,C不再是个常数,而是通过如下公式计算:
C
*
1A0As
UK
*
3-34
U其中,
~
SijSijijij,the mean ij是 is ijij2ijkk,ijijijkk,
rate-of –rotation tensor viewed in a rotating reference frame with the angular velocity k。模型常数A04.04,As
6cos,而:
arccos(6W),式中W=
1
3
SijSjkSkj~1ujui
SSS,,Sij() ijij2xixjS
我们可以发现,C是平均应变率与旋度的函数。在平衡边界层惯性底层,可以得到C
=0.09,与标准κ-ε模型中采用底常数一样。
双方程模型中,无论是标准κ-ε模型、重整化群κ-ε模型还是可实现κ-ε模型,三个模型有类似的形式,即都有κ和ε的输运方程,1,计算湍流粘性的方法不同;2,控制湍流扩散的湍流Prandtl数不同;3,εGk关系不同。但都包含了相同的表示由于平均速度梯度引起的湍动能产生Gk,用于浮力影响引起的湍动能产生Gb;可压速湍流脉动膨胀对总的耗散率的影响YM。
湍动能产生项
xiT
3-36 Gbgit
Prtxi
式中,Prt是能量的湍流普特朗数,对于可实现κ-ε模型,默认设置值为0.85。对于重整化群κ-ε模型,Prt1/,1/Prk/Cp。热膨胀系数于理想气体,浮力引起的湍动能产生项变为:
Gkuiuj
uj
3-35
1,对Tp
Gbgi
t
3-37
Prtxi
在FLUENT程序中,如果有重力作用,并且流场里有密度或者温度的梯度,浮力对湍动能的影响都是存在的。浮力对耗散率的影响不是很清楚,因此,默认设置中,耗散率方程中的浮力影响不被考虑。如果要考虑浮力对耗散率的影响,用“粘性模型”面板来控制。浮力对耗散率影响是用C3来体现。但C3并不是常数,而是如下的函数形式:
C3v
3-38 u
v是平行于重力方向的速度分量;u是垂直于重力方向的速度分量。如果流动速度与重力方向相同的剪切流动,C3=1,对于流动方向与重力方向垂直的剪切流,C3=0。
对于高马赫数的流动问题,可压速性对湍流影响在YM中体现。
YM2Mt2
其中,Mt是马赫数,定义为:Mt
k
(aRT是声速)。 2a
默认设置中,只要选择可压速理想气体,可压速效应都是考虑的。
在上述的双方程模型中,对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。能量方程形式为:
T
uiEpkeffuj(ij)eff(E)Sh 3-39 txixixi
式中,E是总的能量,keff是有效导热系数;(ij)eff是偏应力张量,定义为:
ujui2
effuiij 3-40 (ij)effeff
xxiixj3
(ij)eff表示的是粘性加热,耦合求解时总是计算。如果不是耦合求解时候,作为默认
设置,并不求解该量。如果有需要,需在“粘性模型”面板中设置。
对于重整化群κ-ε模型,有效导热系数为:
keffcpeff 3-41
α用(3-28)计算,式中,01/Prk/Cp。事实上,随着mol/eff的变化而变化,这是重整化群κ-ε模型的一个优点,因为实验中证明,湍流普朗特数随分子普朗特数及湍流而变化。
湍流质量输运处理过程与能量输运过程类似。对于标准κ-ε模型和可实现的κ-ε模型,默认的Schmidt数是0.7,重整化群模型中,是通过方程3-28来计算的,其中,01/Sc,Sc是分子Schimidt数。
3.3.5 雷诺应力模型(RSM)
雷诺应力模型是求解雷诺应力张量的各个分量的输运方程。具体形式为:
(uiuj)(Ukuiuj) txk
对流项Cij
uiujukp(kjuiikuj)xkxk
uiuj xk
TL
湍流扩撒项Dij 分子扩散Dij
UjUi
uuuugiujgjui jkikxxkk
应力产生项Pij 浮力产生项目Gij
uujuj
2ui pi
xxkxkjxi
压力应变项ij 耗散项ij
2kujumikmuiumjkm 3-42
系统旋转产生项Fij
LT
上面方程中,Cij,Dij,Pij,Fij不需要模拟,而Dij,Gij,ij,ij需要模拟以封
闭方程。下面简单对几个需要模拟项的模拟。
T
可以用Delay and Harlow [L38]的梯度扩散模型来模拟,即: Dij
DCs
xk
Tij
kukuluiuj 3-43
xl
但这个模型会导致数值不稳定,因此FLUENT程序中采用标量湍流扩散模型:
3-44
k2
式中,湍流粘性系数用tC来计算,根据Lien and Leschziner [L98],D
xk
Tij
tuiuj
kxk
k0.82,这和标准κ-ε模型中选取1.0有所不同。
根据Gibson and Launder [L58], Fu [L55], Launder [L88,L89], 压力应变项ij可以分解为三项,即:
w
3-45 ijij,1ij,2ij
ij,1,ij,2和wij分别是慢速项,快速项和壁面反射项。
2
,常数C11.8。 uukijijk3
12
ij,2C2PijFijGijCijijPGC,C20.60,PPkk,
23
11
GGkk,CCkk。
22
ij,1C1
壁面反射项用于重新分布近壁的雷诺正应力分布,主要是减少垂直于壁面的雷诺正应力,增加平行于壁面的雷诺正应力。该项模拟为:
w
ij
33k3/2ukumnknmijuiuknjnkujuknink C1
k22Cld
3/2
33k
km,2nknmijik,2njnkjk,2nink 3-46 C2
22Cdl
=0.5,C20.3,nk是xk在垂直于壁面方向上的单位分量,d是到壁面的式中,C1
距离;ClC
3/4
/k,C0.09,k=0.41。
wij
默认设置时候,FLUENT不计算
。如果需要计算时候,在“粘性模型”面板中设
置。
线性压力应变模型
对于小雷诺数流动,特别是用双层模型求解近壁流动问题时,FLUENT中通过改进模
和C2来改进压力应变项Launder [L91]。这一过程只有在选择双层流型常数C1,C2,C1
模型时候,在“粘性模型”面板上调节。
(0.0067Ret)2C112.58AA21e
C20.A
2
C1C11.67
3
12
C2maxC2,0 C2
其中,Retk2/(),参数A和张量不变量A2, A3定义为:
9A1(A2A3)
8A2ikki A3ikkjji
式中,ij是雷诺应力张量各向异性部分,定义为:
2
uiujkij
3ij
k
3-47
二阶压力应变模型
二阶压力应变模型由Spezible {L157}等人提出。
1*
ij(C1C1*P)bijC2(bikbkjbmnbmnij)(C3C3ijbij)kSij
3
2
C4k(bikSjkbjkSikbmnSmnij)C5k(bikjkbjkik) 3-48
3
式中,bij是雷诺应力各向异性张量,定义为:
3-49
u1ujui;ij1jui;模型常数平均应变率Sij定义为: Sij
22xixjxixj*
C13.4,C1*1.8,C24.2,C30.8,C31.3,C41.25,C50.4。
二阶压力应变模型不需要考虑壁面反射影响去模拟对数区湍流边界层过程。
2
uiujkij
3 bij
2k
FLUENT讲义:第三章 湍流模型
浮力对湍流的影响
浮力引起的产生项模拟为:
TT 3-50 Gijgigj
Prtxjxi
其中,Prt是能量的湍流普朗特数,默认设置值为0.85。
对于理想气体,把热膨胀系数的定义代入上式,得:
t
Gij
Prt
t
gi
3-51 gj
xjxi
耗散项ij的模拟
耗散张量ij模拟为:
2
3
式中,YM2Mt2,Mt是马赫数;标量耗散率用标准k-模型中的采用的耗散率输
ijij(YM) 3-52
运方程求解。
雷诺应力模型的边界条件
在流场进口,雷诺应力模型需要各个雷诺应力分量和湍动能耗散率的值。这些值可以直接输入,也可以湍流强度和特征长度来计算。
在壁面,雷诺应力模型通过壁面函数,给出各个雷诺应力分量和耗散率的值。
雷诺应力模型的能量与质量输运方程
在雷诺应力模型中,对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。能量方程形式为:
cptT
uiEp(E)(k)uj(ij)eff
txixiPrtxi
式中,E是总的能量;(ij)eff是偏应力张量,定义为:
Sh 3-53
(ij)eff
ujui
eff
x
ixj2
effuiij 3-54 3xi
(ij)eff表示的是粘性加热,耦合求解时总是计算。如果不是耦合求解时候,作为默认
设置,并不求解该量,并且Prt=0.85。如果有需要,需在“粘性模型”面板中设置。
3.3.6 大涡模拟(LES)
湍流中包含了不同时间与长度尺度的涡旋。最大长度尺度通常为平均流动的特征长度尺度。最小尺度为Komogrov尺度。
1,动量、能量、质量及其它标量主要由大涡输运;2,流动的几何和边界条件决定了大涡的特性,而流动特性主要在大涡中体现;3,小尺度涡旋受几何和边界条件影响较小,并且各向同性;大涡模拟过程中,直接求解大涡,小尺度涡旋模拟,从而使得网格要求比DNS低。
3.3.6.1大涡模拟的控制方程
LES的控制方程是对Navier-Stokes方程在波数空间或者物理空间进行过滤得到的。过滤的过程是去掉比过滤宽度或者给定物理宽度小的涡旋,从而得到大涡旋的控制方程。
过滤变量(上横线)定义为:
(x)D(x)G(x,x)dx 3-55
其中,D表示流体区域;G是决定涡旋大小的过滤函数。 在FLUENT中,有限控制体离散本身暗中包括了过滤运算,
(x)
1
V(x)dx, xV 3-56 V
其中V是计算控制体体积,过滤函数为:
xV1/V
3-57 G(x,x)
0xV
目前,大涡模拟对不可压流动问题得到较多应用,但在可压缩问题中的应用还很少,
因此这里涉及的理论都是针对不可压流动的大涡模拟方法。在FLUENT中,大涡模拟只能针对不可压流体(当然并非说是密度是常数)的流动。
过滤不可压的Navier-Stokes方程后,可以得到LES控制方程:
ui
0 3-58 txi
upij
3-59 (ui)(uiuj)(i)
txjxjxjxixj
其中,ij为亚网格应力,定义为:
ijuiujuiuj 3-60
很明显,上述方程与雷诺平均方程很相似,只不过大涡模拟中的变量是过滤过的量,而非时间平均量,并且湍流应力也不同。
3.3.6.2 亚网格模型
由于LES中亚网格应力项是未知的,并且需要模拟以封闭方程。目前,采用比较多的亚网格模型为涡旋粘性模型,形式为:
ijkkij2tSij 3-61
式中,t是亚网格湍流粘性系数;Sij是求解尺度下的应变率张量,定义为:
13
uuj1
Siji
2xjxi
3-62
求解亚网格湍流粘性系数t时,FLUENT提供了两种方法。第一,Smagorinsky-Lilly模型;第二,基于重整化群的亚网格模型。
最基本的亚网格模型是Smagorinsky [L145]最早提出的,Lilly [L99]把它进行了改善,这就是今天的Smagorinsky-Lilly模型。该模型的涡粘性计算方程为:
tL2sS 3-63
式中,Ls是亚网格的混合长度;S度Ls可以用下式计算。
2SijSij。Cs是Smagorinsky常数,则亚网格混合长
Lsmin(kd,CsV1/3) 3-64
其中,k=0.42,d是到最近壁面的距离,V是计算控制体体积。
Lilly通过对均匀各向同性湍流惯性子区湍流分析,得到了Cs=0.23。但是研究中发现,对于有平均剪切或者过渡流动中,该系数过高估计了大尺度涡旋的阻尼作用。因此,对于比较多的流动问题,Cs=0.1有比较好的模拟结果,该值是FLUENT的默认设置值。 我们再来看看基于重整化群思想的亚网格模型。人们用重整化群理论推导出了亚网格涡旋粘性系数[L182],该方法得到的是亚网格有效粘性系数,efft,而
eff
s2eff1HC3
1/32
1/3
3-65
式中,s(CrngV
)
2SijSij,H(x)是Heaviside函数,
x0x
3-66 H(x)
x00
V是计算控制体体积;重整化群常数Crng0.157,而常数C=100。
对于高雷诺数流动(t ),eff
t,基于重整化群理论的亚网格模型就与
Smagorinsky-Lilly模型相同,只是模型常数有区别。在流动场的低雷诺数区域,上面的函数就小于零,从而只有分子粘性起作用。所以,基于重整化群理论的亚网格模型对流动转捩和近壁流动问题有较好模拟效果。
3.3.6.3 大涡模拟的边界条件
对于给定进口速度边界条件,速度等于各个方向分量与随机脉动量的和,即:
uiuiIu
FLUENT讲义:第三章 湍流模型
如果网格足够密并可以求解层流底层的流动的话,壁面切应力采用线性应力应变关系,即:
uyu
3-67 u
如果网格不够细,则假定与壁面邻近网格质心落在边界层对数区内,则:
uyu1
lnE() 3-68 uk
其中,k=0.418,E=9.793。
表3-1 雷诺平均模型的比较
FLUENT讲义:第三章 湍流模型
第四节 湍流模型算例及其设置
湍流模型设置命令:Define-model-viscous
无粘,层流和湍流
湍流模型选项 近壁处理方法选择 附加湍流选项
算例一:有换热的腔道流动问题
adiabatic wall
cold air T = 0°F
constant temperature wall T = 100°F
10 ft
步骤:
1
, 检查是否湍流:ReDh5980
2, 选择低雷诺数湍流模型-RNG k模型;壁面处理用非平衡壁面函数,考虑压力
梯度影响 3, 网格划分:(1)四边形网格;(2)由于在靠近水平板处,垂直方向梯度较大,则近
壁网格加密,并保证第一个控制体在对数区内;(3)变化流向网格间距,用于捕捉边界层发展;(4),根据计算结果,自动调节网格,用于进一步计算温度梯度。
计算结果:
Velocity contours
BLs on upper & lower surfaces accelerate the core flow
Temperature contours
Important that thermal BL was accurately resolved as well
P
算例二:圆柱绕流
Compute drag coefficient of the cylinder
步骤:
1, 确定雷诺数,ReD24600
2, 钝体绕流,后面有不稳定的涡旋脱落。采用RNG k模型,壁面处理是双层区
模型; 3, 网格处理:近壁网格加密,由于是双层区模型,需要网格划分到y1
计算圆柱绕流的涡旋脱落过程