第3章fluent湍流模型1.docx
- 文档编号:16004346
- 上传时间:2023-07-09
- 格式:DOCX
- 页数:23
- 大小:209.57KB
第3章fluent湍流模型1.docx
《第3章fluent湍流模型1.docx》由会员分享,可在线阅读,更多相关《第3章fluent湍流模型1.docx(23页珍藏版)》请在冰点文库上搜索。
第3章fluent湍流模型1
第三章,湍流模型
第一节,前言
湍流流动模型很多,但大致可以归纳为以下三类:
第一类是湍流输运系数模型,是Boussinesq于1877年针对二维流动提出的,将速度脉动
的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。
即:
Ui
U1U2t3—1
X2
推广到三维问题,若用笛卡儿张量表示,即有:
UiUj2(…
UiUjtkj3—2
xjxi3
模型的任务就是给出计算湍流粘性系数t的方法。
根据建立模型所需要的微分方程的数
目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。
前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。
大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes方程,得
到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。
选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT提供的湍流模型包括:
单方程(Spalart-Allmaras)模型、双方程模型(标准k-
&模型、重整化群k-&模型、可实现(Realizable)k-&模型)及雷诺应力模型和大涡模拟。
DirectNumericalSimulation
湍流模型种类示意图
第二节,平均量输运方程
雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。
对于速
度,有:
UiUiUi3-3
其中,U?
和Ui分别是平均速度和脉动速度(i=1,2,3)类似地,对于压力等其它标量,我们也有:
3-4
其中,表示标量,如压力、能量、组分浓度等。
把上面的表达式代入瞬时的连续与动量方程,并取平均(去掉平均速度Ui上的横线),我
们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式:
(Ui)03-5
tXi
Dui
P
UiUj2
Ui
ij
UiUj
Xj
3-6
Dt
XiXj
Xjxi3
X|
上面两个方程称为雷诺平均的
Navier-Stokes
(RANS)
方程。
他们和瞬时
Navier-Stokes方
程有相同的形式,只是速度或其它求解变量变成了时间平均量。
额外多出来的项UiUj是雷
诺应力,表示湍流的影响。
如果要求解该方程,必须模拟该项以封闭方程。
如果密度是变化的流动过程如燃烧问题,我们可以用法夫雷(Favre)平均。
这样才可以
求解有密度变化的流动问题。
法夫雷平均就是出了压力和密度本身以外,所有变量都用密度加权平均。
变量的密度加权平均定义为:
〜
/3-7
符号〜表示密度加权平均;对应于密度加权平均值的脉动值用表示,即有:
〜。
很显然,这种脉动值的简单平均值不为零,但它的密度加权平均值等于零,即:
—o,—0
Boussinesq近似与雷诺应力输运模型
为了封闭方程,必须对额外项雷诺应力UiUj进行模拟。
一个通常的方法是应用
Boussinesq假设,认为雷诺应力与平均速度梯度成正比,即:
Ui
Uj
2
Ui
UiUj
t
(k
ti)
Xj
Xi
3
Xi
BoUssinesq假设被用于Spalart-Allmaras单方程模型和k双方程模型。
Boussinesq近似的好处是与求解湍流粘性系数有关的计算时间比较少,例如在Spalart-Allmaras单方程模型中,
只多求解一个表示湍流粘性的输运方程;在k双方程模型中,只需多求解湍动能k和耗散
率&两个方程,湍流粘性系数用湍动能k和耗散率&的函数。
Boussinesq假设的缺点是认为湍
流粘性系数t是各向同性标量,对一些复杂流动该条件并不是严格成立,所以具有其应用限制性。
另外的方法是求解雷诺应力各分量的输运方程。
这也需要额外再求解一个标量方程,通常是耗散率£方程。
这就意味着对于二维湍流流动问题,需要多求解4个输运方程,而三维湍流
问题需要多求解7个方程,需要比较多的计算时间,对计算机内存也有更高要求。
在许多问题中,Boussinesq近似方法可以得到比较好的结果,并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。
但是,如果湍流场各向异性很明显,如强旋流动以及应力驱动的二次流等流动中,求解雷诺应力分量输运方程无疑可以得到更好的结果。
第三节,湍流模型
3.3.1单方程(Spalart-Allmaras)模型
Spalart-Allmaras模型的求解变量是~,表征出了近壁(粘性影响)区域以外的湍流运动粘性系数。
~的输运方程为:
其中,f1是粘性阻尼函数,定义为:
湍流粘性产生项,G用如下公式模拟:
3—10
GCb1S
其中,SSf2,而f21-
kd1
壁面的距离;S,2_j_。
j定义为:
1uj5
ij2XiXj
由于平均应变率对湍流产生也起到很大作用,
SijCprodmin(0,Sjj)
其中,
C
prod
2.0
ij
Jijij,
Sij
Sij
1
uj
ui
2
Xi
Xj
—。
其中,Cb1和k是常数,d是计算点到
'1
3—11
FLUENT处理过程中,定义S为:
3—12
2SijSj,平均应变率Sj定义为:
3—13
在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。
这适合涡流靠近涡旋中心的
区域,那里只有“单纯”的旋转,湍流受到抑止。
包含应变张量的影响更能体现旋转对湍流的影响。
忽略了平均应变,估计的涡旋粘性系数产生项偏高。
湍流粘性系数减少项Y为:
其中,Cw1,Cw2,Cw3是常数,SSf2。
在上式中,包括了平均应变率对S
kd
的影响,因而也影响用S计算出来的r。
上面的模型常数在FLUENT中默认值为:
Cb10.1335,Cb20.622,-2/3,
2
C17.1,Cw1Cb1/k(1Cb2)/~,Cw20.3,Cw32・°,k0.41。
壁面条件
在壁面,湍流运动粘性~设置为零。
当计算网格足够细,可以计算层流底层时,壁面切应
力用层流应力一应变关系求解,即:
uuy
3—18
u
如果网格粗错不能用来求解层流底层,则假设与壁面近邻的网格质心落在边界层的对数
区,则根据壁面法则:
3—19
其中,k=0.419,E=9.793。
对流传热传质模型
在FLUENT中,用雷诺相似湍流输运的概念来模拟热输运过程。
给出的能量方程为:
CptT
(E)——[w(Ep)]—kUj(j洽Sh3—20
tXiXiPrtXi
式中,E是总能量,(Jeff是偏应力张量,定义为:
其中,(ij)eff表示粘性加热,耦合求解。
如果默认为分开求解,FLUENT不求解处(ij)eff。
但是可以通过变化“粘性模型”面板上的湍流普朗特数(Prt),其默认值为0.85。
湍流质量输运与热输运类似,默认的Schmidt数是0.7,该值同样也可以在“粘性模型”
面板上调节。
标量的壁面处理与动量壁面处理类似,分别选用合适的壁面法则。
综上所述,Spalart-Allmaras模型是相对简单的单方程模型,只需求解湍流粘性的输运方程,并不需要求解当地剪切层厚度的长度尺度。
该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果,在透平机械湍流模拟方面也有较好结果。
Spalart-Allmaras模型的初始形式属于对低雷诺数湍流模型,这必须很好解决边界层的粘性
影响区求解冋题。
在FLUENT中,当网格不是很细时,米用壁面函数来解决这一冋题。
当网
格比较粗糙时,网格不满足精确的湍流计算要求,用壁面函数也许是最好的解决方案。
另外,该模型中的输运变量在近壁处的梯度要比k中的小,这使得该模型对网格粗糙带来数值误
差不太敏感。
但是,Spalart-Allmaras模型不能预测均匀各向同性湍流的耗散。
并且,单方程模型没有考虑长度尺度的变化,这对一些流动尺度变换比较大的流动问题不太适合。
比如,平板射流问题,
从有壁面影响流动突然变化到自由剪切流,流场尺度变化明显。
3.3.2标准k模型
标准k模型需要求解湍动能及其耗散率方程。
湍动能输运方程是通过精确的方程推导
得到,但耗散率方程是通过物理推理,数学上模拟相似原形方程得到的。
该模型假设流动为完全湍流,分子粘性的影响可以忽略。
因此,标准k模型只适合完全湍流的流动过程模拟。
标准k模型的湍动能k和耗散率&方程为如下形式:
Dk
tk
kXi
GkGbYm
3—22
Dt
Xi
D
t
2
C1「(GkC3Gb)C2
3—23
Dt
Xi
kXi
k
k
在上述方程中,
Gk表示由于平均速度梯度引起的湍动能产生,
Gb是用于浮力影响引起
的湍动能产生;YM可压速湍流脉动膨胀对总的耗散率的影响。
湍流粘性系数
t
cX
在FLUENT中,作为默认值常数,C1=1.44,C2=1.92,C0.09,湍动能k与耗散
率&的湍流普朗特数分别为k=1.0,=1.3。
可以通过调节“粘性模型”面板来调节这些
常数值。
重整化群k-&模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模
型。
模型中的常数与标准k-&模型不同,而且方程中也出现了新的函数或者项。
其湍动能与
耗散率方程与标准k
-£模型有相似的形式:
Dk
k
keff
GkGb
Ym
3—24
Dtxi
Xi
D
2
eff
C1(Gk
C3Gb)C2R
3—25
DtXi
Xi
k
k
Gk表示由于平均速度梯度引起的湍动能产生,
Gb是用于浮力影响引起的湍动能产生;
Ym可
压速湍流脉动膨胀对总的耗散率的影响,这些参数与标准K-£模型中相同。
k和分别是
湍动能k和耗散率的有效湍流普朗特数的倒数。
湍流粘性系数计算公式为:
2k~
d——1.72d~3—26
厂厂31C
其中,~eff/,C100
对上面方程积分,可以精确得到有效雷诺数(涡旋尺度)对湍流输运的影响,这有助于处理低雷诺数和近壁流动问题的模拟。
k2
对于高雷诺数,上面方程可以给出:
tC,C0.0845。
这个结果非常有意
思,和标准k-£模型的半经验推导给出的常数C0.09非常近似。
在FLUENT中,如果是默认设置,用重整化群k-&模型时候是针对的高雷诺数流动问题。
如果对低雷诺数问题进行数值模拟,必须进行相应的设置。
重整化群k-&模型有旋修正
通常,平均运动有旋时候对湍流有重要影响。
FLUENT中重整化群k-&模型通过修正湍流
粘性系数来考虑了这类影响。
湍流粘性的修正形式为:
k
tt°f(s,,—)3—27
其中,t0是不考虑有旋计算出来的湍流粘性系数;Q是FLUENT计算出来的特征旋流数;s是
旋流常数,不同值表示有旋流动的强度不同。
流动可以是强旋或者中等旋度的。
FLUENT默认设
置s=0.05,针对中等旋度的流动问题,对于强旋流动,可以选择较大的值。
湍动能及其耗散率的有效湍流普朗特数倒数的计算公式为:
1.3929
0.6321
2.3929
01.3929
02.3929
0.3679
mol
3-28
式中,0=1,在高雷诺数流动问题中,
湍流耗散率方程右边的R为:
C3(1/0)2
k
04.38,
R对耗散率的影响,
1
其中,Sk/
为了更清楚体现
D
Dt
Xi
则:
C2
C2
eff
Xi
3(1/
3~
Ci
0)
eff
/1
moleffk
1.393。
0.012。
我们把耗散率输运方程重写为:
2
C2
-29
3-30
-31
*
3,C2
C2大于C2。
以对数区为例,
接近。
因此,对于弱旋和中等旋度的流动问题,重整化群K-&模型给出的结果比标准K-&模型的结果要大。
重整化群模型中,C11.42,C21.68。
在
和标准K
0的区域,R的贡献为正;
-£模型中给出的c2=1.92
2.0,这
3.3.4可实现k-&模型
Dk
tk
GkGb
Ym
Dt
Xj
kXj
D
t
2
C1S
C2
C1C3G
Dt
Xj
tXj
k'一
k
可实现K-&模型的湍动能及其耗散率输运方程为:
3-32
3-33
其中,C1
max0.43,
Sk/
在上述方程中,Gk表示由于平均速度梯度引起的湍动能产生,
Gb是用于浮力影响引起
的湍动能产生;Ym可压速湍流脉动膨胀对总的耗散率的影响。
C2和G是常数;k,分
别是湍动能及其耗散率的湍流普朗特数。
在FLUENT中,作为默认值常数,
G=1.44,C2=1.9,
k=1.0,=1.2。
可实现K-&模型的湍动能的输运方程与标准K-&模型和重整化群K-&模型
有相同的形式,只是模型参数不同。
但耗散率方程有较大不同。
首先耗散率产生项(方
程右边第二项)不包含湍动能产生项Gk,现在的形式更能体现能量在谱空间的传输。
另外的特色在于耗散率减少项中,不具有奇异性。
并不象标准K-£模型模型那样把K放在分母上。
该模型适合的流动类型比较广泛,包括有旋均匀剪切流,自由流(射流和混合层),腔道流动和边界层流动。
对以上流动过程模拟结果都比标准K-&模型的结果好,特别是可实现K-£模型对圆口射流和平板射流模拟中,能给出较好的射流扩张角。
-34
双方程模型中,无论是标准K-£模型、重整化群K-£模型还是可实现K-£模型,三个
模型有类似的形式,即都有尺和£的输运方程,它们的区别在于:
1,计算湍流粘性的方法不
同;2,控制湍流扩散的湍流Prandtl数不同;3,£方程中的产生项和G关系不同。
但都包含了相同的表示由于平均速度梯度引起的湍动能产生Gk,用于浮力影响引起的湍动能产生Gb;
可压速湍流脉动膨胀对总的耗散率的影响Ym。
湍动能产生项
在FLUENT程序中,如果有重力作用,并且流场里有密度或者温度的梯度,浮力对湍动能的影响都是存在的。
浮力对耗散率的影响不是很清楚,因此,默认设置中,耗散率方程中的浮力影响不被考虑。
如果要考虑浮力对耗散率的影响,用“粘性模型”面板来控制。
浮力对耗散率影响是用C3来体现。
但C3并不是常数,而是如下的函数形式:
v
C3tanh—3-38
u
v是平行于重力方向的速度分量;u是垂直于重力方向的速度分量。
如果流动速度与重力
方向相同的剪切流动,C3=1,对于流动方向与重力方向垂直的剪切流,C3=0。
对于高马赫数的流动问题,可压速性对湍流影响在ym中体现。
2
Ym2Mt
其中,Mt是马赫数,定义为:
Mt;(a.RT是声速)。
默认设置中,只要选择可压速理想气体,可压速效应都是考虑的。
[(E)匚uEp
—keff
xXi
Uj(ij)eff
Sh
3-39
在上述的双方程模型中,对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。
能量方程形式为:
式中,E是总的能量,keff是有效导热系数;(ij)eff是偏应力张量,定义为:
(、ujUi2Ui
(ij)effeff二effij3—40
xixj3xi
(ij)eff表示的是粘性加热,耦合求解时总是计算。
如果不是耦合求解时候,作为默认设置,
并不求解该量。
如果有需要,需在“粘性模型”面板中设置。
对于重整化群K-&模型,有效导热系数为:
keffcpeff3—41
%用(3—28)计算,式中,01/Prk/Cp。
事实上,随着mol/eff的变化而变化,这是重整化群K-&模型的一个优点,因为实验中证明,湍流普朗特数随分子普朗特数及湍流而变化。
湍流质量输运处理过程与能量输运过程类似。
对于标准K-&模型和可实现的K-£模型,
默认的Schmidt数是0.7,重整化群模型中,是通过方程3—28来计算的,其中,01/Sc,
Sc是分子Schimidt数。
335雷诺应力模型(RSM雷诺应力模型是求解雷诺应力张量的各个分量的输运方程。
具体形式为:
(UiUj)——(UkUiUj)tXk
对流项Cij
这和标准K-&模型中选取1.0有所不同。
根据GibsonandLaUnder[L58],Fu[L55],LaUnder[L88,L89],项,即:
w
ijij,1ij,2ij
压力应变项ij可以分解为三
3—45
tkuM5uj
DijCs
XkXl
但这个模型会导致数值不稳定,因此
上面方程中,Cij,Dj,Pij,Fij不需要模拟,而dT,Gij,ij,ij需要模拟以封闭方程。
下面简单对几个需要模拟项的模拟。
DT可以用DelayandHarlow[L38]的梯度扩散模型来模拟,即:
3—43
FLUENT程序中采用标量湍流扩散模型:
壁面反射项用于重新分布近壁的雷诺正应力分布,主要是减少垂直于壁面的雷诺正应力,增加平行于壁面的雷诺正应力。
该项模拟为:
线性压力应变模型
对于小雷诺数流动,特别是用双层模型求解近壁流动问题时,FLUENT中通过改进模型常
数G,C2,C1和C2来改进压力应变项Launder[L91]。
这一过程只有在选择双层流模型时候,在“粘性模型”面板上调节。
C112.58A「A1e(00067Ret)2
C20.75、A
2
C1C11.67
C2
max
tC2
C2
0
3
其中,Retk2/(),参数a和张量不变量A2,A3定义为:
9
A1-(A2A3)
8
A2ikki
ikkjji
式中,j是雷诺应力张量各向异性部分,定义为:
ij
2.
UiUjkj
3
k
3—47
二阶压力应变模型
二阶压力应变模型由Spezible{L157}等人提出。
*1
(C1C1P)bjC2(bikbkj3bmnbmn
ij
j)
(C3
C3Jbqbj)
kSj
C4
k(bikSjkbSmnij)C5k(bik
jk
bjk
ik)
3—48
式中,
bj
bj是雷诺应力各向异性张量,定义为:
——2
UiUjkj
3
2k
3—49
平均应变率Sj定义为:
Sj
Uj
Ui
Xi
G3.4,C,1.8,C24.2,
Uj
Ui
模型常数
ij
Xj
*
C3
Xi
Xj
1.3
0.8,
二阶压力应变模型不需要考虑壁面反射影响去模拟对数区湍流边界层过程。
C3
C41.25,
C5
0.4。
浮力对湍流的影响
浮力引起的产生项模拟为:
Gj
Prt
T
gi
Xj
Tgj
Xi
3—50
其中,Prt是能量的湍流普朗特数,默认设置值为0.85。
对于理想气体,把热膨胀系数的定义代入上式,得:
Gij
3—51
耗散项j的模拟
耗散张量ij模拟为:
2
ij3ij(Ym)3-52
式中,Ym2Mt2,Mt是马赫数;标量耗散率用标准k-模型中的采用的耗散率输运方
程求解。
雷诺应力模型的边界条件
这些值可以直接
在流场进口,雷诺应力模型需要各个雷诺应力分量和湍动能耗散率的值。
输入,也可以湍流强度和特征长度来计算。
在壁面,雷诺应力模型通过壁面函数,给出各个雷诺应力分量和耗散率的值。
耦合求解时总是计算。
如果不是耦合求解时候,
作为默认设置,
(ij)eff表示的是粘性加热,
雷诺应力模型的能量与质量输运方程
在雷诺应力模型中,
对流传热传质模型都是通过雷诺相似湍流动量输运方程得到的。
能量
方程形式为:
/E)Ui
tx
Ep(k
Xi
cpt)Prt)
TS
Uj(ij)effSh
Xi
3—53
式中,E是总的能量;(
ij)eff是偏应力张量,
定义为:
并不求解该量,并且Prt=0.85。
如果有需要,需在“粘性模型”面板中设置。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- fluent 湍流 模型