大气中传输扩散计算中的数值方法（Ⅱ）：模式研究.pdf

第五届全国风工程及工业空气动力学学术会论文集 大气中传输扩散计算中的数值方法(II):模式 桑建国 林官明 张伯寅 地球物理系 环境科学中心 湍流重点实验室 北京大学，北京，100871 王景书 中国辐射防护研究院，山西太原 030006 引 言 一般在空气污染气象学的观念中，只要地形是非平坦、非均匀的，都可称为是复杂 地形，例如山区、城市、沿海等。从这个意义上，在大气环境研究中所涉及的地形，几乎 都可归属于复杂地形。在复杂地形上由于流场和rv流场空间上的不均匀，时间上不定常， 使得用解析方法得到的评价模式 (例如正态模式)计算出的浓度分布和实际状态之间，有 较大的偏差。地形愈复杂，这类偏差愈大。我国是个多山国家，又有较长的海岸线，发展 工业、能源和交通都必然要涉及到复杂地形上的大气环境问题。因此解决复杂地形上的扩 散问题，也是我国空气污染气象学工作者面临的重要任务。发展复杂地形上的扩散模式， 关键是如何准确地计算出由于地形的动力和动力作用形成流场.目前计算流场的模式可分 为:1)诊断模式，即用实测及资料内插、调整，得到计算格点上的流场.z)动力模式， 即用大气动力一热力方程组数值解得出流场.例如Pielke的Rams模III和Yamada的 HOTPIC模式lzl 本文介绍的 “北京大学大气环境模式”是一动力流场一浓度场模式，它吸收了同类先 进模式的长处。可在微机上运行。经过多年使用，已证明可用于复杂地形所形成的各类大 气环流的模拟.以及在这些环流中的浓度场的计算。其中包括地形热力强迫形成的环流， 如山谷风、海陆风、海风热力内边界层、冷泄流 (下坡风)、城市热岛及内边界层;地形 动力作用引起的气流阻塞、流场扰动、背风涡旋、背风波和地形绕流等等。以及上述环流 中的污染物浓度分布。参见文献-311’!. 所模拟的环流尺度，由水平宽度只有几百米的山丘的扰动流场，到水平范围为千公 里的高原形成的大型山谷风。尺度覆盖为10勺Osm，包括气象学中的微尺度和中尺度 (micor-。、mcso-Y和meso-0), 数值模式除了较真实地反映了复杂地形上流场浓度场状态外，它的采用还有助于对 复杂地形上的大气现象的深入研究.对大气探测提出新的要求。促进空气污染气象学和中 .368. 尺度气象学的进一步发展. 模式的基本结构 北京大学大气环境模式简称PUMA(PekingUniversityModelofAtmospheicenvironment) 主要构成部分如图1所示， 其中准静力动力学流场式构成模式的主要部分.它愉人观侧到的气象场资料 (如原 始风场和温度场)，愉出为模式各个格点上的风速和温度分布. 边界层湍流参数化模式和风场模式联立，根据湍流动能方程仃KE)计算出湍流动能， 得出边界层中动量、热量湍流通量，并反馈给流场校式. 地表能量平衡方程模拟地表温度变化的物理过程.在模拟热力强迫产生的环流中， 它反映了扰动源的状态. 在已知网格点上流场和湍流场状态的基础上，采用平流扩散方程的数值解或质点随 机行走模式，根据输人的污染源资料，计算出污染物浓度的空间分布和时间演变. 下面将介绍模式的各个部分. 1.1中尺度准静力动力学模式 在下述条件下，大气运动可假定是准豁力的，即:(1)对流活动较弱，因而垂直运动 速度较小;(2)地形坡度较小 (例如小于45“); (3)地形特征尺度 例《如山谷的水平宽 度)不小于几公里.合乎上述条件的大气运动可采用准静力模式模拟. 为了反映下边界地表的起伏，模式采用了地形追随的垂直坐标变换， :=H-- (]) H一z# 其中Z，二ZF(X..Y)为地形高度，z为笛卡尔坐标系((c,Y.z)中的垂直坐标，Z为变换后的垂 直坐标，H为模式的顶部高度. 经过变换后的三维大气动力一热力学方程组包括水平动量方程、准静力方程、连续 方程