As a consequence of the rapid growth of the globally installed offshore wind energy capacity, the size of individual wind farms is increasing. This poses a challenge to models that predict energy production. For instance, the current generation of wake models has mostly been calibrated on existing wind farms of much smaller size. This work analyzes annual energy production and wake losses for future, multi-gigawatt wind farms with atmospheric large-eddy simulation. To that end, 1 year of actual weather has been simulated for a suite of hypothetical 4 GW offshore wind farm scenarios. The scenarios differ in terms of applied turbine type, installed capacity density, and layout. The results suggest that production numbers increase significantly when the rated power of the individual turbines is larger while keeping the total installed capacity the same. Even for turbine types with similar rated power but slightly different power curves, significant differences in production were found. Although wind speed was identified as the most dominant factor determining the aerodynamic losses, a clear impact of atmospheric stability and boundary layer height has been identified. By analyzing losses of the first-row turbines, the yearly average global-blockage effect is estimated to between 2 and 3 %, but it can reach levels over 10 % for stably stratified conditions and wind speeds around 8 m s−1. Using a high-fidelity modeling technique, the present work provides insights into the performance of future, multi-gigawatt wind farms for a full year of realistic weather conditions.@en

This study investigates how wind shear and momentum fluxes in the surface- and boundary layer vary across wind and cloud regimes. We use a 9-year-long data set from the Cabauw observatory complemented by (8.2 × 8.2 (Formula presented.)) daily Large Eddy Simulation (LES) hindcasts. An automated algorithm classifies observed and simulated days into different cloud regimes: (a) clear-sky days, (b) days with shallow convective clouds rooted in the surface layer, with two ranges of cloud cover, and (c) non-convective cloud days. Categorized days in observations and LES do not always match, particularly the number of non-convective cloud days are underestimated in the LES, which likes to develop convection. However, the climatology and diurnal cycle of winds for each regime are very similar in LES and observations, strengthening our confidence in LES’ skill to reproduce certain clouds for certain atmospheric states. Along-wind momentum flux profiles are similar across all regimes, but large cloud cover (convective and non-convective) days have larger total momentum flux distributed over a deeper layer, with up to 30% of the surface flux still present near cloud base. The clear-sky and especially shallow cumulus regime with low cloud cover have notably larger crosswind momentum fluxes in the boundary layer. Surface-layer wind shear at daytime is smallest in the shallow cumulus regimes, having deeper boundary layers and a steady increase in surface layer wind speed during daytime. Compared to clear-sky days at a similar stability, convective cloud regimes have smaller surface-layer wind shear and larger surface friction than estimated by Monin-Obukhov Similarity Theory.

@en

This study investigated the surface temperature, air temperature and mean radiant temperature inside an idealized 2D street geometry during daytime. The goal was to unravel the relative impact of radiative transfer, heat conduction and ventilation to the urban heat budget. A building-resolving simulation model has been used, which represents these processes at a 1 m spatial resolution. Different combinations of the canyon height to width ratio (H/W) and physical mechanisms were investigated. Shortwave radiation is the main source of energy, and for small H/W can be higher at the canyon ground level compared to flat terrain due to multiple reflections. The longwave trapping effect has the second largest contribution and becomes relatively more important with increasing H/W ratio. The influence of the interior building temperature is small. Surface temperature and mean radiant temperature are closely related, since both are largely controlled by radiative properties. No straightforward relation was found between surface temperature and air temperature, since air temperature is dependent on the competing mechanisms of forced and natural convection. A small increase in air temperature inside the canyon was observed compared to the ambient temperature above roof level. The inclusion of all key physical processes in high detail resulted in large computational requirements. If multiple reflections by the building facades are small, the more traditional, yet much simpler view factor approach will strongly reduce the computational costs as compared to the Monte Carlo technique. The influence of using the view factors on the results must be investigated.

@en

High-resolution large-eddy simulations of the Antarctic very stable boundary layer reveal a mechanism for systematic and periodic intermittent bursting. A nonbursting state with a boundary layer height of just 3 m is alternated by a bursting state with a height of ≈5 m. The bursts result from unstable wave growth triggered by a shear-generated Kelvin–Helmholtz instability, as confirmed by linear stability analysis. The shear at the top of the boundary layer is built up by two processes. The upper, quasi-laminar layer accelerates due to the combined effect of the pressure force and rotation by the Coriolis force, while the lower layer decelerates by turbulent friction. During the burst, this shear is eroded and the initial cause of the instability is removed. Subsequently, the interfacial shear builds up again, causing the entire sequence to repeat itself with a time scale of ≈10 min. Despite the clear intermittent bursting, the overall change of the mean wind profile is remarkably small during the cycle. This enables such a fast erosion and recovery of the shear. This mechanism for cyclic bursting is remarkably similar to the mechanism hypothesized by Businger in 1973, with one key difference. Whereas Businger proposes that the flow acceleration in the upper layer results from downward turbulent transfer of high-momentum flow, the current results indicate no turbulent activity in the upper layer, hence requiring another source of momentum. Finally, it would be interesting to construct a climatology of shear-generated intermittency in relation to large-scale conditions to assess the generality of this Businger mechanism.@en

