PoisFFT - A Free Parallel Fast Poisson Solver

A fast Poisson solver software package PoisFFT is presented. It is available as a free software licensed under the GNU GPL license version 3. The package uses the fast Fourier transform to directly solve the Poisson equation on a uniform orthogonal grid. It can solve the pseudo-spectral approximation and the second order finite difference approximation of the continuous solution. The paper reviews the mathematical methods for the fast Poisson solver and discusses the software implementation and parallelization. The use of PoisFFT in an incompressible flow solver is also demonstrated

Radiative Transfer Model 3.0 integrated into the PALM model system 6.0

The Radiative Transfer Model (RTM) is an explicitly resolved three-dimensional multi-reflection radiation model integrated into the PALM modelling system. It is responsible for modelling complex radiative interactions within the urban canopy. It represents a key component in modelling energy transfer inside the urban layer and consequently PALM's ability to provide explicit simulations of the urban canopy at metre-scale resolution. This paper presents RTM version 3.0, which is integrated into the PALM modelling system version 6.0. This version of RTM has been substantially improved over previous versions. A more realistic representation is enabled by the newly simulated processes, e.g. the interaction of longwave radiation with the plant canopy, evapotranspiration and latent heat flux, calculation of mean radiant temperature, and bidirectional interaction with the radiation forcing model. The new version also features novel discretization schemes and algorithms, namely the angular discretization and the azimuthal ray tracing, which offer significantly improved scalability and computational efficiency, enabling larger parallel simulations. It has been successfully tested on a realistic urban scenario with a horizontal size of over 6 million grid points using 8192 parallel processes. © 2021 Pavel Krč et al

Sensitivity analysis of the PALM model system 6.0 in the urban environment

Sensitivity of the PALM model 6.0 with respect to land-surface and building properties is tested in a real urban environment in the vicinity of a typical crossroads in a densely built-up residential area in Prague, Czech Republic. The turbulence-resolving PALM is able to simulate the urban boundary layer flow for realistic setups. Besides an accurate representation of the relevant physical processes, the model performance also depends on the input data describing the urban setup, namely the building and land-surface properties. Two types of scenario are employed. The first one is the synthetic scenarios altering mainly surface and material parameters such as albedo, emissivity or wall conductivity, testing sensitivity of the model simulations to potentially erroneous input data. Second, urbanistic-type scenarios are analysed, in which commonly considered urban heat island mitigation measures such as greening of the streets or changing surface materials are applied in order to assess the limits of the effects of a particular type of scenario. For the synthetic scenarios, surface parameters used in radiation balance equations are found to be the most sensitive overall followed by the volumetric heat capacity and thermal conductivity of walls. Other parameters show a limited average effect; however, some can still be significant during some parts of the day, such as surface roughness in the morning hours. The second type, the urbanistic scenarios, shows urban vegetation to be the most effective measure, especially when considering both physical and biophysical temperature indicators. The influence of both types of scenario was also tested for air quality, specifically PM2.5 dispersion, which generally shows opposite behaviour to that of thermal indicators; i.e. improved thermal comfort brings deterioration of PM2.5 concentrations. © 2021 Michal Belda et al

Evaluation of fast atmospheric dispersion models in a regular street network

