瞬时点源的扩散.pdfVIP

这一章内容紧随第一章，研究瞬时点源释放后所演化出来的高斯（Gaussian ）浓度剖面。强 调一些重要的高斯分布特征（以及它们的应用）。讨论点源在3 维方向的扩展，提示沿任意 直线的云团横断面浓度剖面将具有高斯分布形式。动画显示了具有各向异性扩散率的瞬时点 源扩展，突出高斯分布的长度大小与均方根D 的函数关系。 3 ．瞬时点源的扩散 质量守恒方程也被认为是输运方程，因为它描述了流体系统中标量物质种类的输运。在这一 节和接下来的一节中，我们考虑输运方程的解析解，它描述了在不同边界条件和初始条件下， 质量的运动走向。对于具有各向同性和均匀扩散系数的不可压缩流动，我们从以下的输运方 程形式开始。这里我们也将忽略掉源/汇项。 ∂C ∂C ∂C ∂C ∂ C ∂ C ∂ C 1 ∂t u ∂x v ∂y w ∂z D∂x ∂y ∂z 第一个解描述一小片质量M 在t ＝0 时刻，瞬时释放到静止流体中（u,v,w ＝0 ）的演化过程。 我们考虑沿 x 方向的一维系统，也就是系统在 y 和 z 方向上是均匀的，这样∂/ ∂y 0， ∂/ ∂z 0 。这样一个系统的模型如下所示。释放一小片质量（灰色）到一个狭长的管道中。 质量在y －z 平面内是均匀分布的，并可以忽略x 方向的尺度，所以初时浓度C（x ）＝Mδ （x ）， 其中δ （）是狄拉克δ函数（Dirac delta function ）。此系统的输运方程由（2 ）式给出。 ∂C ∂ C 2 D , 在 t 0 时刻，x 0 处，释放M ∂t ∂x 解C （x ，t ）可以按如下方式求得。从（2 ）式中我们预测浓度是M ，D[L2T-1]，x[L]和t[T] 的函数，方括号中的单位根据质量[M]，长度[L]和时间[T]表示。在某个一维系统中，浓度 具有单位长度上质量的单位，C[ML-1] 。由量纲分析{见白金汉π定律（Buckingham pi theorem ） 和重复变量方法（Method of Repeating Variables ）}，以上各参数可形成以下的无量纲组合形 式。 C x 3 M/√4πDt f √4Dt 为了方便加入因子 4π和 4 ，因为它们出现在最终的解中。最初这些因子是通过反复的试验 得到的。定义 ＝x/√4Dt，接着将（3 ）式代入（2 ）式， d df 4 dη dη 2fη 0 设定括号中的项为零，得到（4 ）的一个解， df 5 2fη 0 dη （5 ）式所具有的解是 exp η 6 f A 其中A0 是一个任意常数。（6 ）式给了我们（3 ）式的一个特定形式，这就是（2 ）式的解答。 A x exp 7 C x, t √4πDt 4Dt 因为总质量在整个时间上是守恒的，对于所有时间 8 M Cdx 所以A0=M 。总之我们有以下的解。 一维，瞬时点源 在t ＝0 时刻，x ＝0 处释放M 。 初始条件，C （x ）＝Mδ （x ）。 M x 9 C x, t √4πDt exp 4Dt M/L 注意到以上的解对于一维数学系统是正确的，其中 C[M/L] 。在真实（三维）系统中，浓度 一定具有M/L3 的单位。为了将（9 ）中给出的数学解答转换到真实空间中，除以忽略掉的尺 度，这就是系统在y －z 平面中的横截面面积Ayz 。 M x