Wind environment around buildings

Update: June 1, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/incompressible/simpleFoam/windAroundBuildings

Summary

We will analyze the wind flow around buildings. The fluid flows in at (10, 0, 0) m/s from the region "inlet" (X minimum plane) and out from the region "outlet" (X maximum plane). No-slip conditions are set for the ground and building surfaces. The other walls (region "symmetry") are set to symmetry by including "caseDicts/setConstraintTypes ($WM_PROJECT_DIR/etc/caseDicts/setConstraintTypes)".

The standard k-ε model is used as the turbulence model.

Model geometry Model geometry

The meshes are as follows, and the number of mesh is 334666.

Mesh Meshes
Mesh Mesh (XY-plane at Z=1.5m)

The calculation result is as follows.

Flow velocity in XY-plane at Z=1.5m (U) Flow velocity in XY-plane at Z=1.5m (U)
Pressure in XY-plane at Z=1.5m (p) Pressure in XY-plane at Z=1.5m (p)
Turbulence energy in XY-plane at Z=1.5m (k) Turbulence energy in XY-plane at Z=1.5m (k)
Turbulence dissipation rate in XY-plane at Z=1.5m (epsilon) Turbulence dissipation rate in XY-plane at Z=1.5m (epsilon)

Commands

cp -r $FOAM_TUTORIALS/incompressible/simpleFoam/windAroundBuildings windAroundBuildings
cd windAroundBuildings

surfaceFeatureExtract
blockMesh
snappyHexMesh -overwrite

simpleFoam

paraFoam

Calculation time

6 minutes 18.35 seconds *Single, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz

Reference