Modelling Radiative and Convective Thermal Exchanges over a European City Center and Their Effects on Atmospheric Dispersion

Micro-meteorological studies of urban flow and pollution dispersion often assume a neutral atmosphere and often the three-dimensional variation in temperature fields and flow around buildings is neglected in most building energy balance models. The aim of this work is to present the results of development and validation of a three-dimensional tool coupling thermal energy balance of the buildings and modelling of the atmospheric flow and dispersion in urban areas. To do so, a 3D microscale atmospheric radiative scheme has been developed in the atmospheric module of the computational fluid dynamics (CFD) code Code_Saturne adapted to detailed building geometries. The full coupling of the radiative transfer and fluid dynamics models has been validated with idealized cases. In this paper, our focus is to simulate and compare with measurements the diurnal evolution of the brightness surface temperatures and the momentum and energy fluxes for a neighborhood in the city center of Toulouse, in the southwest part of France. This is performed by taking into account the 3D effects of the flow around the buildings and all thermal exchanges, in real meteorological conditions, and compare them to aircraft infrared images and in situ measurements on a meteorological mast. The calculation mesh developed for the city center and the simulation conditions for the selected day of the field campaign are presented. The results are evaluated with the measurements from the Canopy and Aerosol Particles Interactions in TOulouse Urban Layer experiment (CAPITOUL). In addition, the second purpose of this work is to investigate a hypothetical release of passive pollutant dispersion in the same area of Toulouse under different thermal transfer conditions for the street and the buildings surfaces: neutral and 3D radiative transfer heating. The presence of heat transfer continually modifies the airflow field while the airflow in the neutral case reaches a stationary state. Compared to the neutral case, taking into account the thermal transfer enhances the turbulence kinetic energy and vertical velocity (especially at the roof level) due to buoyancy forces. The simulation results also show that the thermal effects considerably alter the plume shape.