湍流理论学习
1.层流和湍流
粘性流体的运动存在着两种完全不同的流动状态:层流状态和湍流状态。雷诺首先于1883年通过做圆管内流动实验观察到层流与湍流现象。当圆管中流动速度较小时,管中的流线之间层次分明,互不掺混,这样的流动称为层流。当流速增大后,流体作复杂、无规律、
'随机的不定常运动,称为湍流。流动状态与雷诺数R e 、下临界雷诺数R ec 和上临界雷诺数R ec
''时,有关。当R e ≤R ec 时,流动为层流;当R ec ≤R e ≤R ec 流动为不稳定过渡状态;当R e >R ec
时,流动为湍流。
湍流是在连续介质范畴内流体的不规则运动,它有别于物质分子的不规则运动。具体来说,在极不规则的湍流中,流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度。因此湍流运动产生的质量和能量的输运将远远大于分子热运动产生的宏观输运。 2.湍流的平均化、雷诺粘性应力
经典的湍流理论认为,湍流是一种完全不规则的随机运动,湍流场中的物理量在时间和空间上呈随机分布,不同的瞬时有不同的值,关注某个瞬时的值是没有意义的。因此,雷诺首创用统计平均方法来描述湍流的随机运动,即对各瞬时量进行平均得到有意义的平均值。
从N-S 方程出发,利用平均化运算的法则推导平均物理量满足的方程组。只考虑不可压缩流体情形,假设体力可以忽略,此时,N-S 方程具有下列形式
∂u ∂u ∂u 1∂p ⎧∂u +u +v +w =-+υ∆u ⎪∂t ∂x ∂y ∂z ρ∂x ⎪
∂v ∂v ∂v 1∂p ⎪∂v +u +v +w =-+υ∆v ⎪∂t ∂x ∂y ∂z ρ∂y ⎪
(1) ⎨
∂w ∂w ∂w ∂w 1∂p ⎪+u +v +w =-+υ∆w ⎪∂t ∂x ∂y ∂z ρ∂z ⎪
∂u ∂v ∂w ⎪++=0
⎪∂x ∂y ∂z ⎩
运用(1)式中的连续性方程,运动方程可改写为
⎧∂u ∂u 2∂uv ∂uw 1∂p
+++=-+υ∆u ⎪∂t ∂x ∂y ∂z ρ∂x ⎪
⎪∂v ∂uv ∂v 2∂vw 1∂p
+++=-+υ∆v ⎪
∂y ∂z ρ∂y ⎪∂t ∂x
(2) ⎨2
⎪∂w +∂uw +∂vw +∂w =-1∂p +υ∆w ⎪∂t ∂x ∂y ∂z ρ∂z ⎪
∂u ∂v ∂w ⎪++=0⎪∂x ∂y ∂z ⎩
对方程组(2)中各式两边进行平均化运算,并利用平均化运算法则得到
____
⎧22
'∂∂∂∂∂u ∂∂'1∂⎪++++++=-+υ∆⎪∂t ∂x ∂y ∂z ∂x ∂y ∂z ρ∂x ⎪____
⎪∂∂∂2∂∂∂v '2∂'1∂⎪++++++=-+υ∆⎪∂t ∂x ∂y ∂z ∂x ∂y ∂z ρ∂y (3) ⎨
____
⎪22
⎪∂+∂+∂+∂+∂'+∂'+∂w '=-1∂+υ∆⎪∂t ∂x ∂y ∂z ∂x ∂y ∂z ρ∂z ⎪
∂∂∂⎪++=0
⎪∂x ∂y ∂z ⎩
考虑到方程组(3)的第四式,方程组(3)中的头三个方程可改写成另一种形式,把脉动项
移到右边,得到
____
⎧⎛2⎫'∂-ρu ⎪ ⎪
∂(-ρ)∂(-ρ')⎛⎫∂∂∂∂∂⎝⎭⎪ρ+++⎪=-+μ∆+++
⎪ ∂t ∂x ∂y ∂z ∂x ∂x ∂y ∂z ⎝⎭⎪____
⎛⎪2⎫'∂-ρv ⎪⎪⎛∂∂(-ρ)∂(-ρ')⎫∂∂∂∂⎝⎭⎪+++⎪=-+μ∆+++⎪ρ
∂t ∂x ∂y ∂z ∂y ∂x ∂y ∂z ⎨⎝⎭⎪____
⎛⎫⎪∂ -ρw '2⎪
⎪⎛∂∂(-ρ')∂(-ρ')∂∂∂⎫∂⎭+++++⎝⎪ρ ⎪=-+μ∆+
∂x ∂y ∂z ⎭∂z ∂x ∂y ∂z ⎪⎝∂t
⎪∂∂∂⎪++=0
∂x ∂y ∂z ⎪⎩
(4)
将上式和应力形式的运动方程对比
ρ
其中是应力张量,则有
d V
=div dt
=-p I +2μ+P '
其中I 是单位张量。是平均运动变形速度张量,P '是对称的二阶张量,它的表达式
如下
___________⎤⎡2-ρu '-ρu 'v '-ρu 'w '⎥'⎤⎢τxz
___________⎥⎥⎢2
''-ρv '⎢τ'-ρv 'w '⎥ (5) yz ⎥=-ρu v ⎢___⎥________'⎥τzz
⎦⎢-ρu 'w '-ρv 'w '-ρw '2⎥⎢⎥⎣⎦
'τxy '⎡τxx
⎢P '=⎢τ'τ'yx yy
⎢τzx '⎣'τzy
由(5)式得到在湍流运动中除了平均运动的粘性应力外,还多了一项由于脉动所引起的
应力,称为雷诺粘性应力。平均流动元所满足的方程(4)称为雷诺方程。上述雷诺方程是不封闭的,方程个数只有4个,而未知函数有10个,即、及六个雷诺粘性应力分量。为了使方程组封闭,必须在雷诺应力即平均速度之间建立补充关系式。 3.湍流模型
所谓湍流模型理论,就是根据理论和经验,对雷诺平均运动方程和雷诺应力项建立表达
式或方程,然后对雷诺应力方程的某些项提出尽可能合理的模型和假设,使方程组封闭求解的理论。
3.1 零方程模型
所谓零方程模型,就是在雷诺方程和连续方程以外,不需要另外再加任何方程式来使方程组封闭。即雷诺粘性应力能直接用某些物理量和物理常数表达出来,没有引进高阶微分方程,又称代数模型。
(1) Boussinesq涡粘性系数法
一般的方法利用Boussinesq 假设把雷诺压力和平均速度梯度联系起来:
⎛∂u i ∂u j ⎫2⎛∂u ⎫
-ρu i 'u 'j =μt +- ρk +μt i ⎪δij (6) ⎪ ∂x ⎪∂x i ⎭⎝j ∂x i ⎭3⎝
_____
(2) Cebeci-Smith模型
该模型通常写成涡粘性系数形式,在内层有
⎡y ⎤∂U
εm =(κy ) 2⎢1-exp(-) ⎥() (7)
l D ⎦∂y ⎣
式中l D 是衰减长度系数,当壁面不存在质量引射时为常数。在外层有
2
εm =0.0168U δ*γ (8)
式中δ为边界层位移厚度,U 是外部自由流速度,γ为间隙因子。
3.2 一方程模型
(1) The Spalart-Allmaras模型
Spalart-Allmaras 模型的变量中是湍流动粘滞率除了近壁区域,方程是:
*
∂∂1(ρ) +(ρi ) =G v +∂t ∂x i σ⎧⎪∂
⎨⎪⎩∂x j ⎡∂⎤∂2⎫⎪
⎢(μ+ρ) ⎥+C b 2ρ() ⎬-Y v +S (9)
∂x j ⎥∂x j ⎪⎢⎣⎦⎭
这里G v 是湍流粘度生成的,Y v 是被湍流粘度消去,发生在近壁区域。
对于解决动力漩涡粘性,S-A 模型是相对简单的方程,不必要去计算和剪应力层厚度相
关的长度尺度。
(2) 简化的二方程模型
由于ε方程难模化和精度低,可以在保留k 方程的基础上舍弃ε方程,将ε表示为
k 3/2ε= (9)
l
式中,l 可视为混合长度,必须根据具体流场情况另外给出。 c 二方程模型
(1) 标准k -ε模型
标准k -ε模型是个半经验公式,主要是基于湍流动能和扩散率。k 方程是个精确方程,ε方程是个由经验公式导出的方程。
湍流动能方程k 方程为
∂∂∂(ρk ) +(ρku i ) =∂t ∂x i ∂x j
扩散方程ε方程为
⎡μt ∂k ⎤
⎢(μ+) ⎥+G k +G b -ρε-Y M +S k (10)
σ∂x ⎢k j ⎥⎣⎦
⎡μt ∂ε⎤εε2
⎢(μ+) ⎥+C 1ε(G k +C 3εG b ) -C 2ερ+S ε(11)
σ∂x k k ⎢εj ⎥⎣⎦
方程中G k 表示由层流速度梯度而产生的湍流动能,G b 是由浮力产生的湍流动能,Y M 由∂∂∂(ρε) +(ρεu i ) =∂t ∂x i ∂x j
于在可压缩湍流中,过渡的扩散产生的波动。 (2) 标准k -ω模型
标准k -ω模型是一种经验模型,是基于湍流能量方程和扩散速率方程。
⎛∂k ⎫Γ+G k -Y k +S k (12) k ∂x ⎪⎪
j ⎭⎝
∂∂∂⎛∂ω⎫(ρω) +(ρωu i ) =Γω+G ω-Y ω+S ω (13) ⎪ ⎪∂t ∂x i ∂x j ⎝∂x j ⎭
在方程中,G k 是由层流速度梯度而产生的湍流动能。G ω是由ω方程产生的,Γk 和Γω表明了k 和ω的扩散率。Y k 和Y ω由于扩散产生的湍流。
∂∂∂(ρk ) +(ρku i ) =∂t ∂x i ∂x j
4.雷诺应力模式
雷诺应力模式简称为RSM ,由于雷诺平均运动方程所包含的雷诺应力为一未知项,所以最直接的办法是建立雷诺应力v i 'v 'j 的方程。 (1) 雷诺应力方程和湍流动能方程
得到v i 'v 'j 的方程步骤为:把N-S 方程化成
_____
_____
D (v i v j ) Dt
的形式,将瞬时速度表示成平均速
度与脉动速度之和,即v i =i +v i '、v j =j +v '将其代入j ,平均运动的形式
D (v i v j ) Dt
形式方程,减去
D (v i v j ) Dt
D (v i v j ) Dt
_____
,最终得到雷诺应力v i 'v 'j 方程为
_____
D (v i 'v 'j ) Dt
_____
=
∂v i 'v 'j ∂t
_____
+k
∂v i 'v 'j ∂x k
_____
⎛__________⎫_____
∂v 'j p '_______∂v i 'v 'j ⎪∂v i 'p '∂
'-υ=-+δik +v i 'v 'j v k δjk ⎪
∂x k ρρ∂x k ⎪
⎪
湍流扩散项分子扩散项⎝⎭
__⎫__⎛_____⎛⎫_____''∂v ∂v ∂v '''∂v ∂v ∂v p j j j ⎪''i ⎪-2υi - v i 'v k +v 'j v k + i +
∂x k ∂x k ⎪∂x k ∂x k ρ ⎪ ∂x j ∂x i ⎪⎪ ⎝⎭⎝⎭ 耗散项
产生项
压力变形项
__
__
_____________________
____
(14)
在模式理论中要用到湍动能k =v '/2,在式(14)中令i=j得到
___
2i
Dk ∂k ∂k
=+j Dt ∂t ∂x j
=-
∂
∂x j
⎛1⎫1
'''''''''''∂v p +v v v -2υv s i i j k i ij ⎪-v i v k s ij -2υs ij s ij
2⎝ρ⎭
_____
_______
_____
_____
_____
(15)
这就是湍动能方程。
(2) 经模化后的雷诺应力方程和湍动能方程
'、v i 'p '等,因此,需要在雷诺应力方程(14)和湍动能方程(15)中存在未知项,如v i 'v 'j v k
建立这些未知量的表达式,这一过程称为湍流模式建立的过程。最终经模化后的雷诺应力方
程和湍动能方程为:
____________
__________
⎛⎫2
D v i 'v 'j ∂v i 'v 'j ⎪∂ k ∂v i 'v 'j 2
=C k +υ+P ij -δij εDt ∂x l ε∂x l ∂x l ⎪3 ⎪ (16) ⎝⎭
ε⎛_____22⎫⎛⎫-C 1 v i 'v 'j -δij k ⎪-C 2 P ij -δij P k ⎪
k ⎝33⎭⎝⎭
_____
Dk ∂⎛k 2∂k ∂k ⎫=C +υ k ⎪+P k -ε (17) Dt ∂x l ⎝ε∂x l ∂x l ⎭
其中
⎛P ij =- '∂v __
_______'j
_____⎫__'∂v 'i ⎪ v i 'v k +v 'j v , P _____∂v 'k =-v i 'v l 'i ⎝∂x k k ∂x k ⎪⎪⎭
∂x l 式中经验常数C k =0.09~0.11, C 1=1.5~2.2, C 2=0.4~0.5。
以上方程中出项的湍流耗散率ε仍是未知项,建立ε方程,进行模化得到
D ε2
__
=∂⎛_____
Dt ∂x ⎝C k ∂ε∂ε⎫ε∂v i 'ε2
εε∂x +υ⎪-C ε1v i 'v l '
-C ε2 l l ∂x l ⎭k ∂x l k
式中经验常数C ε=0.07~0.09, C ε1=1.41~1.45,C ε2=1.9~1.92。
联立封闭微分方程组,这就是雷诺应力封闭模式。
(18) (19)
湍流理论学习
1.层流和湍流
粘性流体的运动存在着两种完全不同的流动状态:层流状态和湍流状态。雷诺首先于1883年通过做圆管内流动实验观察到层流与湍流现象。当圆管中流动速度较小时,管中的流线之间层次分明,互不掺混,这样的流动称为层流。当流速增大后,流体作复杂、无规律、
'随机的不定常运动,称为湍流。流动状态与雷诺数R e 、下临界雷诺数R ec 和上临界雷诺数R ec
''时,有关。当R e ≤R ec 时,流动为层流;当R ec ≤R e ≤R ec 流动为不稳定过渡状态;当R e >R ec
时,流动为湍流。
湍流是在连续介质范畴内流体的不规则运动,它有别于物质分子的不规则运动。具体来说,在极不规则的湍流中,流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度。因此湍流运动产生的质量和能量的输运将远远大于分子热运动产生的宏观输运。 2.湍流的平均化、雷诺粘性应力
经典的湍流理论认为,湍流是一种完全不规则的随机运动,湍流场中的物理量在时间和空间上呈随机分布,不同的瞬时有不同的值,关注某个瞬时的值是没有意义的。因此,雷诺首创用统计平均方法来描述湍流的随机运动,即对各瞬时量进行平均得到有意义的平均值。
从N-S 方程出发,利用平均化运算的法则推导平均物理量满足的方程组。只考虑不可压缩流体情形,假设体力可以忽略,此时,N-S 方程具有下列形式
∂u ∂u ∂u 1∂p ⎧∂u +u +v +w =-+υ∆u ⎪∂t ∂x ∂y ∂z ρ∂x ⎪
∂v ∂v ∂v 1∂p ⎪∂v +u +v +w =-+υ∆v ⎪∂t ∂x ∂y ∂z ρ∂y ⎪
(1) ⎨
∂w ∂w ∂w ∂w 1∂p ⎪+u +v +w =-+υ∆w ⎪∂t ∂x ∂y ∂z ρ∂z ⎪
∂u ∂v ∂w ⎪++=0
⎪∂x ∂y ∂z ⎩
运用(1)式中的连续性方程,运动方程可改写为
⎧∂u ∂u 2∂uv ∂uw 1∂p
+++=-+υ∆u ⎪∂t ∂x ∂y ∂z ρ∂x ⎪
⎪∂v ∂uv ∂v 2∂vw 1∂p
+++=-+υ∆v ⎪
∂y ∂z ρ∂y ⎪∂t ∂x
(2) ⎨2
⎪∂w +∂uw +∂vw +∂w =-1∂p +υ∆w ⎪∂t ∂x ∂y ∂z ρ∂z ⎪
∂u ∂v ∂w ⎪++=0⎪∂x ∂y ∂z ⎩
对方程组(2)中各式两边进行平均化运算,并利用平均化运算法则得到
____
⎧22
'∂∂∂∂∂u ∂∂'1∂⎪++++++=-+υ∆⎪∂t ∂x ∂y ∂z ∂x ∂y ∂z ρ∂x ⎪____
⎪∂∂∂2∂∂∂v '2∂'1∂⎪++++++=-+υ∆⎪∂t ∂x ∂y ∂z ∂x ∂y ∂z ρ∂y (3) ⎨
____
⎪22
⎪∂+∂+∂+∂+∂'+∂'+∂w '=-1∂+υ∆⎪∂t ∂x ∂y ∂z ∂x ∂y ∂z ρ∂z ⎪
∂∂∂⎪++=0
⎪∂x ∂y ∂z ⎩
考虑到方程组(3)的第四式,方程组(3)中的头三个方程可改写成另一种形式,把脉动项
移到右边,得到
____
⎧⎛2⎫'∂-ρu ⎪ ⎪
∂(-ρ)∂(-ρ')⎛⎫∂∂∂∂∂⎝⎭⎪ρ+++⎪=-+μ∆+++
⎪ ∂t ∂x ∂y ∂z ∂x ∂x ∂y ∂z ⎝⎭⎪____
⎛⎪2⎫'∂-ρv ⎪⎪⎛∂∂(-ρ)∂(-ρ')⎫∂∂∂∂⎝⎭⎪+++⎪=-+μ∆+++⎪ρ
∂t ∂x ∂y ∂z ∂y ∂x ∂y ∂z ⎨⎝⎭⎪____
⎛⎫⎪∂ -ρw '2⎪
⎪⎛∂∂(-ρ')∂(-ρ')∂∂∂⎫∂⎭+++++⎝⎪ρ ⎪=-+μ∆+
∂x ∂y ∂z ⎭∂z ∂x ∂y ∂z ⎪⎝∂t
⎪∂∂∂⎪++=0
∂x ∂y ∂z ⎪⎩
(4)
将上式和应力形式的运动方程对比
ρ
其中是应力张量,则有
d V
=div dt
=-p I +2μ+P '
其中I 是单位张量。是平均运动变形速度张量,P '是对称的二阶张量,它的表达式
如下
___________⎤⎡2-ρu '-ρu 'v '-ρu 'w '⎥'⎤⎢τxz
___________⎥⎥⎢2
''-ρv '⎢τ'-ρv 'w '⎥ (5) yz ⎥=-ρu v ⎢___⎥________'⎥τzz
⎦⎢-ρu 'w '-ρv 'w '-ρw '2⎥⎢⎥⎣⎦
'τxy '⎡τxx
⎢P '=⎢τ'τ'yx yy
⎢τzx '⎣'τzy
由(5)式得到在湍流运动中除了平均运动的粘性应力外,还多了一项由于脉动所引起的
应力,称为雷诺粘性应力。平均流动元所满足的方程(4)称为雷诺方程。上述雷诺方程是不封闭的,方程个数只有4个,而未知函数有10个,即、及六个雷诺粘性应力分量。为了使方程组封闭,必须在雷诺应力即平均速度之间建立补充关系式。 3.湍流模型
所谓湍流模型理论,就是根据理论和经验,对雷诺平均运动方程和雷诺应力项建立表达
式或方程,然后对雷诺应力方程的某些项提出尽可能合理的模型和假设,使方程组封闭求解的理论。
3.1 零方程模型
所谓零方程模型,就是在雷诺方程和连续方程以外,不需要另外再加任何方程式来使方程组封闭。即雷诺粘性应力能直接用某些物理量和物理常数表达出来,没有引进高阶微分方程,又称代数模型。
(1) Boussinesq涡粘性系数法
一般的方法利用Boussinesq 假设把雷诺压力和平均速度梯度联系起来:
⎛∂u i ∂u j ⎫2⎛∂u ⎫
-ρu i 'u 'j =μt +- ρk +μt i ⎪δij (6) ⎪ ∂x ⎪∂x i ⎭⎝j ∂x i ⎭3⎝
_____
(2) Cebeci-Smith模型
该模型通常写成涡粘性系数形式,在内层有
⎡y ⎤∂U
εm =(κy ) 2⎢1-exp(-) ⎥() (7)
l D ⎦∂y ⎣
式中l D 是衰减长度系数,当壁面不存在质量引射时为常数。在外层有
2
εm =0.0168U δ*γ (8)
式中δ为边界层位移厚度,U 是外部自由流速度,γ为间隙因子。
3.2 一方程模型
(1) The Spalart-Allmaras模型
Spalart-Allmaras 模型的变量中是湍流动粘滞率除了近壁区域,方程是:
*
∂∂1(ρ) +(ρi ) =G v +∂t ∂x i σ⎧⎪∂
⎨⎪⎩∂x j ⎡∂⎤∂2⎫⎪
⎢(μ+ρ) ⎥+C b 2ρ() ⎬-Y v +S (9)
∂x j ⎥∂x j ⎪⎢⎣⎦⎭
这里G v 是湍流粘度生成的,Y v 是被湍流粘度消去,发生在近壁区域。
对于解决动力漩涡粘性,S-A 模型是相对简单的方程,不必要去计算和剪应力层厚度相
关的长度尺度。
(2) 简化的二方程模型
由于ε方程难模化和精度低,可以在保留k 方程的基础上舍弃ε方程,将ε表示为
k 3/2ε= (9)
l
式中,l 可视为混合长度,必须根据具体流场情况另外给出。 c 二方程模型
(1) 标准k -ε模型
标准k -ε模型是个半经验公式,主要是基于湍流动能和扩散率。k 方程是个精确方程,ε方程是个由经验公式导出的方程。
湍流动能方程k 方程为
∂∂∂(ρk ) +(ρku i ) =∂t ∂x i ∂x j
扩散方程ε方程为
⎡μt ∂k ⎤
⎢(μ+) ⎥+G k +G b -ρε-Y M +S k (10)
σ∂x ⎢k j ⎥⎣⎦
⎡μt ∂ε⎤εε2
⎢(μ+) ⎥+C 1ε(G k +C 3εG b ) -C 2ερ+S ε(11)
σ∂x k k ⎢εj ⎥⎣⎦
方程中G k 表示由层流速度梯度而产生的湍流动能,G b 是由浮力产生的湍流动能,Y M 由∂∂∂(ρε) +(ρεu i ) =∂t ∂x i ∂x j
于在可压缩湍流中,过渡的扩散产生的波动。 (2) 标准k -ω模型
标准k -ω模型是一种经验模型,是基于湍流能量方程和扩散速率方程。
⎛∂k ⎫Γ+G k -Y k +S k (12) k ∂x ⎪⎪
j ⎭⎝
∂∂∂⎛∂ω⎫(ρω) +(ρωu i ) =Γω+G ω-Y ω+S ω (13) ⎪ ⎪∂t ∂x i ∂x j ⎝∂x j ⎭
在方程中,G k 是由层流速度梯度而产生的湍流动能。G ω是由ω方程产生的,Γk 和Γω表明了k 和ω的扩散率。Y k 和Y ω由于扩散产生的湍流。
∂∂∂(ρk ) +(ρku i ) =∂t ∂x i ∂x j
4.雷诺应力模式
雷诺应力模式简称为RSM ,由于雷诺平均运动方程所包含的雷诺应力为一未知项,所以最直接的办法是建立雷诺应力v i 'v 'j 的方程。 (1) 雷诺应力方程和湍流动能方程
得到v i 'v 'j 的方程步骤为:把N-S 方程化成
_____
_____
D (v i v j ) Dt
的形式,将瞬时速度表示成平均速
度与脉动速度之和,即v i =i +v i '、v j =j +v '将其代入j ,平均运动的形式
D (v i v j ) Dt
形式方程,减去
D (v i v j ) Dt
D (v i v j ) Dt
_____
,最终得到雷诺应力v i 'v 'j 方程为
_____
D (v i 'v 'j ) Dt
_____
=
∂v i 'v 'j ∂t
_____
+k
∂v i 'v 'j ∂x k
_____
⎛__________⎫_____
∂v 'j p '_______∂v i 'v 'j ⎪∂v i 'p '∂
'-υ=-+δik +v i 'v 'j v k δjk ⎪
∂x k ρρ∂x k ⎪
⎪
湍流扩散项分子扩散项⎝⎭
__⎫__⎛_____⎛⎫_____''∂v ∂v ∂v '''∂v ∂v ∂v p j j j ⎪''i ⎪-2υi - v i 'v k +v 'j v k + i +
∂x k ∂x k ⎪∂x k ∂x k ρ ⎪ ∂x j ∂x i ⎪⎪ ⎝⎭⎝⎭ 耗散项
产生项
压力变形项
__
__
_____________________
____
(14)
在模式理论中要用到湍动能k =v '/2,在式(14)中令i=j得到
___
2i
Dk ∂k ∂k
=+j Dt ∂t ∂x j
=-
∂
∂x j
⎛1⎫1
'''''''''''∂v p +v v v -2υv s i i j k i ij ⎪-v i v k s ij -2υs ij s ij
2⎝ρ⎭
_____
_______
_____
_____
_____
(15)
这就是湍动能方程。
(2) 经模化后的雷诺应力方程和湍动能方程
'、v i 'p '等,因此,需要在雷诺应力方程(14)和湍动能方程(15)中存在未知项,如v i 'v 'j v k
建立这些未知量的表达式,这一过程称为湍流模式建立的过程。最终经模化后的雷诺应力方
程和湍动能方程为:
____________
__________
⎛⎫2
D v i 'v 'j ∂v i 'v 'j ⎪∂ k ∂v i 'v 'j 2
=C k +υ+P ij -δij εDt ∂x l ε∂x l ∂x l ⎪3 ⎪ (16) ⎝⎭
ε⎛_____22⎫⎛⎫-C 1 v i 'v 'j -δij k ⎪-C 2 P ij -δij P k ⎪
k ⎝33⎭⎝⎭
_____
Dk ∂⎛k 2∂k ∂k ⎫=C +υ k ⎪+P k -ε (17) Dt ∂x l ⎝ε∂x l ∂x l ⎭
其中
⎛P ij =- '∂v __
_______'j
_____⎫__'∂v 'i ⎪ v i 'v k +v 'j v , P _____∂v 'k =-v i 'v l 'i ⎝∂x k k ∂x k ⎪⎪⎭
∂x l 式中经验常数C k =0.09~0.11, C 1=1.5~2.2, C 2=0.4~0.5。
以上方程中出项的湍流耗散率ε仍是未知项,建立ε方程,进行模化得到
D ε2
__
=∂⎛_____
Dt ∂x ⎝C k ∂ε∂ε⎫ε∂v i 'ε2
εε∂x +υ⎪-C ε1v i 'v l '
-C ε2 l l ∂x l ⎭k ∂x l k
式中经验常数C ε=0.07~0.09, C ε1=1.41~1.45,C ε2=1.9~1.92。
联立封闭微分方程组,这就是雷诺应力封闭模式。
(18) (19)