The authors regret to inform the readers that a programming error was found in the numerical code described and used in the original submission [1]. Longwave trapping was not correctly scaled with the local cell size and is underpredicted for deep canyons. Calculations have been repeated and this corrigendum provides an overview of the new results for a deep canyon. All new results will be published in a forthcoming Phd thesis [2] that will made publicly available through the TU Delft website (https://repository.tudelft.nl/). The authors would like to apologise for any inconvenience caused.

@en

Observations of two typical contrasting weakly stable and very stable boundary layers from the winter at Dome C station, Antarctica, are used as a benchmark for two centimetre-scale-resolution large-eddy simulations. By taking the Antarctic winter, the effects of the diurnal cycle are eliminated, enabling the study of the long-lived steady stable boundary layer. With its homogeneous, flat snow surface, and extreme stabilities, the location is a natural laboratory for studies on the long-lived stable boundary layer. The two simulations differ only in the imposed geostrophic wind speed, which is identified as the main deciding factor for the resulting regime. In general, a good correspondence is found between the observed and simulated profiles of mean wind speed and temperature. Discrepancies in the temperature profiles are likely due to the exclusion of radiative transfer in the current simulations. The extreme stabilities result in a considerable contrast between the stable boundary layer at the Dome C site and that found at typical mid-latitudes. The boundary-layer height is found to range from approximately 50m to just 5m in the most extreme case. Remarkably, heating of the boundary layer by subsidence may result in thermal equilibrium of the boundary layer in which the associated heating is balanced by the turbulent cooling towards the surface. Using centimetre-scale resolutions, accurate large-eddy simulations of the extreme stabilities encountered in Antarctica appear to be possible. However, future simulations should aim to include radiative transfer and sub-surface heat transport to increase the degree of realism of these types of simulations.@en

Accurate short-term power forecasts are crucial for the reliable and efficient integration of wind energy in power systems and electricity markets. Typically, forecasts for hours to days ahead are based on the output of numerical weather prediction models, and with the advance of computing power, the spatial and temporal resolutions of these models have increased substantially. However, high-resolution forecasts often exhibit spatial and/or temporal displacement errors, and when regarding typical average performance metrics, they often perform worse than smoother forecasts from lower-resolution models. Recent computational advances have enabled the use of large-eddy simulations (LESs) in the context of operational weather forecasting, yielding turbulence-resolving weather forecasts with a spatial resolution of 100 m or finer and a temporal resolution of 30 seconds or less. This paper is a proof-of-concept study on the prospect of leveraging these ultra high-resolution weather models for operational forecasting at Horns Rev I in Denmark. It is shown that temporal smoothing of the forecasts clearly improves their skill, even for the benchmark resolution forecast, although potentially valuable high-frequency information is lost. Therefore, a statistical post-processing approach is explored on the basis of smoothing and feature engineering from the high-frequency signal. The results indicate that for wind farm forecasting, using information content from both the standard and LES resolution models improves the forecast accuracy, especially with a feature selection stage, compared with using the information content solely from either source.

@en

Turbulent entrainment is a process of primary importance in the atmospheric boundary layer; however despite several decades of intense study much remains to be understood. Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES) have a tremendous potential to improve the understanding of turbulent entrainment, particularly if combined with theory. We discuss a recently developed framework for turbulent jets and plumes to decompose turbulent entrainment in various physical processes, and modify it for use in a stably stratified shear driven (nocturnal) boundary layer. The decomposition shows that inner layer processes become negligible as time progresses and that the entrainment coefficient is determined by turbulence production in the outer layer only.@en

Simultaneous particle-image velocimetry and laser-induced fluorescence combined with large-eddy simulations are used to investigate the flow and pollutant dispersion behaviour in a rural-to-urban roughness transition. The urban roughness is characterized by an array of cubical obstacles in an aligned arrangement. A plane fence is added one obstacle height h upstream of the urban roughness elements, with three different fence heights considered. A smooth-wall turbulent boundary layer with a depth of 10h is used as the approaching flow, and a passive tracer is released from a uniform line source 1h upstream of the fence. A shear layer is formed at the top of the fence, which increases in strength for the higher fence cases, resulting in a deeper internal boundary layer (IBL). It is found that the mean flow for the rural-to-urban transition can be described by means of a mixing-length model provided that the transitional effects are accounted for. The mixing-length formulation for sparse urban canopies, as found in the literature, is extended to take into account the blockage effect in dense canopies. Additionally, the average mean concentration field is found to scale with the IBL depth and the bulk velocity in the IBL.

@en

