ansys气体扩散仿真污染扩散

一、项目背景与科学问题

工业烟囱与事故泄漏排放的有害气体在大气中的扩散直接关系到厂界浓度达标与环境风险防控。预测下风向污染物浓度的空间分布，是环境影响评价、应急预案与排放源强反演的基础。高斯烟羽模型是应用最广泛的稳态扩散模型，其将湍流扩散近似为高斯分布，以解析形式给出浓度场。本项目的科学问题是：连续点源（烟囱）排放的污染物在下风向中心线上的浓度如何随距离衰减？垂直于风向的浓度横向分布具有怎样的形态？基于ANSYS Fluent的组分输运仿真可以数值求解对流扩散方程，获得比解析模型更精细的浓度场。

二、计算方法与理论背景

高斯烟羽模型认为污染物在水平与垂直方向均按高斯分布扩散，地面浓度表达式为C=Q/(2·pi·u·sigma_y·sigma_z)·exp(-y²/(2·sigma_y²))·[exp(-(z-H)²/(2·sigma_z²))+exp(-(z+H)²/(2·sigma_z²))]，其中Q为源强、u为平均风速、sigma_y与sigma_z为横向与垂直扩散参数、H为有效源高。扩散参数随下风向距离按幂律增长，反映湍流扩散随时间（即随风向距离）增强。中心线浓度在下风向先受源高抬升影响而较低，随sigma_z增大地面浓度逐渐上升至峰值，随后随距离增大而按近似1/x²衰减。CFD方法通过求解Reynolds平均Navier-Stokes方程与组分输运方程，可考虑地形与障碍物影响，是对高斯模型的重要补充。

三、计算设置与模型构建

以10克每秒的源强、3米每秒的平均风速为基准工况，烟囱有效源高取10米，大气稳定度取中性（D类）条件，按帕斯奎尔公式确定sigma_y与sigma_z的扩散参数。高斯模型直接按解析式编程计算下风向20米至1000米范围内的地面中心线浓度与横向剖面。CFD部分在ANSYS Fluent中建立包含烟囱入口的计算域，采用标准k-epsilon湍流模型求解流场，再以组分输运方程求解污染物浓度场，入流边界给定风速廓线，出口采用压力出口，地面设为无滑移壁面。将CFD结果与高斯模型对比，评估两种方法的一致性。

四、结果分析与案例图解读

图13的左侧面板（a）给出地面中心线浓度随下风向距离的双对数曲线。浓度在下风向约20米处较低，随后上升并在约80米处达到峰值，此后随距离增大而稳定衰减，从数百米到一千米范围近似呈线性下降趋势，对应高斯模型中远场浓度与距离平方近似成反比的规律。峰值出现的原因在于烟羽中心线尚未触及地面时，地面浓度随扩散增强而增大；当烟羽完全落地后，稀释效应占主导使浓度持续下降。这一”先升后降”的形态是抬升源高斯扩散的典型特征，也是确定最大落地浓度及其出现距离的依据。

图13的右侧面板（b）给出三个下风向距离处（100米、300米、600米）的横向浓度剖面。各剖面均呈以烟羽轴线为中心的高斯钟形分布：距离越近峰值越高、剖面越窄；距离越远峰值越低、剖面越宽。100米处剖面尖锐，600米处剖面扁平且横向铺展明显，直观反映了湍流扩散随距离对烟羽的横向稀释作用。横向高斯形态说明中性层结下污染物在横风方向的扩散接近正态分布，可用sigma_y统一刻画。三组剖面的对比定量展示了”峰值降低、宽度增大”的扩散演化规律，可用于确定厂界敏感点的浓度贡献。

两面板结合，左面板给出了浓度沿风向的衰减规律，右面板给出了浓度横风方向的扩散形态，共同完成了ANSYS气体扩散仿真的完整分析，为环境风险评价提供了方法与数据支撑。

五、结论与展望

本项目基于高斯烟羽模型并结合CFD方法分析了连续点源气体扩散，中心线浓度先升后降、远场近似1/x²衰减，横向剖面呈高斯分布并随距离展宽。案例图以中心线浓度曲线与横向剖面双面板呈现。后续可拓展：(1)引入非中性层结与不同风速廓线，分析扩散参数的敏感性；(2)考虑建筑物绕流与复杂地形，用CFD精细刻画局地扩散；(3)耦合非稳态泄漏情景，开展事故扩散的瞬态风险评估。