Abstract: The correct prediction of air pollutants dispersed in urban areas is of paramount importance to safety, public health and a sustainable environment. Vehicular traffic is one of the main sources of nitrogen oxides (NO ) and particulate matter (PM), strongly related to human morbidity and mortality. In this study, the pollutant level and distribution in a section of one of the main road arteries of Antwerp (Belgium, Europe) are analyzed. The assessment is performed through computational fluid dynamics (CFD), acknowledged as a powerful tool to predict and study dispersion phenomena in complex atmospheric environments. The two main traffic lanes are modeled as emitting sources and the surrounding area is explicitly depicted. A Reynolds-averaged Navier–Stokes (RANS) approach specific for Atmospheric Boundary Layer (ABL) simulations is employed. After a validation on a wind tunnel urban canyon test case, the dispersion within the canopy of two relevant urban pollutants, nitrogen dioxide (NO) and particulate matter with an aerodynamic diameter smaller than 10 m (PM10), is studied. An experimental field campaign led to the availability of wind velocity and direction data, as well as PM10 concentrations in some key locations within the urban canyon. To accurately predict the concentration field, a relevant dispersion parameter, the turbulent Schmidt number, , is prescribed as a locally variable quantity. The pollutant distributions in the area of interest – exhibiting strong heterogeneity – are finally demonstrated, considering one of the most frequent and concerning wind directions. Possible local remedial measures are conceptualized, investigated and implemented and their outcomes are directly compared. A major goal is, by realistically reproducing the district of interest, to identify the locations inside this intricate urban canyon where the pollutants are stagnating and to analyze which solution acts as best mitigation measure. It is demonstrated that removal by electrostatic precipitation (ESP), an active measure, and by enhancing the dilution process through wind catchers, a passive measure, are effective for local pollutant removal in a realistic urban canyon. It is also demonstrated that the applied ABL methodology resolves some well known problems in ABL dispersion modeling.

Abstract: We report on a novel approach for the Reynolds-averaged Navier-Stokes (RANS) modelling of the neutral atmospheric boundary layer (ABL), using the standard k− turbulence model. A new inlet condition for turbulent kinetic energy is analytically derived from the solution of the k− model transport equations, resulting in a consistent set of fully developed inlet conditions for the neutral ABL. A modification of the standard k− model is also employed to ensure consistency between the inlet conditions and the turbulence model. In particular, the turbulence model constant C μ is generalized as a location-dependent parameter, and a source term is introduced in the transport equation for the turbulent dissipation rate. The application of the proposed methodology to cases involving obstacles in the flow is made possible through the implementation of an algorithm, which automatically switches the turbulence model formulation when going from the region where the ABL is undisturbed to the region directly affected by the building. Finally, the model is completed with a slightly modified version of the Richards and Hoxey rough-wall boundary condition. The methodology is implemented and tested in the commercial code Ansys Fluent 12.1. Results are presented for a neutral boundary layer over flat terrain and for the flow around a single building immersed in an ABL.

Abstract: When considering the modelling of small particle dispersion in the lower part of the Atmospheric Boundary Layer (ABL) using Reynolds Averaged Navier Stokes simulations, the particle paths depend on the velocity profile and on the turbulence kinetic energy, from which the fluctuating velocity components are derived to predict turbulent dispersion. It is therefore important to correctly reproduce the ABL, both for the velocity profile and the turbulence kinetic energy profile. For RANS simulations with the standard kå model, Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the kå turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 4647, 145153.) proposed a set of boundary conditions which result in horizontally homogeneous profiles. The drawback of this method is that it assumes a constant profile of turbulence kinetic energy, which is not always consistent with field or wind tunnel measurements. Therefore, a method was developed which allows the modelling of a horizontally homogeneous turbulence kinetic energy profile that is varying with height. By comparing simulations performed with the proposed method to simulations performed with the boundary conditions described by Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the kå turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 4647, 145153.), the influence of the turbulence kinetic energy on the dispersion of small particles over flat terrain is quantified.