1.3 边界条件
计算域入口采用梯度风速度入口,来流风与街谷垂直,出口边界设为自由出流;计算域上部边界及侧面设为对称边界条件,即为零梯度边界;在建筑物的近壁面区和近地面区,雷诺数较低,湍流应力小,导致湍流不能充分发展,因此在建筑物的壁面、地面,引入了标准壁面函数进行近壁面修正。所有壁面应用无滑移墙边界条件。
1.4 计算方程和收敛情况
本文采用ANSYS FLUENT 19.2 进行计算。来流风通过用户自定义函数(UDF)方法进行编译;对流项的离散方法选用二阶迎风方案;梯度项插值选用格林-高斯节点方案;实验计算法则采用SIMPLE算法;压力插值方法使用Standard算法;采用重整化k-ε模型;计算方程为连续性方程,动量方程,能量方程,组分输运方程以及RNG k-ε方程。能量方程的收敛因子设置为 1×10-8,其它变量收敛因子设置为 1×10-5。