Abstract. In this paper we present the current version of the Parallelized Large-Eddy Simulation Model (PALM) whose core has been developed at the Institute of Meteorology and Climatology at Leibniz Universität Hannover (Germany). PALM is a Fortran 95-based code with some Fortran 2003 extensions and has been applied for the simulation of a variety of atmospheric and oceanic boundary layers for more than 15 years. PALM is optimized for use on massively parallel computer architectures and was recently ported to general-purpose graphics processing units. In the present paper we give a detailed description of the current version of the model and its features, such as an embedded Lagrangian cloud model and the possibility to use Cartesian topography. Moreover, we discuss recent model developments and future perspectives for LES applications.

No abstract

No abstract

A large-eddy simulation (LES) study is presented that investigates the spatial variability of temporal eddy covariance fluxes and the systematic underestimation of representative fluxes linked to them. It extends a prior numerical study by performing high resolution simulations that allow for virtual measurements down to 20 m in a convective boundary layer, so that conditions for small tower measurement sites can be analysed. It accounts for different convective regimes as the wind speed and the near-surface heat flux are varied. Moreover, it is the first LES imbalance study that extends to the stable boundary layer. It reveals shortcomings of single site measurements and the necessity of using horizontally-distributed observation networks. The imbalances in the convective case are attributed to a locally non-vanishing mean vertical advection due to turbulent organised structures (TOS). The strength of the TOS and thus the imbalance magnitude depends on height, the horizontal mean wind and the convection type. Contrary to the results of a prior study, TOS cannot generally be responsible for large energy imbalances: at low observation heights (corresponding to small towers and near-surface energy balance stations) the TOS related imbalances are generally about one order of magnitude smaller than those in field experiments. However, TOS may cause large imbalances at large towers not only in the case of cellular convection and low wind speeds, as found in the previous study, but also in the case of roll convection at large wind speeds.In the stably stratified boundary layer for all observation heights neither TOS nor significant imbalances are observed. Boundary-Layer Meteorol (2007) 123:77-98 Attempting to reduce imbalances in convective situations by applying the conventional linear detrending method increases the systematic flux underestimation. Thus, a new filter method is proposed.

Mesoscale circulations induced by differential boundary layer heating due to surface inhomogeneities on scales of 5 km and more can significantly change the average properties and the structure of the convective boundary layer (CBL) as well as trigger off temporal oscillations. The results of one of the first numerical case studies using large eddy simulation (LES) on the mesoscale suggest that mesoscale circulations exhibit a considerably larger average kinetic energy than convection under homogeneous conditions. This affects turbulent transport processes and should be accounted for in larger-scale models even if their turbulence parameterizations rely on homogeneous control runs of high-resolution models. This case study uses the Hannover parallelized large eddy simulation model (PALM) with prescribed 1D sinusoidal surface heat flux variations on wavelengths from 2.5 to 40 km. The resulting mesoscale circulations are analyzed by means of domain-averaged cross sections, time averaged and normalized with the boundary layer height, as well as power spectra and domain-averaged time series. The simulated mesoscale circulations were periodic. Vertical profiles and time series demonstrate that the onset of the mesoscale circulation triggers off a temporal boundary layer oscillation, whose period and amplitude depend on the surface heat flux perturbation wavelength and amplitude and on the background wind component perpendicular to the surface inhomogeneity orientation. Spectral analysis shows that the mesoscale circulations damp convection equally in all directions. A hypothesis of the oscillation mechanism is briefly discussed.

Recent observational studies have described the non-closure of the energy balance when the eddy covariance (EC) method is used for the measurements. We investigated this problem using a numerical simulation of a heterogeneous surface region. A typical daytime boundary layer was simulated, using the large eddy simulation (LES) method in which horizontal heterogeneity was imposed on the ground surface heating as a one-dimensional sinusoidal variation. This horizontal heterogeneity is expected to produce a mesoscale circulation.We decomposed the total vertical heat flux into the EC turbulent flux, the heat flux due to a mesoscale circulation (hereafter, mesoscale flux), and the ''residual flux''. The sum of the mesoscale flux, and residual flux accounts for the energy imbalance if we estimate the total flux only from the EC method.The numerical results demonstrated that larger amplitude of surface heating caused larger mesoscale flux, but smaller residual flux. As a result, the energy imbalance became minima at some weak amplitude of surface heating.The residual flux was caused by the turbulent organized structure (hereafter, TOS), which is a cluster of thermals moving, with a larger time scale than that of individual plumes. The larger surface heating amplitude weakened the TOS due to the following two mechanisms; (1) the TOS is organized in roll due to the strong horizontal pressure gradient, (2) the higher horizontal wind speed, parallel to the mesoscale circulation, advects the TOS faster then the ergodicity works better.The other cases with gepstrophic winds, resulted in the decrease of the energy imbalance with increasing wind velocity.

Mesoscale circulations induced by differential boundary layer heating due to surface inhomogeneities on scales of 5 km and more can significantly change the average properties and the structure of the convective boundary layer (CBL) as well as trigger off temporal oscillations. The results of one of the first numerical case studies using large eddy simulation (LES) on the mesoscale suggest that mesoscale circulations exhibit a considerably larger average kinetic energy than convection under homogeneous conditions. This affects turbulent transport processes and should be accounted for in larger-scale models even if their turbulence parameterizations rely on homogeneous control runs of high-resolution models. This case study uses the Hannover parallelized large eddy simulation model (PALM) with prescribed 1D sinusoidal surface heat flux variations on wavelengths from 2.5 to 40 km. The resulting mesoscale circulations are analyzed by means of domain-averaged cross sections, time averaged and normalized with the boundary layer height, as well as power spectra and domain-averaged time series. The simulated mesoscale circulations were periodic. Vertical profiles and time series demonstrate that the onset of the mesoscale circulation triggers off a temporal boundary layer oscillation, whose period and amplitude depend on the surface heat flux perturbation wavelength and amplitude and on the background wind component perpendicular to the surface inhomogeneity orientation. Spectral analysis shows that the mesoscale circulations damp convection equally in all directions. A hypothesis of the oscillation mechanism is briefly discussed.

Thermal effects on turbulent flow and dispersion in and above an idealized street canyon with a street aspect ratio of 1 are numerically investigated using the parallelized large-eddy simulation model (''PALM''). Each of upwind building wall, street bottom, and downwind building wall is heated, and passive scalars are emitted from the street bottom. When compared with the neutral (no heating) case, the heating of the upwind building wall or street bottom strengthens a primary vortex in the street canyon and the heating of the downwind building wall induces a shrunken primary vortex and a winding flow between the vortex and the downwind building wall. Heating also induces higher turbulent kinetic energy and stronger turbulent fluxes at the rooftop height. In the neutral case, turbulent eddies generated by shear instability dominate mixing at the rooftop height and appear as band-shaped perturbations in the time-space plots of turbulent momentum and scalar fluxes. In all of the heating cases, buoyancy-generated turbulent eddies as well as shear-generated turbulent eddies contribute to turbulent momentum and scalar fluxes and band-shaped or lump-shaped perturbations appear at the rooftop height. A quadrant analysis shows that at the rooftop height, in the neutral case and in the case with upwind building-wall heating, sweep events are less frequent but contribute more to turbulent momentum flux than do ejection events. By contrast, in the case with street-bottom and downwind building-wall heating, the frequency of sweep events is similar to that of ejection events and the contribution of ejection events to turbulent momentum flux is comparable to that of sweep events.