高等流体力学-第五讲ppt课件.ppt

《高等流体力学-第五讲ppt课件.ppt》由会员分享，可在线阅读，更多相关《高等流体力学-第五讲ppt课件.ppt（51页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论本讲主要内容本讲主要内容一、一、基本概念基本概念二、二、描述扩散运动的基本运动方程描述扩散运动的基本运动方程三、三、扩散系数及其分析确定方法扩散系数及其分析确定方法四、四、扩散运动的解析解扩散运动的解析解五、五、岸边排放与中心排放污染带的计算岸边排放与中心排放污染带的计算2北京工业大学市政学科部马长明 高等流体(水）力学讲稿一、基本概念一、基本概念1、扩散现象、扩散现象烟囱排烟；河流排污；水面蒸发；食糖与食盐的溶解等。烟囱排烟；河流排污；水面蒸发；食糖与食盐的溶解等。2、传输过程、传输过程流体中所含有物质（如各

2、种污染物，也包括动量、能量和热量）在流场流体中所含有物质（如各种污染物，也包括动量、能量和热量）在流场中某一处到另一处转移的过程。中某一处到另一处转移的过程。3、扩散（、扩散（Diffusion）是一类传输过程，指物质由含量高处向含量低处的传输过程。是一类传输过程，指物质由含量高处向含量低处的传输过程。（1）按扩散的机制可分为）按扩散的机制可分为：分子扩散、对流传输扩散及紊动扩散。：分子扩散、对流传输扩散及紊动扩散。1）分子扩散（）分子扩散（Molecular Diffusion）：由分子运动产生。）：由分子运动产生。2）紊动扩散（）紊动扩散（Turbulent Diffusion）：由流体质

3、点的紊动产生的扩散。）：由流体质点的紊动产生的扩散。3）移流传输（）移流传输（Advection）：扩散物随同流体质点的时均运动而转移。）：扩散物随同流体质点的时均运动而转移。第五讲第五讲 扩散理论扩散理论3北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（2）按研究问题的类型分：）按研究问题的类型分：1）剪切流中的离散（）剪切流中的离散（Dispersion）：由于剪切流中速度分布不均匀）：由于剪切流中速度分布不均匀产生含有随流散开的作用，也称弥散。产生含有随流散开的作用，也称弥散。离散中包含有移流扩散和紊动扩散。离散中包含有移流扩散和紊动扩散。2）射流扩散

4、：指从各种排泄口喷出流入周围另一流体域内运动的）射流扩散：指从各种排泄口喷出流入周围另一流体域内运动的一股流体。包括移流扩散和紊动掺混扩散。一股流体。包括移流扩散和紊动掺混扩散。其中射入同种性质的流体内称其中射入同种性质的流体内称淹没射流淹没射流，射入不同性质的流体内为，射入不同性质的流体内为非淹没射流非淹没射流。按射流的原动力还可分为：动量射流、浮力羽流、浮射流。按射流的原动力还可分为：动量射流、浮力羽流、浮射流。动量射流（动量射流（Jets）：射流以出流的动量（：射流以出流的动量（Momentum）为原动力，）为原动力，该动量对射流运动起主要作用。该动量对射流运动起主要作用。浮力羽流（浮力

5、羽流（Plumes）：浮力（：浮力（Bouyancy Forces）是原动力，产生）是原动力，产生的运动形态呈羽毛状。如烟气，水体中泄入污染液体后的运动等。的运动形态呈羽毛状。如烟气，水体中泄入污染液体后的运动等。浮射流（浮射流（Buoyant Jets）：原动力既有动量又有浮力。：原动力既有动量又有浮力。3）分层流（）分层流（Stratified Flowing）：在重力场中密度不均匀的流体）：在重力场中密度不均匀的流体形成有层次的流动。形成有层次的流动。4北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论4、污水泄入河道中问题研究的阶段划分、污水泄入河道中问

6、题研究的阶段划分第一阶段：污水离开排污口与周围水体掺混，以射流和浮力羽流形式扩散，一般第一阶段：污水离开排污口与周围水体掺混，以射流和浮力羽流形式扩散，一般按三维运动问题处理；按三维运动问题处理；第二阶段：污水还没有扩展到河流全断面，污水随河水一起运动，按二维紊动离第二阶段：污水还没有扩展到河流全断面，污水随河水一起运动，按二维紊动离散问题研究；散问题研究；第三阶段：污水已扩展到河流全断面，并且全断面完全混合，污水沿纵向继续离第三阶段：污水已扩展到河流全断面，并且全断面完全混合，污水沿纵向继续离散，可按一维纵向离散问题分析。散，可按一维纵向离散问题分析。5、扩散理论的基本假设、扩散理论的基本假

7、设（1）扩散质的存在不改变流体质点的流动特性。即将）扩散质的存在不改变流体质点的流动特性。即将扩散质扩散质视为视为标志物质标志物质或称或称示踪剂示踪剂（Tracer）。）。（2）流体质点上所带的）流体质点上所带的扩散质扩散质在运动过程中保持不变，即流体质点间不发生扩在运动过程中保持不变，即流体质点间不发生扩散质的转移（不计分子扩散），扩散质的扩散完全是由于带有扩散质的流体质点发生散质的转移（不计分子扩散），扩散质的扩散完全是由于带有扩散质的流体质点发生掺混的结果。掺混的结果。（3）对不可压流体，携带扩散质流体质点的总体积在扩散过程中保持不变，扩）对不可压流体，携带扩散质流体质点的总体积在扩散过

8、程中保持不变，扩散结果反映在携带扩散质的流体质点所占空间位置和轮廓随时间而变化。散结果反映在携带扩散质的流体质点所占空间位置和轮廓随时间而变化。5北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论二、描述扩散运动的基本方程二、描述扩散运动的基本方程1、分子扩散的费克定律、分子扩散的费克定律 （第一定律）（第一定律）费克（费克（Adolph FickAdolph Fick，18551855）提出假设：盐分在其溶液中扩散的物理定）提出假设：盐分在其溶液中扩散的物理定律应等同于傅立叶（律应等同于傅立叶（FourierFourier，18221822）提出的热传导定律相

9、同，即：）提出的热传导定律相同，即：其中：其中：费克定律说明：在扩散溶液浓度场中的时空点上，单位时间内通过单位费克定律说明：在扩散溶液浓度场中的时空点上，单位时间内通过单位面积的扩散质的质量与该点处扩散溶液浓度的梯度成正比，比例系数为该种面积的扩散质的质量与该点处扩散溶液浓度的梯度成正比，比例系数为该种扩散溶液的分子扩散系数；方向与浓度梯度方向相反。扩散溶液的分子扩散系数；方向与浓度梯度方向相反。分子扩散系数由扩散质及溶解质的种类、温度、压强决定，与溶液的运分子扩散系数由扩散质及溶解质的种类、温度、压强决定，与溶液的运动形态无关，是物性参数。动形态无关，是物性参数。imxCDqC扩散液的浓度；

10、量纲：扩散液的浓度；量纲：3MLq在在xi方向的单位面积的扩散质量通量；量纲：方向的单位面积的扩散质量通量；量纲：12TMLmD分子扩散系数；量纲：分子扩散系数；量纲：12TL常见扩散质在水中的扩散系数表常见扩散质在水中的扩散系数表对三维情况，以矢量表示：对三维情况，以矢量表示：CDqm6北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2、分子扩散方程（费克第二定律）、分子扩散方程（费克第二定律）分析分子扩散应满足的控制方程，以一维情况为例，流体静止。分析分子扩散应满足的控制方程，以一维情况为例，流体静止。取控制体如图，由质量守恒定律，可得：取控制体如图，由质

11、量守恒定律，可得：dtxxqqdtqdttxC1)(1) 1(0 xqtC整理得：整理得：将费克第一定律代入，可得：将费克第一定律代入，可得：22xCDtCm22xqDtqm或或CDzCyCxCDtCmm2222222对三维情况：对三维情况：上述方程称之为上述方程称之为分子扩散方程分子扩散方程（Diffusion Equations）。）。7北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论3、移流扩散方程（、移流扩散方程（Advective Diffusion Equations）取控制体如图，以取控制体如图，以x1方向为例。方向为例。假设：层流运动时溶液的扩

12、散与流体静止假设：层流运动时溶液的扩散与流体静止 时的分子扩散相同。时的分子扩散相同。由质量守恒定律，可得：由质量守恒定律，可得：dtdxdxqqdtdxqdxdttdxdxCdx3232321)()(1xxqq11xCDCuqm及及21211)(xCDxCutCm整理可得：整理可得：对三维流动：对三维流动：232222212332211)()()(xCxCxCDxCuxCuxCutCm8北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论4、紊动扩散方程、紊动扩散方程对紊动扩散，瞬时浓度和流速都可分解时均值与脉动值之和，即：对紊动扩散，瞬时浓度和流速都可分解时均

13、值与脉动值之和，即：代入移流扩散方程后取时均值，有：代入移流扩散方程后取时均值，有：仿照分子扩散系数的表示形式，引入紊动扩散系数仿照分子扩散系数的表示形式，引入紊动扩散系数Dij，令：令：紊动扩散方程可表示为（考虑源、汇项后）：紊动扩散方程可表示为（考虑源、汇项后）：iimiiiuCxCDxxCutCuuuCCC和和iuC 其中：其中：是由紊动脉动量产生的沿是由紊动脉动量产生的沿i方向的质量扩散通量。方向的质量扩散通量。jijixCDuCcjijimiiiFxCDxCDxxCutC9北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论紊动扩散质量通量与紊动扩散系数

14、紊动扩散质量通量与紊动扩散系数Dij（Turbulent Diffusion Coefficients）可用矩阵表示为：）可用矩阵表示为：Dij应是空间坐标的函数，当选择坐标使其与二阶张量的主轴方向一致应是空间坐标的函数，当选择坐标使其与二阶张量的主轴方向一致时，九个量中仅有三个主值，即：时，九个量中仅有三个主值，即：D11，D22，D33不为零。当满足各向同性不为零。当满足各向同性条件下，有：条件下，有：321333231232221131211321xCxCxCDDDDDDDDDuCuCuCTDDDD332211iiTiixxCDxCutC2mTDD 因为因为，分子扩散可忽略不计，紊动扩散

15、方程为：，分子扩散可忽略不计，紊动扩散方程为：10北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论三、扩散系数及其分析确定三、扩散系数及其分析确定1、分子扩散系数与概率统计量间的关系、分子扩散系数与概率统计量间的关系（1）分子扩散方程的基本解）分子扩散方程的基本解问题问题： （考虑一维问题）在考虑一维问题）在t=0时刻，坐标原点处（时刻，坐标原点处（x=0）放置质量为放置质量为M的扩散质，确定浓度沿的扩散质，确定浓度沿x轴的扩散过程。轴的扩散过程。22xCDtCm基本方程：基本方程：求解方法：求解方法：1）量纲分析相似解法）量纲分析相似解法 ； 2）数理方程解

16、法）数理方程解法定解条件：定解条件：由质量守恒，在任何时刻，有：由质量守恒，在任何时刻，有： MCdx或或)(),(xMxC011北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论量纲分析相似解法量纲分析相似解法量纲分析：量纲分析：C（x，t）必然与）必然与M、Dm和和t有关，有量纲关系有关，有量纲关系代入扩散方程，得：代入扩散方程，得：得出：得出： 分离变量法结果：分离变量法结果：)4(4)(4),(tDxftDMftDMtxCmmm 02fddf 20)( ecf再由定解条件再由定解条件 确定系数确定系数c0。 MCdx)4exp(4),(2tDxtDMtx

17、Cmm dtDxtDCtxCmm4)(exp4)0 ,(),(212北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（2）分子扩散的随机游走分析）分子扩散的随机游走分析确定分子游走的概率值：分子的自由程确定分子游走的概率值：分子的自由程l，经过经过N次运动后行走次运动后行走X距离距离的概率。设在的概率。设在N次行走中，有次行走中，有p次沿次沿x的正向，的正向，q次沿次沿x的反向行走；每次行走的反向行走；每次行走是相互独立的。是相互独立的。有：有：p + q = N；令；令: p q = S ，故有：，故有：X = Sl沿沿x的正向行走的正向行走X距离的概率距离

18、的概率P为：为：利用利用Sterling公式：公式：可推出：可推出：)!1 (2)!1 (22!2) !/(!NSNNSNNqpNPNN)2ln(21ln)21(!ln nnnnNSNP2exp22 13北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（3）分子扩散与分子随机游走概率密度间的关系）分子扩散与分子随机游走概率密度间的关系令令a为分子运动速度，为分子运动速度，t为分子运动为分子运动N次所用时间，则有：次所用时间，则有：随机游走概率可表示为：随机游走概率可表示为：如令如令概率值概率值分子落于分子落于X，X+dX间的概率为：间的概率为：atlXatlP

19、2exp22 lXSlatN/ ,/tNlalDm2212tDXtDlPmm4exp2 lXtDXtDlPmm24exp2 14北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论可见概率密度：可见概率密度：为扩散方程（为扩散方程（M=1）的解。）的解。可见分子扩散方程的解满足标准差为可见分子扩散方程的解满足标准差为 的正态分布。的正态分布。分析结论：分析结论：1）分子扩散的浓度可以用概率密度表示；）分子扩散的浓度可以用概率密度表示；2）分子扩散系数可由概率分布的方差确定，即：）分子扩散系数可由概率分布的方差确定，即：或积分得或积分得tDXtDXPmm4exp41

20、2 tDm2 dtdDm221 )(2122122ttDm 15北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（4）分子扩散系数与方差的关系证明）分子扩散系数与方差的关系证明1）浓度分布各阶矩的定义）浓度分布各阶矩的定义2）浓度分布各阶矩与统计量间的关系）浓度分布各阶矩与统计量间的关系dxtxxCM),(1一阶矩：一阶矩：dxtxCxM),(22二阶矩：二阶矩：dxtxCxMnn),(n阶矩：阶矩：dxtxCM),(0零阶矩：零阶矩：扩散质总质量：扩散质总质量浓度分布质量中心坐标浓度分布质量中心坐标均值：均值：01MM/浓度分布方差：浓度分布方差：20220

21、21MMdxtxCxM/),()(取：取： 0 , 0 ; 110 MM即即可证明可证明mDdtd22 16北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论3）分子扩散系数与方差的关系证明）分子扩散系数与方差的关系证明扩散方程两边同乘扩散方程两边同乘x2后从后从-到到+积分积分dxxCDxdxtCxm2222222dxxCxxCxDdxCxdtdmmmDdxCxCxCxDdtd22222 上是说明：确定扩散系数可以从研究浓度分布的方差着手。上是说明：确定扩散系数可以从研究浓度分布的方差着手。17北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲

22、扩散理论扩散理论2、紊动扩散系数、紊动扩散系数Dt的分析的分析假设紊动场均匀、各向同性，以一维为例阐述。假设紊动场均匀、各向同性，以一维为例阐述。（1）基本概念）基本概念 1）紊动扩散现象分析）紊动扩散现象分析观测实验结果说明：观测实验结果说明： 浓度中心位置；浓度中心位置； 浓度扩散分析；浓度扩散分析； 系综平均结果。系综平均结果。 2）扩散系数的物理解释）扩散系数的物理解释两质点的相对扩散（系综平均）；两质点的相对扩散（系综平均）；单个质点的扩散（系综平均）。单个质点的扩散（系综平均）。 3）系综平均与时间平均）系综平均与时间平均各态遍历特性各态遍历特性扩散观测扩散观测扩散现象图示扩散现象

23、图示两质点的相对扩散图示两质点的相对扩散图示单个质点的扩散图示单个质点的扩散图示各态遍历图示各态遍历图示18北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（2）单个质点的紊动扩散）单个质点的紊动扩散泰勒扩散理论泰勒扩散理论泰勒（泰勒（G. I. Taylor 1921）研究了单个质点在平稳、各向同性的紊流场中）研究了单个质点在平稳、各向同性的紊流场中的扩散问题。的扩散问题。考虑一维情况下，设紊动速度为考虑一维情况下，设紊动速度为u，平均速度，平均速度=0；在；在t时刻，紊动质点的坐时刻，紊动质点的坐标位置为标位置为X(t)，X(t)的时间平均值（系综平均值）

24、：的时间平均值（系综平均值）：令开始记录的时间为令开始记录的时间为t0，经过，经过t时刻，质点所移动的距离为时刻，质点所移动的距离为其时间平均平方值可表示为：其时间平均平方值可表示为：tt dtutX0)()(0)(tXtt dttuttX000)()(00000000022)()(1)(1)(dttdttutdttuTdtttXTtXTttT t tTtdtddtttuttuT0 00000)()(119北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论)()()()()(200 Rutututtuttu ttt tt dt dt dt d000 02 ttt

25、dtdttuttutX0 0002)()()(积分限变换：积分限变换：由均匀和各向同性条件：由均匀和各向同性条件：ttdRtdutX0022)(2)( 代入时均方差表示式，有：代入时均方差表示式，有：分部积分后分部积分后tdRtutX022)()(2)( 20北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论讨论两种极端情况：扩散时间很短，扩散时间很长。讨论两种极端情况：扩散时间很短，扩散时间很长。1）扩散时间很短，）扩散时间很短，R1，积分，积分即：即：扩散系数：扩散系数：2）扩散时间很长，当）扩散时间很长，当tTL后，后，R（）= = 0 0其中：其中：扩散

26、系数：扩散系数：220022)(2)()(2)(tudtdRtutXtt tutX22)(tudttXddtdDt222)(2121 consttTudRtutXL20222)()(2)( LtTudttXddtdD222)(2121 LTdR0)( 21北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论有紊动理论知，积分时间比尺有紊动理论知，积分时间比尺TL与扩散长度比尺与扩散长度比尺满足：满足：所以有：所以有：实验资料说明：定常均匀紊动的流场在实验资料说明：定常均匀紊动的流场在tTtTL L后是接近正态分布的，紊后是接近正态分布的，紊动质点运动已成为随机运动

27、，紊动扩散和分子扩散遵循相同的规律。动质点运动已成为随机运动，紊动扩散和分子扩散遵循相同的规律。022)( dRuTuL222)(21uTudttXdDLt22北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论例题例题(pp.119-121)：设在一均匀紊流场内，在原点投入许多示踪：设在一均匀紊流场内，在原点投入许多示踪质粒子，量测不同时刻粒子的横向位移质粒子，量测不同时刻粒子的横向位移X，X2的系综平均值的系综平均值 及通过原点后的时间及通过原点后的时间t的数值如表所示。式绘出的数值如表所示。式绘出 关系曲线，关系曲线，推算出紊动扩散系数推算出紊动扩散系数Dt

28、，计算脉动速度均方值，计算脉动速度均方值 及扩散长度及扩散长度比比 。 2XtX 22u数据表格数据表格所作曲线图所作曲线图smXXtXDt/.242702012100370012121tuX22smtXu/.024020480222uDtmuDt2421025102401003.解：解：23北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（3）两质点的相对扩散）两质点的相对扩散1）基本定义）基本定义以一维运动为例，对均匀、各向同性的平稳紊动场，设两质点以一维运动为例，对均匀、各向同性的平稳紊动场，设两质点a、b的的脉动速度速度为脉动速度速度为ua(t), u

29、b(t)，脉动速度差为，脉动速度差为w(t)：两点位移差为：两点位移差为：与单质点定义类似，可定义相对扩散均方值：与单质点定义类似，可定义相对扩散均方值：相对扩散系数：相对扩散系数：定义相对扩散速度：定义相对扩散速度：)()()(tututwabdtsdDr221tabtdtwztXtXtz0)()0()()()(z(t)zs2202ssdtdur24北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2）基本运算公式）基本运算公式两质点的相对扩散的自变量为两质点的相对扩散的自变量为s0、和、和t，其它量都是，其它量都是s0、t的函的函数，所以数，所以),(),(

30、),(21),(00020tszdtdtszdttssdtsDrt ttrt ddswtswst dtsDss0 0002000202),(),(2),(2 dtswswt000),(),(t trt ddswtswdtdssdtdtsu0210002020),(),(2),( 25北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论3）主要结论）主要结论 扩散时间扩散时间t t很短很短constswsw)0 ,()0 ,(00tAdswswtsDtr10000)0 ,()0 ,(),( 212000202),(2tAstdtsDsstr12020),(Assdt

31、dtsur)()()()()(2tXtXtXtXtzzsababXXXXXXaabb221022),(tDtDdtdtsurrr),(),(00tsztsz 、 扩散时间很长扩散时间很长 相互独立相互独立为单个质点位移均方值的二倍，单个质点的扩散系数为常数。为单个质点位移均方值的二倍，单个质点的扩散系数为常数。26北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论 扩散时间扩散时间t t不短不短随着随着t t的增加，初始状态的增加，初始状态s s0 0的影响消失，但紊动尺度在惯性区，影的影响消失，但紊动尺度在惯性区，影响因素仅有能耗响因素仅有能耗 量纲：量纲：L

32、 L2 2T T-3-3 ，由量纲分析，可得定性关系：，由量纲分析，可得定性关系：所以有：所以有：相对扩散系数与相对扩散距离的均方值存在以下关系相对扩散系数与相对扩散距离的均方值存在以下关系tAtwtw 22)()(220)()()(tAdtwtwtDtr 320232)(2)(tAt dtDtstr tAssdtdtur 220223)(322312)(49)(sAtDr 27北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论四、扩散运动的解析解四、扩散运动的解析解1、扩散方程特性、扩散方程特性1）当扩散系数为常数时，浓度扩散方程为线性的二阶偏微分）当扩散系数

33、为常数时，浓度扩散方程为线性的二阶偏微分方程，解可以叠加。方程，解可以叠加。2）将扩散质假设为示踪剂后，浓度与流场的耦合问题不需考）将扩散质假设为示踪剂后，浓度与流场的耦合问题不需考虑，计算步骤是先求出流场中各运动参数，而后用浓度扩散方程求虑，计算步骤是先求出流场中各运动参数，而后用浓度扩散方程求浓度分布。浓度分布。3）对简单的边界条件和初始值问题，扩散方程有解析解。）对简单的边界条件和初始值问题，扩散方程有解析解。4）对复杂问题可通过数值解或实验测定方法确定浓度分布问）对复杂问题可通过数值解或实验测定方法确定浓度分布问题。题。332332222211211332211xxCDxxCDxxCD

34、xCuxCuxCutC28北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2、扩散问题的解析解、扩散问题的解析解静止流体中瞬时源的扩散静止流体中瞬时源的扩散（1）集中投入的情况）集中投入的情况1）一维问题）一维问题问题描述：在问题描述：在t=0时刻，在原点时刻，在原点x=0处集中投放质量为处集中投放质量为M的扩的扩散质。要求确定散质。要求确定t0时刻在一维空间点时刻在一维空间点x处的浓度值（扩散系数已处的浓度值（扩散系数已知）。知）。基本方程：基本方程：定解条件：定解条件：基本解：基本解：在在x=投放，解为：投放，解为：22xCDtCmdxtxCM),(0)(

35、)0,(xMxC 或或tDxtDMtxCmm4exp4),(2 tDxtDMtxCmm4)(exp4),(2 29北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2、扩散问题的解析解、扩散问题的解析解静止流体中瞬时源的扩散静止流体中瞬时源的扩散（1）集中投入的情况）集中投入的情况2）二维问题）二维问题问题描述：在问题描述：在t=0时刻，在原点时刻，在原点x=0，y=0处集中投放质量为处集中投放质量为M的的扩散质。要求确定扩散质。要求确定t0时刻在二维空间点（时刻在二维空间点（x，y）处的浓度值（扩散）处的浓度值（扩散系数系数Dx、Dy已知）。已知）。基本方程：

36、基本方程：定解条件：定解条件：解的表示：解的表示：在在x=，y=投放，解如何表示？三维问题的解如何表示？投放，解如何表示？三维问题的解如何表示？222221yCDxCDtCdyCdxM 0)()()0,(yxMyxC 或或tDytDxDDtMtyxCyxyx44exp4),(22 30北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2、扩散问题的解析解、扩散问题的解析解静止流体中瞬时源的扩散静止流体中瞬时源的扩散（2）空间分布的情况）空间分布的情况1）一维问题）一维问题问题描述：在问题描述：在t=0时刻，在时刻，在x=a至至x=b段分布有质量密度函数段分布有质

37、量密度函数为为f(x)的扩散质。要求确定的扩散质。要求确定t0时刻在一维空间点时刻在一维空间点x处的浓度值（扩处的浓度值（扩散系数散系数Dx已知）。已知）。基本方程：基本方程：定解条件：定解条件：解的表示：解的表示：22xCDtCxbxaxfxC其它 )(0)0,( dtDxDtftxCbax4)(exp4)(),(2瞬时分布源瞬时分布源31北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2、扩散问题的解析解、扩散问题的解析解静止流体中瞬时源的扩散静止流体中瞬时源的扩散（2）空间分布的情况）空间分布的情况2）阶跃函数问题）阶跃函数问题问题描述：在问题描述：在t

38、=0时刻，在时刻，在x=0至至x=+段分布有质量密度函段分布有质量密度函数为数为f(x)=c0的扩散质。要求确定的扩散质。要求确定t0时刻在一维空间点时刻在一维空间点x处的浓度处的浓度值（扩散系数值（扩散系数Dx已知）。已知）。基本方程：基本方程：定解条件：定解条件：解的表示：解的表示：22xCDtCxxxcxC00 0)0,(0 dtDxDtctxCx0204)(exp4),(tDxerfcx4120浓度分布浓度分布32北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论误差函数误差函数定义：定义：性质：性质：其中：其中：erfc(z)成为补误差函数。成为补误差

39、函数。erf(z)及及erfc(z)的值可查数学手册。的值可查数学手册。zdezerf022)( 1)(erf0) 0(erf)()(zerfzerf)()(1zerfczerf33北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论3、扩散问题的解析解、扩散问题的解析解静止流体中连续源的扩散静止流体中连续源的扩散（1）等强度连续点源的扩散）等强度连续点源的扩散问题描述：从问题描述：从t=0时刻开始，在时刻开始，在x=0放置一等强度的连续点源，即在放置一等强度的连续点源，即在x=0处有一浓度为处有一浓度为C0的浓度源。要求确定的浓度源。要求确定t0时刻在一维空间点

40、时刻在一维空间点x处的浓度处的浓度值（扩散系数值（扩散系数Dx已知）。已知）。基本方程：基本方程：定解条件：定解条件：基本解：基本解：等强度连续源放置在等强度连续源放置在x=处其解如何表示？处其解如何表示？22xCDtCxtDxerfcCtDxerfCtxCxx441),(000 ),0(0tCtC等强度连续扩散结果图示等强度连续扩散结果图示34北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论3、扩散问题的解析解、扩散问题的解析解静止流体中连续源的扩散静止流体中连续源的扩散（2）变强度连续点源的扩散）变强度连续点源的扩散问题描述：从问题描述：从t=0时刻开始，

41、在时刻开始，在x=0放置一变强度的连续点源，即在放置一变强度的连续点源，即在x=0处有一浓度为处有一浓度为f(t)的浓度源。要求确定的浓度源。要求确定t0时刻在一维空间点时刻在一维空间点x处的浓度处的浓度值（扩散系数值（扩散系数Dx已知）。已知）。基本方程：基本方程：定解条件：定解条件：解的表示：解的表示：变强度连续源放置在变强度连续源放置在x=处其解如何表示？处其解如何表示？22xCDtCxtxxdtDxtDftxC02)(4exp)(4)(),( 0 )(), 0(ttftC变强度连续源图示变强度连续源图示35北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理

42、论3、扩散问题的解析解、扩散问题的解析解静止流体中连续源的扩散静止流体中连续源的扩散（3）变强度连续分布源的扩散）变强度连续分布源的扩散问题描述：从问题描述：从t=0时刻开始，在时刻开始，在x=a至至x=b放置变强度的连续分布源，放置变强度的连续分布源，即从即从x=a至至x=b上有浓度为上有浓度为f(x,t)的分布浓度源。要求确定的分布浓度源。要求确定t0时刻在一维空时刻在一维空间点间点x处的浓度值（扩散系数处的浓度值（扩散系数Dx已知）。已知）。基本方程：基本方程：定解条件：定解条件：解的表示：解的表示：22xCDtCx0 , ),(0),(tbxatxftxC其它 t baxxddtDxt

43、DftxC02)(4)(exp)(4),(),( 36北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论4、均匀紊流中的瞬时点源扩散问题（沿、均匀紊流中的瞬时点源扩散问题（沿x方向的均匀流方向的均匀流U）基本方程：基本方程：处理方法处理方法引入平行移动坐标系引入平行移动坐标系基本方程化简为：基本方程化简为：222222zCDyCDxCDxCUtCzzyyxxzzyyUtxx111 , ,txUtzzyyxx1111 , , ,), , ,(),(1111tzyUtxCtzyxC212212212zCDyCDxCDtCzzyyxx37北京工业大学市政学科部马长明

44、高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论4、均匀紊流中的瞬时点源扩散问题（沿、均匀紊流中的瞬时点源扩散问题（沿x方向的均匀流方向的均匀流U）（1）一维问题）一维问题方程：方程： 解的表示：解的表示：（2）三维问题）三维问题方程：方程：解的表示：解的表示：212xCDtCxx212212212zCDyCDxCDtCzzyyxxtDUtxtDMtxCxxxx4)(exp4),(2 tDztDytDUtxDDDtMtzyxCzzyyxxzzyyxx444)(exp)4(),(22223 38北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论5、有限边界问题

45、、有限边界问题边界反射扩散（以一维瞬时投入为例）边界反射扩散（以一维瞬时投入为例）（1）单边界图示单边界图示（2）双边界图示双边界图示（3）求解方法）求解方法叠加法叠加法由于边界的制约，在边界处质量通量由于边界的制约，在边界处质量通量q=0，即有：，即有：在边界两侧浓度分布对称，由此确定叠加点源位置与强度。在边界两侧浓度分布对称，由此确定叠加点源位置与强度。（4）解的表示）解的表示1）单边界）单边界2）双边界）双边界0 xCDtLxDtxDttxC4)2(exp4exp41),(22 nnDtnLxDttxC4)2(exp41),(2 39北京工业大学市政学科部马长明 高等流体(水）力学讲稿第

46、五讲第五讲 扩散理论扩散理论五、岸边排放与中心排放污染带的计算五、岸边排放与中心排放污染带的计算1 1、基本定义、基本定义（1 1）中心排放：排放口的位置位）中心排放：排放口的位置位于河流中央，如图所示。于河流中央，如图所示。（2 2）岸边排放：排放口的位置位）岸边排放：排放口的位置位于河流岸边，如图所示。于河流岸边，如图所示。（3 3）污染带宽度）污染带宽度B Bb b：在同一过流：在同一过流断面上，边缘点的浓度为最大浓断面上，边缘点的浓度为最大浓度值的度值的5%5%所对应点的宽度。所对应点的宽度。（4 4）污染带均匀混合长度）污染带均匀混合长度L Lm m：当：当污染带宽度上最大浓度值与最

47、小污染带宽度上最大浓度值与最小浓度值之差不超过浓度值之差不超过5%5%时称污染带时称污染带达到均匀混合状态。达到均匀混合状态。40北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论2 2、扩散问题类型、扩散问题类型 （1 1）坐标系的选择）坐标系的选择 沿水流流向（纵向）为沿水流流向（纵向）为x x轴，河轴，河宽方向（横向）为宽方向（横向）为y y轴，轴，z z轴垂直于水轴垂直于水面，坐标原点建在源处，如图所示。面，坐标原点建在源处，如图所示。Byy yDBuxx2Byy, x000相对点源坐标：相对点源坐标：无量纲坐标系：无量纲坐标系：41北京工业大学市政学科

48、部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论 起始断面平均浓度：起始断面平均浓度：hBuMQMC0 （2 2）扩散问题类型）扩散问题类型2 2）扩散空间：对宽而浅河流，认为在）扩散空间：对宽而浅河流，认为在z z（水深）方向为均匀混合（水深）方向为均匀混合态，故可简化为二维问题；单位线源强度（点源）为：态，故可简化为二维问题；单位线源强度（点源）为：hMM 单位时间内污染物的排放量，单位：单位时间内污染物的排放量，单位：m m3 3/s/s；h （x x，y y）处的水深，单位：）处的水深，单位：m m。3 3）当横向扩散（沿）当横向扩散（沿y y向）到达河岸后，须考虑两河岸

49、边界的反向）到达河岸后，须考虑两河岸边界的反射。射。1 1）扩散源：等强度点源移流扩散；移流断面平均流速度）扩散源：等强度点源移流扩散；移流断面平均流速度 ；u42北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论 3 3、计算的基本任务、计算的基本任务 当扩散处于稳定状态下（当扩散处于稳定状态下（tt）时，确定：）时，确定：xDyuuxDuhMx,yCyy4exp42（1 1）污染带的浓度分布：）污染带的浓度分布：xyxCC4exp4120以相对值表示：以相对值表示：考虑两岸边界反射下的浓度分布：考虑两岸边界反射下的浓度分布：nxynyxynyxCC42exp

50、42exp412020043北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论050max.xCx,bC（3 3）从排放点扩散到河流全断面需要的距离）从排放点扩散到河流全断面需要的距离L LB B：（2 2）污染带的宽度）污染带的宽度B Bb b（x x）：）：中心排放：中心排放：bxBb2岸边排放：岸边排放：bxBb中心排放：中心排放：0502max.LC,B/LCBB岸边排放：岸边排放：050max.xC,BLCB44北京工业大学市政学科部马长明 高等流体(水）力学讲稿第五讲第五讲 扩散理论扩散理论（4 4）从排放点到河流全断面均匀混合所需要的距离）从排放点