第七章第五节湍流谱理论
14
谱密度函数E(k)的说明
可见，在各向同性湍流中，由于湍流能量在 三维波数空间中分布的球对称性．可将三维空 间简化为以球面半径k为自变量的一维空间，则 E(k) 可理解为在这一维波数空间上对湍流能量 的贡献密度。它描述了湍流能量在各个波数 上，也即在各个长度尺度上的分布情况。故称 函数 E(k) 为能谱函数，也像纵向二元速度关联 函数f(r)一样。也是湍流统计理论中的一个主要 研究对象。 E(k) 是速度起伏能量在波数空间的 三维谱密度。 Φii(k) 表示波数空间内单位体积的能量， 15 E(k)为单位波数的能量。
2
E1 (k1 ) + 2 E2 (k1 ) =
1
π
∫
∞
−∞
e − ik1r1 Bii (r1 )dr1
24
E1(k)、 E2(k)、E (k)之间的关系
两边对k1微分，并利用(7.5.18)得
2 ∂ ∂ − i ∞ −ik1r1 E1 (k1 ) + 2 E 2 ( k1 ) = r1e Bii (r1 )dr1 = − E (k1 ) π ∫−∞ ∂k1 ∂k1 k1
23
5 E1(k)、 E2(k)、E(k)之间的关系
由 f(r) 与 g（r） 之间的关系能够得到 E1(k) 和E2(k)之间的关系:
∂E 2 (k1 ) 1 E 2 (k1 ) = [ E1 (k1 ) − k1 ] ∂k1 2
1 ∞ Bi ,i (r1 ) = u [ f ( x1 ) + 2 g ( x1 )] = ∫ [ E1 (k1 ) + 2 E 2 (k1 )]e ik1r1 dk1 2 −∞
v v v ui ( x )u j ( x + r )
ui (t )u j (τ + t )
4
概
述:欧拉时间相关和 拉格朗日时间相关
RL(τ)
假定湍流是准定常的，则这个统计平均值应 只是两时刻的间隔τ的偶函数。定义无量纲的 欧拉时间关联系数 RE(τ) 和拉格朗日时间相关
u 2 (t )u 2 (t + τ ) u2
2 2
u2
∞
22
4、f (r)和g(r)的Fourier变换
设 f (r) 和 g(r) 的 Fourier 变换分别为 E1(k) 和 E2(k),
1 ∞ ik1r1 u f ( r1 ) = ∫ e E1 (k1 ) dk1 2 −∞ u 2 ∞ −ik1r1 E1 (k1 ) = ∫−∞ e f (r1 )dr1 π
复 习
相关，各向同性 零阶矩，一阶矩，二阶矩，三阶矩
v Ai = A1 r
v Bij (r ) = B1 (ri ri )ri r j + B2 (ri ri )δ i j
Bi , j (r ) = [ B LL (r ) − B NN (r )]
ri r j r
2
+ B NN (r )δ ij
′ Bij ,k ≡ u i u j u k = A2 ri r j rk + B2 rk δ ij + C 2 r j δ ik + D2 ri δ jk
∂E 2 (k1 ) 1 ] E 2 (k1 ) = [ E1 (k1 ) − k1 2 ∂k1
1 2 ∂2 1 ∂E ( k ) E ( k 1 ) = k1 E1 (k1 ) − k1 1 1 2 ∂k1 2 2 ∂k1
25
E1(k)、 E2(k)、E (k)之间的关系
应用常数变易法解上面的微分方程得到
11
引入三维能谱函数E(k)
引入一个新的函数：三维能谱函数 E(k),
E (k ) = −4πk A(k )
4
v E (k ) 2 Φ ij (k ) = [ k δ ij − k i k j ] 4 4πk
1 E (k ) Φ ii ( k ) = 2 4πk 2
12Biblioteka Φii(k)Φii(k)是波数空间单位体积内的能量。 Φii(k)只是波矢量的模k的函数，具有球对
RE (τ ) =
u1 (t )u1 (t + τ ) u
2
RL (τ ) =
这些相关可以反映湍流的结构或尺度，但并 不能直接估计出湍涡能量的大小。
5
概述：湍涡的能量
某个大小的湍涡有自己的能量 ：
E ( k ), E ( n), E (ω )
k = 2π / λ
n = u /λ
ω = 2πn
6
概述：湍涡的能量
由相关矩和谱之间的变换及连续性方程，有
vv v v v v v v Bij ( r ) = ui ( x )u j ( r + x ) = ∫∫∫ e ik ⋅r Φ ij ( k ) dk
vv vv v v v v ∂ ∂ v ik ⋅r ik ⋅r Bij (r ) = ∫∫∫ e Φ ij (k )dk = ∫∫∫ik j e Φ ij (k )dk = 0 ∂r j ∂r j v k j Φ ij (k ) = 0
vv v v v v v v i k ⋅r Bij ( r ) = u i ( x )u j ( r + x ) = ∫∫∫ e Φ ij ( k ) dk
v Φ ij (k ) =
1 (2π )
2
3
∫∫∫ e
vv − i k ⋅r
v v Bij (r )dr
v v v Bij (0) = u i ( x ) = ∫∫∫ Φ ij (k )dk
2
2 ∞ ∞
1 ∞ ik1r1 u g ( r1 ) = ∫ e E 2 ( k1 ) dk1 2 −∞ u 2 ∞ −ik1r1 E2 (k1 ) = ∫ e g (r1 )dr1 π −∞
2
u = ∫ E1 (k1 )dk1 = ∫ E 2 (k1 )dk1
0 0
u = u1
2
2
1 2 = u ii 3
Bij ,k
1 (u 2 ) 3 / 2 2 2 ′ ≡ ui u j u k = {(k − h − 2 g )ri r j rk + hr rk δ ij + gr (r j δ ik + ri δ jk )} 3 r
§7．5湍流谱理论
一、概述 二、均匀各向同性湍流谱理论 三、准平衡区湍流谱特征 四、实验个例 五、总结和作业
2
一、概述
湍流是由大小不同的湍涡组成的，这些大大 小小的湍涡之间的关系是湍流的重要因素。 空间某一点同一方向速度起伏之间的相关反 映湍流强度或湍流能量； 空间某一点不同方向的速度起伏之间的相关 反映Reynolds应力或湍流动量通量；
u i u i = 3u
2
τ ij = u i u j
3
概
述
不同空间点或不同时间点之间的速度相关， 则反映湍涡空间尺度大小和时间尺度大小。例 如小于两点间距离的小湍涡，将使两点速度有 较大的差别，即速度相关小；而大于两点间距 离的大湍涡，将使两点间的速度比较一致（因 为处于同一个湍涡中），即两点的速度相关较 大。
∞
26
解释说明
积分下限的物理 意义是只有对应于 2π/k≤2π/k1 ， 即 k≥k1 的三维谱的那 些波长 2π/k 才对一 维 谱 的 波 长 2π/k1 的值有贡献。
2、 E(k)与二元相关矩B(r)关系
根据(6.5.20)
Bij (r ) = 4π ∫
∞ 0
v 2 sin( kr ) Φ ij ( k )k dk kr
E (k ) 1 Φ ii ( k ) = 2 4πk 2
Bii ( r ) = 4π ∫
∞
∞
0
v 2 sin( kr ) Φ ii (k )k dk kr
湍谱密度E(k)，总的湍流能量：
3u = ∫ E (k )dk
2 0 ∞
u = u1
2
2
1 2 = u ii 3
u = ∫ E (n)dn = ∫ E (ω )dω
2 0 0
∞
∞
7
概述：相关矩和谱
相关矩Bij和能谱Φij张量之间的变换：
vv v v v u i ( x ) = ∫∫∫ e ik ⋅x dZ i ( k )
E1 ( k1 ) = k [ ∫
2 1 k1 k1 E (k ) E (k ) dk 3 + C ] − ∫ dk +D 0 k k
0
k1 →∞
lim E1 (k1 ) = 0
∞ 2 1 2
E (k ) k (1 − ) E1 (k1 ) = ∫ dk k1 k k
k12 1 E (k ) (1 + 2 ) E 2 (k1 ) = ∫ dk 2 k1 k k
称性。如让它在半径为k的球面S(k)上积分
E (k ) 1 Φ ii (k )dS (k ) = ∫∫ dS (k ) = E (k ) 2 ∫∫) 2 S (k S ( k ) 4πk
可见E(k)就代表在三维波数空间中，在 半径为k的球面上的能量密度的总和。
13
湍流总能量
湍流总能量为
v 1 E (k ) v ∫∫∫ Φ ii (k )dk = ∫∫∫ 4πk 2 dk 2 E (k ) 2 k sin ϕ sin θdϕdθdk = ∫∫∫ 2 4πk 1 2 = ∫∫∫ E (k )dk = uii 2
π
0
16
E(k)和Bii(k)
E(k)和Bii(k)互为Fourier正弦变换 ：
Bii (r ) = 2 ∫
E (k ) =
∞
0
sin( kr ) v E (k )dk kr
∫ π
1
∞
0
kr sin krBii ( r )dr
高等传热学Chap7
此处，k1是一个常数。如果把k1合并到l中来定义混合长度 混合长度， 混合长度 且采用相同的符号，考虑到τt和(du/dy)有相同的正负号，则 有相同的正负号 有
τ t = ρl 2 du / dy du / dy
ε m = l du / dy
2
普朗特混合长度理论
§7-1 基本概念
一般情况下，混合长度是一个既与速度无关 混合长度是一个既与速度无关、 混合长度是一个既与速度无关、也与流 体的物性无关的量， 体的物性无关的量，只取决于流体微团脉动的距离，普朗 只取决于流体微团脉动的距离 普朗 特假定它与离壁面的法向距离y成正比， 成正比，即
§7-1 基本概念
4. 湍流热流密度
在湍流对流换热时，流体的温度脉动也要引起附加的湍流 热流密度。
q = qi + qt = −λ∂t / ∂y + ρc p v′t ′
分子热扩散 采用类似的方法有
湍流热扩散
qt = ρc p v′t ′ = − ρc p ε h ∂t / ∂y
q = ql + qt = − ρc p (α + ε h )∂t / ∂y
混合长度
L. Prandtl (1875-1953)
§7-1 基本概念
速度u的脉动是上述两个速度差的平均值， 的脉动是上述两个速度差的平均值，它具有以下 的定性关系： 的定性关系：u’~l(du/dy)。另一方面，u’是由v’引起，可以 预见两者具有相同的数量级，由此可得：
τ t = − ρ u ′v′ = k1 ρl 2 (du / dy )2
式中，εh是湍流热扩散率。 是湍流热扩散率 εm和εh的比值称为湍流普朗特数， 的比值称为湍流普朗特数
Prt = ε m / ε h
湍流理论介绍
湍流理论湍流理论theory of turbulence研究湍流的起因和特性的理论，包括两类基本问题：①湍流的起因，即平滑的层流如何过渡到湍流；②充分发展的湍流的特性。
湍流的起因层流过渡为湍流的主要原因是不稳定性。
在多数情况下，剪切流中的扰动会逐渐增长，使流动失去稳定性而形成湍流斑，扰动继续增强，最后导致湍流。
这一类湍流称为剪切湍流。
两平板间的流体受下板面加热或由上板面冷却达到一定程度，也会形成流态失稳，猝发许多小尺度的对流；上下板间的温差继续加大，就会形成充分发展的湍流。
这一类湍流称热湍流或对流湍流。
边界层、射流以及管道中的湍流属于前一类；夏天地球大气受下垫面加热后产生的流动属于后一类。
为了弄清湍流过渡的机制，科学家们开展了关于流动稳定性理论(见流体运动稳定性)分岔(bifurcation)理论和混沌(chaos)理论的研究,还进行了大量实验研究（见湍流实验）。
对于从下加热流层而向湍流过渡的问题，原来倾向于下述观点：随着流层温差的逐渐增加，在发生第一不稳定后，出现分岔流态；继而发生第二不稳定，流态进一步分岔；然后第三、第四以及许多更高程度的不稳定接连发生；这种复杂的流动称为湍流。
实验结果支持这一论点。
但是，这一运动过程在理论上得不出带有连续谱的无序运动，而与实验中观察到的连续谱相违。
最近，对不稳定系统的理论分析提出了另一种观点：在发生第一、第二不稳定之后，第三不稳定就直接导致一个可解释为湍流的无序运动。
这一观点也得到实验的支持。
剪切流中湍流的发生情况更为复杂。
实验发现，平滑剪切流向湍流过渡常会伴有突然发生的、作奇特波状运动的湍流斑或称过渡斑。
可以设想，许多逐渐形成的过渡斑，由于一再出现的新的突然扰动而互相作用和衰减，使混乱得以维持。
把过渡斑作为一种孤立的非线性波动现象来研究，有可能对湍流过渡现象取得较深刻的理解。
因此，存在着不止一条通向湍流的途径。
过去认为，一个机械系统发生无序行为往往是外部干扰或外部噪声影响的结果。
(完整word版)湍流模型理论
湍流模型理论§3.1 引言自然界中的实际流动绝大部分是三维的湍流流动，如河流,血液流动等。
湍流是流体粘性运动最复杂的形式，湍流流动的核心特征是其在物理上近乎于无穷多的尺度和数学上强烈的非线性，这使得人们无论是通过理论分析、实验研究还是计算机模拟来彻底认识湍流都非常困难。
回顾计算流体力学的发展,特别是活跃的80年代，不仅提出和发展了一大批高精度、高分辨率的计算格式,从主控方程看相当成功地解决了Euler方程的数值模拟,可以说Euler方程数值模拟方法的精度已接近于它有效使用范围的极限；同时还发展了一大批有效的网格生成技术及相应的软件，具体实现了工程计算所需要的复杂外形的计算网格；且随着计算机的发展，无论从计算时间还是从计算费用考虑，Euler方程都已能适用于各种实践所需。
在此基础上，80年代还进行了求解可压缩雷诺平均方程及其三维定态粘流流动的模拟。
90年代又开始一个非定常粘流流场模拟的新局面，这里所说的粘流流场具有高雷诺数、非定常、不稳定、剧烈分离流动的特点，显然需要继续探求更高精度的计算方法和更实用可靠的网格生成技术.但更为重要的关键性的决策将是，研究湍流机理,建立相应的模式，并进行适当的模拟仍是解决湍流问题的重要途径。
要反映湍流流场的真实情况,目前数值模拟主要有三种方法：1。
平均N-S方程的求解，2。
大涡模拟（LES），3。
直接数值模拟（DNS)。
但是由于叶轮机械内部结构的复杂性以及目前计算机运算速度较慢，大涡模拟和直接数值模拟还很少用于叶轮机械内部湍流场的计算,更多的是通过求解平均N-S方程来进行数值模拟。
因为平均N-S方程的不封闭性，人们引入了湍流模型来封闭方程组,所以模拟结果的好坏很大程度上取决于湍流模型的准确度。
自70年代以来，湍流模型的研究发展迅速，建立了一系列的零方程、一方程、两方程模型和二阶矩模型，已经能够十分成功的模拟边界层和剪切层流动。
但是,对于复杂的工业流动，比如航空发动机中的压气机动静叶相互干扰问题，大曲率绕流，激波与边界层相互干扰，流动分离，高速旋转以及其他一些原因，常常会改变湍流的结构，使那些能够预测简单流动的湍流模型失效，所以完善现有湍流模型和寻找新的湍流模型在实际工作中显得尤为重要。
turbulent flows 第8版 pope 译注
turbulent flows 第8版 pope 译注Turbulent Flows 是一本由 Stephen B. Pope 所著的著名流体力学领域的经典教材，本文对该书的第8版进行了详细的译注。
以下是对该书部分章节内容的梳理和解析。
第一章简介本书是关于湍流流体力学的权威教材，主要介绍了湍流的基本概念和理论，以及相关的实验和数值模拟方法。
通过深入分析湍流现象和其背后的数学模型，读者将能够更好地理解和预测湍流的行为。
第二章湍流的描述湍流是一种复杂且难以捉摸的流动现象。
书中介绍了湍流的统计描述，包括涡旋相关、相关时间和长度尺度等基本概念。
此外，还详细阐述了湍流的能谱和相关的能量传递机制，为后续的章节打下基础。
第三章计算流体力学中的湍流模型计算流体力学（CFD）作为一种重要的湍流研究方法，被广泛应用于工程和科学领域。
本章介绍了常用的湍流模型，包括雷诺平均湍流模型（RANS），大涡模拟（LES）和直接数值模拟（DNS）。
通过比较不同模型的优劣，读者将能够选择适合自己研究对象的湍流模型。
第四章湍流的数值模拟本章主要介绍了湍流的数值模拟方法，包括有限体积法、有限元法和谱方法等。
通过数值模拟，可以更加深入地研究湍流的特性和行为。
同时，书中还涵盖了一些常见的湍流模拟技巧，如网格生成和边界条件设定等，帮助读者掌握湍流模拟的实质。
第五章湍流的统计理论湍流的统计理论是湍流研究的重要组成部分。
本章详细介绍了湍流的统计特性，包括湍流的概率密度函数、相关函数和湍流湍度等。
此外，还阐述了重要的湍流统计理论模型，如湍流统计平衡理论和尺度相似理论等，为读者进行湍流统计的研究提供了重要参考。
第六章湍流的实验技术湍流的实验研究是湍流研究的基础性工作之一。
本章介绍了一些常见的湍流实验技术，包括激励湍流、热线湍流和粒子图像测速法等。
通过实验手段，可以直接观测湍流的各种特性和行为，为湍流理论的验证提供了有力支撑。
第七章壁湍流壁湍流是湍流研究的重要分支，也是工程流体力学中的核心问题之一。
湍流理论学习
湍流理论学习1．层流和湍流粘性流体的运动存在着两种完全不同的流动状态：层流状态和湍流状态。
雷诺首先于1883年通过做圆管内流动实验观察到层流与湍流现象。
当圆管中流动速度较小时，管中的流线之间层次分明，互不掺混，这样的流动称为层流。
当流速增大后，流体作复杂、无规律、随机的不定常运动，称为湍流。
流动状态与雷诺数e R 、下临界雷诺数ec R 和上临界雷诺数ecR '有关。
当e ec R R ≤时，流动为层流；当ec e ec R R R '≤≤时，流动为不稳定过渡状态；当e ecR R '>时，流动为湍流。
湍流是在连续介质范畴内流体的不规则运动，它有别于物质分子的不规则运动。
具体来说，在极不规则的湍流中，流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度。
因此湍流运动产生的质量和能量的输运将远远大于分子热运动产生的宏观输运。
2．湍流的平均化、雷诺粘性应力经典的湍流理论认为，湍流是一种完全不规则的随机运动，湍流场中的物理量在时间和空间上呈随机分布，不同的瞬时有不同的值，关注某个瞬时的值是没有意义的。
因此，雷诺首创用统计平均方法来描述湍流的随机运动，即对各瞬时量进行平均得到有意义的平均值。
从N-S 方程出发，利用平均化运算的法则推导平均物理量满足的方程组。
只考虑不可压缩流体情形，假设体力可以忽略，此时，N-S 方程具有下列形式1110u u u u p u v w u t x y z x v v v v p u v w v t x y z y w w w w p u v w w t x y z z u v w x y z υρυρυρ∂∂∂∂∂⎧+++=-+∆⎪∂∂∂∂∂⎪∂∂∂∂∂⎪+++=-+∆⎪∂∂∂∂∂⎪⎨∂∂∂∂∂⎪+++=-+∆⎪∂∂∂∂∂⎪∂∂∂⎪++=⎪∂∂∂⎩(1) 运用(1)式中的连续性方程，运动方程可改写为2221110u u uv uw pu t x y z x v uv v vw pv t xy z y w uw vw w p w t x y z z u v w x y z υρυρυρ⎧∂∂∂∂∂+++=-+∆⎪∂∂∂∂∂⎪⎪∂∂∂∂∂+++=-+∆⎪∂∂∂∂∂⎪⎨∂∂∂∂∂⎪+++=-+∆⎪∂∂∂∂∂⎪∂∂∂⎪++=⎪∂∂∂⎩(2) 对方程组(2)中各式两边进行平均化运算，并利用平均化运算法则得到____22____22____221110u u uv uw u u v u w p ut x y z x y z x v uv v vw u v v v w p v t x y z x y z y w uw vw w u w v w w p w t x y z x y z z u v w x y z υρυρυρ⎧'''''∂∂∂∂∂∂∂∂++++++=-+∆∂∂∂∂∂∂∂∂'''''∂∂∂∂∂∂∂∂++++++=-+∆∂∂∂∂∂∂∂∂'''''∂∂∂∂∂∂∂∂++++++=-+∆∂∂∂∂∂∂∂∂∂∂∂++=∂∂∂⎪⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎪⎩(3) 考虑到方程组(3)的第四式，方程组(3)中的头三个方程可改写成另一种形式，把脉动项移到右边，得到()()()()____2____2u u v u w u u u u p u v w u t x y z x x y z v u v v w v v v v p u v w v t x y z y x y z w w w w p u v w w tx y z z ρρρρμρρρρμρρμ⎛⎫'∂- ⎪''''∂-∂-⎛⎫∂∂∂∂∂⎝⎭+++=-+∆+++⎪∂∂∂∂∂∂∂∂⎝⎭⎛⎫'∂- ⎪''''∂-∂-⎛⎫∂∂∂∂∂⎝⎭+++=-+∆+++⎪∂∂∂∂∂∂∂∂⎝⎭'∂-⎛⎫∂∂∂∂∂+++=-+∆+⎪∂∂∂∂∂⎝⎭()()____20w u w v w x y z u v w x y z ρρ⎧⎪⎪⎪⎪⎪⎪⎪⎪⎨⎪⎛⎫⎪'∂- ⎪⎪'''∂-⎝⎭++⎪∂∂∂⎪⎪∂∂∂⎪++=∂∂∂⎪⎩(4)将上式和应力形式的运动方程对比d div dtρ=VP 其中P 是应力张量，则有2p μ'=-++P I S P其中I 是单位张量。
湍流理论和湍流模型(博士课程课件)
工程上，将下临界雷诺数作为流态的判断依据。
1.1 湍流的不规则性
湍流速度场是时间、空间坐标、实验次数的不规则函数
ui ui(x,t, ~)
在不规则湍流中，流动的最小时间尺度和最小空间尺度都远远大于分 子热运动的相应尺度，因此湍流运动产生质量和能量的输运远大于分 子热运动产生的宏观输运，所以湍流场中质量和能量的平均扩散远大 于层流扩散。
1.5 湍流脉动的测量原理
湍流脉动的时间序列具有宽频带，测量仪器准确、响应特性好。
测量点的脉动速度的时间序列测量方法：热丝风速计法、激光多普勒测速法
脉动场的脉动速度的时间序列测量方法：统称为粒子图像测速法(PIV, paticle
数据采集的要求
image velocimetry)
(1) 测量精度：仪器精度+电子系统的高信噪比和宽频带的频率响应特性
脉动速度频谱
Su(u )2 1 R u(u )ex i p)d (
其逆变换为
R u(u ) Su(u )exip)d (
时间相关函数与频谱是一一对应的，
它们是统计量在时域和频域之间的转换。
当τ=0时 u2 Suu()d , Suu(ω)表示湍动能在频带中的分布，它在所有
频段上的积分等于湍动能的系综平均或时间平均值。
湍流：流体作复杂的、无规则的、随机的非定常运动，也称紊流;
上临界流速：层流变湍流 下临界流速：湍流变层流
ReRce(232) 0
Vc',
Rc'
Vc'd
Vc'd
Vc,
Rc
VcdVcd
流动为层流
Vc Vc'
RceReRc'e(138)0流0动为不稳定的过渡状态
湍动能能谱
湍动能能谱一、引言湍流，作为自然界和工程领域中广泛存在的复杂流动现象，一直是科学家们研究的热点。
在湍流的研究中，湍动能能谱是一个关键的指标，用于描述湍流运动的能量分布。
本文将详细探讨湍动能能谱的基本概念、测量方法、应用以及未来展望。
二、湍动能能谱的基本概念湍动能能谱，或湍动能谱，是描述湍流运动中能量分布的数学工具。
它反映了湍流中不同尺度运动的能量，可以用来解析湍流的结构和特性。
湍动能能谱通常被定义为湍流场中不同频率下的湍动能的分布。
根据不同的频率或波数，湍动能能谱可以分为低波数区和高波数区。
低波数区对应于大尺度湍流结构，而高波数区则对应于小尺度湍流结构。
不同的波数区间对应着不同的湍流运动状态和特性。
因此，对湍动能能谱的研究有助于深入理解湍流的复杂性和本质。
三、湍动能能谱的测量方法测量湍动能能谱的方法有多种，主要分为直接测量法和模型反演法两大类。
1.直接测量法：利用各种传感器直接测量湍流场中的速度或压力脉动，通过统计方法得到湍动能能谱。
这种方法需要精密的测量设备和复杂的实验条件，但结果较为准确。
2.模型反演法：利用湍流模型对数值模拟或实验室实验得到的数据进行反演计算，得到湍动能能谱。
这种方法不需要高精度的测量设备，但结果受模型精度和数值方法的影响较大。
在实际应用中，研究者通常根据具体的研究条件和需求选择合适的测量方法。
四、湍动能能谱的应用1.工程流动模拟与控制：在流体动力学、航空航天、环境工程等领域，了解湍流场的能量分布对于流动模拟和优化控制至关重要。
通过研究湍动能能谱，可以深入了解流动特性，为工程设计和优化提供依据。
2.地球科学：在气象学、海洋学等领域，湍动能能谱被用于研究大气和海洋中的湍流运动。
例如，在气象预测中，了解大气中的能量分布有助于提高天气预报的准确性。
3.生物医学：在生物流体力学中，湍动能能谱可用于研究血液流动、肺部呼吸等生理过程中的湍流现象。
这有助于深入了解人体循环系统和呼吸系统的功能和机制。
第七章 湍流 流体力学课件
t x y z x y z
x
将上式展开，利用平均化的连续方程，进行简化，可 以得到：
u u u v u w u 1 p 2 u uu uv uw
t x y z x
x y z
u(u v w ) 0 x y z
这就是 x 方向的平均运动方程（雷诺方程）。
Chen Haishan NIM NUIST
同理，可以得到 y ，z 方向的平均运动方程，最终得到形式如
下的平均运动（雷诺）方程：
(
u t
u
u x
v
u y
w
u) z
p x
2
( uu) x
( uv) y
( uw) z
(
v t
u
v x
v
v y
w
v) z
p y
2 v
( vu) x
( vv) y
( vw) z
(
w
u
w
v
w
w
w)
p
2 w
( wu)
如何判断流体运动的属性？确定湍流发生的条 件－－湍流判据问题。
以下简单介绍相关的 雷诺实验 在次基础上过给出确定湍流发生的判据－－临界 雷诺数及其在湍流研究中的应用。
Chen Haishan NIM NUIST
雷诺试验（1883年） 有色液体
流体
流速V V
管道直径d 流体的粘性
d
层流
过渡流
湍流
Chen Haishan NIM NUIST
p
pyx pyy pyz pzx pzy pzz
vu vv vw wu wv ww
Chen Haishan NIM NUIST
湍流的理论与分析
湍流的理论与分析湍流是一种复杂的流动形式，并且广泛存在于自然界和工程实践中。
对湍流的理论研究和分析不仅有助于深入理解流体现象，还可以为湍流控制和能源利用等方面提供支持。
本文将从湍流的定义、产生机理、湍流统计理论和湍流模拟等方面进行探讨。
一、湍流的定义湍流是指一种相对瞬态的流体运动状态，其中流体的速度和方向发生剧烈变化，造成流体的混合和扰动，呈现出随机不规则的涡动结构。
与层流（稳态流动）相比，湍流的运动特征更加复杂，无法用简单的数学公式描述。
湍流的主要特征为不规则、随机、涡动等。
二、湍流的产生机理湍流的产生机理复杂，其中包括传统的机械湍流、自然湍流、边界层失稳等多种因素。
机械湍流是由于固体物体运动时与周围介质相互作用产生的湍流现象，如风力机翼片和涡轮机叶片的湍流。
自然湍流是由于自然界中各种复杂流动引起的，如河流、海洋和大气的运动等。
边界层失稳是当涡旋从高速的流动区进入低速的流动区时产生的，例如水流从管道进入膨胀段时发生的湍流现象。
三、湍流统计理论湍流统计理论是对湍流运动规律的理论分析，是研究湍流基本性质和湍流现象的一种方法。
湍流统计理论中有两个重要的概念，一个是湍流的集成时间，另一个是湍流脉动，这两个概念分别给出了湍流时间与空间扰动中的统计特征。
其中湍流的集成时间是指机械能向湍流能转化和湍流能转化为机械能时所需的时间因子，而脉动是指在一个给定点的流动路径上，流体参数波动的相对不稳定性。
四、湍流模拟湍流模拟是一种基于数值计算的湍流研究方法，主要有两种方式：直接数值模拟（DNS）和大涡模拟（LES）。
直接数值模拟是对湍流运动的一种高精度的数值计算方法，它通过离散化流动中的微小物理尺度，运用数值方法以求解流场运动方程，得到高精度的湍流场数据。
但DNS需要的计算量庞大，计算成本高昂。
大涡模拟是在保留湍流中大尺度涡旋信息的同时，模拟和模拟所得的速度与涡旋脉动能谱于实验结果的吻合程度。
而LES所需要的计算量较之DNS低，同时保留的流场尺度也比DNS更大，能够得到更加直观的湍流现象展示。
高等流体-第八讲,湍流理论
U y u y t
U x u x
U y u y x
(U y u y )
2、涡粘性
涡粘性是流体流动的性质，而不像分子粘性那样是 流体的物理性质。涡粘性是指由于涡旋的运动引起的
类似于分子运动导致的粘性效应。湍流中小尺度的涡
旋对大尺度的涡旋的作用是一种扩散作用。涡粘性系 数起到削弱湍流平均运动速度的效果 。 3、级串 cascade 是指大小涡旋的能量传递，在级串过程中，第一级 大涡的能量一般来自外界。大涡失后产生第二级的小 涡，小涡失稳后产生更小的涡旋。 4、标度律 充分发展了的湍流运动，应服从标度律。
uy x i (x t ) u i ( y ) e y x
( y )ei (x t )
2 u u u u u d Ux y y x x x 代入： ( ) Ux ( ) uy t y x x y x dy2 uy 1 p 2 u x v ( ) x y x
u 主流流动（平均流动）Ux、 Uz P +小扰动流动 、U y x , u y , uz , p
u x U x u x u y U y u y z u z U z u
p P p
为了简单起见，考虑不可压缩流体，二维平行定常流动+二维非定常流
uy
1 P p Ux v 2u y t x y uy
由于主流满足： 1 P 2U x v y 2 x P 0 y u u dU x 1 p x x Ux uy v 2u x t x dy x uy uy 1 p Ux v 2u y t x y