扩散模型解析

作者：金良（[email protected]） csdn博客：http://blog.csdn.net/u012176591

PDF版本下载地址下载链接

1 若干数学工具

1.1 三维傅里叶变换的性质

1. 傅里叶变换
   

   (1)

2. 傅里叶逆变换
   

   (2)

3. 线性性质
   

   (3)

4. 微分性质
   

   (4)

1.2 一阶线性微分求解通式

形如 +PC=Q的公式叫做一阶线性微分方程，其通解形式用Maple求解如图1

图 1: 一阶线性偏微分方程求解

从图1可以看到其通解形式为

(5)

1.3 有限差分

现有连续性函数C=f(t),取相邻两点的距离均为h的横坐标点t_i-2,t_i-1,t_i,t_i+1，将这些点的纵坐标分别标记为C_i-2,C_i-1,C_i,C_i+1在i点附近把函数C=f(t)展开为泰勒级数：

(6)

其中，_i,_i,分别表示各阶导数在t=t_i处的值

      在点i-1处有     t=t_i-1 =x_i-h

      在点i+ 1处有     t=t_i+1 =x_i +h

把上述两式分别带入函数C泰勒级数展开式中，求得C在点t_i-1和t_i+1处的值分别为

(7)

如果h足够小可以忽略h³及更高项，根据上述两式解出_i,_i的值：

(8)

1.4 三维矩阵的理解

可以把三维矩阵的第一项指标看成“行”,第二项指标看成是“列”，第三项指标看成是“页”

1.5 通量,散度和高斯公式

1.5.1 通量

设有向量场

(9)

其中函数PQR均有一阶连续偏导数，Σ是场内的一片有向曲面，n是Σ在点(x,y,z)处的单位法向量，则积分

(10)

称为向量场A通过曲面Σ向着制定侧的通量（或流量）

有两类曲面积分的关系，通量又可表达为

(11)

1.5.2 散度

对于向量场A(x,y,z) = P(x,y,z)i+Q(x,y,z)j+R(x,y,z)k， + +叫做向量场A的散度，记作divA，即

(12)

利用向量微分算子∇，A的散度divA也可以表达为∇⋅A，即

(13)

1.5.3 高斯公式

设空间闭区域Ω是由分片光滑的闭曲面Σ所围成，函数P(x,y,z),Q(x,y,z),R(x,y,z)在Ω上具有一阶连续偏导数，则有

(14)

或

(15)

利用向量场的通量和散度，高斯公式可以写成如下的向量形式

(16)

或

(17)

2 模型

2.1 条件假设与符号约定

1. 河流的雷诺数已经达到一定的阈值（>2000），河水流动的方式为紊流而不是层流;

   流体力学中，雷诺数（Reynolds number）是流体惯性力与黏性力比值的量度，它是一个无量纲量

   雷诺数较小时，黏滞力对流场的影响大于惯性力，流场中流速的扰动会因黏滞力而衰减，流体流动稳定，为层流；反之，若雷诺数较大时，惯性力对流场的影响大于黏滞力，流体流动较不稳定，流速的微小变化容易发展增强，形成紊乱不规则的紊流流场

   层流与紊流两种流动形态传递动量热量和质量的方式不同：层流通过分子间相互作用，紊流主要通过质点间的混掺紊流的传递速率远大于层流水利工 ____________________________________________________________________________

   程所涉及的流动，一般为紊流

   ____________________________________________________________________________

2. 任取河流上一闭曲面S所围成水域为Ω，C(x,y,z,t)表示t时刻(x,y,z)位置点的污染物浓度;
3. 污染物因为河水流动和分子自由运动而发生扩散，D_x,D_y,D_z分别表示x,y,z方向的扩散系数;
4. 河水中的微生物分解等作用使河水具有自我净化能力，设自我降解系数为K(0< K < 1);
5. θ(x,y,z,t)是(x,y,z)位置点t时刻单位体积单位时间污染物的排放量;
6. 河水沿x,y,z三方向的流动速度u_x,u_y,u_z在一定时间范围内是恒定的

2.2 模型构建

为了使得公式看起来不那么臃肿，下面用C,θ分别代表C(x,y,z,t)和θ(x,y,z,t)

1. 由多重积分的意义可知，单位时间通过闭曲面S扩散进来的Ω的污染物的量为

   (18)

   由高斯定理可知：

   (19)

2. 由于河流是不断流动的，单位时间通过闭曲面S流出Ω的污染物的量为

   (20)

   同理，由高斯定理知：

   (21)

3. 由于河水自我降解，单位时间Ω水域的污染物减少量为

   (22)

4. 单位时间Ω水域的污染物的排放量为

   (23)

5. 单位时间Ω水域的污染物浓度变化引起的污染物的增量为

   (24)

由质量守恒定理得

(25)

将式19,21,22,23,24带入式25,得到河流水质污染的4D数学模型即四维水质模型为

(26)

3 模型求解

3.1 四维水质一般模型的求解

1. 对四维水质的一般模型，即式26左右两边同时进行傅里叶变换得

   (27)

2. 求解线性微分方程27得

   (28)

   用Maple求解过程如图2

   
   

   图 2: 一般模型的偏微分方程求解

3. 对式28两边同时进行傅里叶逆变换得

   (29)

式29中的是关于ω_x,ω_y,ω_z的函数，但是没有具体的表达式，所以式29不能求得解析解下面分析一般模型的几个特殊情况，也就是一般模型的简化模型

3.2 瞬时污染点源四维扩散模型

