No abstract

Relationships between parameters of convective entrainment into a shear-free, linearly stratified atmosphere predicted by the zero-order jump and general-structure bulk models of entrainment are reexamined using data from large eddy simulations (LESs). Relevant data from other numerical simulations, water tank experiments, and atmospheric measurements are also incorporated in the analysis. Simulations have been performed for 10 values of the buoyancy gradient in the free atmosphere covering a typical atmospheric stability range. The entrainment parameters derived from LES and relationships between them are found to be sensitive to the model framework employed for their interpretation. Methods of determining bulk model entrainment parameters from the LES output are proposed and discussed. Within the range of investigated free-atmosphere stratifications, the LES predictions of the inversion height and buoyancy increment across the inversion are found to be close to the analytical solutions for the equilibrium entrainment regime, which is realized when the rate of time change of the CBL-mean turbulence kinetic energy and the energy drain from the CBL top are both negligibly small. The zero-order model entrainment ratio of about 0.2 for this regime is generally supported by the LES data. However, the zero-order parameterization of the entrainment layer thickness is found insufficient. A set of relationships between the general-structure entrainment parameters for typical atmospheric stability conditions is retrieved from the LES. Dimensionless constants in these relationships are estimated from the LES and laboratory data. Power-law approximations for relationships between the entrainment parameters in the zero-order jump and general-structure bulk models are evaluated based on the conducted LES. In the regime of equilibrium entrainment, the stratification parameter of the entrainment layer, which is the ratio of the buoyancy gradient in the free atmosphere to the overall buoyancy gradient across the entrainment layer, appears to be a constant of about 1.2.

The reported study examines the dynamics of entrainment and its effects on the evolution of the dry atmospheric convective boundary layer (CBL) when wind shear is present. The sheared CBL can be studied by means of direct measurements in the atmosphere, laboratory studies, and numerical techniques. The advantages and disadvantages of each technique are discussed in the present paper, which also describes the methodological background for studying the dynamics of entrainment in sheared CBLs. For the reported study, large-eddy simulation (LES) was chosen as the primary method of convective entrainment investigation. Twenty-four LES runs were conducted for CBLs growing under varying conditions of surface buoyancy flux, free-atmospheric stratification, and wind shear. The simulations were divided into three categories: CBL with no mean wind (NS), CBL with a height-constant geostrophic wind of 20 m s Ϫ1 (GC), and CBL with geostrophic wind shear (GS). In the simulated cases, the sheared CBLs grew fastest, relative to the NS CBLs, when the surface buoyancy flux was weak and the atmospheric stratification was moderate or weak.Three fundamental findings resulted from the investigated CBL cases: (i) the entrainment zone shear is much more important than the surface shear in enhancing CBL entrainment, although entrainment zone shear is indirectly affected by surface shear; (ii) the sheared entrainment zone features a sublayer of nearly constant flux Richardson number, which points to a balance between shear production and buoyancy consumption of turbulence kinetic energy (TKE) that regulates entrainment; and (iii) the fraction of entrainment zone shear-generated TKE spent on the entrainment is lower than suggested by earlier studies.

A large‐eddy simulation framework, dubbed as the Virtual Wind Simulator (VWiS), for simulating turbulent flow over wind turbines and wind farms in complex terrain is developed and validated. The wind turbines are parameterized using the actuator line model. The complex terrain is represented by the curvilinear immersed boundary method. The predictive capability of the present method is evaluated by simulating two available wind tunnel experimental cases: the flow over a stand‐alone turbine and an aligned wind turbine array. Systematic grid refinement studies are carried out, for both single turbine and multi‐turbine array cases, and the accuracy of the computed results is assessed through detailed comparisons with wind tunnel experiments. The model is further applied to simulate the flow over an operational utility‐scale wind farm. The inflow velocities for this case are interpolated from a mesoscale simulation using a Weather Research and Forecasting (WRF) model with and without adding synthetic turbulence to the WRF‐computed velocity fields. Improvements on power predictions are obtained when synthetic turbulence is added at the inlet. Finally the VWiS is applied to simulate a yet undeveloped wind farm at a complex terrain site where wind resource measurements have already been obtained. Good agreement with field measurements is obtained in terms of the time‐averaged streamwise velocity profiles. To demonstrate the ability of the model to simulate the interactions of terrain‐induced turbulence with wind turbines, eight hypothetical turbines are placed in this area. The computed extracted power underscores the significant effect of site‐specific topography on turbine performance. Copyright © 2014 John Wiley & Sons, Ltd.