环境水力学-射流、羽流及浮射流

1、3/22/20221环境水力学射流、羽流及浮射流射流、羽流及浮射流环境工程教研室郑天柱3/22/20222射流、羽流及浮射流基本概念基本概念平面淹没紊动射流平面淹没紊动射流圆形淹没紊动射流圆形淹没紊动射流静止液体中的浮力羽流静止液体中的浮力羽流静止均质及线性密度分层环境中圆形浮射流静止均质及线性密度分层环境中圆形浮射流流动环境中的紊动射流流动环境中的紊动射流3/22/202231、基本概念 定义：定义： 射流（射流（Jet）：一股流体从几何尺寸远小于接纳）：一股流体从几何尺寸远小于接纳流体所占空间尺寸的喷口注入接纳流体，并与流体所占空间尺寸的喷口注入接纳流体，并与之混合的流动状态，叫做射流。之

2、混合的流动状态，叫做射流。 实际问题多为紊动射流。 羽流（羽流（Plume）：射流的初始动量很小，靠浮力）：射流的初始动量很小，靠浮力的作用使其混合，其扩散云团形似羽毛漂浮在的作用使其混合，其扩散云团形似羽毛漂浮在空中，故称羽流。如废水泄入含盐量大的海水；空中，故称羽流。如废水泄入含盐量大的海水；废热水排入河流。废热水排入河流。 浮射流（浮射流（Buoyancy jet）：是兼受动量和浮力两）：是兼受动量和浮力两种作用而运动的射流。种作用而运动的射流。3/22/20224 从水力学的观点来看，浮射流产生的稀释扩从水力学的观点来看，浮射流产生的稀释扩散与环境水流被动的紊流扩散规律有一定的散与环境

3、水流被动的紊流扩散规律有一定的差别。差别。 因为后者的流速场与浓度场无关，并且许多问因为后者的流速场与浓度场无关，并且许多问题是与固体中热传导规律相似。浮射流排放问题是与固体中热传导规律相似。浮射流排放问题，则必须把流速场和浓度场耦合起来联立求题，则必须把流速场和浓度场耦合起来联立求解。解。3/22/20225基本概念 分类分类 按射流周围环境边界条件来分：按射流周围环境边界条件来分： 自由射流：射流射入无限空间时，称为自由射流。自由射流：射流射入无限空间时，称为自由射流。 有限空间射流：射流射入有限空间时，称为非自由射有限空间射流：射流射入有限空间时，称为非自由射流。其中：流。其中： 若射流

4、是沿着固壁发展的，叫做附壁射流；若射流是沿着固壁发展的，叫做附壁射流； 沿水体自由水表面发展的，叫做表面射流。沿水体自由水表面发展的，叫做表面射流。 按射流进入的流体介质来分：按射流进入的流体介质来分： 淹没射流：射流射入性质一样的同种流体中。淹没射流：射流射入性质一样的同种流体中。 非淹没射流非淹没射流3/22/20226 另外：若喷口是一很长的狭缝，则可以处理成平面射流；喷口是一圆孔，则为孔口射流。 喷口法线方向沿铅直向的射流，叫做铅直射流；沿水平向的，叫做水平射流；与水平向呈一角度的，叫做倾斜射流。 3/22/202272、平面淹没紊动射流 平面淹没紊动射流又称二维紊动射流平面淹没紊动射

5、流又称二维紊动射流 假设一个起始流速假设一个起始流速u0为的长孔紊动射流，长孔宽为的长孔紊动射流，长孔宽度为度为B，浓度为，浓度为C0的污染物射入同密度的静止环的污染物射入同密度的静止环境水体中，如图所示。境水体中，如图所示。静水中平面射流紊动扩散3/22/202282、平面淹没紊动射流BC0 u0umbuy初始段主体段3/22/202292、平面淹没紊动射流 通过试验观测可得平面淹没紊动射流的一般特性：通过试验观测可得平面淹没紊动射流的一般特性： 射流断面上的静止压力近似为常数且可假设等于外部射流断面上的静止压力近似为常数且可假设等于外部压力。压力。 靠近排出口，靠近排出口，x5.2B的范围

