CFD

大气边界层的数值方法实现:以OpenFOAM为例(一)

前言

基本概念

地球表面通过地面的摩擦对空气水平运动产生阻力,从而使气流速度减慢,该阻力对气流的作用随高度的增加而减弱。当超过特定高度 (zG)(z_G) 时,这一作用即可被忽视,气流将沿等压线以梯度风速流动,则将此高度以下的区域称作大气边界层Atmosphere Boundary Layer, ABL)。一般的,大气边界层内近地层的气流是湍流,湍流掺混使地表阻力的影响扩展到大气边界层的整个区域。而在边界层高度以上的自由大气中的风流动属于层流。(1)

风速剖面

对数律

平均风剖面是微气象学研究风速变化的一种主要方法。目前,气象学家认为用对数律表示大气底层强风风速廓线比较理想,其表达式为(2,3)

u = \frac{u^*}{\kappa} \ln \left( \frac{z - d + z_0}{z_0} \right)\tag{1}

式中

uu — 高度 zz 处的平均风速(m/s)

uu^* — 剪切风速(m/s)

κ\kappa — 卡曼常数,一般取 0.4

dd — 零平均位移(m)

z0z_0 — 粗糙高度(m),取值可见下表(1)

地面类型z0/mz_0/\mathrm{m}地面类型z0/mz_0/\mathrm{m}
砂地0.0001~0.001矮棕榈0.10~0.30
雪地0.001~0.006松树林0.90~1.00
割过的草地(~0.01m)0.001~0.01稀疏建成市郊0.20~0.40
矮草地、空旷草原0.01~0.04密集建成市郊、市区0.80~1.20
休耕地0.02~0.03大城市中心2.00~3.00
高草地0.04~0.10

城市地区的修正下次再写。

指数律

懒得写了。

湍流特征

参考文献

1.
Huang B. Principles and Application of Wind-Resistant Analysis of Structures. Tongji University Press; 2001.
2.
Yang Y, Gu M, Chen S, Jin X. New inflow boundary conditions for modelling the neutral equilibrium atmospheric boundary layer in computational wind engineering. Journal of Wind Engineering and Industrial Aerodynamics. 2009;97(2):88–95.
3.
Richards PJ, Hoxey RP. Appropriate boundary conditions for computational wind engineering models using the k-epsilon turbulence model. In: Computational Wind Engineering 1. Elsevier; 1993. p. 145–53.