Update: June 1, 2017

$FOAM_TUTORIALS/incompressible/simpleFoam/windAroundBuildings

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

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

Meshes

Mesh (XY-plane at Z=1.5m)

The calculation result is as follows.

Flow velocity in XY-plane at Z=1.5m (U)

Pressure in XY-plane at Z=1.5m (p)

Turbulence energy in XY-plane at Z=1.5m (k)

Turbulence dissipation rate in XY-plane at Z=1.5m (epsilon)

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

cd windAroundBuildings

surfaceFeatureExtract

blockMesh

snappyHexMesh -overwrite

simpleFoam

paraFoam

cd windAroundBuildings

surfaceFeatureExtract

blockMesh

snappyHexMesh -overwrite

simpleFoam

paraFoam

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