ap Size (42m, 0m) (x, y)
**************************

.. html metadata
.. meta::
   :description: Gap Size (x, y) = (42, 0) m
   :keywords: Digital Wind Tunnel, Validation, Gap 42 x 0

.. Public Project
.. ==============
.. .. link to insight page
.. `AeroSim Web App`_

.. .. _AeroSim Web App: https://insight.aerosim.io/cp/TODEFINE/


.. animations, surfaces and contours
.. figure_expand:: /_static/adj/x3-y0/visuals_1.png
   :alt: Gap Size (x, y) = (42, 0) m
   :figwidth: 100%
   :align: center

1. Project Description
======================

The current study numerically reproduces the experiment `P114_A114_case19`_ from the Tokyo Polytechnic University (TPU) Aerodynamic Database with AeroSim's CFD solver.
The experimental results of this case are analyzed and discussed in :footcite:t:`hui2012mutual` and :footcite:t:`hui2013pressure`.

.. _P114_A114_case19: https://db.wind.arch.t-kougei.ac.jp/aerodynamic/case/P114_A114_case19/
.. Cite experimental database or paper

.. footbibliography::

A high-rise building with :math:`B \times H` dimensions:

.. Comprehensive description of the geometries and dimensions
.. figure_expand:: /_static/adj/geometry.svg
   :alt: Building Geometry
   :max-width: 640px
   :align: center

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

.. Probes location and wind direction orientation
.. figure_expand:: /_static/adj/probes.svg
   :alt: Probes Layout
   :max-width: 640px
   :align: center

The neighbour building of same size is placed as:

.. figure_expand:: /_static/adj/x3-y0/position.svg
   :alt: Building Position
   :max-width: 640px
   :align: center

The wind directions chosen to be simulated were:

.. list of analyzed cases
.. csv-table:: Wind Directions
   :file: /_static/adj/x3-y0/cases.csv
   :widths: 20, 16, 16, 16, 16, 16
   :header-rows: 0


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 :math:`63B` from inlet, and 6 grid refinement levels (:math:`lvl\,0` to :math:`lvl\,5`) were adopted:

.. Top and lateral views, indicating the box and body refinements
.. figure_expand:: /_static/adj/x3-y0/grid.svg
   :alt: Computational Grid
   :figwidth: 100%
   :align: center

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

.. constant simulation parameters
.. csv-table:: Dimensionless Parameters
   :file: /_static/adj/x3-y0/params.csv
   :widths: 80, 20
   :header-rows: 0

The equivalent parameters in full scale are:

.. full scale parameters
.. csv-table:: Full-scale Parameters
   :file: /_static/adj/x3-y0/full_scale.csv
   :widths: 80, 20
   :header-rows: 0

The computational resources required were:

.. hardware and execution aspects
.. csv-table:: Computational Resources
   :file: /_static/adj/x3-y0/perf.csv
   :widths: 30, 14, 14, 14, 14, 14
   :header-rows: 0

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
-------------

.. profiles of velocity, turbulence and integral lenght scale
.. figure_expand:: /_static/adj/inflow/inflow_profile.svg
      :alt: Inflow Profile
      :max-width: 640px
      :align: center


Wind Spectra
------------
The power spectral density of the velocity components at height :math:`H` are compared with theoretical Von Kármán curves to validate the atmospheric flow.

.. spectra of velocity components
.. figure_expand:: /_static/adj/inflow/spectrum.svg
      :alt: Inflow Spectrum
      :figwidth: 100%
      :align: center

.. change spectrum legend to top left, over bottom left


4. Results: Local Statistics
============================

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

  - A sample of size (:math:`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):

.. math:: 
  \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}}|

.. math::
  \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

.. figure_expand:: /_static/adj/x3-y0/scatter_000.svg
  :alt: Scatter 000
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/scatter_045.svg
  :alt: Scatter 045
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/scatter_090.svg
  :alt: Scatter 090
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/scatter_135.svg
  :alt: Scatter 135
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/scatter_180.svg
  :alt: Scatter 180
  :figwidth: 100%
  :align: center
  

Mean and Peak Pressures
-----------------------

.. figure_expand:: /_static/adj/x3-y0/mean_peak_000.svg
  :alt: Mean and Peaks 000
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/mean_peak_045.svg
  :alt: Mean and Peaks 045
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/mean_peak_090.svg
  :alt: Mean and Peaks 090
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/mean_peak_135.svg
  :alt: Mean and Peaks 135
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/mean_peak_180.svg
  :alt: Mean and Peaks 180
  :figwidth: 100%
  :align: center


RMS Pressures
-------------

.. figure_expand:: /_static/adj/x3-y0/rms_000.svg
  :alt: RMS 000
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/rms_045.svg
  :alt: RMS 045
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/rms_090.svg
  :alt: RMS 090
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/rms_135.svg
  :alt: RMS 135
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/rms_180.svg
  :alt: RMS 180
  :figwidth: 100%
  :align: center


Skewness and Kurtosis
---------------------

.. figure_expand:: /_static/adj/x3-y0/highstats_000.svg
  :alt: Skew and Kurt 000
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/highstats_045.svg
  :alt: Skew and Kurt 045
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/highstats_090.svg
  :alt: Skew and Kurt 090
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/highstats_135.svg
  :alt: Skew and Kurt 135
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/highstats_180.svg
  :alt: Skew and Kurt 135
  :figwidth: 100%
  :align: center


Pressure Spectrum
-----------------

.. figure_expand:: /_static/adj/x3-y0/spectra_000.svg
  :alt: Cp Spectrum 000
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/spectra_045.svg
  :alt: Cp Spectrum 045
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/spectra_090.svg
  :alt: Cp Spectrum 090
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/spectra_135.svg
  :alt: Cp Spectrum 135
  :figwidth: 100%
  :align: center

.. figure_expand:: /_static/adj/x3-y0/spectra_180.svg
  :alt: Cp Spectrum 135
  :figwidth: 100%
  :align: center

Execution Notes
===============

.. csv-table:: Execution Notes
   :file: /_static/adj/x3-y0/exec.csv
   :widths: 50, 50
   :header-rows: 0

Changelog
=========

* **30 Oct 2024**: Added scattering plots