环境流体力学第二章分子扩散说课材料.ppt

《环境流体力学第二章分子扩散说课材料.ppt》由会员分享，可在线阅读，更多相关《环境流体力学第二章分子扩散说课材料.ppt（62页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、环境流体力学第二章分子扩散费克定律：费克定律：1855年德国生理学家费克（年德国生理学家费克（Fick）提出）提出静水中的污染物由于分子扩散作用，在单位时间内按一定方向通过单位面静水中的污染物由于分子扩散作用，在单位时间内按一定方向通过单位面积的扩散输送的物质与该方向的浓度梯度成正比。各向同性的介质。积的扩散输送的物质与该方向的浓度梯度成正比。各向同性的介质。式中：式中：Q是单位时间通过单位面积的扩散物质，也称为通量；是单位时间通过单位面积的扩散物质，也称为通量；C是扩散物质的浓度。是扩散物质的浓度。:x方向的浓度梯度。方向的浓度梯度。D是比例系数，称为分子扩散系数，量纲为是比例系数，称为分子

2、扩散系数，量纲为L2T-1 一般约为一般约为10-610-5cm2s-1 。x用等号用等号一维费克扩散示意图一维费克扩散示意图对一维扩散，费克定律可表示为：对一维扩散，费克定律可表示为：费克定律第一定律费克定律第一定律第三节第三节 费克定律费克定律公式中的负号公式中的负号三维的费克定律三维的费克定律:哈密顿算子哈密顿算子说明：只要存在浓度梯度，必然产生物质的扩散说明：只要存在浓度梯度，必然产生物质的扩散费克定律第二定律费克定律第二定律一滴红墨水在玻璃杯中的扩散一滴红墨水在玻璃杯中的扩散分子的扩散系数分子的扩散系数D与介质与物质本身的特性有关，又与温与介质与物质本身的特性有关，又与温度和压力有关

3、。度和压力有关。第三节第三节 费克定律费克定律某些物质在水中的分子扩散系数（某些物质在水中的分子扩散系数（cm2s-1，水温为，水温为20）物物 质质扩散系数扩散系数D物物 质质扩散系数扩散系数D氧氧1.8010-5醋酸醋酸0.8810-5二氧化碳二氧化碳1.5010-5甲醇甲醇1.2810-5一氧化氮一氧化氮1.5110-5乙醇乙醇1.0010-5氨氨1.7610-5酚酚0.8410-5氯氯1.2210-5甘汕甘汕0.7210-5氢氢5.1310-5尿素尿素1.0610-5氮氮1.6410-5葡萄糖葡萄糖0.6010-5氯化氢氯化氢2.6410-5蔗糖蔗糖0.4510-5硫化氢硫化氢1.80