6、内是流动主体段。该段的紊动已完的范围内是流动主体段。该段的紊动已完全透入射流核心区，流速和浓度分布都是自相似全透入射流核心区，流速和浓度分布都是自相似的，可表示为的，可表示为 式中，式中，um,Cm分别表示射流带轴线的流速和浓度值。分别表示射流带轴线的流速和浓度值。2b为射流带宽度；为射流带宽度；为浓度与速度宽度之比，一般为常数。为浓度与速度宽度之比，一般为常数。 射流的横向尺度与纵向尺度相比是较小的，射流的横向尺度与纵向尺度相比是较小的，b(x)x，并且，射流在横向上的扩展是线性的，并且，射流在横向上的扩展是线性的，即射流边界扩展角度为常数。即射流边界扩展角度为常数。 BxbyxCyxCby

7、xuyxumm2 . 5,)(exp)(),()(exp)(),(223/22/202212断面流速分布 射流轴线流速射流轴线流速um沿程变化可据动量守恒得到，沿程变化可据动量守恒得到，将上式代入下式将上式代入下式 可得半经验公式：可得半经验公式： 断面流速分布公式：断面流速分布公式：mdyuudmM200228. 2uxbum)(exp228. 2200byuxbu3/22/202213断面浓度分布 同样，根据物质守恒原理，可得断面浓度分同样，根据物质守恒原理，可得断面浓度分布公式：布公式：)(exp234. 2200byCxbC3/22/2022143、圆断面紊动射流 二维紊动射流的许多试

8、验特征也适用于圆形二维紊动射流的许多试验特征也适用于圆形孔口的射流情况孔口的射流情况 ： 扩散厚度呈线性扩展；扩散厚度呈线性扩展； 静水压力近似为常数；静水压力近似为常数； 对此种轴对称情况，主流核心区长度一般为对此种轴对称情况，主流核心区长度一般为6.2D，D为圆孔直径；为圆孔直径； 3/22/202215 流速和浓度分布为流速和浓度分布为 流动形成段，流动形成段，x6.2D 式中，式中，Cm为轴心浓度，它与射流出口浓度有为轴心浓度，它与射流出口浓度有如下关系。如下关系。DxbrxCrxCbrxurxumm2 . 6,)(exp)(),()(exp)(),(22xdCCm49. 503/22

9、/202217静止均质环境淹没紊动射流主要特性表参数二维平面射流轴对称射流断面上最大流速断面上最大流速um断面上最大浓度断面上最大浓度Cm射流半宽度射流半宽度b轴线稀释度轴线稀释度断面平均稀释度断面平均稀释度浓度分布与速度分浓度分布与速度分布宽度比布宽度比卷吸系数卷吸系数断面流速分布断面流速分布断面浓度分布断面浓度分布21028. 2BxuumxDuum2 . 6021034. 2BxCCmxb154. 0 xDCCm59. 50 xb114. 02143.0BxSmDxSm18. 021062. 0BxqqSDxQQS32. 0041. 112. 1069. 0056. 02)(exp),(

10、brurxum2)(exp),(byuyxum2)(exp),(brCrxCm2)(exp),(byCyxCm3/22/2022184、静止环境水体中的浮射流 实际的扩散质排放出流是介于纯射流和纯羽实际的扩散质排放出流是介于纯射流和纯羽流之间，其密度佛汝德数流之间，其密度佛汝德数Fd为有限值，这种为有限值，这种浮射流在排放口附近动量的影响占优势，而浮射流在排放口附近动量的影响占优势，而在远离排放口的地方，出流的动能已经耗散，在远离排放口的地方，出流的动能已经耗散，于是完全混合的水流就像纯羽流一样了。于是完全混合的水流就像纯羽流一样了。 这里密度佛汝德数这里密度佛汝德数Fd定义为定义为gLuFa

11、d静水中平面浮射流紊动扩散3/22/202219 式中，式中，u为射流的特征流速，为射流的特征流速， L为特征长度，为特征长度， 为射流密度为射流密度, a为环境流体密度。为环境流体密度。 密度佛汝德数密度佛汝德数Fd反映了作用于射流的惯性力与反映了作用于射流的惯性力与浮力之比。当浮力之比。当Fd很大时，表明射流是由动量起很大时，表明射流是由动量起支配作用，当支配作用，当Fd很小时，则由浮力起支配作用。很小时，则由浮力起支配作用。若以当若以当Fd0代表射流出口处密度佛汝德数，显然，代表射流出口处密度佛汝德数，显然，是由动量起支配作用，当是由动量起支配作用，当Fd0 0时，属于浮力时，属于浮力羽

12、流。若羽流。若Fd0 ，浮力作用趋于零，为纯射流。，浮力作用趋于零，为纯射流。若若Fd0处于两者之间则为浮射流。处于两者之间则为浮射流。3/22/202220 浮射流的轨迹线如图所示。由于浮射流取决于起始出浮射流的轨迹线如图所示。由于浮射流取决于起始出流角度、浮力、惯性力和紊动剪切力，所以轨迹曲线流角度、浮力、惯性力和紊动剪切力，所以轨迹曲线是未知的。在局部的是未知的。在局部的“自然自然”坐标系统（坐标系统（s,n）中，）中，可写出连续性方程、动量方程和质量守恒方程，这里可写出连续性方程、动量方程和质量守恒方程，这里s是沿浮射流带轴线在水流方向的坐标，而是沿浮射流带轴线在水流方向的坐标，而n是

13、垂直于是垂直于s方向的坐标。方向的坐标。 XYbbDSu(S,n)3/22/2022214、静止环境水体中的浮射流 假定掺混假说沿弯曲的轨迹线有效，可得下假定掺混假说沿弯曲的轨迹线有效，可得下述几种分布：述几种分布： )(exp)(),(2bnsunsum)(exp)(),(2bnsCnsCm)(exp)(),(2bnsnsm3/22/202222点源羽流基本方程 考虑一种简化模式，认为羽流的源是出自一点，考虑一种简化模式，认为羽流的源是出自一点，周围环境为无限空间静止液体，由于铅垂方向周围环境为无限空间静止液体，由于铅垂方向的浮力作用而形成了流体的上升运动。由于紊的浮力作用而形成了流体的上升

14、运动。由于紊流作用，不断卷吸周围液体，羽流断面逐渐扩流作用，不断卷吸周围液体，羽流断面逐渐扩大。因为周围液体阻力，在横断面上流速分布大。因为周围液体阻力，在横断面上流速分布不均匀，沿轴心线上流速最大，然后向边缘部不均匀，沿轴心线上流速最大，然后向边缘部分逐渐减小。分逐渐减小。 连续性方程连续性方程 bubudsdmm2)(23/22/202223 x方向动量方程方向动量方程 y方向动量方程方向动量方程 相对密度守恒方程相对密度守恒方程 0)cos2(22budsdm2222)sin2(bgbudsdamm0)(2budsdmm3/22/202224 污染物浓度方程污染物浓度方程 浮射流的轨迹定

15、义为浮射流的轨迹定义为 0)(2buCdsdmmsin,cosdsdydsdx3/22/202225静止均质环境中羽流主要特性表参数二维平面羽流轴对称羽流断面上最大流速断面上最大流速um断面上最大浓度断面上最大浓度Cm单位起始浮力通量单位起始浮力通量B0任意断面体积流量任意断面体积流量Q羽流半宽度羽流半宽度b浓度分布与速度分布宽浓度分布与速度分布宽度比度比卷吸系数卷吸系数单位动量通量单位动量通量M任意断面密度佛汝德数任意断面密度佛汝德数F31005. 2Bum000gqBaxBQ310535.0 xb147. 024. 116. 1130. 0085. 0 xBM320774. 048. 3F

16、3131074. 4xBum1310004 . 2xBCQCm353100017.11xBCQCm000gQBa35310156.0 xBQ xb102. 045. 4F3432037. 0 xBM 3/22/202226例题1 某排污管将生活污水排至湖泊，其出口为狭长矩形，某排污管将生活污水排至湖泊，其出口为狭长矩形，孔宽孔宽0.2米，污水出流方向垂直向上，初始射流流米，污水出流方向垂直向上，初始射流流速速u0=4.0米米/秒，出口平面位于湖面下秒，出口平面位于湖面下24米，排泄污米，排泄污水浓度为水浓度为1200ppm，设污水与湖水密度基本相同，设污水与湖水密度基本相同，试求到达湖面时的最

17、大流速、最大浓度及平均稀释试求到达湖面时的最大流速、最大浓度及平均稀释度。度。3/22/202227 解：不计污水与湖水的密度差，可按二元淹解：不计污水与湖水的密度差，可按二元淹 没射流计算。到达水面时的没射流计算。到达水面时的 1.到达湖面时的最大流速到达湖面时的最大流速 2.到达湖面时的最大浓度到达湖面时的最大浓度 3.到达湖面时的平均稀释度到达湖面时的平均稀释度1202 . 024Bx210)(28. 2Bxuumsmum/83. 04)120(28. 2210913. 034. 2)(34. 2210BxCCmppmCCm2560913. 034. 208 . 6)120(62. 0)

