Solution of the Advection-Diffusion Equation in Two Dimensions using the Deposition Velocity

技术摘要：利用沉降速度求解二维平流扩散方程

问题陈述
本研究探讨了在不稳定大气条件下，二维稳态平流扩散方程（ADE）的解析解。主要目标是模拟污染物的扩散，特别侧重于中心线上的横风向积分归一化浓度。该模型在地面层引入了沉降速度边界条件，并利用幂函数分别表示垂直涡扩散系数和风速随垂直高度的变化关系。本研究旨在通过将该解析方法与丹麦哥本哈根野外实验中观测到的六氟化硫（SF6）浓度数据进行对比验证。

方法论
作者采用分离变量法来求解控制方程中的 ADE。其数学框架定义如下：

• 控制方程： 二维 ADE 的构建中，风速 $u(z)$ 和垂直涡扩散系数 $k_z(z)$ 被表示为垂直高度的幂函数（分别为 $z^p$ 和 $z^q$）。
• 边界条件：
  1. 顶层边界（混合层高度 $h$）： 零通量条件（$k_z \frac{\partial C}{\partial z} = 0$）。
  2. 地面边界（$z=0$）： 沉降通量条件，即通量等于沉降速度（$V_d$）与浓度之积（$k_z \frac{\partial C}{\partial z} = -V_d C$）。
  3. 源项： 在烟囱高度 $h_s$ 处设有的连续点源，由狄拉克 $\delta$ 函数表示。
• 解析推导：
  • 方程被分解为空间分量 $X(x)$ 和 $Z(z)$。
  • 垂直分量导致了一个具有变系数的 Sturm-Liouville 型微分方程。
  • 通过假设特定的指数值（例如 $q=0.85$）并利用变量代换，该方程被简化为广义贝塞尔方程。
  • 垂直浓度分布的解利用第一类和第二类贝塞尔函数导出，对于半阶贝塞尔函数，这些函数被简化为三角形式（正弦和余弦）。
  • 通过应用边界条件确定特征值和特征函数，从而得到浓度的级数解 $C_y(x, z)$。
  • 系数通过基于风速分布的正交性条件进行计算。

主要贡献

• 解析框架： 本文提出了不稳定条件下二维 ADE 的闭合形式解析解，明确纳入了沉降速度边界条件，而这通常是扩散建模中的一个复杂因素。
• 幂律参数化： 模型利用幂律表示风速和涡扩散系数，从而能够灵活应对不同的稳定度等级（由指数 $p$ 和 $q$ 定义）。
• 验证： 该解析模型针对著名的哥本哈根 SF6 数据集进行了严格测试，涵盖了多种稳定度等级（A、B、C、D）及下风向距离。

结果
研究将预测的横风向积分归一化浓度与哥本哈根实验中观测到的 SF6 浓度进行了对比。

• 数据对比： 表 2 和图 1 展示了在不同下风向距离下，多个运行组（Run 1 至 Run 9）观测浓度与预测浓度的直接对比。
• 统计性能： 模型性能采用标准统计指标（Hanna, 1989）进行评估：
  • 归一化均方误差 (NMSE)： 0.07（表明误差较低）。
  • 分数偏差 (FB)： 0.003（表明偏差可以忽略不计）。
  • 相关系数 (COR)： 0.82。
  • 两倍因子内比例 (FAC2)： 0.997。
• 视觉一致性： 图 1 说明了观测浓度剖面与预测浓度剖面之间高度的一致性。图 2 证实了预测值落在观测数据的两倍因子范围内。

意义与主张
作者声称所提出的解析模型实现了与观测野外数据的高度闭合一致。具体而言，论文指出，由于 FAC2 值为 0.997 且偏差接近于零，预测的中心线浓度实现了与观测浓度 100% 的吻合。研究结论认为，结合风速和扩散系数的幂律参数化，分离变量法提供了一种强大的解析工具，用于预测不稳定大气条件下的污染物扩散，而无需进行复杂的数值模拟。这项工作强化了基于变换和基于分离的解析解在捕捉不稳定边界层中污染物传输基本物理过程方面的效用。