Wind environment around buildings
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
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)
Commands
cp -r $FOAM_TUTORIALS/incompressible/simpleFoam/windAroundBuildings windAroundBuildings
cd windAroundBuildings
surfaceFeatureExtract
blockMesh
snappyHexMesh -overwrite
simpleFoam
paraFoam
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