第3章-均匀各项同性湍流
第七章第三节均匀各向同性湍流
二、两点脉动速度相关矩
以Bij表示M点脉动速度和M’点脉动速度间 的相关矩,即
Bij = u i u ′j
Bi , j (r ) = [ B LL (r ) − B NN (r )]
BLL = u f (r )
2
ξ iξ j
r
2
+ B NN (r )δ ij
若以f(r)和g(r)分布表示纵向和横向相关系数,则
∂ri ∂r ri ∂ri = = δ ij =3 ∂ri r ∂ri ∂r j
u 2 r df r df 2 Bi , j (r ) = 2 [− ri r j + ( f + )r δ ij ] 2 dr 2 dr r
r1,2,3为r在(x1,x2,x3)在直角坐标系的分量
9
三、两点速度三元相关矩
3
Keller与Friedman的思路
Keller与Friedman的思路
Reynolds平均 一点平均
N-S方程
Reynolds方程
− (u i′u ′j )
方程不封闭
怎么办?
R=0,τ=0 两点相关矩 v v v v v Bij ( x , t ; r ,τ ) = u i′ ( x , t )u ′j ( x + r , t + τ )
B p, j v r = B1 (r ) r
B p , j = pu ′j = B1 (r )
rj r
= pu ′ ρ
rj r
18
各向同性湍流场 压力速度相关矩
由各向同性的定义
⎧ puθ = − puθ ′ ⎪ ′ ⎨ ′ ⎪ puϕ = − puϕ ⎩ ′
⎧ p ′uθ = − p ′uθ ⎪ ⎨ ⎪ p ′uϕ = − p ′uϕ ⎩
03-均匀各向同性湍流
质疑的核心:湍动能耗散率不是一个常数
浙江大学航空航天学院流体工程研究所
24
局部各向同性湍流的结构函数(5)
标度律
Kolmogorov的修正
层次结构标度律-佘振苏
IFE , Zhejiang University
浙江大学航空航天学院流体工程研究所
25
局部各向同性湍流的结构函数(6)
随机函数的增量
大尺度的脉动相互抵消
脉动速度增量的统计矩
IFE , Zhejiang University
二阶 结构 函数
三阶 函数
浙江大学航空航天学院流体工程研究所
18
局部各向同性湍流的结构函数(补充)
IFE , Zhejiang University
浙江大学航空航天学院流体工程研究所
19
不可压缩均匀各向同性湍 流
− − − − 动力学方程 卡门-霍华斯方程 卡门-霍华斯方程应用 能量传输链
IFE , Zhejiang University
浙江大学航空航天学院流体工程研究所
2
均匀各向同性的相关函数谱张量(1)
均匀各向同性湍流
均匀湍流
均匀各向同性湍流
IFE , Zhejiang University
浙江大学航空航天学院流体工程研究所
14
不可压缩均匀各向同性湍流(8)
能量传输链
能量谱方程
记
IFE , Zhejiang University
浙江大学航空航天学院流体工程研究所
15
不可压缩均匀各向同性湍流(9)
能量传输链
特征尺度
惯性区
IFE , Zhejiang University
各向同性湍流Karman-Howarth方程的精确解
a. for all of these four solutions: a1 > 0,
b. only for the secondary solution:σ > 0 ; c. only for the third solution: 0 < σ < 5 .
2
The use of these above exact solutions would be presented in our further publications.
2a1
2
when
k
=1−σ
,
f
(ξ
)
=
− a1 ξ 2
e4
F⎜⎛σ ,2,
a1
ξ
2
⎟⎞
⎝ 4⎠
when
σ
=
1,
f
(ξ
)
=
− a1 ξ 2
e4
F ⎜⎛
1
,1,
a1
ξ
2
⎟⎞
⎝2 4 ⎠
The detailed analysis could be seen in appendix 1. The convergent conditions at infinite must be satisfied
Key words: isotropic turbulence, KH equations, exact solution
Introduction
The ideas of similarity and self-preservation have played an important role in the development of turbulence theory for more than a half-century. The traditional approach to search for similarity solutions in turbulence has been to assume the existence of a single length and velocity scale, then ask whether and what conditions the averaged equations of motion admit to such solutions. In the development of this ideas, Karman (von Karman,1938) firstly gave the analysis solutions of Karman-Howarth equation for the decaying turbulence. Sedov (Sedov,1944,1951) gave anther method to obtain the special solutions of KH equations. In this paper, we found that a complete set of exact solutions exists. Following the methods adopted by Sedov, new solutions have been given, both for three-dimensional and two-dimensional isotropic turbulence.
第三章_湍流模型
第三章 湍流模型第一节 前言湍流流动模型很多,但大致可以归纳为以下三类:第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。
即:2121x u u u t ∂∂=''-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有:ij ijj i t j i k x u xu u u δρμρ32-⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。
根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。
前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。
大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。
选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
湍流模型种类示意图Direct Numerical Simulation包含更多 物理机理每次迭代 计算量增加提的模型选RANS-based models第二节 平均量输运方程雷诺平均就是把Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分。
对于速度,有:i i i u u u '+= 3-3其中,i u 和i u '分别是平均速度和脉动速度(i=1,2,3)类似地,对于压力等其它标量,我们也有:φφφ'+= 3-4 其中,φ表示标量,如压力、能量、组分浓度等。
湍流的数学模型简介精心整理版共88页
一般认为,无论湍流流动多么复杂,非稳态的连续性方 程和N-S方程(动量方程)仍然适用于湍流的瞬时流动。
第1章 湍流导论
1.3、湍流的基本方程(不可压) ❖ N-S方程
ui ui ui'
将非稳态N-S方程对时间作平均,即把湍流的运动看成是时间平均
流动与瞬间脉动流动的叠加:
'
及 t的概念,直接建立以雷诺应力为因变量的微分方程,然
后作适当假设使之封闭。这种模型也称为二阶封闭模型。
代数应力方程模型(Algebraic Stress Model,ASM)
主要思想是设法将应力的微分方程简化为代数表达式, 以减少RSM模型过分复杂的弱点,同时保留湍流各项异性 的基本特点。
3.2 湍流模型具体介绍
第2章 湍流的数值模拟方法简介
2.2 模型比较
湍流模型方法 (RANS方法)
大涡模拟方法 (LES方法)
给出了时间平均的流动信息,易于工程应用
抹去了流动的瞬态特性及细观结构,适合高雷 诺数,不具普适性
介于RANS与DNS之间,非常成功的应用于RANS
不能满足要求的高端应用,如燃烧、混合、外部空 气动力学。
、 k-g 模型等 。其中,应用最普遍的是 k-ε模型。
针对k-ε模型不足,许多学者对标准的模型进行了修正。
▪ 重整化群k-ε模型(renormalization group,RNG model) ▪ 可实现k-ε模型(realizable k-ε model) ▪ 多尺度k-ε模型(multiscale model of turbulence)
Contents
1
湍流导论
2
湍流的数学模型简介
3
第3章 均匀各项同性湍流
上面 称为洛强斯基不变量,它表示各向同性湍流场 衰减过程中某种统计特征量的守恒性。 3.均匀各向同性湍流的后期衰变 各向同性湍流的衰减后期,湍流脉动量很小,这时3阶相关 和2阶相关相比是高阶小量,因而可略去不计,于是不可压缩 各向同性湍流的后期衰变方程简化为:
3.4 不可压缩均匀各向同性湍流动力学的若干性质
3.4 不可压缩均匀各向同性湍流动力学的若干性质
3.4.1 不可压缩均匀湍流的 阶速度相关动力学方程和谱张量动力学方程 不可压缩均匀湍流的2阶速度相关动力学方程和谱张量动力学方程 动力学方程
谱张量动力学方程
3.4 不可压缩均匀各向同性湍流动力学的若干性质
3.4.2 不可压缩均匀各向同性湍流的 不可压缩均匀各向同性湍流的Karman-Howarth方程 方程
3.3 不可压缩均匀各向同性湍流的动力学方程
3.3.4 均匀湍流中的湍动能传输链
在粘性作用下,脉动速度逐渐衰减,而且小尺度的成分衰减 得最快,于是在耗散过程中大尺度脉动成分占更多份额。 由于惯性在速度脉动的各个尺度间进行动量输运,它将大尺 度脉动的动能传输给小尺度脉动。 于是在粘性和惯性的联合作用下,湍流脉动场形成一种能量 传输链:大尺度湍流脉动通过惯性作用向小尺度湍流脉动不断输 送能量,这股能量在小尺度湍流脉动中耗散殆尽。 在能量传输过程中,压强在各个脉动分量间起调节作用,如 果物理空间中初始脉动场的动能在各个分量间分配不均匀,压强 梯度将使它们逐渐均分。
3.5 不可压缩均匀各向同性湍流中的湍动能传输
1 含能波数和含能尺度 谱最大值的波数定义为含能波数 kin ,它的倒数定义为含能尺度L,即
可以估计含能尺度的量级等于
在含能尺度范围内(又称含能区),湍动能通过惯性传输能量, 而湍动能耗散几乎可以忽略,也就是说,含能尺度范围内,惯性主 宰湍流运动,因此含能尺度范围又称惯性区。
[理学]第七章第四节局地均匀各向同性湍流
六、Ko1mogorov的修正部分
几十年来,Ko1mogorov理论已获得很大成 功,在大气中,在海洋中,在实验室中,人们 几乎到处都观测到了这个理论所预测的2/3湍 谱。它不仅有重要的理论意义,是湍流发展史 中一个具有里程碑意义的成就,而且在许多工 程应用中有直接应用价值。
人们在大量实验观测中看到了Ko1mogorov 理论所预测的结果。
υ3
(
) −α / 4
1 − 3α
=υ2 4 ε
r 1 +α
24
α
ε
α=2/3
Dll (r) = Cε 2 / 3r 2 / 3 Dnn (r) = C ′ε 2 / 3r 2 / 3
粘性子区湍流速度结构函数
对于粘性子区,r<<l0,将β(r/l0)展开得
β ( r ) = β (0) +
l0
∂β
有
Dpp (r) ~ ρ υx1 εx2 x3
x1=?,x2=?,x3=?
压力结构函数Dpp(r)
根据量纲理论,得到
M 2 L−2T −4 = ( ML −3 ) x1 ( L2T −1 ) x2 ( L2T −3 ) x3
⎧ ⎪⎨−
2
=
2 = x1 −3x1 + 2x2
+
2x3
⎪⎩ − 4 = −x2 − 3x3
Dln n (r)
=
(r
∂ ∂r
+ 1)Blll
(r)
结构函数的微分方程
从 Karman-Howarth 方 程 出 发 可 以 得 到 相 应的支配结构函数的微分方程
(d dr
+
4 r )[Dlll
内燃机缸内的流动
∫ lI′ =
∞
g ( x)dx
0
(3-1b)
它相应于普朗特理论中的混合长度,(定义:流体微团从一层跳入另一层,经过一段不与 其它流体微团相碰撞的距离),可以证明
lI′ = 0.5lI
(3-1c)
与此相似,可定义湍流积分时间尺度
∫ τ I
=
∞
f (t)dt
0
(3-2)
其中f(t)是同一空间点(x0),不同时间脉动速度的欧拉时间自相关系数
空间相关系数能较好地反映涡团的平均尺度。于是可引入湍流长度积分尺度或简称湍流尺度:
∞
∫ lI = 0 f (x)dx
(3-1a)
-3-
第 3 章 缸内气流运动 式中,f(x)为湍流纵向自相关系数,其定义为
f (x) = u(x0 )u(x0 + x) u2 (x0 ) u2 (x0 + x)
其中, x0 和 x0 + x 分别为 A、B 两点的坐标,u(x)是与两点连线平行的脉动速度分量(图
宋金瓯 内燃机中的流体运动
第三章 缸内气流运动
在内燃机整个工作循环中,缸内气体充量始终进行着极其复杂而又强烈瞬变的湍流流动。 这种湍流运动决定了各种量在缸内的输运及其空间分布,它对可燃混和气形成、火焰传播、 燃烧品质、缸壁传热及污染物形成等都具有直接的、本质的影响。组织良好的缸内空气运动 可以提高汽油机的火焰传播速率、降低燃烧循环变动、适应稀燃和层燃;同样可以提高柴油 机的燃油空气混合速率,提高燃烧速率,促进燃烧过程中空气与未燃燃料的混合(热混合)。 但是,内燃机缸内流动极其复杂,它受到进气状态、工况和燃烧室结构等多种因素制约,因 而不存在对各种发动机都通用的流动规律,甚至不同研究者所得结果不乏相互矛盾之处,这 更加表明了深入研究缸内流动的必要性。
第三章-湍流模型
第三章,湍流模型第一节, 前言湍流流动模型很多,但大致可以归纳为以下三类:第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。
即:2121x u u u t ∂∂=''-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有:ij ijj i t j i k x u xu u u δρμρ32-⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。
根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。
前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。
大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。
选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
湍流模型种类示意图第二节,平均量输运方程包含更多 物理机理每次迭代 计算量增加提的模型选RANS-based models雷诺平均就是把Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分。
对于速度,有:i i i u u u '+= 3-3其中,i u 和i u '分别是平均速度和脉动速度(i=1,2,3)类似地,对于压力等其它标量,我们也有:φφφ'+= 3-4 其中,φ表示标量,如压力、能量、组分浓度等。
化工流体力学 第三章
u x u y u x u z xy yx xz zx y z x x
法向应力与应变率的关系如下:
u y u z yz zy z y
概
适用,对于这些问题,必须考虑流体的黏性。
述
前面所述理想流体,作为真实流体的近似模型,在解决阻力、能耗及扩散等问题是并不
解决问题思路:两种方法-(1)建立运动基本方程-作出合理简化-求得速度分布进而计算 其他量。 (2)所研究问题进行物理分析-适当模型简化-建立数学方程。 例如:对于非等温流动,考虑能量方程再结合初始条件和边界条件就可以确定系统中温度、 速度、浓度随位置及时间的变化。相应的,对于有化学反应时,考虑化学反应的速率方程再 结合初始条件和边界条件就可以求解。
不可压缩流体连续性方程可写为:
u x u y u z 0 x y z
将Stokes流方程组分别对x、y、z求导,然后相加,根据不可压缩流体连续性方程,可得:
2 2 2 2 2 2 u u u p p p ux ux ux u u u y y y z z z 2 2 2 2 2 2 2 2 2 2 2 2 x x x y z x x y z y y z z y z u x u y u z x 2 y 2 z 2 x y z 0 2 2 2 2 2 2
忽略变位加速度的非线性项而求得的线性解,常见的流动类型有库特流、管流、驻点流和旋转流等。
归纳出的利用奈维-斯托克斯方程解决问题的步骤如下: 1、选择适当的坐标系;2、作出合理简化;3、删去N-S方程中的可以忽略的项;4-给出合理边界条件; 5、求解方程;6、所求结果与实验或者观测值进行比较。 2. 近似解 为解决更多的工程实际问题,发展了许多近似解法,如:摄动法、润滑近似、准定常近似等,其中, 对N-S方程最为重要的近似处理是依照雷诺数的大小,分成低雷诺数流(Re≤1)近似和高雷诺数(
《环境流体力学》第七章 各向同性均匀湍流
形式(7-3-9)不变,2 阶张量的一般形式可以直接写作
Aij f (xi xi )xi x j g(xi xi )ij C(xi xi )ijk rk
(7-3-14)
7.4 各向同性湍流的相关张量函数及其性质
根据各向同性湍流场的定义,各向同性湍流场中 n 阶相关的表达式必为
Ri1i2 in Ri1i2 in (ξ1, ξ2 , , ξn1) 。 称 ξ i 为相关向量。应用张量函数表达式,可以导出各向同性湍流场的各阶相关函数张量的
1)两点速度相关张量具有反对称性
Ri, j (ξ) ui (x)u j (x ξ) ui (x'ξ)u j (x') Rj,i (ξ)
2)一点自相关函数总是大于等于两点自相关函数
(7-2-1)
Rii (ξ) Rii (0)
(7-2-2)
对一般 2 阶互相关 Rij (ξ) ui (x)u j (x ξ) 应用 Schwartz 不等式,有
定义 1:沿相对向量方向的脉动速度分量的 2 阶相关称作两点纵向相关 R(ll ξ)。
定义 2:垂直于相对向量方向脉动速度分量的 2 阶相关称作两点横向相关 Rn(n ξ)。 由(7-4-1),有 R(ll ξ) ξ2 f (ξ) g(ξ), Rnn (ξ) g(ξ) 。
可以解出 f (ξ) (Rll (ξ) Rnn (ξ)) / ξ2 和 g(ξ) Rn(n ξ),它们分别称为纵向相关系数和横向
我们还要利用第二个重要原理:如果左边 Bi C j 是多重线性的,则右边也应是多重线
性的。这样我们就可以排除 Bi Bi 等高次幂函数, BiCi 必然出现为线性乘子。也就是说,函
数式中只能包含 n 个向量的线性积。所以(7-3-2)成为
(2021年整理)湍流的数值模拟综述
湍流的数值模拟综述编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(湍流的数值模拟综述)的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为湍流的数值模拟综述的全部内容。
湍流的数值模拟一、引语流体的流动形态分为湍流与层流。
而层流是流体的最简单的一种流动状态。
流体在管内流动时,其质点沿着与管轴平行的方向作平滑直线运动。
此种流动称为层流或滞流,亦有称为直线流动的.流体的流速在管中心处最大,其近壁处最小。
管内流体的平均流速与最大流速之比等于0.5,根据雷诺实验,当雷诺准数引Re<2320时,流体的流动状态为层流.当雷诺数Re〉2320时,流体流动状态开始向湍流态转变,湍流是一种很复杂的流动状态,是流体力学中公认的难题。
自从19世纪末O.Reynolds提出湍流的统计理论以来,已经有一个多世纪了,经过几代科学家的努力,湍流研究取得很大进展,但是仍然不能满足工程应用的需要,以至于经常有悲观的论调侵袭湍流研究。
为什么湍流问题没有圆满地解决会受到如此关注呢?因为湍流是自然界和工程中十分普遍的流动现豫,对于湍流问题的正确认识和模化直接影响到对自然环境的预测和工程的质量.例如,当前影响航天器气动力和气动热预测准确度的主要障碍是缺乏可靠的湍流模型。
和其他一些自然科学的准题不同,解决湍流问题具有迫切性。
湍流运动的最主要特征是不规则性,这是大家公认的。
对于湍流不规则性的深入认识,是一百多年来湍流研究的上要成就之一。
早期的科学家认为,像分子运动一样,湍流是完全不规则运动。
类似于分子运动产生黏性,湍流的耗散可以用涡黏系数来表述。
20世纪初,一些杰出的流体力学家,相继对涡黏系数提出各种流体力学的模型,如Taylor(1921年)的涡模型,Praudtl(1925年)的混合长模型和von Karman (1930年)相似模型等。
湍流力学复习题(附答案)
1、 湍流定义及其主要特点?湍流定义(Definition of Turbulence )① When the flow reaches the certain conditions, generally Re numbersover certain value, flow pattern will transit from laminar to turbulence. ② For turbulent flow, flow pattern shows irregular and random motion inspace and time. A large number of eddies of different size consist the flow.③ It shows the characteristics of dissipation.湍流是一种流动的状态,当雷诺数超过某一临界值时流动会从层流转变为湍流。
这种流动状态显示出空间和时间上的不规则性和随机性,具有大量的不同尺度的涡结构,具有特殊的耗散特性。
湍流的特点:① 湍流被大含能涡主导,具有很强湍流扩散能力和较大湍流应力; ② 湍流涡存在使得其特性参数不能由当地参数确定,其特性与大涡历史行程有较大关系;③ 旋涡之间发生能量传递,其耗散速度受小涡从大涡接受能量速率控制。
2、 湍流随机性特性产生原因是什么?① 在任何湍流流动中不可避免存在小的摄动,主要存在于初始条件、边界条件和物性参数里。
② 湍流场特性表现为对这些摄动很强的敏感性。
在高雷诺数下流动体系高度敏感于小的摄动。
3、 概率密度函数(Probability density function )如何定义,其表征的物理意义是什么? ① PDF 可以由累积函数给出:任何一个事件的概率可以由累积概率(CDF)计算,其定义为F(V)=P{U<V}② 表征在样本空间内单位长度上随机变量概率的大小,其量纲为随机变量量纲的倒数。
fluent教程 第三章,湍流模拟
ui u j 2 ul ij uiu j x j xi 3 xl x j
上面两个方程称为雷诺平均的Navier-Stokes(RANS)方程。
uiu j
如果要求解该方程,必须模拟该项以封闭方程。
计算流体与传热传质
湍流模型
湍流模拟
D1
热科学与能源工程系
计算流体与传热传质
湍流模型
什么是湍流?
湍流:非定常,非周期性的三维速度脉动、 强化物质、动量和 能量的输运. 瞬时速度分解为平均速度和脉动速度: Ui(t) Ui + ui(t)
ui(t) U i (t) Ui
Time
压力、温度、组分浓度值具有类似的脉动
D18 热科学与能源工程系
计算流体与传热传质
湍流模型
目前采用的标准k-e模型方程为
t e e ( e ) ( e u j ) [( ) ] (Ce1 P Ce2 e ) t x j x j e x j k
u P uiuj i x j
其中
gTL3 Ra
D9
热科学与能源工程系
计算流体与传热传质
湍流模型
湍流特点
额外应变率
流向曲率 测向分离 加速或减速 有旋 回流 (或分离) 二次流
3D振荡流动 Transpiration (吹风/吸气) 自由湍流 剪切层相互作用
D10
热科学与能源工程系
i 1
D7 热科学与能源工程系
3
计算流体与传热传质
湍流模型
湍流应力
若以主对角线上的三个分量作为对称轴,则对称的两个切应力分 量是相等的,很显然,这是一个对称的二阶张量。很容易可以证 明,在各向同性湍流中,湍流正应力的三个分量相等,即
第3章,fluent湍流模型-1
第三章,湍流模型第一节, 前言湍流流动模型很多,但大致可以归纳为以下三类:第一类是湍流输运系数模型,是Boussinesq 于1877年针对二维流动提出的,将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。
即:2121x u u u t ∂∂=''-μρ 3-1 推广到三维问题,若用笛卡儿张量表示,即有:ij ijj i t j i k x u xu u u δρμρ32-⎪⎪⎭⎫ ⎝⎛∂∂+∂∂=''- 3-2 模型的任务就是给出计算湍流粘性系数t μ的方法。
根据建立模型所需要的微分方程的数目,可以分为零方程模型(代数方程模型),单方程模型和双方程模型。
第二类是抛弃了湍流输运系数的概念,直接建立湍流应力和其它二阶关联量的输运方程。
第三类是大涡模拟。
前两类是以湍流的统计结构为基础,对所有涡旋进行统计平均。
大涡模拟把湍流分成大尺度湍流和小尺度湍流,通过求解三维经过修正的Navier-Stokes 方程,得到大涡旋的运动特性,而对小涡旋运动还采用上述的模型。
实际求解中,选用什么模型要根据具体问题的特点来决定。
选择的一般原则是精度要高,应用简单,节省计算时间,同时也具有通用性。
FLUENT 提供的湍流模型包括:单方程(Spalart-Allmaras )模型、双方程模型(标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型)及雷诺应力模型和大涡模拟。
湍流模型种类示意图包含更多 物理机理每次迭代 计算量增加提的模型选RANS-based models第二节,平均量输运方程雷诺平均就是把Navier-Stokes 方程中的瞬时变量分解成平均量和脉动量两部分。
对于速度,有:i i i u u u '+= 3-3其中,i u 和i u '分别是平均速度和脉动速度(i=1,2,3)类似地,对于压力等其它标量,我们也有:φφφ'+= 3-4 其中,φ表示标量,如压力、能量、组分浓度等。
湍流实验
流动观察
流动观察
图2湍流边界层底层内的湍斑流动流动观察这是直接获得湍流的各种流动图案和大尺度涡旋的形成、发展和 衰变过程的直观方法。通常观察气流可以用纹影法和干涉法光学技术(见风洞测试仪器),也可以用烟迹法;观 察液体可以用染色法和氢气泡技术。近年来又发展了激光干涉与全息等技术。虽然流动观察是古老而又简单的实 验技术,但在湍流结构上的不少重要发现主要是用此方法得到的。例如二维混合层流动中的布朗-罗什科大尺度 涡旋结构及其合并过程(图1),湍流边界层底层内的湍斑流动及其猝发过程(图2),就是从实验上首先观察到 的。
谢谢观看
湍流实验测量技术原则上都是由感受部分和数据处理部分组成的。用于湍流实验的感受部分的传感元件有压 电元件、热敏元件、热丝和热膜、等离子束和激光束等。但迄今用得最广泛的还是热丝探头。它是一条长度远大 于直径的细铂丝或钨丝,其直径一般在1~5微米之间。测量原理是把热丝置于待测流体介质中,用电加热热丝, 使其温度高于流体介质温度。由热丝与流体之间的热传递的变化引起热丝两端的电压发生变化,从而可以测量流 速的平均值和脉动值。近年来又把激光技术应用到湍流测量中。
2.LianQixiang,Observationson Structures ofTurbulentBoundaryLayer in a Flow with Strong Adverse Pressure Gradiant, Proceedings of Second AsianCongressofFluidMechanics, pp.70~ 77,Science Press,Beijing, China, 1983.
对某些简单的均匀时均流场,如果湍流脉动是均匀的、各向课件
强烈的化学反应、气流的旋转流动、颗粒的 存在以及浮力或电磁场的作用,都会影响当地 的湍流结构。
精
3
2、湍流的基本方程
无论湍流运动多么复杂,非稳态的连续方程和N-S方程对
于湍流的瞬时运动仍然是适用的。在此,考虑不可压流 动,使用笛卡尔坐标系,速度矢量在x、y和z方向的分量 分别为u、v和w,写出湍流瞬时控制方程如下:
精
5
时均量与脉动量的关系
'
(4)
现在用平均值和脉动值之和代替流动变量,即:
u u u';u u u';v v v';w w w'; p p p'
(5)
将(5)代入瞬时状态下的连续性方程(1)和动量方 程(2),并对时间取平均,得到湍流时均流动的控制 方程如下:
u'
x
'
v'
y
'
w' '
z
S
(8)
精
8
密度脉动的影响
以上是假设流体密度为常数; 但是在实际流动中,密度可能是变化的。 Bradshaw等指出,细微的密度变动并不对流动造成明显
的影响 在此,忽略密度脉动的影响,但考虑平均密度的变化,
写出可压湍流平均流动的控制方程如下 注意,为方便起见,除脉动值的时均值外,下式中去掉
湍流是流体力学中的难题
对某些简单的均匀时均流场,如果湍流脉动是 均匀的、各向同性的,可以用经典的统计理论进 行分析。 但实际上,湍流是不均匀的。
均匀各向同性湍流中颗粒所见被动标量的数值研究
(. 东 电 网公 司 电 力 科 学 研 究 院 ,广 州 5 0 8 ;2 1广 10 0 .华 中 科 技 大 学 煤 燃 烧 国家 重 点 实 验 室 ,武 汉 4 0 7 ) 304
摘
要 :对均匀各 向同性湍流 中惯性 颗粒所见被 动标量 的统计特 性进行 了数值模拟研 究 ,探讨 了颗粒惯 性的差异对
根 据 概 率 密 度 函数 其 状 态空 间所 包 含 的 物理 量
的不 同 , DF 方 法 可 以划 分 为两 类 模 型 : P 一类 是 状态
空 间仅 包 含如 颗粒 速 度 和颗粒 温 度等 颗 粒 自身参 数 ,
5562 NE T0 —7 8) 基 金 项 目 : 国家 自然 科 学 基 金 资 助 项 目 ( 0 7 0 7);教 育 部 新 世 纪 优 秀 人 才 计 划 资 助 项 目 ( C -40 0
气 固两相 流广 泛存 在 于 工程 应 用 中 , 如煤 粉 锅炉
相 关 思路 , 如 颗粒 湍 动能 、 粒雷 诺应 力 等概 念 , 例 颗 在 数 值求 解 中也表 现 出一 定 的优 越性 , 这 种类 比在 物 但 理 含 义上 是 不 够严 谨 的 , 以 , 所 从微 观 的颗粒 拉 氏动 力 学 方程 出发 导 出宏 观 的统 计矩 方 程 的 P DF 方 法 , 因其 物 理含 义更 为清 晰 而 成 为 多相 流 统 观 模 型 发 展 的一个 重要 方 向.
项 计算 中产 生 的混淆误 差.
时 间推进 格 式 采 用 了兼 顾 效 率 和精 度 的二 阶显 式 R n— t u gKut 法 , 算 中 C ua t数 ( ) 取 为 a方 计 o rn C选
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
上面三式构成了不可压湍流的基本方程
3.3 不可压缩均匀各向同性湍流的动力学方程
下面考察不可压缩均匀湍流场中的湍动能和雷诺应力的演化。对于不 可压缩均匀湍流,一点的统计相关量的空间导数等于零,因此它的湍动能 利雷诺应力方程可以简化为:
由上式可见,均匀湍流场中湍动能总是耗散的
3.4.2 不可压缩均匀各向同性湍流的Karman-Howarth方程
脉动压强和速度相关项的作用使湍流脉动速度各向同性化, 一旦湍流场达到各向同性状态,压强—速度相关项就不再有任 何作用。
上式即为不可压缩各向同性湍流的2阶纵向速度相关方程。最早 由Karman和 Howarth(1938)导出,称为卡门--霍华斯方程。Karman-Howarth方程是线性偏微分方程,较之原始变量的N—S方程要简单 得多,不过,Karmnn--Howarth方程仍然是不封闭的。
3.3 不可压缩均匀各向同性湍流的动力学方程
3.3.4 均匀湍流中的湍动能传输链
在粘性作用下,脉动速度逐渐衰减,而且小尺度的成分衰减 得最快,于是在耗散过程中大尺度脉动成分占更多份额。
由于惯性在速度脉动的各个尺度间进行动量输运,它将大尺 度脉动的动能传输给小尺度脉动。
于是在粘性和惯性的联合作用下,湍流脉动场形成一种能量 传输链:大尺度湍流脉动通过惯性作用向小尺度湍流脉动不断输 送能量,这股能量在小尺度湍流脉动中耗散殆尽。
,或者说湍流脉
动总是衰减的,初始的湍动能在演化过程中将耗散殆尽。
3.3 不可压缩均匀各向同性湍流的动力学方程
3.3.2 不可压缩均匀湍流的谱理论 对湍流速度场和压强场进行傅立叶展开,经过一系列变换,
可以得到谱空间中湍流脉动的演化方程:
忽略非线性作用,单独考察粘性耗散,有如下积分式:
粘性耗散使脉动以指数函数衰减,衰减指数和波数的平方成 正比,因此高波数成分迅速衰减,单纯衰减过程中各波段间相位 关系不变。
在能量传输过程中,压强在各个脉动分量间起调节作用,如 果物理空间中初始脉动场的动能在各个分量间分配不均匀,压强 梯度将使它们逐渐均分。
3.4 不可压缩均匀各向同性湍流动力学的若干性质
3.4.1 不可压缩均匀湍流的2阶速度相关动力学方程和谱张量动力学方程 动力学方程
谱张量动力学方程
3.4 不可压缩均匀各向同性湍流动力学的若干性质
根据各向同性湍流的定义,要求相关张量函数在坐标系转动时 具有不变性。应用前面知识,就可以导出各向同性湍流场中各阶相 关函数的表达式。
3.2 均匀各向同性湍流场的相关函数和谱张量
1.各向同性湍流场中两点2阶速度相关函数的表达式
各向同性湍流场中2阶速度相关函数必有以下的函数式:
选定两个特定几何构形的相关:纵向相关和横向相关, 它们的定义如图3-3所示。
3.3 不可压缩均匀各向同性湍流的动力学方程
3.3.3 不可压缩均匀湍流中的湍动能输运过程 谱空间中脉动动量输运的动力学方程:
湍流各个脉动成分间的相互作用有双重含义:一种是各个波段 间速度分量相互作用;另一种是同一波数下各个速度分量间的相互 作用,压强梯度和惯性分别担当这两种作用。 (1)压强在不同脉动分量间重新分配能量,而不改变给定尺度(波数) 的湍动能 (2)惯性作用产生各波段间动量传输,但不改变物理空间中各脉动分量 的平均能量
3.2 均匀各向同性湍流场的相关函数和谱张量
按照纵向相关和横向相关的定义带入,各向同性湍流场中 2阶速度相关的一般表达式,就有:
定义:沿相对位移方向的脉动速度分量的2阶相关称做两点纵向 相关 定义:垂直于相对位移方向的脉动速度分量的2阶相关称做两点 横向相关
3.2 均匀各向同性湍流场的相关函数和谱张量
4. 旋转系统中均匀湍流的两点速度相关
在无界旋转系统山的湍流场仍然可能是均匀的,但是由 于旋转系统中的哥氏惯性力,湍流场不再保持各向同性。即 使是无界湍流场是均匀的,这时2阶速度相关函数将是两点 相对向量 和角速度的函数,并可以写为:
坐标系转动无关的2阶张量函数表达式应是:
3.2 均匀各向同性湍流场的相关函数和谱张量
图3-7是各向同性湍流的典型传输谱 T k , 图3-8是典型的能谱 Ek
宣 传 解 释 工 作,使
上图表示均匀湍流的定义(只表示4阶相关),在均匀湍流场中4点 统计相关函数只和4点间3个相对向量有关
3.1 均匀湍流流场的相关函数和谱张量
均匀各向同性
定义:如果任意n点空间几何构形在空间中平移和转动时,脉动速 度的任意n阶统计相关函数的值不变(或任意n点联合概率密度不变), 则称该湍流场是均匀各向同性的。
3.4 不可压缩均匀各向同性湍流动力学的若干性质
3.4.3 Karman--Howarth方程的应用
1.不可压缩均匀各向同性湍流场中湍动能耗散方程和Taylor微尺度 湍动能耗散率直接和Taylor微尺度直接相关
或者
或者
上式表明Taylor微尺度是各向同性湍流中的湍动能耗散特征尺度
3.4 不可压缩均匀各向同性湍流动力学的若干性质
张量函数的坐标不变性: 任何物理定律或定理,不论它足用标量表示还是用张量表示,
都应当和坐标系的刚体转动无关。
3.2 均匀各向同性湍流场的相关函数和谱张量
3.2.2 各向同性湍流的相关张量函数及其性质
N点几何构形由n-1个相对向量 1,2,....,n1完全确定,因
此各向同性湍流流场中n阶相关的表达式为:
2 洛强斯基(1939)不变量
上面 称为洛强斯基不变量,它表示各向同性湍流场 衰减过程中某种统计特征量的守恒性。
3.均匀各向同性湍流的后期衰变 各向同性湍流的衰减后期,湍流脉动量很小,这时3阶相关 和2阶相关相比是高阶小量,因而可略去不计,于是不可压缩 各向同性湍流的后期衰变方程简化为:
3.4 不可压缩均匀各向同性湍流动力学的若干性质
通过分析表明,不可压缩各向同性湍流的后期衰变过程中,湍动 能随时间以(-2.5)次方的幂函数衰减;Taylor微尺度则随时间以0.5次 方增长;以上结果得到Batchelor和Townsend(1948)实验的很好验证, 见图3-6。
3.5 不可压缩均匀各向同性湍流中的湍动能传输
3.5.1 不可压缩均匀各向同性湍流中的湍动能输运方程 积分的谱空间湍动能方程可写为:
第三章 均匀各项同性湍流
3.1 均匀湍流流场的相关函数和谱张量 3.2 均匀各向同性湍流场的相关函数和谱张量 3.3 不可压缩均匀各向同性湍流的动力学方程 3.4 不可压缩均匀各向同性湍流动力学的若干性质 3.5 不可压缩均匀各向同性湍流中的湍动能传输链 3.6 局部各向同性湍流的结构函数 3.7 各向同性湍流相关方程的封闭
3.1 均匀湍流流场的相关函数和谱张量
均匀湍流场
定义:如果任意n点空间几何构形在空间中平移时,脉动速度任意 n阶统计相关函数的值不变(或任意n点联合概率密度不变),则称该 湍流场是均匀的。
述 职 报 告 从内 容上分 ,可以 分为综 合性述 职报告 、专题 性述职 报告、 单项工 作 述 职 报 告 三种, 你知道 它们有 什么区 别吗? 你知道 写述职 报告吗 ?如下 是中国 人 才 网 给 大 家整理 的销售 业务员 述职报 告范文 ,希望 对大家 有所作 用。 销 售 业 务 员 述 职 报告范 文篇【 一】 尊 敬 的 各位 领导、 同事们 : 今 年 以 来 ,作 为 销 区 经 理 ,我能 够认真 履行职 责,团 结带领 ××销 区全体 人员, 在厂部 总体工 作 思路指 引下, 在×厂 长和销 售部各 位经理 的正确 领导下 ,积极 进取、 扎实工 作, 完 成 了 全 年 目标任 务,总 销量达 到箱, 营销工 作取得 了可喜 的成绩 。下面 ,根据 领 导 要 求 , 我进行 述职, 不妥之 处,敬 请领导 和同事 们批评 指正。 一 、 加 强 宣 传 促 销 力 度,较 好地完 成了全 年目标 任务。 今 年 我 们 在销售 工作中 遇到了 一 定 的 困 难 ,特别 是因为 我厂面 临兼并 重组, 各种不 实传闻 使商业 公司对 我厂产 品 的 信 心 不 足,不 少零售 户甚至 不卖我 厂产品 ,面对 不利局 面,我 们在销 售部统 一 指 挥 和 安 排部署 下,发 挥全体 人员的 聪明才 智,进 一步加 强宣传 促销力 度。首 先 , 以 我 厂 产品进 入行业 优等品 为契机 ,迅速 传播信 息。通 过拜访 商业公 司、走 访 零 售 户 、 及时分 送《× ×企业 报》、 《宣传 画报》 等企业 宣传品 的方式 ,做好
因此它的2阶速度谱张量函数必有以下形式:
同样,不可压缩各向同性湍流场中的3阶速度谱张量 的一般表达式为:
3.3 不可压缩均匀各向同性湍流的动力学方程
3.3.1 不可压缩均匀湍流的基本方程 湍流满足Navier-Stokes方程。对于不可压缩流体的均匀湍流场,
它的平均流速是常向量,因此不失其一般性的讨论,可以令平均速度等 于零。于是,均匀湍流的脉动场(取消脉动量的上标“-”’号)满足NavierStokes方程:
(7)不可压缩均匀湍流场巾脉动涡量的2阶谱张量
(8)均匀不可压缩湍流场中的湍动能耗散 湍动能耗散率:
湍流耗散张量和湍动能耗散的积分表达式:
3.2 均匀各向同性湍流场的相关函数和谱张量
3.2.1 张量的不变量和张量函数
张量不变量概念:
有限阶张量只有有限个独立不变量,应用Coyley-Hamiton定 理可证明:2阶张量只有三个独立不变量,即2阶张量的迹以及它 的平方和立方的迹是三个独立不变量,2阶张量的高阶幂函数的 迹可以由以上三个独立不变量算出。
(4)不可压缩均匀湍流场中2阶速度谱张量有以下等式:
(5)均匀湍流场中2阶速度谱张量是Hermit张量,2阶速度 谱张量等于它的转置张量的复共轭:
(6)均匀湍流场中脉动涡量的2阶相关函数和谱张量,脉动涡量 是脉动速度的旋度