城市小区空气流动与交通污染大涡数值模拟.pdfVIP

自避．科誊遗展 第19卷 第2期 2009年2月 城市小区空气流动与交通污染 的大涡数值模拟* 史瑞丰 崔桂香一 王志石 许春晓 张兆顺 1．清华大学工程力学系 ，北京 i00084；2．澳 门大学科技学院，澳门 摘要 城市小区内的大气环境属于微环境流动．复杂的建筑物下垫面决定 了城市小 区大气环境的 精细数值模拟的实现难度 ．文 中采用大涡模拟方法，针对微环境流动的特点研究边界条件的处理方 法和亚格子模 式的合理选用，实现 了小区徽环境风场和污染物浓度场的并行计算．将澳门荷兰园小 区进行几何建模，针对该小区的地形特征，采用浸没边界法和分布 的阻力元法相结合，实现 了复 杂建筑物下垫面的有效处理，计算 了荷兰园小区从早 7：00至晚 17：OO的风场和交通污染物浓度 的 分布 ，取得了与风洞实验和监测数据一致的结果．通过 比较不同位置处的污染物浓度，分析 了建筑 物几何结构对污染物扩散的影响；通过考察不同街道 内污染物浓度的分布情况，得到了城市小区 内污染物扩散的细节特征． 关键词 微环境流动 大涡模拟 复杂下垫面 污染物扩散 浸没边界 目前大气环境问题引起 了全世界范围内的，泛 出建筑物几何结构 的数值方法．由于城市小区包含 关注 ，城市化进程的加快和机动车数量的急刚增长 了若干条街道和众多形状各异 的建筑物 ，其几何边 更使得城市空气质量成为人们关注的焦点问题．目 界的复杂程度远远超过单个街 区，这决定 了实现城 前 ，城市空气质量的研究手段主要为现场监测、模 市小区大气环境 CFD模拟的复杂性．然而，城市交 型实验和数值模拟．由于数值模拟具有灵活方便、 通污染物排放和有毒有害气体泄漏等问题大都 出现 时间周期短和费用低廉等特点 ，随着计算机资源的 在城市冠层以下，其扩散范围也超出单个街区，这 迅速 改 善，计算 流 体力 学 (computationalfluid 就对城市小区微环境 问题的研究提出了迫切要求． dynamics，CFD)方法已经成为研究大气环境流动问 由于实际城市小 区的建筑物下垫面和气象条件 题经济而有效的手段． 都 比较复杂 ，采用 CFD方法研究城 市小 区大气环 城市建筑群间的空气流动受建筑物的阻碍 ，不 境时，常将小区内建筑物的几何构型和空间分布进 断伴随着旋涡的产生 、分离和破碎 ，形成复杂的湍 行简化 (如将建筑物简化为立方体)，研究模型建筑 流风场．机动车排 出的尾气在气流的卷吸携带下， 群 内的流动和污染物扩散 ，而能够分辨实际小 区内 扩散于建筑物之问，由于建筑物的扰动作用 ，可能 建筑物几何结构的 CFD模拟研究 尚处于起步阶段． 使得局部的污染物浓度改变一个量级 以上l1j．因此 Eichhorn等提出的MISKAM 数值模式 ]考虑 了建 若要准确预测城市小区内污染物浓度 的分布，不能 筑物几何结构对城市小区内空气流动的影响 ，该方 像中尺度数值模式 (如 mm5)那样将建筑物下垫面作 法采用 Reynolds统计平 均模式 (Reynoldsaverage 平均化处理来给定地表粗糙度，需要采用能够分辨 navier—stokes，RANS)和壁 函数来模拟平均风场， 2008—0305收稿 ，200806—20收修改稿 * 国家 自然科学基金 (批准号 和澳门科技发展基金 (批 准号：FDCT一022／2006／A)资助项 目 **通信作者 ，Email：cgx@mail．tsinghua．edu．CFI 自鞋科乎建展 第19卷 第2期 2009年2月 213 在湍流模式中如何选取合适的经验参数 以适于城市 方程 (1)和 (2)分别