The need to balance computational speed and simulation accuracy is a key challenge in designing atmospheric dispersion models that can be used in scenarios where near real-time hazard predictions are needed. This challenge is aggravated in cities, where models need to have some degree of building-awareness, alongside the ability to capture effects of dominant urban flow processes. We use a combination of high-resolution large-eddy simulation (LES) and wind-tunnel data of flow and dispersion in an idealised, equal-height urban canopy to highlight important dispersion processes and evaluate how these are reproduced by representatives of the most prevalent modelling approaches: (i) a Gaussian plume model, (ii) a Lagrangian stochastic model and (iii) street-network dispersion models. Concentration data from the LES, validated against the wind-tunnel data, were averaged over the volumes of streets in order to provide a high-fidelity reference suitable for evaluating the different models on the same footing. For the particular combination of forcing wind direction and source location studied here, the strongest deviations from the LES reference were associated with mean over-predictions of concentrations by approximately a factor of 2 and with a relative scatter larger than a factor of 4 of the mean, corresponding to cases where the mean plume centreline also deviated significantly from the LES. This was linked to low accuracy of the underlying flow models/parameters that resulted in a misrepresentation of pollutant channelling along streets and of the uneven plume branching observed in intersections. The agreement of model predictions with the LES (which explicitly resolves the turbulent flow and dispersion processes) greatly improved by increasing the accuracy of building-induced modifications of the driving flow field. When provided with a limited set of representative velocity parameters, the comparatively simple street-network models performed equally well or better compared to the Lagrangian model run on full 3D wind fields. The study showed that street-network models capture the dominant building-induced dispersion processes in the canopy layer through parametrisations of horizontal advection and vertical exchange processes at scales of practical interest. At the same time, computational costs and computing times associated with the network approach are ideally suited for emergency-response applications

A high resolution model of air-flow

Katedra fyziky atmosféryDepartment of Atmospheric PhysicsFaculty of Mathematics and PhysicsMatematicko-fyzikální fakult

A high resolution model of air-flow

Katedra fyziky atmosféryDepartment of Atmospheric PhysicsFaculty of Mathematics and PhysicsMatematicko-fyzikální fakult

Mathematical modelling of air-flow in geometrically complicated areas

The Charles University Large-eddy Microscale Model (CLMM) and its application are presented. It is a numerical model for simulation of turbulent flow and dispersion in the planetary boundary layer. CLMM solves the incompressible Navier-Stokes equations in the Boussinesq approximation and describes turbulence using the large eddy simulation. Three applications of the model are presented. In the first case, the model is applied to the stable boundary layer over a flat terrain. The second case presents the simulations of stably stratified flows over obstacles. The last case deals with the dispersion of a hazardous material within an urban canopy. It was performed in the frame of the COST Action ES1006 and uses the Michelstadt flow and dispersion dataset for model validation

Matematické modelování proudění v geometricky složitých prostředích

Tato práce prezentuje počítačový model Charles University Large-eddy Microscale Model (CLMM) a příklady jeho využití. Je to numerický model pro výpočty proudění a rozptylu příměsí v mezní vrstvě atmosféry. CLMM řeší Navierovy-Stokesovy rovnice nestlačitelného proudění v Boussinesqově aproximaci a popisuje turbulentní proudění pomocí metody simulace velkých vírů. Dále jsou představeny tři příklady použití modelu. V prvním případě je model použit pro simulaci stabilní mezní vrstvy nad plochým terénem. Druhý příklad ukazuje použití modelu pro řešení proudění stabilně teplotně zvrstveného vzduchu přes překážky. Poslední příklad ukazuje výpočet proudění a rozptylu nebezpečné látky v městské zástavbě. Tato poslední studie byla provedena v rámci akce COST ES1006 a používá soubor experimentů "Michelstadt" pro validaci modelu.The Charles University Large-eddy Microscale Model (CLMM) and its application are presented. It is a numerical model for simulation of turbulent flow and dispersion in the planetary boundary layer. CLMM solves the incompressible Navier-Stokes equations in the Boussinesq approximation and describes turbulence using the large eddy simulation. Three applications of the model are presented. In the first case, the model is applied to the stable boundary layer over a flat terrain. The second case presents the simulations of stably stratified flows over obstacles. The last case deals with the dispersion of a hazardous material within an urban canopy. It was performed in the frame of the COST Action ES1006 and uses the Michelstadt flow and dispersion dataset for model validation.Department of Atmospheric PhysicsKatedra fyziky atmosféryFaculty of Mathematics and PhysicsMatematicko-fyzikální fakult

LES of flow in the street canyon

Results of computer simulation of flow over a series of street canyons are presented in this paper. The setup is adapted from an experimental study by [4] with two different shapes of buildings. The problem is simulated by an LES model CLMM (Charles University Large Eddy Microscale Model) and results are analysed using proper orthogonal decomposition and spectral analysis. The results in the channel (layout from the experiment) are compared with results with a free top boundary