Large-eddy simulation (LES) models are widely used to study atmospheric turbulence. The effects of small-scale motions that cannot be resolved need to be modeled by a subfilter-scale (SFS) model. The SFS contribution to the turbulent fluxes is typically significant in the surface layer. This study presents analytical solutions of the classical Smagorinsky SFS turbulent kinetic energy (TKE) model including a buoyancy flux contribution. Both a constant length scale and a stability-dependent one as proposed by Deardorff are considered. Analytical expressions for the mixing functions are derived and Monin-Obukhov similarity relations that are implicitly imposed by the SFS TKE model are diagnosed. For neutral and weakly stable conditions, observations indicate that the turbulent Prandtl number (PrT) is close to unity. However, based on observations in the convective boundary layer, a lower value for PrT is often applied in LES models. As a lower Prandtl number promotes a stronger mixing of heat, this may cause excessive mixing, which is quantified from a direct comparison of the mixing function as imposed by the SFS TKE model with empirical fits from field observations. For a strong stability, the diagnosed mixing functions for both momentum and heat are larger than observed. The problem of excessive mixing will be enhanced for anisotropic grids. The findings are also relevant for high-resolution numerical weather prediction models that use a Smagorinsky-type TKE closure.

@en

Both large-eddy simulations (LES) and water-tunnel experiments, using simultaneous stereoscopic particle image velocimetry and laser-induced fluorescence, have been used to investigate pollutant dispersion mechanisms in regions where the surface changes from rural to urban roughness. The urban roughness was characterized by an array of rectangular obstacles in an in-line arrangement. The streamwise length scale of the roughness was kept constant, while the spanwise length scale was varied by varying the obstacle aspect ratio l / h between 1 and 8, where l is the spanwise dimension of the obstacles and h is the height of the obstacles. Additionally, the case of two-dimensional roughness (riblets) was considered in LES. A smooth-wall turbulent boundary layer of depth 10h was used as the approaching flow, and a line source of passive tracer was placed 2h upstream of the urban canopy. The experimental and numerical results show good agreement, while minor discrepancies are readily explained. It is found that for (Formula presented.) the drag induced by the urban canopy is largest of all considered cases, and is caused by a large-scale secondary flow. In addition, due to the roughness transition the vertical advective pollutant flux is the main ventilation mechanism in the first three streets. Furthermore, by means of linear stochastic estimation the mean flow structure is identified that is responsible for street-canyon ventilation for the sixth street and onwards. Moreover, it is shown that the vertical length scale of this structure increases with increasing aspect ratio of the obstacles in the canopy, while the streamwise length scale does not show a similar trend.

@en

This study investigates the effect of different high-albedo adaptation strategies on air temperature, mean radiant temperature and the Universal Temperature Climate Index (UTCI) for an idealized 2D street canyon. The used numerical model computes the heat transport in the canyon, and specifically takes into account the effect of multiple scattering of radiation. In general the mean radiant temperature has a much larger impact on the UTCI than the air temperature. Moreover, the mean radiant temperature exhibits strong spatial variations in the canyon due to its sensitivity to shading. The impact of albedo-differences on the UTCI is thus relatively small compared to the large shading effects. The best strategy to minimize the UTCI for the outdoor environment with building height to width ratio H/W = 0.5 is found to be a uniform albedo of 0.2. For H/W = 1.0, an albedo gradient from high at the bottom part to low at the top of the vertical walls showed the lowest UTCI. Although using high-albedo materials can mitigate the atmospheric urban heat island effect, it is very likely to increase pedestrian heat stress@en

The predictability horizon of convective boundary layers is investigated in this study. Large-eddy simulation (LES) and direct numerical simulation (DNS) techniques are employed to probe the evolution of perturbations in identical twin simulations of a growing dry convective boundary layer. Error growth typical of chaotic systems is observed, marked by two phases. The first comprises an exponential error growth as , with δ0 as the initial error, δ(t) as the error at time t, and Λ as the Lyapunov exponent. This phase is independent of the perturbation wavenumber, and the perturbation energy grows following a self-similar spectral shape dominated by higher wavenumbers. The nondimensional error growth rate in this phase shows a strong dependence on the Reynolds number (Re). The second phase involves saturation of the error. Here, the error growth follows Lorenz dynamics with a slower saturation of successively larger scales. An analysis of the spectral decorrelation times reveals two regimes: an Re-independent regime for scales larger than the boundary layer height and an Re-dependent regime for scales smaller than , which are found to decorrelate substantially faster for increasing Reynolds numbers@en

We revisit the classical entrainment experiments for gravity currents on inclined slopes (Ellison and Turner, J. Fluid Mech. 6, 423-448, 1959). We derive an entrainment relation that couples the entrainment rate E to the production of turbulence kinetic energy, the net effect of buoyancy and inner layer. Using direct numerical simulations that are run for durations long enough for the flow to reach universal self-similarity, we show that the net effect of inner layer processes on entrainment is very small and that buoyancy has an almost negligible effect on E. It is demonstrated that the dominant process causing entrainment is turbulence production due to shear. Second, we observe that for all simulations the eddy diffusivity and dissipation rate can be parameterised using the turbulence kinetic energy and shear parameter. This information can be used to derive an entrainment law which is in good agreement with the Direct Numerical Simulation (DNS) results. We discuss the potential reasons for why this result is significantly different from experiments and the classical entrainment law introduced by Ellison and Turner.@en