环境水利学第3章 随流扩散与紊动扩散 (5)

1、第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度 l脉动性：各种流动参量如流速、压力等的值呈现强烈的脉动性：各种流动参量如流速、压力等的值呈现强烈的脉动现象，具有一定的随机性脉动现象，具有一定的随机性l不规则性：流体质点做极不规则的运动不规则性：流体质点做极不规则的运动l扩散性：流体的各项特性如动量、能量、温度和含有物扩散性：流体的各项特性如动量、能量、温度和含有物质的浓度等通过紊动向各方向传递质的浓度等通过紊动向各方向传递l三维有涡性：紊流是有涡运动，而且总具有三维的特性三维有涡性：紊流是有涡运动，而且总具有三维的特性l大雷诺数：流体的雷诺数超过某个临界值后，流动不稳大雷诺数：流体的雷

2、诺数超过某个临界值后，流动不稳定，扰动才能发展形成紊流。定，扰动才能发展形成紊流。紊流的特性紊流的特性第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度v因为紊流和紊动扩散是随机过程，在描述它的运动时，因为紊流和紊动扩散是随机过程，在描述它的运动时，常用统计平均方法。常用统计平均方法。v统计平均方法通常有时间平均法（简称时均法）、空间统计平均方法通常有时间平均法（简称时均法）、空间平均法和总体平均法共三种平均法。平均法和总体平均法共三种平均法。dtxTxT 01dVxVxV 1dqqxPxp )(第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度v各态历经：一个随机过程在重复多次试验出

3、现的所有各态历经：一个随机过程在重复多次试验出现的所有样本，亦将在一次试验的相当长时间或相当大的范围样本，亦将在一次试验的相当长时间或相当大的范围内出现，并且出现的概率相同。内出现，并且出现的概率相同。v因为紊流和紊动扩散是随机过程，在描述它的运动时，因为紊流和紊动扩散是随机过程，在描述它的运动时，常用统计平均方法。常用统计平均方法。v统计平均方法通常有时间平均法（简称时均法）、空间统计平均方法通常有时间平均法（简称时均法）、空间平均法和总体平均法共三种平均法。平均法和总体平均法共三种平均法。v如果随机过程是各态历经的，则时均值、空间平均值、如果随机过程是各态历经的，则时均值、空间平均值、总体

4、平均值三者是互等的。总体平均值三者是互等的。一、紊流的分类一、紊流的分类v紊流按其流动特点可分为可分两大类：均匀紊流按其流动特点可分为可分两大类：均匀各向同性紊各向同性紊流流和和剪切紊流剪切紊流。v在均匀紊流中，各种物理量的统计平均值当坐标平移时在均匀紊流中，各种物理量的统计平均值当坐标平移时, ,均保持不变，例如有均保持不变，例如有: :式中：式中：u u1 、u u2 、u u3 分别为沿三条直角坐标的脉动流速分量；分别为沿三条直角坐标的脉动流速分量； 字母上方的字母上方的“”示取统计平均示取统计平均( (例如取时间平均例如取时间平均);); C1、C2 、C3均为常量。均为常量。3232

5、22121,CCC u uu uu u第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度一、紊流的分类一、紊流的分类 3 2 1vvv232221u u u u u u v在均匀紊流中，如果各种物理量的统计平均值还与方在均匀紊流中，如果各种物理量的统计平均值还与方向无关，亦即当坐标轴作任何旋转或镜射时，各种物向无关，亦即当坐标轴作任何旋转或镜射时，各种物理量的统计平均值仍保持不变，例如有理量的统计平均值仍保持不变，例如有 则称这种紊流为均匀各向同性紊流，或简称为各向同则称这种紊流为均匀各向同性紊流，或简称为各向同性紊流。性紊流。v各向同性紊流只是一种理想化的最简单的紊流各向同性紊流只是一种

6、理想化的最简单的紊流232221vvv 第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度v凡不满足均匀性要求的紊流（当然也不满足各向同性）凡不满足均匀性要求的紊流（当然也不满足各向同性）, ,称为剪切紊流。称为剪切紊流。v当紊流中存在切应力时，就有流速梯度，导致各处的紊流当紊流中存在切应力时，就有流速梯度，导致各处的紊流统计量不相同，从而破坏了紊流的均匀性和各向同性。这统计量不相同，从而破坏了紊流的均匀性和各向同性。这种紊流是最常见的，它比各向同性紊流复杂得多。种紊流是最常见的，它比各向同性紊流复杂得多。v在剪切紊流中，存在着尺度由大到小的一系列涡体。研究在剪切紊流中，存在着尺度由大到小

7、的一系列涡体。研究证实，大涡区和中涡区受外界条件的明显影响，不是各向证实，大涡区和中涡区受外界条件的明显影响，不是各向同性的，但小涡区不受外界条件的直接影响，常近似地具同性的，但小涡区不受外界条件的直接影响，常近似地具有各向同性的性质，这称为局部各向同性。有各向同性的性质，这称为局部各向同性。第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度 在紊流中，常要分析两个脉动流速分量的相关矩在紊流中，常要分析两个脉动流速分量的相关矩( (即协方即协方差差) )，它们表征着紊动的重要性质。，它们表征着紊动的重要性质。二、欧拉相关和紊流尺度二、欧拉相关和紊流尺度)()( u u u u aiai第三

8、节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度1 1、欧拉空间相关、欧拉空间相关 定义为：定义为：图图 脉动流速分量示意图脉动流速分量示意图 其中其中: 为在为在i方向上的脉动流速分量；方向上的脉动流速分量； a 为为a点的某一方向的坐标，例如取为点的某一方向的坐标，例如取为xa或或ya ; a + 为为另一点的同一方向的坐标，相应为另一点的同一方向的坐标，相应为xa + D Dx 或或 ya + D Dy。iu uxxa + D Dxxaiu u)()( u u u u aiai第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度1 1、欧拉空间相关、欧拉空间相关 定义为：定义为：指同

9、一瞬时、不同两点的指同一瞬时、不同两点的同一方向脉动流速分量的同一方向脉动流速分量的乘积的统计平均值。乘积的统计平均值。图图 脉动流速分量示意图脉动流速分量示意图 xxa + D Dxxaiu u相应的相关系数为相应的相关系数为: : 对均匀紊流有对均匀紊流有: :均匀紊流的欧拉空间相关系数为均匀紊流的欧拉空间相关系数为: : (3-3-2)(3-3-2)(3-3-1)(3-3-1)222iaiaiu u u u u u )()(2iaiaiiRu u u u u u )()()()()()()(),(22 u u u u u u u u aiaiaiaiaiR当当等于零时，等于零时，Ri()

10、应等于应等于1；愈大，愈大，Ri()愈小愈小;当当超过一定的值，超过一定的值，Ri()渐趋于零渐趋于零(两点分别位于不同的涡体）。两点分别位于不同的涡体）。第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度图图 欧拉空间相关系数欧拉空间相关系数 Ri(x)第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度2iaiaiiRu u u u u u )()()(2 2、欧拉时间相关、欧拉时间相关定义为定义为: : 相应的相关系数为相应的相关系数为: : (3-3-3)(3-3-3)()( u uu u ttii)()()()(),( u uu u u uu u tttttRiiiii22如果

11、紊流在恒定流中发生，紊流场是平稳的，便有如果紊流在恒定流中发生，紊流场是平稳的，便有: :222iiittu uu u u u )()(所以恒定流的欧拉时间相关系数为所以恒定流的欧拉时间相关系数为: :(3-3-4)(3-3-4)2iiiittRu u u uu u )()()(指同一空间点，同一方向的脉指同一空间点，同一方向的脉动流速分量，在不同瞬时动流速分量，在不同瞬时( (相隔相隔时段为时段为 ) )的乘积的统计平均值。的乘积的统计平均值。第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度 3 3、欧拉紊流尺度（比尺）、欧拉紊流尺度（比尺）从紊流统计理论看，其空间点的脉动量可以视为各

12、种不从紊流统计理论看，其空间点的脉动量可以视为各种不同尺度同尺度( (或不同脉动频率或不同脉动频率) )的涡体经过该点所造成的涨落，的涡体经过该点所造成的涨落，较大尺度涡体包含着较小尺度涡体。较大尺度涡体包含着较小尺度涡体。大尺度涡体频率低，小尺度涡体频率高。大尺度涡体频率低，小尺度涡体频率高。由相关系数的概念，引入涡体的平均尺度由相关系数的概念，引入涡体的平均尺度( (积分尺度积分尺度) )。 第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度 3 3、欧拉紊流尺度（比尺）、欧拉紊流尺度（比尺）l对均匀紊流来说，取距离为对均匀紊流来说，取距离为的两点，如果涡体的平均的两点，如果涡体的平均

13、尺度较大，两点处于同一涡体，则空间相关系数尺度较大，两点处于同一涡体，则空间相关系数Ri()就大；就大；l如果涡体平均尺度小，两点分别处于两个涡体中，如果涡体平均尺度小，两点分别处于两个涡体中，Ri()就小。就小。lRi()与涡体平均尺度有密切关系。与涡体平均尺度有密切关系。(3-3-5)(3-3-5) dRLiEi 0)(第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度（1 1）欧拉空间平均尺度）欧拉空间平均尺度涡体空间在涡体空间在i方向上的空间平均尺度定义为：方向上的空间平均尺度定义为：图中假想的矩形面积与图中假想的矩形面积与Ri()曲线下的面积相等。曲线下的面积相等。LEi的意义是

14、体现了涡体尺度在的意义是体现了涡体尺度在i方向上的空间平均值。方向上的空间平均值。图图 欧拉空间平均尺度欧拉空间平均尺度L LEiEi第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度 dRLiEi 0)( (2) (2)欧拉时间平均尺度欧拉时间平均尺度 类似地，当紊流场是平稳的，可以用时间相关系数定义类似地，当紊流场是平稳的，可以用时间相关系数定义时间平均尺度：时间平均尺度：(3-3-6)(3-3-6) dRTiEi0)(第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度 从拉格朗日观点出发，一从拉格朗日观点出发，一个液体质点在运行过程中的脉个液体质点在运行过程中的脉动流速的相关矩定

15、义为：动流速的相关矩定义为： 指跟踪一个质点看，在不同时刻、同一方向的脉动流速分量指跟踪一个质点看，在不同时刻、同一方向的脉动流速分量的乘积的统计平均值。的乘积的统计平均值。 三、拉格朗日相关和紊流尺度三、拉格朗日相关和紊流尺度)()( u uu u ttii)(tiu)( u u ti第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度图图 拉格朗日相关流速拉格朗日相关流速分量示意图分量示意图 相应的相关系数为：相应的相关系数为： 如果紊流场是平稳的，上式变为：如果紊流场是平稳的，上式变为：三、拉格朗日相关和紊流尺度三、拉格朗日相关和紊流尺度)()()()()( u uu u u uu u

16、 ttttRiiiiLi222iiiLittRu u u uu u )()()(3-3-7)(3-3-7)(3-3-8)(3-3-8)第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度)(tiu)( u u ti图图 拉格朗日相关流速拉格朗日相关流速分量示意图分量示意图 在在 i 方向上，拉格朗日时间平均尺度的定义为：方向上，拉格朗日时间平均尺度的定义为：其中假想的以其中假想的以TLi为底的矩形面为底的矩形面积与积与RLi曲线下的面积相等。它曲线下的面积相等。它反映同一质点，不同时刻的随反映同一质点，不同时刻的随机变量之间保持有关所经历的机变量之间保持有关所经历的时间长度。时间长度。 dR

17、TLiLi0)(图图 拉格朗日时间平均尺度拉格朗日时间平均尺度(3-3-9)(3-3-9)第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度在在 i 方向上，拉格朗日时间平均尺度的定义为：方向上，拉格朗日时间平均尺度的定义为：拉格朗日空间平均尺度（或称为扩散平均尺度）定义为：拉格朗日空间平均尺度（或称为扩散平均尺度）定义为：LiiLiTv2 式中：式中： 称为称为i 方向的紊流强度。方向的紊流强度。2iv dRTLiLi 0)(3-3-9)(3-3-9)(3-3-10)(3-3-10)第三节第三节 紊流统计量和紊流尺度紊流统计量和紊流尺度第四节第四节 紊动扩散理论紊动扩散理论 v由于紊流脉

18、动流速的作用使污染物质发生输移的现象称由于紊流脉动流速的作用使污染物质发生输移的现象称为紊动扩散。为紊动扩散。在紊动扩散问题上的距离方差是否也具有与历时成正比的在紊动扩散问题上的距离方差是否也具有与历时成正比的同样规律呢？分子扩散系数同样规律呢？分子扩散系数D是反映由于液体分子运动使是反映由于液体分子运动使污染物质点发生位移时用于表达通量与浓度梯度成比例的污染物质点发生位移时用于表达通量与浓度梯度成比例的一个系数（费克定律），那么，由于紊动而使示踪质点发一个系数（费克定律），那么，由于紊动而使示踪质点发生输移时，是否还要使用类似的系数呢？生输移时，是否还要使用类似的系数呢？ 由于液体分子运动的

19、作用使污染物质点发生随机运动时的由于液体分子运动的作用使污染物质点发生随机运动时的距离方差与扩散历时成正比：距离方差与扩散历时成正比：s sx2 2 =2=2D 。紊流的脉动流速是随机性，因而使污染物质质点的位移也紊流的脉动流速是随机性，因而使污染物质质点的位移也是随机性的，它们的随机性与液体分子运动的随机性是否是随机性的，它们的随机性与液体分子运动的随机性是否相类似？相类似？v液体分子运动的速度和紊流脉动流速的随机特征是不同的。液体分子运动的速度和紊流脉动流速的随机特征是不同的。v液体分子运动的下一步状态，只与当前的状态有关，而与以液体分子运动的下一步状态，只与当前的状态有关，而与以前的运动

20、历史无关。前的运动历史无关。从前面介绍的拉格朗日从前面介绍的拉格朗日脉动流速相关系数脉动流速相关系数可以看到，示踪质可以看到，示踪质点的脉动速度在运动过程的一段时间内都是相关的。点的脉动速度在运动过程的一段时间内都是相关的。只有当经历的时间间隔只有当经历的时间间隔t TLi之后（或经历的位移之后（或经历的位移LLi之之后）才能认为此时的脉动流速与后）才能认为此时的脉动流速与t 之前的脉动流速无关。之前的脉动流速无关。第四节第四节 紊动扩散理论紊动扩散理论一、单个质点的紊动扩散一、单个质点的紊动扩散 TTdtdyydtdy)2()(2 泰勒（泰勒（TaylorTaylor）19211921年首先

21、应用年首先应用统计理论统计理论和和拉格朗日拉格朗日方法方法来研究单个质点的紊动扩散问题。来研究单个质点的紊动扩散问题。(3-4-1)(3-4-1) 不失一般性，讨论水流以时间平均流速不失一般性，讨论水流以时间平均流速 沿沿 x 方向作均方向作均匀流动时，示踪质点由于紊动在匀流动时，示踪质点由于紊动在 y 方向的扩散。如果在速度方向的扩散。如果在速度为为 的动坐标系统上观察，质点由于在的动坐标系统上观察，质点由于在 y 方向上脉动速度方向上脉动速度u u 的作用，在的作用，在 y 方向上的随机位移为方向上的随机位移为y ，则在某一指定时刻，则在某一指定时刻T, y2随时间随时间t 的变化率为的变

22、化率为uuuxyv第四节第四节 紊动扩散理论紊动扩散理论图图 单个质点的紊动扩散单个质点的紊动扩散 )()(TdtdyyTu u 因为：因为：有有将上两式代入式：将上两式代入式：上式是对一个示踪质点而论的上式是对一个示踪质点而论的TTdtdyydtdy)2()(2 dttTyTy)()( u u 0dtTtdtdyyTyT)()(2)(02u uu u 第四节第四节 紊动扩散理论紊动扩散理论TTdtdyydtdy)2()(2 有有22)()()()()(yyyyyyLyttTtRu u u uu uu uu uu u dtTtdtydyTyT)()(2)(02uu如果观察大量示踪质点，然后取总

23、体平均，则有：如果观察大量示踪质点，然后取总体平均，则有：(3-4-2)(3-4-2)设：时间间隔设：时间间隔 =T-t y y(t t) )y y(T T ) )为同一个质点在时间间隔为为同一个质点在时间间隔为 的两个脉动速的两个脉动速度的乘积度的乘积 为拉格朗日自相关矩为拉格朗日自相关矩 在平稳的紊流场中，相应的相关系数为在平稳的紊流场中，相应的相关系数为:)()(Ttyyu uu u 第四节第四节 紊动扩散理论紊动扩散理论式中的式中的t T 。它表示在平稳的紊流场中扩散距离方差与它表示在平稳的紊流场中扩散距离方差与 成正比，也与扩成正比，也与扩散历时散历时t 有关。有关。对其他方向对其他

24、方向 x 和和 z 也有类似公式。也有类似公式。2yu u TyLyTdRdtyd022)(2)( u u故式故式(3-4-2)(3-4-2)可变为：可变为：(3-4-3)(3-4-3)上式的上式的T 也可理解为任一指定时刻也可理解为任一指定时刻t，故有扩散距离的方差：，故有扩散距离的方差： u u s sTyLtyydRdTtyt00222)(2)()(3-4-4)(3-4-4)称为泰勒扩散公式称为泰勒扩散公式第四节第四节 紊动扩散理论紊动扩散理论（1） tTLy（TLy为为y方向上的拉格朗日时间平均尺度）方向上的拉格朗日时间平均尺度） 此时有此时有趋近于趋近于0， RLy 1：( 3-4-

25、5a )或或上式表明，对短的扩散历时（上式表明，对短的扩散历时（ tTLy）距离方差）距离方差i2与与t 2成正比，这是与分子扩散规律不同的，属于成正比，这是与分子扩散规律不同的，属于非费克型扩散非费克型扩散。( 3-4-5b )s sx2与扩散历时与扩散历时t t成正比的扩散成正比的扩散费克型扩散费克型扩散两种特殊情况：两种特殊情况：2222)()(ttytyyu us s ttytyy22)()(u us s 第四节第四节 紊动扩散理论紊动扩散理论 u u s sTyLtyydRdTtyt00222)(2)()(（2 2）t TLy此时有：此时有：LyLytLyTdRdR 00)()(或或

26、( 3-4-6a )( 3-4-6b ) u u s sTyLtyydRdTtyt00222)(2)()(LyyytTtyt2222u us s )()(tTtytLyyy222u us s )()(第四节第四节 紊动扩散理论紊动扩散理论对较大的扩散历时（对较大的扩散历时（tTLy），距离方差与扩散历时成正），距离方差与扩散历时成正比。这样的紊动扩散规律与分子扩散规律相同，符合马尔可比。这样的紊动扩散规律与分子扩散规律相同，符合马尔可夫过程，为费克型扩散。夫过程，为费克型扩散。可以定义一个与分子扩散系数类似的紊动扩散系数可以定义一个与分子扩散系数类似的紊动扩散系数: ),(zyxi 这也是求这

27、也是求Ei的矩法公式。如果的矩法公式。如果Ei是常数，可将上式改写为是常数，可将上式改写为tEii22 s sdtdEii221s s 式（式（3-4-7a3-4-7a）可写为：）可写为：( 3-4-7b )( 3-4-7a )( 3-4-8 )LiiiTvE2 LyyytTtyt2222u us s )()(第四节第四节 紊动扩散理论紊动扩散理论LiiiE 2/12)(u u 上式表明：紊动扩散系数与拉格朗日时间平均尺度成正比。上式表明：紊动扩散系数与拉格朗日时间平均尺度成正比。可以认为紊动扩散系数主要与较大尺度的涡体运动有关。可以认为紊动扩散系数主要与较大尺度的涡体运动有关。 若以拉格朗日

28、空间平均尺度若以拉格朗日空间平均尺度Li表示紊动扩散系数，有表示紊动扩散系数，有:( 3-4-9 )LiiiTvE2 第四节第四节 紊动扩散理论紊动扩散理论分子扩散系数由污染物质的物理性质决定，紊动扩散系分子扩散系数由污染物质的物理性质决定，紊动扩散系数与液体的物理性质、污染物质的物理性质和流场的紊数与液体的物理性质、污染物质的物理性质和流场的紊动结构有关。动结构有关。在河流中，最大涡体的尺寸约与水深相等，由于河流的在河流中，最大涡体的尺寸约与水深相等，由于河流的水深一般都不大，所以水深一般都不大，所以拉格朗日空间平均尺度拉格朗日空间平均尺度也不大，也不大，相应的拉格朗日时间平均尺度也较小，使

29、污水自排污口相应的拉格朗日时间平均尺度也较小，使污水自排污口注入河流之后不久，其扩散历时就远大于注入河流之后不久，其扩散历时就远大于T TL L，因此，河，因此，河流中的紊动扩散一般是属于费克型。流中的紊动扩散一般是属于费克型。？第四节第四节 紊动扩散理论紊动扩散理论表表1 质点紊动扩散的实验数据质点紊动扩散的实验数据 smu/305. 0 例例1 在在x方向有时均速度方向有时均速度 ，量测了许多示踪质点的，量测了许多示踪质点的y 向位移，得相应于不同时间的总体平均值向位移，得相应于不同时间的总体平均值 ，见表，见表1。试对。试对 ，TLy和和Ey进行估计。进行估计。)(2ty2u u st/

30、 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.06 0.23 0.53 0.93 1.44 2.00 2.59 3.19 3.78 4.38 0.245 0.480 0.728 0.964 1.20 1.41 1.61 1.79 1.94 2.09 24210/ )(mty mty2210/)( 第四节第四节 紊动扩散理论紊动扩散理论 图图3-11 3-11 质点紊动扩散的实验数据曲线质点紊动扩散的实验数据曲线第四节第四节 紊动扩散理论紊动扩散理论解：由图解：由图3-113-11可以看到，当可以看到，当 与与t t为线性关系。为线性关系。 根据式根据式(

31、3-4-5b)(3-4-5b)，用，用t t =0.1s=0.1s的值代入计算：的值代入计算：)(1 . 02tyst， )/(0245. 01 . 0/10245. 0/222smty u u22y2y第四节第四节 紊动扩散理论紊动扩散理论图图3-11 3-11 质点紊动扩散的实验数据曲线质点紊动扩散的实验数据曲线根据式根据式(3-4-9)(3-4-9)，有：，有：)/(100 . 67 . 00 . 110)59. 238. 4(2442smdtyd 根据式根据式(3-4-6)(3-4-6)，得：，得： 由图由图3-113-11看到，当看到，当 线性增加，线性增加，用用t t =0.7s=

32、0.7s和和1.0s1.0s的的 值计算：值计算：)(7 . 02tyst， )( 5 . 00245. 02100 . 6212422 sdtydTLy u u)/(1035 . 0)0245. 0(2422 smTELyLy u u22y2y2y第四节第四节 紊动扩散理论紊动扩散理论图图3-11 3-11 质点紊动扩散的实验数据曲线质点紊动扩散的实验数据曲线二、两个质点的相对扩散二、两个质点的相对扩散v污染物云团可以认为是由大量污染物质点组成的。因此，污染物云团可以认为是由大量污染物质点组成的。因此，可以分别研究污染物云团质心的运动和云团的质点相对于可以分别研究污染物云团质心的运动和云团的

33、质点相对于此质心的运动。此质心的运动。v对质心的运动应按上述的单个质点的紊动扩散进行分析。对质心的运动应按上述的单个质点的紊动扩散进行分析。对云团的变形则要研究质心与其他质点的相对运动，也就对云团的变形则要研究质心与其他质点的相对运动，也就是两质点的相对扩散问题。是两质点的相对扩散问题。v白切勒（白切勒（BathelorBathelor）对相对扩散问题的研究做出了较大的）对相对扩散问题的研究做出了较大的贡献。贡献。第四节第四节 紊动扩散理论紊动扩散理论第五节第五节 随流紊动扩散方程随流紊动扩散方程 ( 3-5-1 )( 3-5-2 ) 对对紊流紊流讲，点瞬时速度讲，点瞬时速度i可表为点时均速度

34、可表为点时均速度 与点脉动与点脉动速度速度ui之和，即之和，即iuiiiu uu uu u ccc 在紊流中，由于流体的紊动使污染物质的点瞬时浓度在紊流中，由于流体的紊动使污染物质的点瞬时浓度 c 也具有随机性质。因此，可以假设点瞬时浓度也具有随机性质。因此，可以假设点瞬时浓度 c 是时均浓度是时均浓度 和脉动浓度和脉动浓度c 之和，即：之和，即：ciiiiixxccDxcctcc )()()()(2u uu uiiiiiiiixcxcxcxc )()(u uu uu uu u第五节第五节 随流紊动扩散方程随流紊动扩散方程 ( 3-5-4 )因为脉动流连续方程为因为脉动流连续方程为 ，便有，便

35、有: 0 iixu uiiiixxcDxctc u u 2 三维随流扩散方程是研究随流紊动扩散的基础，将三维三维随流扩散方程是研究随流紊动扩散的基础，将三维随流扩散方程简写为：随流扩散方程简写为：( 3-5-3 )iiiiiixcxxcDxctc )(2u uu u)()(22iiiiiiiiiiiixxcxxcDxcxcxcxctctc u uu uu uu u将上式代入将上式代入(3-5-4)(3-5-4)，有：，有：( 3-5-5 )运用运用雷诺运算法则雷诺运算法则，并注意到，并注意到 ，便得，便得: 0 u u c三维随流紊动扩散方程的初形三维随流紊动扩散方程的初形第五节第五节 随流紊

36、动扩散方程随流紊动扩散方程iiiiiixcxxcDxctc )(2u uu u)()(222222zcycxcDzcycxcutc u u比较比较第五节第五节 随流紊动扩散方程随流紊动扩散方程三维随流紊动扩散方程多了三个紊动通量梯度：三维随流紊动扩散方程多了三个紊动通量梯度： 上式包含了四个未知函数上式包含了四个未知函数 。为求解为求解 , ,需要对需要对 三个未知函数进行处理。三个未知函数进行处理。 iixc /u ucccczyx u uu uu u、cccczyx u uu uu u、式中式中: :E Eijij称为紊动扩散系数张量，是一个二阶张量。称为紊动扩散系数张量，是一个二阶张量。

37、 相对于紊动扩散张量的主轴来说：相对于紊动扩散张量的主轴来说： 当当i=j,i=j,E Eijij=E=Ejiji；当；当ij, ij, E Eijij= =0 0 则则( 3-5-6 )( 3-5-6 )式变为：式变为：jijixcEc u u 最常用的方法是采用经验模式，对最常用的方法是采用经验模式，对t t T TL L的扩散，把紊的扩散，把紊动扩散与分子扩散相比拟，参照动扩散与分子扩散相比拟，参照费克定律费克定律的形式，令的形式，令: :( 3-5-6 )iiiixcEc u u ( 3-5-7 )第五节第五节 随流紊动扩散方程随流紊动扩散方程)(2iiiiiiiixcExxxcDxc

38、tc u u 紊动扩散系数紊动扩散系数E Eijij与液体的物理性质、污染物质的物理性质与液体的物理性质、污染物质的物理性质以及流场的紊动结构有关。因为剪切紊流不是各向同性的，所以及流场的紊动结构有关。因为剪切紊流不是各向同性的，所以它的值在不同方向上是不同的。以它的值在不同方向上是不同的。 各种情形下的各种情形下的E Eiiii值主要靠值主要靠实验实验确定。确定。将紊动扩散系数的定义：将紊动扩散系数的定义： 第五节第五节 随流紊动扩散方程随流紊动扩散方程代入式代入式iiiixcEc u u iiiiiixcxxcDxctc )(2u uu u有：有：)(2iiiiiiiixcExxxcDxc

39、tc u u 在紊流扩散中，一般的紊动扩散系数都要比分子扩散系在紊流扩散中，一般的紊动扩散系数都要比分子扩散系数大几个数量级，故在紊流中可数大几个数量级，故在紊流中可忽略分子扩散项忽略分子扩散项，得：，得： )(iiiiiixcExxctc u u( 3-5-8 )第五节第五节 随流紊动扩散方程随流紊动扩散方程如果紊流场是均匀紊流，则如果紊流场是均匀紊流，则E Eiiii为常数，即不随位置和时为常数，即不随位置和时间而变间而变, , 可将可将E Eiiii简写为简写为E Ei i。22iiiixcExctc u u( 3-5-9a )将上式改用直角坐标表示为：将上式改用直角坐标表示为： 222

40、222zcEycExcEzcycxcutzyx u uu u式中：式中： 分别是点时均流速在分别是点时均流速在x x、y y和和z z方向上的分量。方向上的分量。 u u、u( 3-5-9b )第五节第五节 随流紊动扩散方程随流紊动扩散方程)(iiiiiixcExxctc u uiiiixxcExct 2u uu u)(222222zcycxcEzcycxcut u uu u( 3-5-10a )( 3-5-10b )以上结果均是针对示踪物质而得到的以上结果均是针对示踪物质而得到的如果紊流场是各向同性的，有如果紊流场是各向同性的，有 ，EEEEzyx 第五节第五节 随流紊动扩散方程随流紊动扩散

41、方程随流紊动扩散方程为随流紊动扩散方程为: :BzcEycEycEzcycxcutczyx 222222 u uBzcycxcEzcycxcutc )(222222 u u 如果进一步考虑非守恒物质，则要计算由于化学、物理等如果进一步考虑非守恒物质，则要计算由于化学、物理等各种原因使水体中的非守恒物质发生变化（各种原因使水体中的非守恒物质发生变化（增生或降解增生或降解）。）。 式中：当式中：当B为增生时，为增生时，B是正值（称为源项）；是正值（称为源项）； 当当B为降解时，为降解时，B是负值（称为汇项）。是负值（称为汇项）。( 3-5-11 )( 3-5-12 )第五节第五节 随流紊动扩散方程

42、随流紊动扩散方程 设在单位时间单位水体积内由于上述原因使非守恒物质发设在单位时间单位水体积内由于上述原因使非守恒物质发生的变化量为生的变化量为B( (称为称为源汇项源汇项) )，则随流紊动扩散方程为：，则随流紊动扩散方程为：第六节第六节 紊动扩散系数的确定紊动扩散系数的确定 一、雷诺比拟假说一、雷诺比拟假说 前面讨论的扩散都是指污染物质在水中的扩散；从广义前面讨论的扩散都是指污染物质在水中的扩散；从广义上说，水流本身的某些属性由于上说，水流本身的某些属性由于分子运动分子运动、水流流速水流流速和和水流水流紊动紊动而传递到另一部分水体中去的现象都可看成是扩散现象而传递到另一部分水体中去的现象都可看

43、成是扩散现象, , 例如水流本身具有的动量、动能和热量（温度）的扩散。因例如水流本身具有的动量、动能和热量（温度）的扩散。因此，对它们的扩散过程也可以用类似的扩散方程表示。此，对它们的扩散过程也可以用类似的扩散方程表示。 ( 3-6-1 )iixxTtT 2 iiiixxcDxctc 2u u 例如在随流紊动扩散方程中，以温度例如在随流紊动扩散方程中，以温度T T 代替浓度代替浓度c c ，以，以导热系导热系数导热系导热系数代替分子扩散系数代替分子扩散系数D，并有，并有vi=0=0，则有热传，则有热传导方程：导方程：zcEczz u uzxzx u u u u u u 上式是由于紊动上式是由于

44、紊动 z 方向发生动量传递而产生的紊动切应方向发生动量传递而产生的紊动切应力表示式，从而也可以把力表示式，从而也可以把 看作是动量扩散系数。因此，也看作是动量扩散系数。因此，也有人在浓度扩散的计算中采用有人在浓度扩散的计算中采用 值作为值作为Ei值。值。( 3-6-2 )第六节第六节 紊动扩散系数的确定紊动扩散系数的确定在上式中，若以主流方向（在上式中，若以主流方向（x 向）的动量向）的动量 ( 为水的密度为水的密度 )代替代替 ，以，以 u ux代替代替c，以运动粘性系数，以运动粘性系数 代替代替Ez，则有：，则有：xu u c如在式如在式 中，对垂直方向中，对垂直方向( (z z向向) )

45、的扩散有：的扩散有：jijixcEc u u 雷诺比拟假说：雷诺比拟假说：对示踪物质来说对示踪物质来说, 不论是哪种扩散物质，不论是哪种扩散物质，紊动对物质、热量、动能或动量的扩散系数存在着完全紊动对物质、热量、动能或动量的扩散系数存在着完全的比拟关系，扩散系数彼比都是相等的。已在近壁紊流的比拟关系，扩散系数彼比都是相等的。已在近壁紊流中的实验得到证实。中的实验得到证实。第六节第六节 紊动扩散系数的确定紊动扩散系数的确定l在明渠二度均匀流中，取垂向紊动扩散系数在明渠二度均匀流中，取垂向紊动扩散系数Ez= 是符是符合实际的。合实际的。l其他方向的紊动扩散系数，对剪切紊流讲，只能认为其他方向的紊动

46、扩散系数，对剪切紊流讲，只能认为Ei 与与 是同量级的，是同量级的， 成一定的比例关系。成一定的比例关系。一些研究成果表明：一些研究成果表明：二、河渠水流的紊动扩散系数二、河渠水流的紊动扩散系数 1、垂向紊动扩散系数、垂向紊动扩散系数 对二度明渠均匀流的纵向流速分布，对二度明渠均匀流的纵向流速分布， 可采用卡门可采用卡门(Karman)对数型公式：对数型公式：)ln1(hzkuuu 式中：式中： 为点时均流速为点时均流速 在垂线上的平均流速；在垂线上的平均流速； k为卡门常数；为卡门常数； u*为剪切流速；为剪切流速； h为水深，为水深，z 轴的原点在渠底，向上为正。轴的原点在渠底，向上为正。

47、uu( 3-6-3 )表示沿垂线取平均表示沿垂线取平均第六节第六节 紊动扩散系数的确定紊动扩散系数的确定 第六节第六节 紊动扩散系数的确定紊动扩散系数的确定dzuddzud 或或）（1)1(22hzu ）（( 3-6-4 )( 3-6-5 )切应力的两种表达方式：切应力的两种表达方式：dzud ）（dzud / )ln1(hzkuuu zkudzud1* )1(hzzku )1(2hzu 6/1|hungghsu 根据雷诺比拟假说，可得：根据雷诺比拟假说，可得： 式中：式中：k k=0.40=0.400.420.42。式中：式中：s 为河道底坡；为河道底坡；n为河道糙率；为河道糙率；g为重力加

48、速度。为重力加速度。)1(hzzkuEz ( 3-6-6 )( 3-6-7 )第六节第六节 紊动扩散系数的确定紊动扩散系数的确定u u*可由下式计算可由下式计算: :Ez是水深的函数，沿水深各点不同是水深的函数，沿水深各点不同 hukdzhzzkuhdzEhEhhzz6)1(11|00 Ez沿水深取平均沿水深取平均: 取卡门常数取卡门常数k=0.41，有：，有： huEz068.0|( 3-6-8 )第六节第六节 紊动扩散系数的确定紊动扩散系数的确定二度明渠均匀流垂向紊二度明渠均匀流垂向紊动扩散系数计算公式动扩散系数计算公式 对非矩形断面的棱柱体河道均匀流来说，由于垂向混合主对非矩形断面的棱柱

49、体河道均匀流来说，由于垂向混合主要受水深制约，因而认为也可用式（要受水深制约，因而认为也可用式（3-6-63-6-6）求，此时式）求，此时式中的中的h是随横坐标是随横坐标y而变的。而变的。对各种不规则河道的恒定渐变流，河底和边坡都可能沿程对各种不规则河道的恒定渐变流，河底和边坡都可能沿程变化，也可能有弯道，二次流比较显著。但一般认为，这变化，也可能有弯道，二次流比较显著。但一般认为，这些因素对垂向混合的影响不会很大，仍可近似使用式（些因素对垂向混合的影响不会很大，仍可近似使用式（3-3-6-66-6）求，只是其中的）求，只是其中的h是随点（是随点（x，y）而变的。）而变的。 第六节第六节 紊动

50、扩散系数的确定紊动扩散系数的确定)1(hzzkuEz （3-6-63-6-6）2、横向紊动扩散系数、横向紊动扩散系数Ey （横向混合系数（横向混合系数My ） 在二度明渠均匀流中，因为假设不存在横向流速，不在二度明渠均匀流中，因为假设不存在横向流速，不可能建立类似于垂向紊动扩散系数可能建立类似于垂向紊动扩散系数的方程的方程(3-6-6)(3-6-6) 来确定来确定Ey值。值。 对其他非二度水流，由于影响因素复杂，目前也不可对其他非二度水流，由于影响因素复杂，目前也不可能像垂向扩散一样通过沿河宽的横向流速分布来建立能像垂向扩散一样通过沿河宽的横向流速分布来建立Ey的的计算式，而只能根据室内和现场

51、实验的结果大致给出一个计算式，而只能根据室内和现场实验的结果大致给出一个变化范围。变化范围。第六节第六节 紊动扩散系数的确定紊动扩散系数的确定2、横向紊动扩散系数、横向紊动扩散系数Ey （横向混合系数（横向混合系数My ） lEy的意义是反映由于横向脉动流速作用而导致的扩散的强的意义是反映由于横向脉动流速作用而导致的扩散的强弱程度。弱程度。 l在实际的河渠中，即使是平直的渠道均匀流，也必然有二在实际的河渠中，即使是平直的渠道均匀流，也必然有二次环流次环流 (螺旋流螺旋流)，虽然二次环流的，虽然二次环流的 =0，但因，但因 在垂线在垂线上的分布不均匀而导致了分散，使横向扩散大大加强，这上的分布不

52、均匀而导致了分散，使横向扩散大大加强，这是横向扩散的主因。是横向扩散的主因。l由于在实验和分析中都以将这两种成因的结果分开，得到由于在实验和分析中都以将这两种成因的结果分开，得到的系数值是综合结果。故对河渠水流来说，常将的系数值是综合结果。故对河渠水流来说，常将Ey称为横称为横向混合系数向混合系数My。uu第六节第六节 紊动扩散系数的确定紊动扩散系数的确定（1 1）矩法）矩法122122222xxVdxdVMyyyy s ss ss s D D D D iiiiiiWWwWyycyycdycdyyc|22/2/2/2/22s s第六节第六节 紊动扩散系数的确定紊动扩散系数的确定一般只能近似地一

53、般只能近似地用于比较规则的用于比较规则的平直河流。平直河流。式中：式中： V为断面平均流速；为断面平均流速； 、 分别为分别为x1和和x2断面处横向浓度分布的方差。断面处横向浓度分布的方差。 22ys s21ys s式中：式中：W为河宽。为河宽。（1 1）矩法）矩法 对不规则的天然河流，对不规则的天然河流，h和和 有横向变化，有横向变化， ，这些，这些都对横向扩散有重要影响，特别是必须对已包含在实验结果都对横向扩散有重要影响，特别是必须对已包含在实验结果中由于中由于 产生的随流扩散加以扣除，使产生的随流扩散加以扣除，使My主要反映二次环主要反映二次环流和横向脉动流速等作用，否则，得出的值会很大

54、或很小，流和横向脉动流速等作用，否则，得出的值会很大或很小，甚至出现负值。为此，甚至出现负值。为此，Holley等提出了等提出了普遍矩法普遍矩法。 u0 u uu u 在进行示踪实验时，当浓度场达到稳定，可导出普遍矩在进行示踪实验时，当浓度场达到稳定，可导出普遍矩法的基本方程：法的基本方程： wWywWwWdycuhydyychMdycuhydycvhdycuhdyycuhdxd0000002|2|2|第六节第六节 紊动扩散系数的确定紊动扩散系数的确定（2）经验公式）经验公式 在整理实验结果时，一般采用在整理实验结果时，一般采用无量纲数无量纲数：式中：式中： h 和和 u*分别为被研究河段的水

55、深和剪切流速的平均值分别为被研究河段的水深和剪切流速的平均值; 为比例系数，视河渠水流条件而定。为比例系数，视河渠水流条件而定。( 3-6-9 ) huMy第六节第六节 紊动扩散系数的确定紊动扩散系数的确定（2）经验公式）经验公式 费希尔对费希尔对顺直的矩形断面明槽顺直的矩形断面明槽水流收集了约水流收集了约75个实验个实验资料，除了在野外实验得到的资料，除了在野外实验得到的=0.240.25 之外，几乎所有之外，几乎所有的室内实验都得到的室内实验都得到=0.10.2，于是取其平均作为估算，有，于是取其平均作为估算，有估算式估算式: ( 3-6-9 )( 3-6-10 )15. 0 huMy h

56、uMy6/1RnVggRsu ( 3-6-11 )式中：式中：u*为摩阻流速；为摩阻流速；第六节第六节 紊动扩散系数的确定紊动扩散系数的确定 对天然河流的对天然河流的My值虽然做了一些研究，值虽然做了一些研究， 但仍然不够充分。但仍然不够充分。费希尔总结了一些实验成果，认为费希尔总结了一些实验成果，认为河道的弯曲河道的弯曲和和两岸的不规则两岸的不规则都会使都会使 值很少值很少于于0.40.4，如果河道的弯曲是缓慢的，边岸的不规，如果河道的弯曲是缓慢的，边岸的不规则性是中等的，则性是中等的， 值通常为值通常为0.8，因而建议采用估算式：因而建议采用估算式：( 3-6-12 )式

57、中：式中：h 和和 u* 用河段平均值。用河段平均值。%506 . 0* huMy第六节第六节 紊动扩散系数的确定紊动扩散系数的确定如果河段有急弯或边岸有突变，则如果河段有急弯或边岸有突变，则 值将更高。值将更高。张永兴等人在长江宜昌段进行示踪实验，得张永兴等人在长江宜昌段进行示踪实验，得 =0.290.290.950.95。 Lau Lau 等人曾总结分析了等人曾总结分析了1111条天然条天然河流的资料，为了更好地河流的资料，为了更好地反映二次环流对反映二次环流对My的影响，采用无量纲数的影响，采用无量纲数My /(u*W) ，并认为，并认为该数与宽深比该数与宽深比W/h有关有关( (u*

58、、W和和h用河段平均值），得到了图用河段平均值），得到了图3-123-12，图中有，图中有9个点都落于一条平均曲线的附近，这些点的个点都落于一条平均曲线的附近，这些点的My /(u*W) 值在值在1.81.81010-3-310.610.61010-3-3之间，之间， W/h值在值在2727170170之间之间；另有两个点因为河流弯曲较大，使；另有两个点因为河流弯曲较大，使My /(u*W)值较高（分别为值较高（分别为22.422.41010-3-3和和55.855.81010-3-3）。）。图图3-12 3-12 河流的河流的M My y经验曲线经验曲线第六节第六节 紊动扩散系数的确定紊动扩

59、散系数的确定3、纵向紊动扩散系数、纵向紊动扩散系数Ex 由于紊动引起的纵向扩散和横向扩散受边界的约束作由于紊动引起的纵向扩散和横向扩散受边界的约束作用相对垂向来说小很多，所以纵向和横向的紊动扩散系数用相对垂向来说小很多，所以纵向和横向的紊动扩散系数有有相同的数量级相同的数量级，这已为某些实验证实。，这已为某些实验证实。lSayre等用聚乙稀质点在水面上进行观察，发现纵向紊动等用聚乙稀质点在水面上进行观察，发现纵向紊动扩散系数约为横向紊动扩散系数的三倍。扩散系数约为横向紊动扩散系数的三倍。l从实用上看，可以将纵向紊动扩散作用忽略不计。这是从实用上看，可以将纵向紊动扩散作用忽略不计。这是因为在常见

60、的河渠水流中，由于纵向流速分布的不均匀因为在常见的河渠水流中，由于纵向流速分布的不均匀而产生的随流分散作用比纵向紊动扩散作用要大得多，而产生的随流分散作用比纵向紊动扩散作用要大得多，同时这两种作用是混在一起出现，在对实验进行分析时，同时这两种作用是混在一起出现，在对实验进行分析时，难以将两者分开。难以将两者分开。第六节第六节 紊动扩散系数的确定紊动扩散系数的确定3、纵向紊动扩散系数、纵向紊动扩散系数Ex l过去很少对纵向紊动扩散系数进行研究，而将注意力过去很少对纵向紊动扩散系数进行研究，而将注意力放在研究纵向随流分散及其分散系数上。纵向分散系放在研究纵向随流分散及其分散系数上。纵向分散系数是纵