环境水利学第4章剪切流动的分散(6).pptVIP

第四章 剪切流动的分散 在第三章中，当求随流扩散方程和随流紊动扩散方程的解析解时，为了能容易得解答，都假定时均流速是均匀分布的。 由于分散只是一种特殊情况的随流扩散，所以在原则上可以通过随流扩散方程或随流紊动扩散方程求解。前人在这方面已做过许多研究，但只有一些简单的情况才能得到解析解。 第一节 一维纵向剪切流的分散 第二节 圆管均匀流的纵向分散 第三节 明渠二度均匀流的纵向分散 埃尔德(Elder)1959年应用泰勒分析方法对二维明渠均匀流在紊流状态下的纵向分散问题进行了研究。 例：设有两块平行板，间距为h，如图所示。上下板分别以速度U/2向右、向左移动, 其间为层流，速度为 ， 在某断面上瞬时投放均匀线源示踪质，经过较长时间后，推导纵向分散系数K的计算公式。 ( 4-2-25 ) 以上推证中未考虑纵向紊动扩散。 根据泰勒求得的纵向紊动扩散系数为： ( 4-2-26) 将纵向分散系数和纵向紊动扩散系数同时计入，纵向分散系数变为： ( 4-2-27 ) 均匀紊流状态下的纵向分散系数： 第二节 圆管均匀流的纵向分散 式中：x、y、z分别为纵向、横向和竖向坐标。 一维随流紊流扩散方程： ( 4-3-1 ) 对二度明渠流有?/ ?y=0，因为纵向的紊动扩散项远小于纵向随流项，故可将纵向的紊动扩散项忽略，于是上式简化为： 设 u = ub + V 代入式 得： 取坐标变换： ，有： 第三节 明渠二维均匀流的纵向分散 ( 4-3-2 ) ( 4-3-3 ) ( 4-3-4 ) 将上两式代入式(4-3-2)，有： ( 4-3-5 ) 上式是Elder在研究二度均匀流的纵向分散问题时作为一个假设而直接给出的。 当扩散（即使是瞬时源的扩散）经过足够长的时间之后,站在动坐标看，c 随时间的变化很慢，有： 在 情况下： 0 / = ? ? t c ( 4-3-6 ) 近似为零 第三节 明渠二维均匀流的纵向分散 式中: h为明渠的水深; z从水面向下计算。 则 可变为： 令无量纲坐标： ( 4-3-7 ) ( 4-3-8 ) 式（4-3-8）表明了经过足够长的扩散历时之后 , 垂向的紊动扩散作用与纵向的分散作用相平衡。 以 a b C c c + = 代入上式, 因 , 并假设 ，有： 第三节 明渠二维均匀流的纵向分散 上式满足边界条件：当h =0或1，均有 。 对 进行积分, 并将 改写为 , 得： 0 / = ? ? h b c ( 4-3-9 ) 对二维明渠均匀流 ( 4-3-10 ) 将式(4-3-9)代入上式，得： ( 4-3-11 ) 对层流： x ? ? / x ? ? / 第三节 明渠二维均匀流的纵向分散 Elder对二度明渠均匀流采用如下的对数流速分布： ( 4-3-12 ) 式中：k为卡门常数; um为最大流速。 根据雷诺比拟，类似于式 ， 有垂向紊动扩散系数： ( 4-3-13 ) 第三节 明渠二维均匀流的纵向分散 ( 4-3-14 ) 通过流速分布公式计算断面平均流速： ( 4-3-15 ) 流速分布公式 第三节 明渠二维均匀流的纵向分散 对上式右边的积分按γ函数的级数计算，其值约为0.4041；取κ=0.41，便有: 偏离速度： ( 4-3-16 ) ( 4-3-17 ) ( 4-3-18 ) 以上结果是在未考虑纵向紊动扩散的情况下取得的K值。 第三节 明渠二维均匀流的纵向分散 在补充考虑纵向紊动扩散情况下： 由垂向紊动扩散沿水深的平均值 。假设可按各向同性紊动处理，有： * = hu E z 068 . 0 || || 纵向分散系数变为： * * * = + = hu hu hu K 93 . 5 068 . 0 86 . 5 ( 4-3-19 ) ( 4-3-20 ) 以上对Ex的估值可能偏低， Elder在实验室试验得Ey=0.23hu＊, (见第三章第六节)，按各向同性紊动考虑有Ex=Ey，则： ( 4-3-21 ) 第三节 明渠二维均匀流的纵向分散 对一些近似矩形断面的宽浅河道，