TY - JOUR
T1 - A Cartesian grid refinement method for simulating thermally stratified urban environments
AU - Angelidis, Dionysios
AU - Assimakopoulos, Vasiliki D.
AU - Bergeles, George
N1 - Publisher Copyright:
© 2015 Elsevier Ltd.
PY - 2015/7/1
Y1 - 2015/7/1
N2 - The capability to accurately simulate thermally stratified environments highly depends on the numerical model which governs the physical problem, but also on the temporal and spatial discretization resolution. This paper presents a multi-resolution Cartesian mesh refinement method for solving the incompressible Navier-Stokes equations, while the standard k-ε turbulence model is employed to resolve stratified turbulent flows. Hanging nodes are used to facilitate the implementation of high order discretization schemes and the algebraic system of the discretized differential equations is solved with a generic modified Alternating Direction Implicit (m-ADI) method developed to resolve locally refined grids. The overall performance of the new flow solver is initially evaluated against experimental data related to the laminar flow past a square cylinder for Reynolds numbers, defined based on the mean ambient wind speed and the characteristic length scale, of 87≤Re≤250 and neutral thermal stratification. The unsteady character of the flow field is successfully predicted and the grid refinement approach enables accurate low cost and low blockage simulations. The flow solver is applied to simulate the flow over thermally stratified idealized urban areas. Good agreement with experimental measurements, where the bulk Richardson number was defined based on the mean wind speed at the top of the reference street canyon, is obtained for -0.21≤Rb≤0.78 while the grid refinement method allowed simulation of the experimental setup in detail, including all the roughness elements and the idealized street canyons. The numerical calculations indicate that even weak thermal stability can lead to laminarization of the induced flow field; reducing at the same time the potential for air mixing and results to air degradation. Even though the recirculation patterns are influenced under weak thermal stratification, extreme thermal instability within the region of -16≤Rb≤-64 with Re=6800 results to the development of two counter-rotating vortices of equivalent kinetic energy. However, the linear correlation of the average kinetic energy at the roof-level with the Rb number indicates that the buoyancy forces determine the air mixing levels.
AB - The capability to accurately simulate thermally stratified environments highly depends on the numerical model which governs the physical problem, but also on the temporal and spatial discretization resolution. This paper presents a multi-resolution Cartesian mesh refinement method for solving the incompressible Navier-Stokes equations, while the standard k-ε turbulence model is employed to resolve stratified turbulent flows. Hanging nodes are used to facilitate the implementation of high order discretization schemes and the algebraic system of the discretized differential equations is solved with a generic modified Alternating Direction Implicit (m-ADI) method developed to resolve locally refined grids. The overall performance of the new flow solver is initially evaluated against experimental data related to the laminar flow past a square cylinder for Reynolds numbers, defined based on the mean ambient wind speed and the characteristic length scale, of 87≤Re≤250 and neutral thermal stratification. The unsteady character of the flow field is successfully predicted and the grid refinement approach enables accurate low cost and low blockage simulations. The flow solver is applied to simulate the flow over thermally stratified idealized urban areas. Good agreement with experimental measurements, where the bulk Richardson number was defined based on the mean wind speed at the top of the reference street canyon, is obtained for -0.21≤Rb≤0.78 while the grid refinement method allowed simulation of the experimental setup in detail, including all the roughness elements and the idealized street canyons. The numerical calculations indicate that even weak thermal stability can lead to laminarization of the induced flow field; reducing at the same time the potential for air mixing and results to air degradation. Even though the recirculation patterns are influenced under weak thermal stratification, extreme thermal instability within the region of -16≤Rb≤-64 with Re=6800 results to the development of two counter-rotating vortices of equivalent kinetic energy. However, the linear correlation of the average kinetic energy at the roof-level with the Rb number indicates that the buoyancy forces determine the air mixing levels.
KW - Atmospheric stability
KW - Cartesian grid refinement
KW - M-ADI solver
KW - Street canyon
KW - Thermal stratification
KW - Urban environment
UR - http://www.scopus.com/inward/record.url?scp=84928350515&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84928350515&partnerID=8YFLogxK
U2 - 10.1016/j.jweia.2015.03.012
DO - 10.1016/j.jweia.2015.03.012
M3 - Article
AN - SCOPUS:84928350515
SN - 0167-6105
VL - 142
SP - 149
EP - 163
JO - Journal of Wind Engineering and Industrial Aerodynamics
JF - Journal of Wind Engineering and Industrial Aerodynamics
ER -