18、(62. 02121BxS3/22/202228例题2 某排污管将生活污水排至湖泊，其出口为直径某排污管将生活污水排至湖泊，其出口为直径0.2米的圆形喷口，污水出流方向垂直向上，初始射米的圆形喷口，污水出流方向垂直向上，初始射流流速流流速u0=4.0米米/秒，出口平面位于湖面下秒，出口平面位于湖面下24米，排米，排泄污水浓度为泄污水浓度为1200ppm，设污水与湖水密度基本相，设污水与湖水密度基本相同，试求到达湖面时的最大流速、最大浓度及平同，试求到达湖面时的最大流速、最大浓度及平均稀释度。均稀释度。3/22/202229 解：不计污水与湖水的密度差，可按圆断面解：不计污水与湖水的密度差，可按

19、圆断面淹没射流计算。到达水面时的淹没射流计算。到达水面时的 1.到达湖面时的最大流速到达湖面时的最大流速 2.到达湖面时的最大浓度到达湖面时的最大浓度 3.到达湖面时的平均稀释度到达湖面时的平均稀释度1202 . 024DxxDuum2 . 60smum/21. 00 . 412012 . 6xDCCm59. 50ppmCCm9 .55120159. 504 .3832. 0DxS计算表明，圆形喷口比宽度与直径相同的二元喷计算表明，圆形喷口比宽度与直径相同的二元喷口混合和稀释效率要高。口混合和稀释效率要高。3/22/202230流动环境水体中的浮射流 在均匀环境水体流动的情况下，浮射流的特在均

20、匀环境水体流动的情况下，浮射流的特性是很复杂的，射流与环境水流的相互作用，性是很复杂的，射流与环境水流的相互作用，能引起射流发生不同的变化，从流动的轨迹能引起射流发生不同的变化，从流动的轨迹简单地发生弯曲变形到射流结构的全部崩解。简单地发生弯曲变形到射流结构的全部崩解。 3/22/202231流动环境水体中的浮射流 有横向水流的三维浮射流如图所示。浮射流上有横向水流的三维浮射流如图所示。浮射流上升并和周围水体混合时，由于受到两个机械力升并和周围水体混合时，由于受到两个机械力的作用而向下游弯曲，因掺混作用而得到环境的作用而向下游弯曲，因掺混作用而得到环境水流的动量。周围流动的水流在射流下游形成水

21、流的动量。周围流动的水流在射流下游形成一个低压区，产生一拖曳力。一个低压区，产生一拖曳力。 此外，可观测到射流的流速和浓度分布有明显此外，可观测到射流的流速和浓度分布有明显的变形。因为浮射流的上表部分的稀释速率比的变形。因为浮射流的上表部分的稀释速率比内核部分的要快，外表面部分的浮力变得小一内核部分的要快，外表面部分的浮力变得小一些。因此，随身射流转弯，射流外表面部分往些。因此，随身射流转弯，射流外表面部分往往有向内测倾斜的趋势，产生一个马蹄形的分往有向内测倾斜的趋势，产生一个马蹄形的分布。布。 3/22/202232流动环境水体中的浮射流 xy射流轴外边界内边界MM3/22/202233流动环境中的紊动射流 根据很多学者的研究，横流中的紊流射流可根据很多学者的研究，横流中的紊流射流可分为流动特性有明显差别的三部分：分为流动特性有明显差别的三部分： 近区近区，即射流离开喷口后距离不长的一段范围。，即射流离开喷口后距离不长的一段范围。射流离开喷口后，由于紊动卷吸作用，只有中射流离开喷口后，由于紊动卷吸作用，只有中间核心部分保持初始流速，边缘部分