Wind energy, together with other renewable energy sources, are expected to grow substantially in the coming decades and play a key role in mitigating climate change and achieving energy sustainability. One of the main challenges in optimizing the design, operation, control, and grid integration of wind farms is the prediction of their performance, owing to the complex multiscale two-way interactions between wind farms and the turbulent atmospheric boundary layer (ABL). From a fluid mechanical perspective, these interactions are complicated by the high Reynolds number of the ABL flow, its inherent unsteadiness due to the diurnal cycle and synoptic-forcing variability, the ubiquitous nature of thermal effects, and the heterogeneity of the terrain. Particularly important is the effect of ABL turbulence on wind-turbine wake flows and their superposition, as they are responsible for considerable turbine power losses and fatigue loads in wind farms. These flow interactions affect, in turn, the structure of the ABL and the turbulent fluxes of momentum and scalars. This review summarizes recent experimental, computational, and theoretical research efforts that have contributed to improving our understanding and ability to predict the interactions of ABL flow with wind turbines and wind farms.

This work is dedicated to systematically studying and predicting the wake characteristics of a yawed wind turbine immersed in a turbulent boundary layer. To achieve this goal, wind tunnel experiments were performed to characterize the wake of a horizontal-axis wind turbine model. A high-resolution stereoscopic particle image velocimetry system was used to measure the three velocity components in the turbine wake under different yaw angles and tip-speed ratios. Moreover, power and thrust measurements were carried out to analyse the performance of the wind turbine. These detailed wind tunnel measurements were then used to perform a budget study of the continuity and Reynolds-averaged Navier–Stokes equations for the wake of a yawed turbine. This theoretical analysis revealed some notable features of the wakes of yawed turbines, such as the asymmetric distribution of the wake skew angle with respect to the wake centre. Under highly yawed conditions, the formation of a counter-rotating vortex pair in the wake cross-section as well as the vertical displacement of the wake centre were shown and analysed. Finally, this study enabled us to develop general governing equations upon which a simple and computationally inexpensive analytical model was built. The proposed model aims at predicting the wake deflection and the far-wake velocity distribution for yawed turbines. Comparisons of model predictions with the wind tunnel measurements show that this simple model can acceptably predict the velocity distribution in the far wake of a yawed turbine. Apart from the ability of the model to predict wake flows in yawed conditions, it can provide valuable physical insight on the behaviour of turbine wakes in this complex situation.

A scale-dependent dynamic subgrid-scale model for large-eddy simulation of turbulent flows is proposed. Unlike the traditional dynamic model, it does not rely on the assumption that the model coefficient is scale invariant. The model is based on a second test-filtering operation which allows us to determine from the simulation how the coefficient varies with scale. The scale-dependent model is tested in simulations of a neutral atmospheric boundary layer. In this application, near the ground the grid scale is by necessity comparable to the local integral scale (of the order of the distance to the wall). With the grid scale and/or the test-filter scale being outside the inertial range, scale invariance is broken. The results are compared with those from (a) the traditional Smagorinsky model that requires specification of the coefficient and of a wall damping function, and (b) the standard dynamic model that assumes scale invariance of the coefficient. In the near-surface region the traditional Smagorinsky and standard dynamic models are too dissipative and not dissipative enough, respectively. Simulations with the scale-dependent dynamic model yield the expected trends of the coefficient as a function of scale and give improved predictions of velocity spectra at different heights from the ground. Consistent with the improved dissipation characteristics, the scale-dependent model also yields improved mean velocity profiles.

Wind-tunnel experiments were performed to study turbulence in the wake of a model wind turbine placed in a boundary layer developed over rough and smooth surfaces. Hot-wire anemometry was used to characterize the cross-sectional distribution of mean velocity, turbulence intensity and kinematic shear stress at different locations downwind of the turbine for both surface roughness cases. Special emphasis was placed on the spatial distribution of the velocity deficit and the turbulence intensity, which are important factors affecting turbine power generation and fatigue loads in wind energy parks. Non-axisymmetric behaviour of the wake is observed over both roughness types in response to the non-uniform incoming boundary-layer flow and the effect of the surface. Nonetheless, the velocity deficit with respect to the incoming velocity profile is nearly axisymmetric, except near the ground in the far wake where the wake interacts with the surface. It is found that the wind turbine induces a large enhancement of turbulence levels (positive added turbulence intensity) in the upper part of the wake. This is due to the effect of relatively large velocity fluctuations associated with helicoidal tip vortices near the wake edge, where the mean shear is strong. In the lower part of the wake, the mean shear and turbulence intensity are reduced with respect to the incoming flow. The non-axisymmetry of the turbulence intensity distribution of the wake is found to be stronger over the rough surface, where the incoming flow is less uniform at the turbine level. In the far wake the added turbulent intensity, its positive and negative contributions and its local maximum decay as a power law of downwind distance (with an exponent ranging from −0.3 to −0.5 for the rough surface, and with a wider variation for the smooth surface). Nevertheless, the effect of the turbine on the velocity defect and added turbulence intensity is not negligible even in the very far wake, at a distance of fifteen times the rotor diameter.