In this paper, we describe the PALM model system 6.0. PALM (formerly an abbreviation for Parallelized Large-eddy Simulation Model and now an independent name) is a Fortran-based code and has been applied for studying a variety of atmospheric and oceanic boundary layers for about 20 years. The model is optimized for use on massively parallel computer architectures. This is a follow-up paper to the PALM 4.0 model description in Maronga et al. (2015). During the last years, PALM has been significantly improved and now offers a variety of new components. In particular, much effort was made to enhance the model with components needed for applications in urban environments, like fully interactive land surface and radiation schemes, chemistry, and an indoor model. This paper serves as an overview paper of the PALM 6.0 model system and we describe its current model core. The individual components for urban applications, case studies, validation runs, and issues with suitable input data are presented and discussed in a series of companion papers in this special issue.

In this paper, we present a newly developed mesoscale nesting interface for the PALM model system 6.0, which enables PALM to simulate the atmospheric boundary layer under spatially heterogeneous and non-stationary synoptic conditions. The implemented nesting interface, which is currently tailored to the mesoscale model COSMO, consists of two major parts: (i) the preprocessor INIFOR (initialization and forcing), which provides initial and time-dependent boundary conditions from mesoscale model output, and (ii) PALM's internal routines for reading the provided forcing data and superimposing synthetic turbulence to accelerate the transition to a fully developed turbulent atmospheric boundary layer. We describe in detail the conversion between the sets of prognostic variables, transformations between model coordinate systems, as well as data interpolation onto PALM's grid, which are carried out by INIFOR. Furthermore, we describe PALM's internal usage of the provided forcing data, which, besides the temporal interpolation of boundary conditions and removal of any residual divergence, includes the generation of stability-dependent synthetic turbulence at the inflow boundaries in order to accelerate the transition from the turbulence-free mesoscale solution to a resolved turbulent flow. We demonstrate and evaluate the nesting interface by means of a semi-idealized benchmark case. We carried out a large-eddy simulation (LES) of an evolving convective boundary layer on a clear-sky spring day. Besides verifying that changes in the inflow conditions enter into and successively propagate through the PALM domain, we focus our analysis on the effectiveness of the synthetic turbulence generation. By analysing various turbulence statistics, we show that the inflow in the present case is fully adjusted after having propagated for about two to three eddy-turnover times downstream, which corresponds well to other state-of-the-art methods for turbulence generation. Furthermore, we observe that numerical artefacts in the form of grid-scale convective structures in the mesoscale model enter the PALM domain, biasing the location of the turbulent up- and downdrafts in the LES. With these findings presented, we aim to verify the mesoscale nesting approach implemented in PALM, point out specific shortcomings, and build a baseline for future improvements and developments.

Ventilation in cities is crucial for the well being of their inhabitants. Therefore, local governments require air ventilation assessments (AVAs) prior to the construction of new buildings. In a standard AVA, however, only neutral stratification is considered, although diabatic and particularly unstable conditions may be observed more frequently in nature. The results presented here indicate significant changes in ventilation within most of the area of Kowloon City, Hong Kong, included in the study. A new definition for calculating ventilation was introduced, and used to compare the influence of buildings on ventilation under conditions of neutral and unstable stratification. The overall ventilation increased due to enhanced vertical mixing. In the vicinity of exposed buildings, however, ventilation was weaker for unstable stratification than for neutral stratification. The influence on ventilation by building parameters, such as the plan area index, was altered when unstable stratification was considered. Consequently, differences in stratification were shown to have marked effects on ventilation estimates, which should be taken into consideration in future AVAs.

We demonstrate the capability of the PALM model system version 6.0 to simulate neutrally stratified urban boundary layers. Our simulation uses the real-world building configuration of the HafenCity area in Hamburg, Germany. Using PALM's virtual measurement module, we compare simulation results to wind-tunnel measurements of a downscaled replica of the study area. Wind-tunnel measurements of mean wind speed agree within 5 % on average while the wind direction deviates by approximately 4 ^∘. Turbulence statistics similarly agree. However, larger differences between measurements and simulation arise in the vicinity of surfaces where building geometry is insufficiently resolved. We discuss how to minimize these differences by improving the grid layout and give tips for setup preparation. Also, we discuss how existing and upcoming features of PALM like the grid nesting and immersed boundary condition help improve the simulation results.

Courtyards are an omnipresent feature within the urban environment. Residents often use courtyards as recreation areas, which makes them crucial for the physical and psychological comfort of the urban population. However, considering that courtyards represent enclosed cavities, they are often poorly ventilated spaces and pollutants from neighboring traffic, once entrained, can pose a serious threat to human health. Here, we studied the effects of lateral openings on courtyard pollution and ventilation. Therefore, we performed a set of large-eddy simulations for idealized urban environments with different courtyard configurations. While pollutant concentration and ventilation are barely modified by lateral openings for wide courtyards, lateral openings have a significant effect on the mean concentration, the number of high-concentration events and the ventilation within narrower and deeper courtyards. The impacts of lateral openings on air quality within courtyards strongly depend on their orientation with respect to the flow direction, as well as on the upstream flow conditions and upstream building configuration. We show that lateral openings, in most cases, have a negative impact on air quality; nevertheless, we also present configurations where lateral openings positively impact the air quality within courtyards. These outcomes may certainly contribute to improve future urban planning in terms of health protection.