Atmospheric Wind Field Modelling with OpenFOAM for Near-Ground Gas Dispersion

CFD simulations of near-ground gas dispersion depend significantly on the accuracy of the wind field. When simulating wind fields with conventional RANS turbulence models, the velocity and turbulence profiles specified as inlet boundary conditions change rapidly in the approach flow region. As a result, when hazardous materials are released, the extent of hazardous areas is calculated based on an approach flow that differs significantly from the boundary conditions defined. To solve this problem, a turbulence model with consistent boundary conditions was developed to ensure a horizontally homogeneous approach flow. Instead of the logarithmic vertical velocity profile, a power law is used to overcome the problem that with the logarithmic profile, negative velocities would be calculated for heights within the roughness length. With this, the problem that the distance of the wall-adjacent cell midpoint has to be higher than the roughness length is solved, so that a high grid resolution can be ensured even in the near-ground region which is required to simulate gas dispersion. The evaluation of the developed CFD model using the German guideline VDI 3783/9 and wind tunnel experiments with realistic obstacle configurations showed a good agreement between the calculated and the measured values and the ability to achieve a horizontally homogenous approach flow.