TPU Building - Corner Cut#

TPU Building - Corner Cut

1. Project Description#

The current study numerically reproduces the experiment CornerCut_Cat2 from the Tokyo Polytechnic University (TPU) Aerodynamic Database with AeroSim’s CFD solver. The experimental results of this case, along with many other complex high-rise building geometries, are discussed in Tanaka et al.[1].

A high-rise building with \(B \times H\) dimensions:

TPU Building - Corner Cut geometry

is positioned in a wind tunnel setup with 184 pressure probes distributed over its surface:

TPU Building - Corner Cut coordinates

The wind directions chosen to be simulated were:

TPU Building - Corner Cut Cases#

Wind direction

0 \(^\circ\)

15 \(^\circ\)

30 \(^\circ\)

45 \(^\circ\)

2. Simulation Setup#

The Synthetic Eddy Method (SEM) boundary condition is applied at the inlet of the computational domain. Solid fins are distributed across the floor to ensure the desired velocity and turbulence profiles during flow development length. A Neumann boundary condition is applied at the remaining boundaries.

The building is positioned \(48B\) from inlet, and 6 grid refinement levels (\(lvl\,0\) to \(lvl\,5\)) were adopted:

TPU Building - Corner Cut grid

A 1:2 refinement ratio is estabilished between levels, and the simulation parameters at grid \(lvl\,5\) were:

TPU Building Corner Cut - Parameters#

\(\boldsymbol{\Delta x/B}\) (spatial resolution)

1.00e-02

\(\boldsymbol{\Delta t/CTS}\) (temporal resolution)

6.96e-04

exports/\(\boldsymbol{CTS}\) (pressure acquisition frequency)

1.00e+01

\(\boldsymbol{T/CTS}\) (statistical sample size)

8.40e+02

\(\boldsymbol{Re_{B}=U_{H}B/\nu}\)

4.15e+04

The equivalent parameters in full scale are:

TPU Building Corner Cut - Full Scale#

\(\boldsymbol{\Delta x[m]}\)

0.50

\(\boldsymbol{\Delta t[ms]}\)

0.50

\(\boldsymbol{f[Hz]}\)

14.00

\(\boldsymbol{T[s]}\)

600.00

The computational resources required were:

TPU Building Corner Cut - Resources#

Device

NVIDIA RTX A5500

NVIDIA A10G

NVIDIA A10G

NVIDIA RTX A5500

Wind direction

0 \(^\circ\)

15 \(^\circ\)

30 \(^\circ\)

45 \(^\circ\)

Node count (million)

130.73

133.70

134.52

132.49

Allocated memory (Gb)

16.89

17.29

17.40

17.15

Ellapsed time (h)

36.89

48.71

50.18

39.75

3. Inflow#

An empty domain simulation is performed to measure the incident velocity and turbulence profiles. A probe line is placed at the position where the building will be located. The average velocities used for calculating the pressure coefficient and convective time scale are taken from this simulation.

Wind Profiles#

TPU Building - Corner Cut

Wind Spectra#

The power spectral density of the velocity components at height \(H\) are compared with theoretical Von Kármán curves to validate the atmospheric flow.

TPU Building - Corner Cut

4. Results: Local Statistics#

The pressure coefficient is calculated using the velocity at the building height \(H\) and the reference pressure measured from a position far above the building. The same statistical processing of the pressure coefficient \(C_{p}\) at each probe position is performed for both experimental and numerical data:

  • A sample of size (\(T / CTS\)) is evaluated over moving averages of the original signals, with window size of 3s.

  • The processed data is subdivided in 10 intervals, from which the minimum and maximum values will be taken for each.

  • Those min/max values are fitted to a Gumbel distribution.

  • The shape parameters from Gubel distribution are then rescaled to 3600s.

Scatter on Points Statistics#

The deviation between numerical and experimental data is quantified using the mean absolute error (MAE) and the normalized mean absolute error (NMAE):

\[\mathrm{MAE} = \frac{1}{N_{\mathrm{probes}}}\sum_{i=1}^{N_{\mathrm{probes}}}|q^{\left(i\right)}_{\mathrm{EXP}}-q^{\left(i\right)}_{\mathrm{NUM}}|\]
\[\mathrm{NMAE} = \frac{1}{N_{\mathrm{probes}}}\sum_{i=1}^{N_{\mathrm{probes}}}\frac{|q^{\left(i\right)}_{\mathrm{EXP}}-q^{\left(i\right)}_{\mathrm{NUM}}|}{\left[q^{\left(\mathrm{max}\right)}_{\mathrm{EXP}}-q^{\left(\mathrm{min}\right)}_{\mathrm{EXP}}\right]}\times 100\]
TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut

Mean and Peak Pressures#

TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut

RMS Pressures#

TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut

Skewness and Kurtosis#

TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut

Pressure Spectrum#

TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut
TPU Building - Corner Cut

Execution notes#

TPU Building Corner Cut - Execution notes#

Execution Date (YYYY-MM-DD)

2024-10-26

Solver Version

1.6.0a2

Changelog#

  • 30 Oct 2024: Added scattering plots