重气扩散中复杂障碍物的数值处理方法.pdfVIP

第 30 卷 第 10 期 计算机与应用化学 Vol.30, No.10 2013 年 10 月 28 日 Computers and Applied Chemistry October 28, 2013 重气扩散中复杂障碍物的数值处理方法 1 2* 2 1 王正 ，李向阳 ，陈杰 ，赵祥迪 (1. 中国石化安全工程研究院，山东，青岛，266071 ； 2. 中国科学院过程工程研究所，北京，100190) 摘要：CFD 方法逐渐成为国内外研究重气扩散物理过程的重要手段。在泄漏源或扩散过程中存在的不同类型的障碍物改变了大 气的流动，使扩散过程更为复杂，在模拟过程中对障碍物的处理是获得准确模拟的关键问题之一。本文综述了目前文献中处理 复杂障碍物的方法的基本原理，并通过指出各自的优缺点，为自主开发程序进行重气扩散过程 CFD 模拟的研究者提供参考。商 业软件普遍采用的贴体坐标法, 该方法实质上是通过特定的变换方法把物理平面上的不规则区域转换成计算平面上的规则区 域，原则上是与实际障碍物最接近的处理方法。该方法的缺点在于生成网格过程工作量较大，也会使微分方程和数值方法复杂 化，同时自编程序实现难度较大；浸入边界法和镜像流体法类似，都在固定的欧拉坐标系中求解流体运动方程。不同之处在于 浸入边界法是把固体界面对流体的作用看作流体在界面所在位置受到的弹性力作用，在流体运动方程中加入一个外力项来描述 界面作用在流体上的力。镜像流体法则把固体所占用区域格点上的流体力学参数由镜像关系指定，固体的运动通过拉格朗日方 法追踪。这两种方法生成网格过程工作量小，处理流场中复杂固体边界均能获得准确的模拟结果，具有良好的应用前景。近似 台阶法原理上最简单也最易实现，但其应用需要进一步改进以提高其精度。 关键词：重气扩散；数值模拟；障碍物 中图分类号：X169; X937 文献标识码：A 文章编号：1001-4160(2013)10-1087-1090 DOI: 10.11719/com.app.che 法。 2.1 贴体坐标系法 1 引言 贴体坐标系(Body-Fitted Coordinates，BFC)是近些年 根据原劳动部的有关统计数据，在因毒物泄漏造成 来发展起来的一种网格坐标系统，通过构造与不规则计 的人员伤亡中，约有90 %与重气泄漏有关[1] 。在泄漏源 算域相适应的坐标系，使网格边界与物理模型边界相吻 [7]和图 或扩散过程中往往存在不同类型的障碍物(如山坡、建 合，能适用于具有复杂几何形状的积分区域( 图 1 筑、树木等) ，改变了大气的流动，使得有障碍物情况下 2) 。贴体坐标系的生成过程可以看成是一种变换，即把 [2] 物理平面上的不规则区域变换成计算平面上的规则区 的重气扩散过程更为复杂 。 [12] [13] CFD 方法能够较为准确地描述重气扩散的物理过 域。通常有 3 种方式：求解微分方程