粘性流体力学第六章(6-2)

卡门理论似乎比普朗特理论完善一些，但比较繁
琐上，多此用P外r当andtldd理x2u2论21 。0 时， t12 为奇点 ，所以工程
17
试验表明, 在湍流边界层据壁面的某个范围内，
速度与距离x2按对数关系变化，即u1 lnx2,根据卡门 的混合长度计算公式（6-47)可得
lm =k
du1 dx2
2
l12
d 2 u1 dx22
2
2
(6－45）
l1 是附加混合长度，不过由试验确定 l1 很复杂。 d、混合长度理论没有考虑压力脉动对动量传递的影响，
而压力脉动可以跨越lm而传递的，考虑到压力脉动的影 响，流体微团的动量不可能在lm范围内保持不变。
12
Karman相似理论：Karman（1930）提出了一种湍流 局部相似性假设。他认为在自由湍流场中各空间点的湍 流脉动具有几何相似性，也就是说，各点的湍流脉动对 同一个时间尺度和长度尺度只有比例系数的差别，因此 只要用一个时间和速度比尺就能确定湍动结构。对于二
到达点 x2 时，恰巧在 x2 l 微团的左边时，就会产
生碰撞，而产生横向运动 u1 ，这样 u2 ~ u1 。同样, 当向两中个间微 补团 充到 也达 会产x2生点u2时。向相反运动时，周围的微团会
8
图6－7 u2的产生
认为
u
2
常数 l du1 dx2
12
u1u2
lm 2
du1 dx2
du1 dx2
对三维流动
u1u2
t
u1 x2
t
u1 x2
uiu
j
t
Dij
t
ui x j
u j xi
2 t Sij
(6－34） (6－35）
3
应力张量表示
ij pij t Dij
(6－36）
其中： t湍流动力粘性系数，或称涡旋粘性系数，t t 是湍流运动粘性系数。
上述假设所以能提出是基于对湍流脉动引起的动量
23
tr 1 exp G
x1 x1tr
x1 x1tr
1 Ue
dx1
其中 xtr为转捩开始的位置，G为经验系数：
G
1 1200
U
3 e
Re 1.34 xtr
Rextr U e x1tr
BL模型（1978）同样适用于湍流边界层。其主要 特点为：采用分区的涡粘公式；用涡量取代变形率； 对混合长度做近壁修正。具体的表达式为：
(6－39）
9
t
lm2
du1 dx2
(6－40）
根据实验研究可以得到以下几点： a、由试验得到的 lm ，不象假设的那样为流体微团的
尺寸，而是与流动的平均尺度一样的量级。 b 、lm 不 是 空 间 常 数 。 在 边 界 层 中 根 据 尼 古 拉 兹 和
Klebanoff试验，在内层（壁面区）
lm
i
kx2
(6－41）
t
i
kx2v*
式中i表示内层，k=0.40~0.41，v*
阻力。
w
,
w是壁面摩擦
10
在边界层的外层（核心区）：
lm o 1 t o 2v*
(6－42）
式中o表示外层，1 0.075 ~ 0.09 ， 2 0.06 ~ 0.075 ，
为间隙因子。
Van Driest（1959）年提出在内层：
是壁面处的流体运动粘性系数，v* w 是壁面的磨
阻速度。
25
外层的涡粘性系数公式为：
式中，
t o ＝CFwakeFkleb (x2 )
(6－53）
Fwake＝min
x2
max
Fmax
,
Cwk
x2
maxU
2 dif
/ Fmax
Fwake为尾流函数，Fmax和x2max分别为F (x2 ) x2[1 exp
交换与气体分子运动引起的粘性切应力进行简单的类比
的结果。对于 一般在定温下可认为是常数，但 t不
是常量，因为湍流的动量交换取决于湍流的平均运动。
流动只在一个方向上有明确的速度梯度时，可以认
为 t是个标量。在一般情况下，当i=j时
uiui
2k
2 t
ui xi
(6－37）
4
式中k为湍动能(k
1 2
uiui )
lm
i
k
x2
1
exp
x2
25.3v
(6－43）
11
式中k 0.435 ，在外层：
lm o 0.09
(6－44）
c、根据上述公式在管流中 u1
u1 max
时，u1
x2
0
，那么，
t12 0 ，而事实上 t12 0 ，为此Prandtl提出修正：
1
t12
lm
2
d u1 dx1
19
(3)、零方程模型
在上述理论的基础上，一些学者提出了湍流粘 性系数的代数型模型，也称为零方程模型:
Cebci－Smith（1968）（CS）模型， Mellor－Herring（1968）（MH）， Patanka－Spalding（1968）（PS）和 Baldwin－Lomax（BL）等模型。 这些模型的共同点是根据湍流边界层的结构， 对 t 在边界层的内层和外层须用不同的尺度。
速度与混合长度lm的乘积成正比
t u1lm
混合长度lm类比于气体分子运动的自由行程，在lm一段
特征长度之内湍流微团保持自己的动量不变。
6
图6－6 混合长度与脉 动速度
对于图中二维不可压定常湍流：
ux2 l ux2 ux2 l
假设微团从x2 l 或 x2 l运动至 x2，对于 x2 来讲，
脉动速度
u
2
0
或
u
2
0
，
7
u1
x2
l
u1
x2
l
d u1(x2 ) dx2
x x2
... ...
u1
u1 x2
l
u1 x2
l
d u1 dx2
由于u1 ~ u1, 所以微团从 x2 l 运动到 x2 时 u1 0,
u2还 可0，以也认就为是u说2 u~1u和1，u这2 是是异因号为的当。x2 l 处的微团
(6－47b）
18
其中A为衰减长度因子，定义为
A 26 w / 1 2
w为壁面剪切应力。式(6-47b)表明, 当x2很小时，粘性
作用很大；而当x2增大后，粘性作用逐渐消失。将该式中 的指数函数用泰勒级数展开后容易看出，当x2 0，lm x22。所以式(6-47b)综合了lm x22和lm x2两个区域混合 长度的变化，已称为现在许多实用的零方程的基础。
x2 x20
同时
u
2
u1
，并令 l1 kl2
，那么
l1
d u1 dx2
l12 k
d 2 u1 dx22
lm l1 k
du1
dx2 d 2u1 dx22
(6－47）
16
所以雷诺切应力为：
du1
3
t12 u1u2
式中k=0.40。
k2
dx2
d
2u1
dx22
2
du1 dx2
(6－48）
CS模型发展了Van Priest的模型，得到广泛的 应用，其公式为：
20
对于内层：
t
l2
u1 x2
rtr
式中 0 x2 x2c
l kx2 1 exp x2 A
k=0.40,A是衰减因子。
A A v*1
N
1
N
p v
1 exp 11.8v
exp 11.8v
第六章 湍流基本理论
第一节 第二节 第三节 第四节 第五节 第六节 第七节
湍流的基本特征和统计平均方法 湍流连续方程和雷诺方程 湍流能量方程 雷诺平均统计模式 湍流的相关函数和谱分析 拟序结构 湍流大涡数值模拟
1
第四节 雷诺平均统计模式
在雷诺方程中的不封闭量是雷诺应力，因此统计模式的 目标是封闭雷诺平均方程，建立足够的雷诺应力方程组（代 数的、微分的或一般泛函形式的）使得平均运动方程可解。
ui 0 xi
那么 uiui 0
，如果 t 是一个标量，
，而实际上
uiui
2k 3
。
为此Boussinesq修正（6－35），提出对于三维湍流
uiu
j
2 3
kij
t
Dij
(6－38）
Townsend.A.A测得在圆柱尾迹的充分湍流区，
t
0.0164U
0
d
,
(U
是来流速度，d是圆柱直径）。
0
Hinze J.O.在空气的圆截面射流中测得t 0.0116U pd
代数涡粘模式的最大缺点是它的局部性，代数表达 式中雷诺应力之核当地的平均变形率有关。代数模式 完全忽略了湍流通计量之间关系的历史效应，而历史 效应很难做局部的修正，因此发展包含历史效应的模
式是必要的，常用的 k 模式包含部分历史效应，称
为目前工程湍流计算的主要封闭模式。
24
t
t t
i o
, ,
内层的涡粘公式为：
x2 x2c x2 x2c
(6－51）
t i ＝l2
(6－52）
ii 1 2 (i ijk uk x j )是当地的涡量绝对值；l是
考虑壁面修正的混合长度：
l＝kx2 1 exp
x2 / A
式中，k＝0.4,为卡门常数，A＝26，x2＝v*x2 / w。 w