In this study, large-eddy simulation is combined with a turbine model to investigate the influence of atmospheric thermal stability on wind-turbine wakes. The simulation results show that atmospheric stability has a significant effect on the spatial distribution of the mean velocity deficit and turbulence statistics in the wake region as well as the wake meandering characteristics downwind of the turbine. In particular, the enhanced turbulence level associated with positive buoyancy under the convective condition leads to a relatively larger flow entrainment and, thus, a faster wake recovery. For the particular cases considered in this study, the growth rate of the wake is about 2.4 times larger for the convective case than for the stable one. Consistent with this result, for a given distance downwind of the turbine, wake meandering is also stronger under the convective condition compared with the neutral and stable cases. It is also shown that, for all the stability cases, the growth rate of the wake and wake meandering in the vertical direction are smaller compared with the ones in the lateral direction. This is mainly related to the different turbulence levels of the incoming wind in the different directions, together with the anisotropy imposed by the presence of the ground. It is also found that the wake velocity deficit is well characterized by a modified version of a recently proposed analytical model that is based on mass and momentum conservation and the assumption of a self-similar Gaussian distribution of the velocity deficit. Specifically, using a two-dimensional elliptical (instead of axisymmetric) Gaussian distribution allows to account for the different lateral and vertical growth rates, particularly in the convective case where the non-axisymmetry of the wake is stronger. Detailed analysis of the resolved turbulent kinetic energy budget in the wake reveals also that thermal stratification considerably affects the magnitude and spatial distribution of the turbulence production, dissipation and transport terms.

Large-eddy simulation is used to study the influence of free-atmosphere stratification on the structure of atmospheric boundary-layer flow inside and above very large wind farms, as well as the power extracted by the wind turbines. In the simulations, tuning-free Lagrangian scale-dependent dynamic models are used to model the subgrid-scale turbulent fluxes, while the turbine-induced forces are parameterized with an actuator-disk model. It is shown that for a given surface cover (with and without turbines) thermal stratification of the free atmosphere reduces the entrainment from the flow above compared with the unstratified case, leading to lower boundary-layer depth. Due to the fact that in very large wind farms vertical energy transport associated with turbulence is the only source of kinetic energy, lower entrainment leads to lower power production by the wind turbines. In particular, for the wind-turbine arrangements considered in the present work, the power output from the wind farms is reduced by about 35% when the potential temperature lapse rate in the free atmosphere increases from 1 to 10 K/km (within the range of values typically observed in the atmosphere). Moreover, it is shown that the presence of the turbines has significant effect on the growth of the boundary layer. Inspired by the obtained results, a simple one-dimensional model is developed to account for the effect of free-atmosphere stability on the mean flow and the power output from very large wind farms.

No abstract

In this article, a new model is developed to parameterize the effect of wind farms in large-scale atmospheric models such as weather models. In the new model, wind turbines in a wind farm are parameterized as elevated sinks of momentum and sources of turbulence. An analytical approach is used to estimate the turbineinduced forces as well as the turbulent kinetic energy (TKE) generated by the turbines inside the atmospheric boundary layer (ABL). In addition, the proposed model can take into account not only the effect of wind-farm density but also the effect of wind-farm layout and wind direction. The performance of the new model is tested with large-eddy simulations of ABL flows over very large wind farms with different turbine configurations. The results show that the new model is capable to accurately predict the turbine-induced forces as well as the TKE generated by the turbines inside the ABL. V C 2015 AIP Publishing LLC.

Minimum-dissipation models are a simple alternative to the Smagorinsky-type approaches to parametrize the subfilter turbulent fluxes in large-eddy simulation. A recently derived model of this type for subfilter stress tensor is the anisotropic minimum-dissipation (AMD) model [Rozema et al., Phys. Fluids 27, 085107 (2015)], which has many desirable properties. It is more cost effective than the dynamic Smagorinsky model, it appropriately switches off in laminar and transitional flows, and it is consistent with the exact subfilter stress tensor on both isotropic and anisotropic grids. In this study, an extension of this approach to modeling the subfilter scalar flux is proposed. The performance of the AMD model is tested in the simulation of a high-Reynolds-number rough-wall boundary-layer flow with a constant and uniform surface scalar flux. The simulation results obtained from the AMD model show good agreement with well-established empirical correlations and theoretical predictions of the resolved flow statistics. In particular, the AMD model is capable of accurately predicting the expected surface-layer similarity profiles and power spectra for both velocity and scalar concentration.

A generalized form of a recently developed minimum dissipation model for subfilter turbulent fluxes is proposed and implemented in the simulation of thermally stratified atmospheric boundary-layer flows. Compared with the original model, the generalized model includes the contribution of buoyant forces, in addition to shear, to the production or suppression of turbulence, with a number of desirable practical and theoretical properties. Specifically, the model has a low computational complexity, appropriately switches off in laminar and transitional flows, does not require any ad hoc shear and stability corrections, and is consistent with theoretical subfilter turbulent fluxes. The simulation results show remarkable agreement with well-established empirical correlations, theoretical predictions, and field observations in the atmosphere. In addition, the results show very little sensitivity to the grid resolution, demonstrating the robustness of the model in the simulation of the atmospheric boundary layer, even with relatively coarse resolutions.

No abstract

No abstract