Marc Calaf

Large eddy simulation study of fully developed wind-turbine array boundary layers

It is well known that when wind turbines are deployed in large arrays, their efficiency decreases due to complex interactions among themselves and with the atmospheric boundary layer (ABL). For wind farms whose length exceeds the height of the ABL by over an order of magnitude, a “fully developed” flow regime can be established. In this asymptotic regime, changes in the streamwise direction can be neglected and the relevant exchanges occur in the vertical direction. Such a fully developed wind-turbine array boundary layer (WTABL) has not been studied systematically before. A suite of large eddy simulations (LES), in which wind turbines are modeled using the classical “drag disk” concept, is performed for various wind-turbine arrangements, turbine loading factors, and surface roughness values. The results are used to quantify the vertical transport of momentum and kinetic energy across the boundary layer. It is shown that the vertical fluxes of kinetic energy are of the same order of magnitude as the power...

Large eddy simulation study of scalar transport in fully developed wind-turbine array boundary layers

Wind harvesting is fast becoming an important alternative source of energy. As wind farms become larger, they begin to attain scales at which two-way interactions with the atmospheric boundary layer (ABL) must be taken into account. Several studies have shown that there is a quantifiable effect of wind farms on the local meteorology, mainly through changes in the land-atmosphere fluxes of heat and moisture. In particular, the observed trends suggest that wind farms increase fluxes at the surface and this could be due to increased turbulence in the wakes. Conversely, simulations and laboratory experiments show that underneath wind farms, the friction velocity is decreased due to extraction of momentum by the wind turbines, a factor that could decrease scalar fluxes at the surface. In order to study this issue in more detail, a suite of large eddy simulations of an infinite (fully developed) wind turbine array boundary layer, including scalar transport from the ground surface without stratification, is perf...

Coherent structures and the k(-1) spectral behaviour

Here we present unique evidence of a k −1 scaling behaviour in the atmospheric boundary layer and its connection to large scale coherent structures within the boundary layer. Wind lidar measurements were conducted above a lake under cold atmospheric conditions. The large coherent structures could be visually observed over Lake Geneva in Switzerland when cold air met the relatively warm water. Proper orthogonal decomposition of the experimental data acquired with the wind lidar clearly reveals coherent oscillations of both the fluctuating velocity field and the water aerosol field over the surface of the lake. Precise identification of the large coherent structures propagating in the flow allows for detailed analysis of their contribution to the total spectral budget. Additionally, it is shown that the experimental data agree well with recent theoretical predictions.

Controls on the diurnal streamflow cycles in two subbasins of an alpine headwater catchment

In high-altitude alpine catchments, diurnal streamflow cycles are typically dominated by snowmelt or ice melt. Evapotranspiration-induced diurnal streamflow cycles are less observed in these catchments but might happen simultaneously. During a field campaign in the summer 2012 in an alpine catchment in the Swiss Alps (Val Ferret catchment, 20.4 km2, glaciarized area: 2%), we observed a transition in the early season from a snowmelt to an evapotranspiration-induced diurnal streamflow cycle in one of two monitored subbasins. The two different cycles were of comparable amplitudes and the transition happened within a time span of several days. In the second monitored subbasin, we observed an ice melt-dominated diurnal cycle during the entire season due to the presence of a small glacier. Comparisons between ice melt and evapotranspiration cycles showed that the two processes were happening at the same times of day but with a different sign and a different shape. The amplitude of the ice melt cycle decreased exponentially during the season and was larger than the amplitude of the evapotranspiration cycle which was relatively constant during the season. Our study suggests that an evapotranspiration-dominated diurnal streamflow cycle could damp the ice melt-dominated diurnal streamflow cycle. The two types of diurnal streamflow cycles were separated using a method based on the identification of the active riparian area and measurement of evapotranspiration.

Perturbations to the Spatial and Temporal Characteristics of the Diurnally-Varying Atmospheric Boundary Layer Due to an Extensive Wind Farm

The effect of extensive terrestrial wind farms on the spatio-temporal structure of the diurnally-evolving atmospheric boundary layer is explored. High-resolution large-eddy simulations of a realistic diurnal cycle with an embedded wind farm are performed. Simulations are forced by a constant geostrophic velocity with time-varying surface boundary conditions derived from a selected period of the CASES-99 field campaign. Through analysis of the bulk statistics of the flow as a function of height and time, it is shown that extensive wind farms shift the inertial oscillations and the associated nocturnal low-level jet vertically upwards by approximately 200 m; cause a three times stronger stratification between the surface and the rotor-disk region, and as a consequence, delay the formation and growth of the convective boundary layer (CBL) by approximately 2 h. These perturbations are shown to have a direct impact on the potential power output of an extensive wind farm with the displacement of the low-level jet causing lower power output during the night as compared to the day. The low-power regime at night is shown to persist for almost 2 h beyond the morning transition due to the reduced growth of the CBL. It is shown that the wind farm induces a deeper entrainment region with greater entrainment fluxes. Finally, it is found that the diurnally-averaged effective roughness length for wind farms is much lower than the reference value computed theoretically for neutral conditions.

