模型建筑群间风场和污染物浓度场的大涡数值模拟.docVIP

模型建筑群间风场和污染物浓度场的大涡数值模拟 史瑞丰1, 崔桂香1*, 许春晓1, 张兆顺1, 王志石2 1工程力学系,清华大学,北京,100084,中国 2技术学院, 澳门大学,澳门,中国 摘要: 本文采用大涡数值模拟方法研究了建筑群间的风场和污染物浓度场. 计算中采用具有四阶精度的有限体积方法, 并用浸没边界方法实现固壁边界条件. 以一种规则排布的建筑群为例, 数值模拟得到的结果与实验数据相符较好. 关键词: 大涡模拟, 浸没边界方法, 有限体积方法, 建筑群, 污染物 引言 空气质量问题是现代城市可持续发展中需要关注的重要课题. 本文采用大涡数值模拟(LES)来预测居民小区内的风场和污染物浓度场. 居民区内的空气流动具有几何边界复杂,湍涡类型多变,流动非定常等特性, 所以适合采用LES来模拟这类流动问题. 本文采用具有四阶精度的有限体积方法, 并用浸没边界法(IBM)来处理固壁边界条件. 该数值方法已通过一系列复杂流动算例的验证(Xu et al. 2006). 在污染物扩散的大涡数值模拟中, 采用拉格朗日动力模式计算亚格子应力, 并假定湍流普朗特数为常数采用湍涡扩散模型处理浓度控制方程. 以一种规则排布的建筑群为例, 数值模拟得到的结果与实验数据相符较好. 结果表明本文所用LES方法能够预测居民小区内的风场和污染物浓度分布. 数值方法 大涡模拟的控制方程可以写成: (1) (2) (3) 其中?ij和Tj 是亚格子应力和污染物浓度的亚格子通量, 分别采用涡粘模式和湍涡扩散模式处理 (4) , (5) 其中?t, ?t分别是亚格子涡粘系数和湍涡扩散系数, 湍流普朗特数Prt = 0.70. 对于Smagorinsky模式, 比较了Cs = 0.1 ()时单纯Smagorinsky模式和动力模式的结果. 在Meveneau (1996)提出的拉格朗日动力模式中模式系数C = ILM/IMM (6) (7) 其中. 在有限体积离散格式的推导中采用了四阶Padé型紧致插值, 具体过程参见Xu et al. (2006). 时间推进采用四阶龙格-库塔法. 正交网格上的壁面无滑移条件采用Peskin (1972)提出的浸没边界方法来实现. 应用浸没边界方法时在边界附近网格点的控制方程中引入一个附加的力项 (8) 其中附加力项由下式给出 (9) RHS包含了压力梯度项, 对流项和扩散项, 是固壁边界的速度. 结果分析 THE RESULTS 将一组交叉排布的立方体块作为模型建筑群, 如图1所示, 每个立方体的长H, 宽B和高W均为0.12米. 基于来流风速和建筑物高度的雷诺数为4(106. 计算区域和采用的非均匀网格如图2所示. 在建筑群前方10W远处布置了一个污染物排放点源. 图1 建筑群排布 图2 计算区域和网格 边界条件的设定如下: 进口处平均风速给定按对数分布的边界层风廓线, 并添加按高斯分布的随机扰动; 出口采用无反射条件, 展向采用充分发展条件, 分别如下 , at x/B=70 (10) , at y/B=0 and y/B=30