该模型对式26的模型进行了简化，如下

其 中 式 31 是 初 始 条 件 ， 式 32 是 边 界 条 件

1. 对式30进行傅里叶变换，得

   (34)

   对式31进行傅里叶变换，得

   (35)

2. 根据式34和式35求解得到Ĉ表达式

   (36)

   用Maple求解过程如图3

   
   

   图 3: 瞬时污染点模型的偏微分方程求解

3. 对式36进行逆傅里叶变换得到C

   (37)

   用Maple求解过程如图4

   
   

   图 4: 瞬时污染点模型的解析式

3.3 连续污染点源四维扩散模型

这里我首先介绍一下参考文献中的解法，然后讲一下我对这种解法的看法

3.3.1 积分积出来的解

如果污染物释放是以连续污染点源的方式进行，则形成的浓度相当于单位时间内连续释放的瞬时点源的积分，实际上相当于对瞬时污染点模型的解析解进行时间区间上的积分：

(38)

其中，v和c分别表示单位时间内排放的含有污染物的污水量和污染物的浓度

3.3.2 积分方法的正确性论证

从当前时间向前推时间T，然后将该段时间均分成N份，时间点标记为0 = t₀ < t₁ < t₂< t_N-1< t_N =T，且有t_i+1-t_i =h(i = 0,1N-1,h=).假设在每个时间点都将量为hvc的污染物排放到独立的具有相同参数的水流中（由于各个水流是独立的，所以互相不影响），则当前时间各个水流中的污染物的当前浓度可以用瞬时污染点四维扩散模型（式36）来表示为

(39)

则对所有水流的浓度求和，有

(40)

当N→∞时，求和就变成了积分

(41)

可以发现这就是式38

上述分析过程存在的问题在于，它把各个水流看成了独立的，各个水流中的污染物浓度不相互影响，因而最终的浓度变成了各个时间点投放污染物的河流的当前污染物浓度的简单叠加

下面将为什么能简单叠加为了方便，我们重新抄写式26

下面的几个等式的存在说明叠加方法的正确，也说明了式38的正确性

3.4 连续污染点源的扩散稳态模型

3.4.1 连续污染点源的一维扩散稳态模型

假定旨在x方向上存在污染物的浓度梯度，则稳态的一维模型

(46)

求解得到

(47)

用maple求解如图5

图 5: 连续污染点源的一维扩散稳态模型解析式

如果忽略弥散作用，则模型变成了

(48)

求解得到

(49)

用maple求解如图6

图 6: 连续污染点源的一维扩散稳态模型解析式

3.4.2 连续污染点源的三维扩散稳态模型

模型为

(50)

仿照式47结构，我们直接给出式50的解

(51)

用maple验证的过程如图7

图 7: 连续污染点源的三维扩散稳态模型解析式

• Two-Dimensional Fourier Transform
  http://fourier.eng.hmc.edu/e101/lectures/image_processing/node6.html 
• 2D Fourier Transforms
  http://www.cs.toronto.edu/~jepson/csc320/notes/linearFilters2.pdf 
• 基于MATLAB的四维水质模型仿真
  http://mag.ieechina.com/oa/pdfdow.aspx?Sid=201108079 
• 四维水质模型
  MATLAB在数学建模中的应用，卓金武
• 长江水质综合评价与预测的数学模型
  http://wenku.baidu.com/view/9ab43330f111f18583d05a11 
• 菲克定律
  http://zh.wikipedia.org/wiki/%E8%8F%B2%E5%85%8B%E5%AE%9A%E5%BE%8B 
• 高斯散度定理
  http://zh.wikipedia.org/wiki/%E9%AB%98%E6%96%AF%E6%95%A3%E5%BA%A6%E5%AE%9A%E7%90%86 
• 基于MATLAB的四维水质模型仿真
  https://raw.githubusercontent.com/wolfbrother/csdnPICs/master/water_pollute_model/%E5%9F%BA%E4%BA%8EMATLAB%E7%9A%84%E5%9B%9B%E7%BB%B4%E6%B0%B4%E8%B4%A8%E6%A8%A1%E5%9E%8B%E4%BB%BF%E7%9C%9F.pdf
• 大气污染的数学模型
  https://raw.githubusercontent.com/wolfbrother/csdnPICs/master/water_pollute_model/%E5%A4%A7%E6%B0%94%E6%B1%A1%E6%9F%93%E7%9A%84%E6%95%B0%E5%AD%A6%E6%A8%A1%E5%9E%8B.pdf
• 纺织品湿热舒适性影响因素的研讨
  https://raw.githubusercontent.com/wolfbrother/csdnPICs/master/water_pollute_model/%E7%BA%BA%E7%BB%87%E5%93%81%E6%B9%BF%E7%83%AD%E8%88%92%E9%80%82%E6%80%A7%E5%BD%B1%E5%93%8D%E5%9B%A0%E7%B4%A0%E7%9A%84%E7%A0%94%E8%AE%A8.pdf
• 长江水质综合评价与预测的数学模型
  https://raw.githubusercontent.com/wolfbrother/csdnPICs/master/water_pollute_model/%E9%95%BF%E6%B1%9F%E6%B0%B4%E8%B4%A8%E7%BB%BC%E5%90%88%E8%AF%84%E4%BB%B7%E4%B8%8E%E9%A2%84%E6%B5%8B%E7%9A%84%E6%95%B0%E5%AD%A6%E6%A8%A1%E5%9E%8B_%E9%9F%A9%E4%B8%AD%E5%BA%9A.pdf