A New Wall Shear Stress Model for Atmospheric Boundary Layer Simulations

AbstractA new wall shear stress model to be used as a wall boundary condition for large-eddy simulations of the atmospheric boundary layer is proposed. The new model computes the wall shear stress and the vertical derivatives of the streamwise velocity component by means of a modified, instantaneous, and local law-of-the-wall formulation. By formulating a correction for the modeled shear stress, using experimental findings of a logarithmic region in the streamwise turbulent fluctuations, the need for a filter is eliminated. This allows one to model the wall shear stress locally, and at the same time accurately recover the correct average value. The proposed model has been applied to both unique high Reynolds number experimental data and a suite of large-eddy simulations, and compared to previous models. It is shown that the proposed model performs equally well or better than the previous filtered models. A nonfiltered model, such as the one proposed, is an essential first step in developing a universal wa...

Large Wind Farms and the Scalar Flux over an Heterogeneously Rough Land Surface

The influence of surface heterogeneities extends vertically within the atmospheric surface layer to the so-called blending height, causing changes in the fluxes of momentum and scalars. Inside this region the turbulence structure cannot be treated as horizontally homogeneous; it is highly dependent on the local surface roughness, the buoyancy and the horizontal scale of heterogeneity. The present study analyzes the change in scalar flux induced by the presence of a large wind farm installed across a heterogeneously rough surface. The change in the internal atmospheric boundary-layer structure due to the large wind farm is decomposed and the change in the overall surface scalar flux is assessed. The equilibrium length scale characteristic of surface roughness transitions is found to be determined by the relative position of the smooth-to-rough transition and the wind turbines. It is shown that the change induced by large wind farms on the scalar flux is of the same order of magnitude as the adjustment they naturally undergo due to surface patchiness.

Carbon monoxide as a tracer of gas transport in snow and other natural porous media

The movement of air in natural porous media is complex and challenging to measure. Yet gas transport has important implications, for instance, for the evolution of the seasonal snow cover and for water vapor transport in soil. A novel in situmulti-sensor measurement system providing high-resolution observation of gas transport in snow is demonstrated. Carbon monoxide was selected as the tracer gas for having essentially the same density as air, low background concentration, low water solubility, and for being detectable to ? 1 ppmv with small, low-cost, low-power sensors. The plume of 1% CO injections 30 cm below the snow surface was monitored using 28 sensors (4 locations, 7 depths). The CO breakthrough curves obtained at distances of 0.5–1 m were in good agreement with a simple analytical advection-diffusion model. The tracer system appears suitable for a wide range of applications in experimental soil science and hydrology addressing moisture transport and evapotranspiration processes.

Dependence of near‐surface similarity scaling on the anisotropy of atmospheric turbulence

Turbulence data from the CASES‐99 field experiment, over comparatively horizontally homogeneous and flat terrain, are separated based on the anisotropy of the Reynolds stress tensor (into isotropic, two‐component axisymmetric and one‐component turbulence) and flux‐variance similarity scaling relations are tested. Results illustrate that different states of anisotropy correspond to different similarity relations, especially under unstable stratification. Experimental data with close to isotropic turbulence match similarity relationships well. On the other hand, very anisotropic turbulence deviates significantly from the traditional scaling relations. We examine in detail the characteristics of these states of anisotropy, identify conditions in which they occur and connect them with different governing parameters. The governing parameters of turbulence anisotropy are shown to be different for stable and unstable stratification, but are able to delineate clearly the conditions in which each of the anisotropy states occurs.

Anisotropy stress invariants of thermally stratified wind turbine array boundary layers using large eddy simulations

In the interest of understanding the dynamics and energy transfer between the atmospheric boundary layer and large wind turbine arrays, a description of the turbulence anisotropy in the wake region of wind turbines is necessary. Flux of high momentum flow into the wind turbine array by anisotropic turbulence is a dominant factor of mean kinetic energy resupply for the wind farm. Under thermal stratification, the behavior of the turbulence field, and the energy flux, is significantly altered. Here, a thermally stratified wind turbine array boundary layer is studied in detail using a combination of Lumley and barycentric maps together with the recently introduced spheroid and color maps. The atmospheric flow is modelled using a large eddy simulation driven by a constant geostrophic wind and a time-varying ground surface temperature, obtained from a selected period of the Cooperative Atmosphere-Surface Exchange Study-99 field experiment. The wind farm is modelled using the traditional actuator-disk with rota...