4、10-5食盐食盐1.3510-5硫酸硫酸1.7310-5氢氧化钠氢氧化钠1.5110-5D值由实验确定，值由实验确定，D值大，扩散快；反之，扩散慢。值大，扩散快；反之，扩散慢。第三节第三节 费克定律费克定律单位时间进入单位时间进入x面的扩散质通量为：面的扩散质通量为：Q(x,t)从从(x+x)面出去的通量为面出去的通量为:设设c(x,t)是时刻是时刻t位于位于x处上扩处上扩散质（溶质）的浓度。在该控散质（溶质）的浓度。在该控制体积内扩散质对时间的制体积内扩散质对时间的变化率为：变化率为：第二节、分子扩散方程的推导第二节、分子扩散方程的推导(单纯扩散）单纯扩散）一维为例一维为例 第四节第四节 分

5、子扩散方程分子扩散方程一维输移的控制体：两个具有单位面一维输移的控制体：两个具有单位面积的平行面与积的平行面与x轴垂直轴垂直变化量：变化量：根据质量守恒定律有：单位时间流入的污染物质量根据质量守恒定律有：单位时间流入的污染物质量-流出流出的污染物质的污染物质=污染物质量对时间的变化率相等，即污染物质量对时间的变化率相等，即:Fick定律：定律：如将如将Q(x,t)作为热通量（即热流密度），作为热通量（即热流密度），c(x,t)作为热浓度（即温度），以作为热浓度（即温度），以热扩散系数热扩散系数（或导温系数）代替分子扩散系数（或导温系数）代替分子扩散系数D,变为热传导傅里叶方程。变为热传导傅里叶

6、方程。分子扩散与热传导是数学形式相同的两个过程。分子扩散与热传导是数学形式相同的两个过程。二阶线性抛物二阶线性抛物型偏微分方程型偏微分方程第四节第四节 分子扩散方程分子扩散方程推广到三维：推广到三维：故有故有用直角坐标表示用直角坐标表示时变项时变项分子扩散项分子扩散项扩散方程本质上是质量守恒定律在扩散问题上的体现扩散方程本质上是质量守恒定律在扩散问题上的体现第四节第四节 分子扩散方程分子扩散方程Fick定律：定律：在扩散特性各向同性的液体中，在在扩散特性各向同性的液体中，在x、y、z三个方向上，三个方向上，D为常数。为常数。在扩散特性各向异性的液体中在扩散特性各向异性的液体中第三节第三节 一维

7、扩散方程的基本解一维扩散方程的基本解&扩散方程的定解条件（初始条件、边界条件）。扩散方程的定解条件（初始条件、边界条件）。&解的形式：解析解、数值解。解的形式：解析解、数值解。&污染源（按空间）：点源、线源、面源、有限分布源、不存在污染源（按空间）：点源、线源、面源、有限分布源、不存在绝对的点源、无限长线源、无限大面源，只是一种近似处理。绝对的点源、无限长线源、无限大面源，只是一种近似处理。&污染源（按时间）：污染源（按时间）：瞬时源瞬时源、时间连续源（、时间连续源（事故排放事故排放、正常排、正常排放）。放）。&瞬时源是指污染物在瞬时内排放入水域，实际上一种近似，如瞬时源是指污染物在瞬时内排放

8、入水域，实际上一种近似，如热核武器试验的核污染或者油轮事故突然泄漏的油污染。热核武器试验的核污染或者油轮事故突然泄漏的油污染。&连续源又分为恒定和非恒定源。连续源又分为恒定和非恒定源。&污染物扩散：根据水域是几维，对应一维、二维、三维扩散方污染物扩散：根据水域是几维，对应一维、二维、三维扩散方程。程。第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 第三节第三节 一维扩散方程的基本解一维扩散方程的基本解 第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 集中投入的情况，在集中投入的情况，在t=0时刻，在原点瞬时投入质量为时刻，在原点瞬时投入质量为M M的的扩散质，分析以后任意时刻在

9、无界空间中的浓度分布，这扩散质，分析以后任意时刻在无界空间中的浓度分布，这是扩散方程的最基本的解。是扩散方程的最基本的解。是在静止水域中的扩散，而且是瞬时集中源与坐标原点重是在静止水域中的扩散，而且是瞬时集中源与坐标原点重合的一维扩散方程的特解。因为扩散方程是线性的，在线合的一维扩散方程的特解。因为扩散方程是线性的，在线性的边界条件下，可用这个特解式叠加来构造其他定解条性的边界条件下，可用这个特解式叠加来构造其他定解条件下的解。件下的解。0 x-x第三节第三节 一维扩散方程的基本解一维扩散方程的基本解 第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 瞬时瞬时单位单位平面源的扩散平面源的

10、扩散瞬时源：瞬时源：t=0t=0时，在原点瞬时集中投放质量为时，在原点瞬时集中投放质量为M M的扩散质。的扩散质。1 1、一根无限长断面均匀的直水管，截面积是一个、一根无限长断面均匀的直水管，截面积是一个单位单位2 2、垂直管轴，瞬时投入一包含质量、垂直管轴，瞬时投入一包含质量M M的薄片红色染液的薄片红色染液3 3、染液薄片充满了整个断面、染液薄片充满了整个断面4 4、染料只沿长度方向扩散、染料只沿长度方向扩散令染液投入点为坐标原点令染液投入点为坐标原点0 x-x瞬瞬时时点点源源或或称称瞬瞬时时无无限限平平面面源源在在无无界界空空间间的的定定解解条条件件下下的的解析解。定解条件在数学上表达为

11、：解析解。定解条件在数学上表达为：c(x,0)=M(x)狄拉克（狄拉克（Dirac）函数函数 当当t=0时，在通过时，在通过x=0处且与处且与x轴垂直的平面上，污染物质量轴垂直的平面上，污染物质量为为M,它位于它位于x=0处以无限大的浓度强度浓缩在无限小的空间处以无限大的浓度强度浓缩在无限小的空间（2）边界条件：边界条件：c(,t)=0,c(,t)/x=0（1）初始条件：初始条件：一维分子扩散方程：一维分子扩散方程：1.1.定解条件定解条件第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 M(x)表示质量表示质量M集中于微小容积集中于微小容积内。相对概念。例如把一小桶颜色内。相对概念。例

12、如把一小桶颜色水倾注到大河里，可以认为起始浓水倾注到大河里，可以认为起始浓度集中于微小体积内。度集中于微小体积内。物理含义：物理含义：2.解析方法：如拉普拉斯变换、分离变量法和量纲分析法解析方法：如拉普拉斯变换、分离变量法和量纲分析法量纲分析量纲分析，物理方程中各项物理量的量纲之间存在的规律：，物理方程中各项物理量的量纲之间存在的规律：l量纲和谐性，物理方程中各项的量纲应当相同；量纲和谐性，物理方程中各项的量纲应当相同；l任一有量纲的物理方程可以改写为无量纲项组成的方程而任一有量纲的物理方程可以改写为无量纲项组成的方程而不会改变物理过程的规律性；不会改变物理过程的规律性；l物理方程中各物理量之

13、间的规律性以及相应各量纲之间的物理方程中各物理量之间的规律性以及相应各量纲之间的规律性，不会因所选的基本量纲不同而发生改变。规律性，不会因所选的基本量纲不同而发生改变。定律（布金汉定律）：任何一个物理过程，包含有定律（布金汉定律）：任何一个物理过程，包含有k+1个有个有量纲的物理量，如果选择其中量纲的物理量，如果选择其中m个作为基本物理量，那么该物个作为基本物理量，那么该物理过程可以由理过程可以由(k+1)-m个无量纲数所组成的关系来描述。个无量纲数所组成的关系来描述。第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 式中式中:f为待定函数，在上式中写上为待定函数，在上式中写上4和和4，

14、目的是使最终，目的是使最终的解较为简明的解较为简明;M是全部污染物的质量，量纲是是全部污染物的质量，量纲是MM假设有函数：假设有函数：F(c,M,D,x,t)=0 方程线性方程线性 利用利用定律，选定律，选c、D、t为基本变量，可得：为基本变量，可得：从物理概念上分析，浓度从物理概念上分析，浓度c是是M、D、x、t的函数的函数第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 一维一维 扩散中，浓度的量纲扩散中，浓度的量纲 ML-1,浓度浓度c应与应与M除以某一特征长度成除以某一特征长度成正比。正比。是一个合适的特征长度是一个合适的特征长度 进一步令进一步令 ，有，有:。边界条件由原来的边

15、界条件由原来的c(,t)=0,c(,t)/x=0f()=0()=0，df()/()/dh=0h=0以以f的边界条件代入上式得的边界条件代入上式得k1 1=0=0，故上式变为，故上式变为:第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 设变量设变量进一步令进一步令 ，有，有:即即=常数常数k1,因此有：因此有：它的通解为：它的通解为：确定待定函数确定待定函数f为任何时刻源点浓度（坐标为任何时刻源点浓度（坐标原点与源点重合的情况下）原点与源点重合的情况下）根据污染物质的质量守恒定律，有根据污染物质的质量守恒定律，有对上式分别通过求对上式分别通过求t0、x0和和t0（x0）的极限，可）的极限

16、，可得到得到c=和和c=0，这说明了该解也是满足初始条件的。，这说明了该解也是满足初始条件的。此外，上式虽然是对此外，上式虽然是对x0的定解条件求解，但也可用于的定解条件求解，但也可用于x0情形。情形。,推出推出k0=1第五节第五节 一维扩散方程的基本解一维扩散方程的基本解 瞬时点源一维无界空间的浓度分布瞬时点源一维无界空间的浓度分布瞬时点源一维无界空间的浓度瞬时点源一维无界空间的浓度场在任一时刻场在任一时刻t沿沿x轴是正态分布，轴是正态分布，随时间随时间t的增加，浓度的峰值的增加，浓度的峰值Cm变小，而扩散的范围变宽，分变小，而扩散的范围变宽，分布曲线趋于平坦。布曲线趋于平坦。第五节第五节

17、一维扩散方程的基本解一维扩散方程的基本解 浓度分布符合正态分布（即高斯分布）浓度分布符合正态分布（即高斯分布）污染源点和坐标原点重合的情况污染源点和坐标原点重合的情况1、浓度对距离的各阶矩定义浓度对距离的各阶矩定义 零阶矩零阶矩 一阶矩一阶矩 二阶矩二阶矩对原点的任意对原点的任意p阶矩阶矩 对瞬时点源来说，零阶矩对瞬时点源来说，零阶矩 M0 0=全部扩散质的质量，对任全部扩散质的质量，对任意时刻意时刻M0 0是一常数，但一般情况下，矩都是时间的函数。是一常数，但一般情况下，矩都是时间的函数。各式的右端可供各式的右端可供当具有实验资料当具有实验资料时，计算浓度各时，计算浓度各阶矩之用。阶矩之用。

18、第四节第四节 浓度分布的各阶矩浓度分布的各阶矩第六节第六节 浓度分布的各阶矩浓度分布的各阶矩2、浓度分布的统计特征值浓度分布的统计特征值（1）浓度分布的距离均值（数学期望）浓度分布的距离均值（数学期望）表示浓度分布曲线重心距表示浓度分布曲线重心距x坐标原点的水平距离坐标原点的水平距离，当曲线对称于，当曲线对称于c轴时轴时 x=0=0。（2）浓度分布的距离方差浓度分布的距离方差 2 2 表示浓度分布对于平均浓度值的离散程度，表示浓度分布对于平均浓度值的离散程度，2 2值愈大，值愈大，分布曲线愈平坦。分布曲线愈平坦。第六节第六节 浓度分布的各阶矩浓度分布的各阶矩质量中心坐标质量中心坐标 x对于正态

19、分布曲线（标准）有：对于正态分布曲线（标准）有：将瞬时点源的解代入将瞬时点源的解代入M2 2，得距离方差：，得距离方差：当当已已求求得得 ，可可用用上上式式反反求求D。由由于于D是是常常数数，将将上上式式对对t t求导，有：求导，有：称为矩法公式，可以称为矩法公式，可以差分差分代替微分。代替微分。对任何其它分布，只要在无界空间情况下满足边界条件：对任何其它分布，只要在无界空间情况下满足边界条件：或或仍存在仍存在上式表明上式表明 方差与扩散历时方差与扩散历时t t成正比。凡符合这个规律的成正比。凡符合这个规律的扩散，都称为扩散，都称为费克型扩散费克型扩散。第六节第六节 浓度分布的各阶矩浓度分布的

20、各阶矩曲线的分布区间曲线的分布区间-2 ,2 -2 ,2 占总面积的占总面积的95%源与坐标原点不重合源与坐标原点不重合源与坐标原点重合源与坐标原点重合-2 2 ,+2 2 证明此结论证明此结论（3）三阶中心矩三阶中心矩表示曲线偏斜度：表示曲线偏斜度：=0 左右对称左右对称;正态分布正态分布 0左右不对称，长尾伸向正轴方向；左右不对称，长尾伸向正轴方向；0，长尾伸向负轴方向。，长尾伸向负轴方向。=0=0 00 0 0图图 对浓度分布图形的影响对浓度分布图形的影响第六节第六节 浓度分布的各阶矩浓度分布的各阶矩偏态系数偏态系数（4）四阶中心矩四阶中心矩表示曲线峰态或平坦度的一个指标，值愈大表示表示

21、曲线峰态或平坦度的一个指标，值愈大表示峰型愈大。峰型愈大。第六节第六节 浓度分布的各阶矩浓度分布的各阶矩2 2、对静止和动水环境中射流的一些基本原理、基本、对静止和动水环境中射流的一些基本原理、基本规律的研究现状。规律的研究现状。源与坐标原点重合时，浓度源与坐标原点重合时，浓度曲线的分布区间曲线的分布区间-2 2,2,2 范围内，分布曲线与范围内，分布曲线与x x轴所围面积轴所围面积占总面积的占总面积的95%。1、证明此结论、证明此结论作业作业第五节第五节 一维扩散方程的若干定解条件下的解一维扩散方程的若干定解条件下的解设只当设只当t=0时在时在x=处投放污染物质（瞬时点源）处投放污染物质（瞬

22、时点源）初始条件：初始条件：c(x,0)=M(x-)边界条件边界条件：c(,t)=0 第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解如果示综物质如果示综物质M不是集中到一处，而是非均匀地分布在一定不是集中到一处，而是非均匀地分布在一定范围上同时瞬时投放，这就是瞬时投放源，这种情况可考范围上同时瞬时投放，这就是瞬时投放源，这种情况可考虑为若干个瞬时集中源的迭加，按迭加原理求解。虑为若干个瞬时集中源的迭加，按迭加原理求解。现将初始条件改为：现将初始条件改为：c(x,0)=f(x),-x 其中其中f(x)为任意给定的函数，亦即该初始分布是沿无限为任意给定的函数，亦即该

23、初始分布是沿无限长直线上给定的浓度为长直线上给定的浓度为f()，d微小长度上投放示踪质的质微小长度上投放示踪质的质量为量为M=f()d。位于位于处由该微小污染单元的扩散处由该微小污染单元的扩散而导致在时刻而导致在时刻t位于位于x的浓度应为的浓度应为:用一系列质量为用一系列质量为f()d的的团块来求浓度分布团块来求浓度分布第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解下面讨论下面讨论f()为常数的两种特殊情况：为常数的两种特殊情况：单侧阶梯浓度函数的浓度分布单侧阶梯浓度函数的浓度分布1.当当f(x)为阶梯函数：为阶梯函数：该问题的物理模型可认该问题的物理模型可认为

24、是在一条无限长的等为是在一条无限长的等截面管（渠）的静水中，截面管（渠）的静水中，左端（左端（x0）为清水，现闸）为清水，现闸门突然打开，左边的污门突然打开，左边的污染物质向右边扩散。解染物质向右边扩散。解的形式为：的形式为：第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解t=0时时 取变换取变换式中：式中：erf(z)为为z的的误差函数，误差函数，erfc(z)为为z的的余误差函数，即余误差函数，即第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解误差函数的值可查误差函误差函数的值可查误差函数数值表或计算软件得到数数值表或计算软件得到

25、误差函数的定义误差函数的定义:从而有：从而有：余误差函数的定义余误差函数的定义:第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解000误差函数的计算是把被积函数展开为麦克劳林级数，然后逐项积分误差函数的计算是把被积函数展开为麦克劳林级数，然后逐项积分第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解 取变换取变换=x-,有有该问题的物理模型可认为是在一条无限长的等截面渠该问题的物理模型可认为是在一条无限长的等截面渠（管）的静水中，突然发生事故，在渠中出现一段污染（管）的静水中，突然发生事故，在渠中出现一段污染源而向两端扩散的情形。解的形

26、式为：源而向两端扩散的情形。解的形式为：2.当当f(x)为阶梯函数：为阶梯函数：x=0 x=x1x=-x1初始浓度分布图初始浓度分布图第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解再取变换再取变换 ,有有t=0时时双侧阶梯浓度函数的浓度分布双侧阶梯浓度函数的浓度分布第七节第七节 一维扩散方程空间瞬时线源的解析解一维扩散方程空间瞬时线源的解析解随着随着 增大，浓度增大，浓度分布曲线愈平坦化。分布曲线愈平坦化。第六节第六节 一维扩散方程的时间连续源的解析解一维扩散方程的时间连续源的解析解 一、时间连续点源一、时间连续点源 在流场的某一点上，连续不断地投入浓度为在流场

27、的某一点上，连续不断地投入浓度为c c0 0（常数）（常数）的污染物质，即时间连续恒定点源。的污染物质，即时间连续恒定点源。如果一维扩散区域无限长，则可将投放点位置取为坐标原如果一维扩散区域无限长，则可将投放点位置取为坐标原点，初始条件点，初始条件c(0,t)=c0 在在x=0处浓度突然从零增加到处浓度突然从零增加到c0，以后保，以后保持不变，持不变，无限边界条件无限边界条件c(,0)=0 初瞬时初瞬时t=0，沿，沿x轴各处的浓度均为零轴各处的浓度均为零 本问题的解也是一个有用的基本解，可以用来构造其他某本问题的解也是一个有用的基本解，可以用来构造其他某些问题的解。些问题的解。第八节第八节 一

28、维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解1、点源处给定投放浓度、点源处给定投放浓度第八节第八节 一维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解边界条件为边界条件为f(0)=1,f()=0,显然有显然有c(-x,t)=c(x,t),解对称于原点，解对称于原点，只需沿只需沿x正向求解。正向求解。借助借助量纲分析法量纲分析法来求解浓度分布来求解浓度分布c(c(x,t)显然，显然，c c与与c c0 0,D,x和和t有关，利用有关，利用定理定理,选选 c、D和和t为基本变量，可得如下关系为基本变量，可得如下关系式：式：式中：式中：f是某一待确定的函数。令是某一待确定的函

29、数。令 ,有有二阶变系数齐次常微分方程二阶变系数齐次常微分方程代入扩散方程代入扩散方程等强度连续点源的浓度分布等强度连续点源的浓度分布1 1、点源处给定投放浓度问题的解、点源处给定投放浓度问题的解 在任一时刻在任一时刻t总的浓度是总的浓度是t以前全以前全部时段内浓度分布的总和，由上部时段内浓度分布的总和，由上式通过积分可得的解：式通过积分可得的解：用于求解浓度用于求解浓度c0(t)的迭加法的迭加法第八节第八节 一维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解更一般的情形是更一般的情形是c0不是常数，不是常数，即不是时间连续恒定点源。而即不是时间连续恒定点源。而是是c0随时间随时间t

30、 t而变，即而变，即c0=c0(t)t)。则可以看作无数不同强度的。则可以看作无数不同强度的瞬时源产生的扩散在时间上扩散的结果。瞬时源产生的扩散在时间上扩散的结果。当时间增加当时间增加tt，位于，位于x=0处的浓度增量为处的浓度增量为 ，由于浓度增，由于浓度增量的关系，相应的扩散结果可借助上式表示为量的关系，相应的扩散结果可借助上式表示为:由于是时间连续点源，故可得由于是时间连续点源，故可得如果如果 为常数，上式变为为常数，上式变为 第八节第八节 一维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解以下讨论在以下讨论在x=0处，给定单位时间投入的污染物质量的速度处，给定单位时间投入的污

31、染物质量的速度 (简称质量投放率简称质量投放率)，即在，即在tt时间投放质量为时间投放质量为 。此时。此时，根据瞬时点源的解式可得在瞬时，根据瞬时点源的解式可得在瞬时t t投放质量投放质量 的浓度场：的浓度场：2、点源处给定投放质量、点源处给定投放质量第六节第六节 一维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解时间连续点源的浓度分布时间连续点源的浓度分布污染的范围和污染的范围和浓度均随时间浓度均随时间的增加而增大的增加而增大 第八节第八节 一维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解通过数值积通过数值积分进行计算分进行计算二、时间连续线源二、时间连续线源第八节第八

32、节 一维扩散方程时间连续源的解析解一维扩散方程时间连续源的解析解全部分布连续源一维扩散产生的浓度应将上式几分二次成为全部分布连续源一维扩散产生的浓度应将上式几分二次成为如果连续源不集中在原点，而是分布在沿如果连续源不集中在原点，而是分布在沿x轴一定范围如轴一定范围如a x1 b之上，则加入的扩散质最一般的情况是时间和空之上，则加入的扩散质最一般的情况是时间和空间的函数，设间的函数，设f(x,t)为在单位时间内单位体积上投放的污为在单位时间内单位体积上投放的污染物质质量，在染物质质量，在x=x x处处d dx x上，于时刻上，于时刻t t在在d dt t时间内加入扩时间内加入扩散质的量为散质的量

33、为 M=f(x x,t t)d dx xdt t，一维扩散经时间，一维扩散经时间(t-t t)在在x处得处得浓度为：浓度为：第七节第七节 有界一维扩散和叠加方法有界一维扩散和叠加方法 以上讨论的全是无限空间的扩散，实际河渠或水库湖泊以上讨论的全是无限空间的扩散，实际河渠或水库湖泊都有岸和底存在，污染物质在河渠、水库中扩散至边界时，都有岸和底存在，污染物质在河渠、水库中扩散至边界时，有两种可能，一种是扩散物质到达边界后被边界吸收或粘结有两种可能，一种是扩散物质到达边界后被边界吸收或粘结在边界上，称为在边界上，称为完全吸收完全吸收；另外一种情况是遇到边界就反射；另外一种情况是遇到边界就反射回去，称

34、为回去，称为完全反射完全反射。介于两种状态之间的为不完全吸收和。介于两种状态之间的为不完全吸收和不完全反射，这在不完全反射，这在实际中居多。实际中居多。显然，吸收和反射与污染物显然，吸收和反射与污染物性质和边界的性质有关。当然，最不利的情况是发生完全反性质和边界的性质有关。当然，最不利的情况是发生完全反射。如果天然河岸是土壤，就不能看做是最不利的情况。射。如果天然河岸是土壤，就不能看做是最不利的情况。下面仅研究完全反射的情况。下面仅研究完全反射的情况。第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法 讨论最简单的情况：当讨论最简单的情况：当t=0时，在时，在x=0处与处与x轴垂直的单

35、轴垂直的单位面积上，投放的污染物质量为位面积上，投放的污染物质量为M。在正方向的边界为无穷。在正方向的边界为无穷远，但在远，但在x=-L处有一阻止物质扩散的壁存在，并设该壁不吸处有一阻止物质扩散的壁存在，并设该壁不吸收扩散物质（完全反射），则任意时刻通过该岸壁的示踪量收扩散物质（完全反射），则任意时刻通过该岸壁的示踪量的净通量为零。的净通量为零。第七节第七节 有界一维扩散和叠加方法有界一维扩散和叠加方法 对扩散被各种边界所限制的问题，通常运用叠加原理来对扩散被各种边界所限制的问题，通常运用叠加原理来解决。因为扩散方程是线性的，如果边界条件也是线性的，解决。因为扩散方程是线性的，如果边界条件也是

36、线性的，则可以叠加任意数量的单独解，从而构成新的解。则可以叠加任意数量的单独解，从而构成新的解。假设边界为完全反射壁，即不吸收扩散物质。假设边界为完全反射壁，即不吸收扩散物质。一、一边反射的瞬时点源情形一、一边反射的瞬时点源情形边界条件：壁面上的浓度梯度必须是零，由费克定律得到边界条件：壁面上的浓度梯度必须是零，由费克定律得到:初始条件：初始条件：第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法一边侧壁的像源法一边侧壁的像源法像源法：像源法：设想有一平面镜位于固体边界处，在平面镜后面有一反射设想有一平面镜位于固体边界处，在平面

37、镜后面有一反射源（又称像源），反射源到真源的距离源（又称像源），反射源到真源的距离x=-2L，像源的强度与，像源的强度与真源的强度相同，标准差也相同，因而像源在边界上的通量与真源的强度相同，标准差也相同，因而像源在边界上的通量与真源在边界上的通量大小相等，方向相反，故形成边界上扩散真源在边界上的通量大小相等，方向相反，故形成边界上扩散物质的通量为零。物质的通量为零。第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法一边侧壁的像源法一边侧壁的像源法 真源与像源相距为真源与像源相距为2L，在，在x轴上任意点的浓度轴上任意点的浓度应该为由像源和真实源各自产生的浓度之和，即应该为由像源和真实源

38、各自产生的浓度之和，即第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法一边侧壁的像源法一边侧壁的像源法 对于完全反射的边界，对于完全反射的边界，在反射壁边界处的浓度等于在反射壁边界处的浓度等于不存在该壁时的两倍。不存在该壁时的两倍。当当x=-L时时（即固体（即固体边边界上），其界上），其浓浓度度为为二、两边反射的瞬时点源情形二、两边反射的瞬时点源情形 在在x=-L和和x=L均有完全反射壁：均有完全反射壁：无穷多像源各自的浓度分布叠加便得到问题的解：无穷多像源各自的浓度分布叠加便得到问题的解：两面侧壁的像源法两面侧壁的像源法第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法在在

39、x=-2L的的像像源源的的扩扩散散又又会会在在x=L处处产产生生一一个个正正的的浓浓度度梯梯度度，需需要要在在x=4L处放一个像源。处放一个像源。在岸壁处，即在在岸壁处，即在x=L处要求浓度为零，在处要求浓度为零，在x=2L处必须放置强度为负值的瞬时源，在处必须放置强度为负值的瞬时源，在x=4L处处放置强度为正值的瞬时源，如此类推。完整解为放置强度为正值的瞬时源，如此类推。完整解为第九节第九节 有界一维扩散和叠加方法有界一维扩散和叠加方法另一种边界条件另一种边界条件在实际问题中，因河流较宽，在实际问题中，因河流较宽，L比较大，一般只考虑一、比较大，一般只考虑一、两次反射就可以满足实际需要，如两

40、次反射就可以满足实际需要，如n=0,1，2。第八节第八节 二维和三维扩散方程的解析解二维和三维扩散方程的解析解 一、瞬时点源一、瞬时点源二维扩散方程为：二维扩散方程为：上式中的上式中的Dx和和Dy分别为分别为x和和y方向的扩散系数，虽然在方向的扩散系数，虽然在分子扩散中，分子扩散中，Dx=Dy=D，但因为我们将来可以借用该方程，但因为我们将来可以借用该方程的解来解决某具有非各向同性性质的紊流扩散问题，所以的解来解决某具有非各向同性性质的紊流扩散问题，所以在这里以在这里以DxDy进行讨论。进行讨论。第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解c(x,y,0)=M(x

41、)(y)（2）边界条件：）边界条件：c(,y,t)=0,c(,t)/x=0c(x,t)=0,c(,t)/y=0（1）初始条件：）初始条件：上式只有当两个括号内的量分别等于零才能得到满足，从上式只有当两个括号内的量分别等于零才能得到满足，从而得到两个一维扩散方程，它们的瞬时点源无界空间的解均具而得到两个一维扩散方程，它们的瞬时点源无界空间的解均具有扩散方程基本解的形式，将这两个解相乘，就得到解答。有扩散方程基本解的形式，将这两个解相乘，就得到解答。利用利用“乘积法则乘积法则”求解：则本问题的解可以表为求解：则本问题的解可以表为c(x,y,t)=c1(x,t)c2(y,t)式中式中c1不依赖于不依

42、赖于y,c2不依赖于不依赖于x，代入扩散方程，故有，代入扩散方程，故有第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解当当Dx=Dy=D时，上式变为：时，上式变为：在在z轴上单位面积轴上单位面积上的质量上的质量 对于一维扩散，浓度对于一维扩散，浓度c c的单位是单位长度的质量；对于二的单位是单位长度的质量；对于二维扩散它是单位面积上的质量；对于三维扩散则是单位体积内维扩散它是单位面积上的质量；对于三维扩散则是单位体积内的质量。的质量。可将可将“乘积法则乘积法则”求解的方法推广到瞬时点源无界空间的三求解的方法推广到瞬时点源无界空间的三维扩散。三维扩散方程为：维扩散。三维

43、扩散方程为：当当Dx=Dy=Dz=D时，时，r2=x2+y2+z2,有有第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解r是自点源起算的距离。在是自点源起算的距离。在 处上式指数为处上式指数为0.01，即这里的浓，即这里的浓度等于原点瞬时浓度值的百分之一。度等于原点瞬时浓度值的百分之一。初始条件：初始条件：c(x,y,z0)=M(x)(y)(z)二、瞬时无限长线源二、瞬时无限长线源瞬时无限长线源情形瞬时无限长线源情形第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解 一个瞬时无限长线源是沿一无限长直线上的每一个单位长一个瞬时无限长线源是沿一无限

44、长直线上的每一个单位长度瞬时投放质量为度瞬时投放质量为mz z所构成的，所构成的，mz z的量纲为的量纲为MLML-1-1。对于一个。对于一个沿沿z轴分布的无限长线源来讲，根据三维瞬时点源的解可得由轴分布的无限长线源来讲，根据三维瞬时点源的解可得由于于h h处的点源处的点源mz zdh h所产生的所产生的P点（点（x,y,z）处的浓度为）处的浓度为与瞬时点源的二维情形相同与瞬时点源的二维情形相同无限长线源无限长线源h h积分从积分从负无穷到正无穷负无穷到正无穷令令变更上下限变更上下限其解可由瞬时点源的解迭加得到其解可由瞬时点源的解迭加得到单位面积瞬时引入质量为单位面积瞬时引入质量为m,对在对在

45、yz平面上的一个平面源平面上的一个平面源来讲，由位于来讲，由位于x x处沿处沿z方向单方向单位宽度上质量位宽度上质量mz z=md dx x的无限的无限长线源在长线源在P点（点（x,y,z）处产生）处产生的浓度：的浓度：三、瞬时无限平面源三、瞬时无限平面源 第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解于是由无限平面源在于是由无限平面源在P点处产生的浓度为：点处产生的浓度为：当当Dx=Dy=Dz=D时，有时，有瞬时无限平面源情形瞬时无限平面源情形l瞬时无限平面的分子瞬时无限平面的分子扩散只沿与该平面垂扩散只沿与该平面垂直的方向进行，是一直的方向进行，是一维扩散维扩散

46、l一维扩散问题中，点一维扩散问题中，点源就是无限平面源。源就是无限平面源。设在坐标原点处（设在坐标原点处（x=y=z=0=0），单位时间内投放的污染物），单位时间内投放的污染物质量为质量为M（常数）。在（常数）。在d dt t的微小时间内，投放质量为的微小时间内，投放质量为M dt t，将，将每一个每一个M dt t看作是一个瞬时点源，借助瞬时点源的解在瞬时看作是一个瞬时点源，借助瞬时点源的解在瞬时投入质量投入质量M dt t的浓度场的浓度场:四、三维时间连续恒定点源四、三维时间连续恒定点源 第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解由于时间连续源，可对上式对时

47、间进行积分，有由于时间连续源，可对上式对时间进行积分，有当当t t=0=0，有，有当当t t=t，有，有u=当当Dx=Dy=Dz=D时，上式简化为时，上式简化为:第十节第十节 二维和三维扩散方程的某些解析解二维和三维扩散方程的某些解析解第九节第九节 随机游动法随机游动法(Random walk)在以上各节中，是用欧拉（在以上各节中，是用欧拉（Euler）的观点和确定性的数）的观点和确定性的数学方法来确究费克扩散的。但是，也可以采用拉格朗日学方法来确究费克扩散的。但是，也可以采用拉格朗日的观点，跟踪污染物质点的不规则运动，以及采用概率的观点，跟踪污染物质点的不规则运动，以及采用概率统计的方法来进

48、行研究。统计的方法来进行研究。第十一节第十一节 随机游动法随机游动法进一步考虑有大量污染物质点当进一步考虑有大量污染物质点当t=0=0在位置在位置 处释放，于是，处释放，于是，当时刻当时刻t在位置在位置 处的浓度可表示为：处的浓度可表示为：污染物质点由于受到分子运动的作用而作不规则的运动污染物质点由于受到分子运动的作用而作不规则的运动:第十一节第十一节 随机游动法随机游动法质点的位移质点的位移 是一个是一个随机过程，该过程可以通随机过程，该过程可以通过概率密度过概率密度 表达表达:l某一质点当时刻某一质点当时刻t=0时处时处于位置于位置l当时刻当时刻t时则处于位置时则处于位置根据大数定律，概率

49、根据大数定律，概率 可近似地由下式确定：可近似地由下式确定：马尔可夫马尔可夫（Markow）链过程：由于受到液体分子碰撞，质点可能在时刻链过程：由于受到液体分子碰撞，质点可能在时刻t落入某一指定的体积元，也可能不落入该体积元。而在每受一次碰撞之落入某一指定的体积元，也可能不落入该体积元。而在每受一次碰撞之后，下一步走到什么位置，只是与当前的位置有关，而与前一个位置无后，下一步走到什么位置，只是与当前的位置有关，而与前一个位置无关。也可以说，下一步处于某一位置的概率，只与质点位于当前位置的关。也可以说，下一步处于某一位置的概率，只与质点位于当前位置的概率有关，而与质点位于前一个位置的概率无关，亦

50、即质点的运动已失概率有关，而与质点位于前一个位置的概率无关，亦即质点的运动已失去历史的影响。去历史的影响。设设P(x,t)为一个质点当某一时刻为一个质点当某一时刻t位于某一位置位于某一位置x的概率的概率由全概率公式得到质点当时刻由全概率公式得到质点当时刻t位于位于x处的概率处的概率:第十一节第十一节 随机游动法随机游动法质点的一维随机运动质点的一维随机运动第十一节第十一节 随机游动法随机游动法因为质点只有向左或向右运动的两种可能性，故有因为质点只有向左或向右运动的两种可能性，故有p1+p2=1，且，且p1=p2，因，因此上式右边的末项为零。此上式右边的末项为零。进一步对上式取进一步对上式取0和