BSMPKU 与 OpenFOAM 应用在北京某小区的模拟比较研究.docVIP

BSMPKU 与 OpenFOAM 应用在北京某小区的模拟比较研究 1 引 言 随着城市化进程的推进，城市面积不断扩张，城市建筑区向着高层且密集的方向发展；同时城市人口的不断增加使得人类活动产生的各种有毒气体以及气溶胶显著改变了大气的组成成分，影响着空气质量，造成严重的城市环境问题。另一方面，全球恐怖事件频发，各种生化和放射性试剂在人口密集区释放所带来的后果令人担忧。研究表明，当高大的建筑物分布在狭长的街道两侧时，会形成峡谷效应，使得污染物很难扩散出去，并且建筑区内建筑物高度的增加会显著降低街区内的扩散能力。因此，了解污染物在实际城区的扩散特征具有十分重要的现实意义，并对城市规划、应急预警等有一定参考价值。 研究城市街区尺度的流场及扩散问题主要有以下几种方法：外场实验、风洞实验、计算流体力学（Computational Fluid Dynamics, CFD）模式模拟。 外场实验始于 20 世纪 70 年代，通过这种方法可以得到真实大气条件下，城市街区内流场以及扩散情况。DePaul 等测量了不同宽高比街谷中的风场，发现当风速超过一定值时，街谷内的涡旋才会形成。Stein 等监测了阿根廷某城市部分街区 CO 的排放情况，并与模式模拟的结果进行了对比。Macdonald 等观测了污染物在建筑物阵列迎风面释放时，发现污染物的浓度分布是符合高斯扩散模式的。 相较于外场实验，风洞实验的实验条件易于控制且成本较为低廉，自 20 世纪 80 年代后期，开始利用风洞实验对大气及污染物的扩散进行模拟研究。国内外学者通过风洞实验发现，街谷内污染物的扩散情况取决于街谷的几何结构以及排放源的位置。一些研究比较了风洞与外场实验，指出两者的差异主要来自外场实验中真实大气的风向变化以及大尺度湍流的影响。 2 模式与数据介绍 2.1 BSMPKU BSMPKU 分为两个子模式：第一个子模式是街区尺度的风场预报模式，该模式是基于稳态 RANS 方程组耦合标准 k minus;epsilon; 湍流模型的 CFD 模式，该模式在给定适当的边界条件下，用于计算城区的风场与湍流场；第二个子模式是一个欧拉扩散模式，在给定源强数据且 CFD 计算获得风场和湍流场的基础上，通过求解非稳态的雷诺平均输运方程获得浓度场。 RANS 方程组是对 Navier-Stokes 方程做雷诺时间平均分解获得的。一般，大气运动的雷诺数很大，粘性项可以忽略。同时考虑污染物浓度水平较低，不会显著影响空气密度的情况，将流体流动状态假设为不受浓度分布影响。 2.2 OpenFOAM OpenFOAM（Open Source Field Operation and Manipulation）是在 linux 平台下基于 C++的面向对象计算流体力学软件包，它自带一系列求解器可以对连续介质力学问题进行数值计算。关于 OpenFOAM 的详细信息可以参考以下网站：。 本研究中使用的 OpenFOAM 版本为 2.1.1，求解器为 simpleFoam，simpleFoam 是OpenFOAM 中稳态不可压缩流动（层流或者 RAS 湍流）求解器，该求解器求解的基本方程与 BSMPKU 模式中的一致，OpenFOAM 使用有限体积法来对控制方程组进行离散，压力速度耦合和 BSMPKU 一样采用 simple 算法。 2.3 风洞实验数据 上个世纪 90 年代，R.S.Thompson 汇集了一系列在美国环境保护署风洞中进行的长方体建筑物扩散实验的数据。这个数据集系统地涵盖了不同建筑物形状，不同排放源高度和位置的风洞实验。2006 年，丹麦国家环境保护研究所的一个研究组利用 Thompson 数据集评估了一些扩散模式在存在建筑物情况下的模拟能力。这些模式包括 OML（OperationelleMeteorologiske Luftkvalitetsmodeller ）、 AERMOD/PRIME （ American MeteorologicalSociety/Environmental Protection Agency Regulatory Model with Plume Rise ModelEnhancement）和 MISKAM（microscale climatic and dispersion model）模式。 接下来的研究，将取自 2006 年丹麦科学家所用的 Thompson 数据和基于高斯扩散理论的 AERMOD 模拟数据，与 BSMPKU 以及 OpenFOAM 的模拟结果进行比较，验