Effects of Thermal Stability and Incoming Boundary-Layer Flow Characteristics on Wind-Turbine Wakes: A Wind-Tunnel Study

Wind-tunnel experiments were carried out to study turbulence statistics in the wake of a model wind turbine placed in a boundary-layer flow under both neutral and stably stratified conditions. High-resolution velocity and temperature measurements, obtained using a customized triple wire (cross-wire and cold wire) anemometer, were used to characterize the mean velocity, turbulence intensity, turbulent fluxes, and spectra at different locations in the wake. The effect of the wake on the turbulence statistics is found to extend as far as 20 rotor diameters downwind of the turbine. The velocity deficit has a nearly axisymmetric shape, which can be approximated by a Gaussian distribution and a power-law decay with distance. This decay in the near-wake region is found to be faster in the stable case. Turbulence intensity distribution is clearly non-axisymmetric due to the non-uniform distribution of the incoming velocity in the boundary layer. In the neutral case, the maximum turbulence intensity is located above the hub height, around the rotor tip location and at a distance of about 4–5.5 rotor diameters, which are common separations between wind turbines in wind farms. The enhancement of turbulence intensity is associated with strong shear and turbulent kinetic energy production in that region. In the stable case, the stronger shear in the incoming flow leads to a slightly stronger and larger region of enhanced turbulence intensity, which extends between 3 and 6 rotor diameters downwind of the turbine location. Power spectra of the streamwise and vertical velocities show a strong signature of the turbine blade tip vortices at the top tip height up to a distance of about 1–2 rotor diameters. This spectral signature is stronger in the vertical velocity component. At longer downwind distances, tip vortices are not evident and the von Kármán formulation agrees well with the measured velocity spectra.     

The Effect of Wind-Turbine Wakes on Summertime US Midwest Atmospheric Wind Profiles as Observed with Ground-Based Doppler Lidar

We examine the influence of a modern multi-megawatt wind turbine on wind and turbulence profiles three rotor diameters ( $D$ D ) downwind of the turbine. Light detection and ranging (lidar) wind-profile observations were collected during summer 2011 in an operating wind farm in central Iowa at 20-m vertical intervals from 40 to 220 m above the surface. After a calibration period during which two lidars were operated next to each other, one lidar was located approximately $2D$ 2 D directly south of a wind turbine; the other lidar was moved approximately $3D$ 3 D north of the same wind turbine. Data from the two lidars during southerly flow conditions enabled the simultaneous capture of inflow and wake conditions. The inflow wind and turbulence profiles exhibit strong variability with atmospheric stability: daytime profiles are well-mixed with little shear and strong turbulence, while nighttime profiles exhibit minimal turbulence and considerable shear across the rotor disk region and above. Consistent with the observations available from other studies and with wind-tunnel and large-eddy simulation studies, measurable reductions in wake wind-speeds occur at heights spanning the wind turbine rotor (43–117 m), and turbulent quantities increase in the wake. In generalizing these results as a function of inflow wind speed, we find the wind-speed deficit in the wake is largest at hub height or just above, and the maximum deficit occurs when wind speeds are below the rated speed for the turbine. Similarly, the maximum enhancement of turbulence kinetic energy and turbulence intensity occurs at hub height, although observations at the top of the rotor disk do not allow assessment of turbulence in that region. The wind shear below turbine hub height (quantified here with the power-law coefficient) is found to be a useful parameter to identify whether a downwind lidar observes turbine wake or free-flow conditions. These field observations provide data for validating turbine-wake models and wind-tunnel observations, and for guiding assessments of the impacts of wakes on surface turbulent fluxes or surface temperatures downwind of turbines.     

Wind-Turbine Wakes in a Convective Boundary Layer: A Wind-Tunnel Study

Thermal stability changes the properties of the turbulent atmospheric boundary layer, and in turn affects the behaviour of wind-turbine wakes. To better understand the effects of thermal stability on the wind-turbine wake structure, wind-tunnel experiments were carried out with a simulated convective boundary layer (CBL) and a neutral boundary layer. The CBL was generated by cooling the airflow to 12–15 °C and heating up the test section floor to 73–75 °C. The freestream wind speed was set at about 2.5 m s^?1, resulting in a bulk Richardson number of ?0.13. The wake of a horizontal-axis 3-blade wind-turbine model, whose height was within the lowest one third of the boundary layer, was studied using stereoscopic particle image velocimetry (S-PIV) and triple-wire (x-wire/cold-wire) anemometry. Data acquired with the S-PIV were analyzed to characterize the highly three-dimensional turbulent flow in the near wake (0.2–3.2 rotor diameters) as well as to visualize the shedding of tip vortices. Profiles of the mean flow, turbulence intensity, and turbulent momentum and heat fluxes were measured with the triple-wire anemometer at downwind locations from 2–20 rotor diameters in the centre plane of the wake. In comparison with the wake of the same wind turbine in a neutral boundary layer, a smaller velocity deficit (about 15 % at the wake centre) is observed in the CBL, where an enhanced radial momentum transport leads to a more rapid momentum recovery, particularly in the lower part of the wake. The velocity deficit at the wake centre decays following a power law regardless of the thermal stability. While the peak turbulence intensity (and the maximum added turbulence) occurs at the top-tip height at a downwind distance of about three rotor diameters in both cases, the magnitude is about 20 % higher in the CBL than in the neutral boundary layer. Correspondingly, the turbulent heat flux is also enhanced by approximately 25 % in the lower part of the wake, compared to that in the undisturbed CBL inflow. This study represents the first controlled wind-tunnel experiment to study the effects of the CBL on wind-turbine wakes. The results on decreased velocity deficit and increased turbulence in wind-turbine wakes associated with atmospheric thermal stability are important to be taken into account in the design of wind farms, in order to reduce the impact of wakes on power output and fatigue loads on downwind wind turbines.     

Observing and Simulating Wind-Turbine Wakes During the Evening Transition

Wind-turbine-wake evolution during the evening transition introduces variability to wind-farm power production at a time of day typically characterized by high electricity demand. During the evening transition, the atmosphere evolves from an unstable to a stable regime, and vertical stratification of the wind profile develops as the residual planetary boundary layer decouples from the surface layer. The evolution of wind-turbine wakes during the evening transition is examined from two perspectives: wake observations from single turbines, and simulations of multiple turbine wakes using the mesoscale Weather Research and Forecasting (WRF) model. Throughout the evening transition, the wake’s wind-speed deficit and turbulence enhancement are confined within the rotor layer when the atmospheric stability changes from unstable to stable. The height variations of maximum upwind-downwind differences of wind speed and turbulence intensity gradually decrease during the evening transition. After verifying the WRF-model-simulated upwind wind speed, wind direction and turbulent kinetic energy profiles with observations, the wind-farm-scale wake evolution during the evening transition is investigated using the WRF-model wind-farm parametrization scheme. As the evening progresses, due to the presence of the wind farm, the modelled hub-height wind-speed deficit monotonically increases, the relative turbulence enhancement at hub height grows by 50%, and the downwind surface sensible heat flux increases, reducing surface cooling. Overall, the intensifying wakes from upwind turbines respond to the evolving atmospheric boundary layer during the evening transition, and undermine the power production of downwind turbines in the evening.     

Impact of the Diurnal Cycle of the Atmospheric Boundary Layer on Wind-Turbine Wakes: A Numerical Modelling Study

The wake characteristics of a wind turbine for different regimes occurring throughout the diurnal cycle are investigated systematically by means of large-eddy simulation. Idealized diurnal cycle simulations of the atmospheric boundary layer are performed with the geophysical flow solver EULAG over both homogeneous and heterogeneous terrain. Under homogeneous conditions, the diurnal cycle significantly affects the low-level wind shear and atmospheric turbulence. A strong vertical wind shear and veering with height occur in the nocturnal stable boundary layer and in the morning boundary layer, whereas atmospheric turbulence is much larger in the convective boundary layer and in the evening boundary layer. The increased shear under heterogeneous conditions changes these wind characteristics, counteracting the formation of the night-time Ekman spiral. The convective, stable, evening, and morning regimes of the atmospheric boundary layer over a homogeneous surface as well as the convective and stable regimes over a heterogeneous surface are used to study the flow in a wind-turbine wake. Synchronized turbulent inflow data from the idealized atmospheric boundary-layer simulations with periodic horizontal boundary conditions are applied to the wind-turbine simulations with open streamwise boundary conditions. The resulting wake is strongly influenced by the stability of the atmosphere. In both cases, the flow in the wake recovers more rapidly under convective conditions during the day than under stable conditions at night. The simulated wakes produced for the night-time situation completely differ between heterogeneous and homogeneous surface conditions. The wake characteristics of the transitional periods are influenced by the flow regime prior to the transition. Furthermore, there are different wake deflections over the height of the rotor, which reflect the incoming wind direction.     

Fractional Flow Speed-Up from Porous Windbreaks for Enhanced Wind-Turbine Power

The potential for porous windbreaks to enhance wind-turbine power production is studied using linearized theory and wind-tunnel experiments. Results suggest that windbreaks have the potential to substantially increase power production, while lowering mean shear, and leading to negligible changes in turbulence intensity. The fractional increase in turbine power output is found to vary roughly linearly with windbreak height, where a windbreak 10% the height of the turbine hub increases power by around 10%. Wind-tunnel experiments with a windbreak imposed beneath a turbulent boundary layer show the linearized predictions to be in good agreement with particle-image-velocimetry data. Power measurements from a model turbine further corroborate predictions in power increase. Moreover, the wake of the windbreak showed a significant interaction with the turbine wake, which may inform windbreak use in large wind farms. Power measurements from a second turbine downwind of the first with its own windbreak show that the net effect for multiple turbines is dependent on windbreak height.     

Wind profile development above a locally adjusted sea surface

Results are presented from a numerical experiment of wind and shear stress profile development away from a shore line; the water surface is assumed to obey the Charnock-Ellison relation between surface roughness and friction velocity. In typical cases the upwind, land surface is rough relative to the sea and the velocity and shear stress results are qualitatively similar to those for flows from relatively rough to relatively smooth solid surfaces. In the present case, however, the downwind surface roughness and friction velocity vary with position and we find that wind profile development may play a significant role in the relationship between sea surface roughness and fetch.     

Impact of Neutral Boundary-Layer Turbulence on Wind-Turbine Wakes: A Numerical Modelling Study

The wake characteristics of a wind turbine in a turbulent boundary layer under neutral stratification are investigated systematically by means of large-eddy simulations. A methodology to maintain the turbulence of the background flow for simulations with open horizontal boundaries, without the necessity of the permanent import of turbulence data from a precursor simulation, was implemented in the geophysical flow solver EULAG. These requirements are fulfilled by applying the spectral energy distribution of a neutral boundary layer in the wind-turbine simulations. A detailed analysis of the wake response towards different turbulence levels of the background flow results in a more rapid recovery of the wake for a higher level of turbulence. A modified version of the Rankine–Froude actuator disc model and the blade element momentum method are tested as wind-turbine parametrizations resulting in a strong dependence of the near-wake wind field on the parametrization, whereas the far-wake flow is fairly insensitive to it. The wake characteristics are influenced by the two considered airfoils in the blade element momentum method up to a streamwise distance of 14D (D = rotor diameter). In addition, the swirl induced by the rotation has an impact on the velocity field of the wind turbine even in the far wake. Further, a wake response study reveals a considerable effect of different subgrid-scale closure models on the streamwise turbulent intensity.     

Turbulent Flow Properties Around a Staggered Wind Farm

The fundamental properties of turbulent flow around a perfectly staggered wind farm are investigated in a wind tunnel. The wind farm consisted of a series of 10 rows by 2–3 columns of miniature wind turbines spaced 5 and 4 rotor diameters in the streamwise and spanwise directions respectively. It was placed in a boundary-layer flow developed over a smooth surface under thermally neutral conditions. Cross-wire anemometry was used to obtain high resolution measurements of streamwise and vertical velocity components at various locations within and above the wind farm. The results show that the staggered configuration is more efficient in terms of momentum transfer from the background flow to the turbines compared to the case of an aligned wind turbine array under similar turbine separations in the streamwise and spanwise directions. This leads to improved power output of the overall wind farm. A simplified analysis suggests that the difference in power output between the two configurations is on the order of 10%. The maximum levels of turbulence intensity in the staggered wind farm were found to be very similar to that observed in the wake of a single wind turbine, differing substantially with that observed in an aligned configuration with similar spacing. The dramatic changes in momentum and turbulence characteristics in the two configurations show the importance of turbine layout in engineering design. Lateral homogenization of the turbulence statistics above the wind farm allows for the development of simple parametrizations for the adjustment of flow properties, similar to the case of a surface roughness transition. The development of an internal boundary layer was observed at the upper edge of the wind farm within which the flow statistics are affected by the superposition of the ambient flow and the flow disturbance induced by the wind turbines. The adjustment of the flow in this layer is much slower in the staggered situation (with respect to its aligned counterpart), implying a change in the momentum/power available at turbine locations. Additionally, power spectra of the streamwise and vertical velocity components indicate that the signature of each turbine-tip vortex structure persists to locations deep within the wind farm.     

风电场流场特性及风机布局数值模拟研究

针对风电场流场特性研究对风力机工作性能提高的重要意义,采用计算流体力学(CFD)方法在单机风力机模拟验证的基础上,对某风电场单风力机和三种布局条件下的风电场流场特性进行了数值模拟研究.考察了不同布局条件下风电场速度、叶轮表面压力以及湍流涡的分布特性.结果表明:叶轮后方尾流效应明显,速度损失随着相对距离的增加而逐渐减小,...     

An extended Miles' theory for wave generation by wind

Miles' inviscid theory of surface wave generation by wind is (a) modified by replacing the logarithmic shear velocity profile with one which applies right down to the wave surface and which exhibits an explicit dependence on the roughness of the surface, and (b) extended to include the effects of the interaction of wave with air flow turbulence by considering the wave-modified mean flow as the mean of the actual turbulent air flow over water waves and using this in a mixing-length model.The surface pressure is shown to depend significantly on the flow conditions being aerodynamically smooth or rough. Its component in phase with the surface elevation is practically unaffected by the wave-turbulence interaction. However, such interaction tends to increase the rate of energy input ß from wind to waves travelling in the same direction, e.g., the increase is 2gk ² for aerodynamically rough flow, where gk is the Von Karman constant. It also provides damping of waves in an adverse wind which can be about 10% of the growth rate in a favourable wind.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

A Numerical Study of Wind-Turbine Wakes for Three Atmospheric Stability Conditions

The effects of atmospheric stability on wind-turbine wakes are studied via large-eddy simulations. Three stability conditions are considered: stable, neutral, and unstable, with the same geostrophic wind speed aloft and the same Coriolis frequency. Both a single 5-MW turbine and a wind farm of five turbines are studied. The single-turbine wake is strongly correlated with stability, in terms of velocity deficit, turbulence kinetic energy (TKE) and temperature distribution. Because of the Coriolis effect, the wake shape deviates from a Gaussian distribution. For the wind-farm simulations, the separation of the core region and outer region is clear for the stable and neutral cases, but less distinct for the unstable case. The unstable case exhibits strong horizontal variations in wind speed. Local accelerations such as related to aisle jets are also observed, whose features depend on stability. The added TKE in the wind farm increases with stability. The highest power extraction and lowest power deficit are observed for the unstable case.     

基于激光雷达测风数据的复杂地形下风力机尾流特性研究

以张家口某处风电场为实验场地,采用两台多普勒激光雷达(Wind3D 6000和WindMast WP350)分别测量风力机的尾流和来流风速,对全尾流、半尾流和独立尾流3种工况进行研究。结果表明,3种工况下随着尾流发展尾流宽度均不断变大,而尾流深度和速度衰减则不断减小;全尾流和半尾流工况中,上游风力机的存在会增加下游风力机尾流宽度,且全尾流比半尾流的增加量大;全尾流和半尾流工况中,相较于上游风力机,下游风力机尾流深度和速度衰减均较小。     

Simulation of Turbulent Flow Inside and Above Wind Farms: Model Validation and Layout Effects

A recently-developed large-eddy simulation framework is validated and used to investigate turbulent flow within and above wind farms under neutral conditions. Two different layouts are considered, consisting of thirty wind turbines occupying the same total area and arranged in aligned and staggered configurations, respectively. The subgrid-scale (SGS) turbulent stress is parametrized using a tuning-free Lagrangian scale-dependent dynamic SGS model. The turbine-induced forces are modelled using two types of actuator-disk models: (a) the ‘standard’ actuator-disk model (ADM-NR), which calculates only the thrust force based on one-dimensional momentum theory and distributes it uniformly over the rotor area; and (b) the actuator-disk model with rotation (ADM-R), which uses blade-element momentum theory to calculate the lift and drag forces (that produce both thrust and rotation), and distributes them over the rotor disk based on the local blade and flow characteristics. Validation is performed by comparing simulation results with turbulence measurements collected with hot-wire anemometry inside and above an aligned model wind farm placed in a boundary-layer wind tunnel. In general, the ADM-R model yields improved predictions compared with the ADM-NR in the wakes of all the wind turbines, where including turbine-induced flow rotation and accounting for the non-uniformity of the turbine-induced forces in the ADM-R appear to be important. Another advantage of the ADM-R model is that, unlike the ADM-NR, it does not require a priori specification of the thrust coefficient (which varies within a wind farm). Finally, comparison of simulations of flow through both aligned and staggered wind farms shows important effects of farm layout on the flow structure and wind-turbine performance. For the limited-size wind farms considered in this study, the lateral interaction between cumulated wakes is stronger in the staggered case, which results in a farm wake that is more homogeneous in the spanwise direction, thus resembling more an internal boundary layer. Inside the staggered farm, the relatively longer separation between consecutive downwind turbines allows the wakes to recover more, exposing the turbines to higher local wind speeds (leading to higher turbine efficiency) and lower turbulence intensity levels (leading to lower fatigue loads), compared with the aligned farm. Above the wind farms, the area-averaged velocity profile is found to be logarithmic, with an effective wind-farm aerodynamic roughness that is larger for the staggered case.     

基于测风塔实测台风威马逊登陆过程的强风特性分析

利用广东省徐闻县西连镇90 m测风塔在1409号超强台风威马逊登陆期间获取的具备完整的台风代表性的观测数据以及处于台风外围的广东省茂名市博贺镇100 m测风塔的观测数据,对台风威马逊的近地层强风特性进行了分析,西连测风塔结果表明:风速时程曲线呈明显的"M"型分布特征,台风中心经过测风塔前后,风向沿逆时针方向大幅偏转约170°。风速随高度增加而增大,风速廓线较好地符合对数和幂指数律;台风过境前后,各强风区的风速廓线幂指数和粗糙长度呈先减小后增大的特点;粗糙陆地下垫面的风速廓线幂指数和粗糙长度较大。湍流强度和阵风系数在前外围强风区或后外围强风区较大,在前眼壁强风区或后眼壁强风区较小,湍流强度和阵风系数随高度增加而减小,基本符合指数为负值的幂指数律;粗糙下垫面对湍流强度和阵风系数有增大的作用。外围强风区和眼壁强风区的10 min风向变率变化较为平稳,而在眼区变动较为剧烈,在眼区,当风速达到最低值或次低值时,10 min风向变率幅值达到最大值。博贺测风塔结果表明其总体上与西连测风塔台风前外围和前眼壁强风区的情形相似。     

Mean horizontal wind profiles measured in the atmospheric boundary layer about a simulated block building

Instrumented wind towers are used to measure the three components of wind about a simulated block building. The mean horizontal wind profiles over the building are compared with wind profiles measured in the absence of the building and the wind speed deficit in the wake of the building is correlated.Horizontal mean wind speeds measured in the natural atmospheric boundary layer with and without the presence of a simulated building show excellent reproducibility and agreement with fundamental concepts of fluid mechanics. The data possess all the characteristic features reported from wind-tunnel studies of building flows. In the present study the turbulence intensity is of the order of 20% in the undisturbed flow whereas the free stream turbulence intensity of wind-tunnel studies is generally not more than 5%. The effect of smaller averaging periods and the structure of the turbulence will be reported at a later time.The velocity profiles measured in the undisturbed flow zones support the representation of a neutrally stable atmospheric boundary layer with a logarithmic wind profile.     

A wind-tunnel boundary-layer study of the effects of a surface roughness change: Rough to smooth

The changes imposed on mean velocities and turbulence statistics in the lower atmosphere by an abrupt change in surface roughness, from very rough to smooth, were modelled in a wind tunnel. The influence of a change in the effective surface level, which often accompanies such a variation in surface roughness, was also studied. A deep, turbulent flow was generated upstream of the change, which had a logarithmic mean velocity profile and constant shear-stress for approximately 200 mm above the floor, except for a region near the surface which was influenced by the three-dimensional nature of the random rough surface.When the surface roughness change coincided with a change in surface level, the downstream flow close to the surface was in the wake of the upstream roughness elements, and measured Reynolds shear-stress values were lower than those obtained when the downstream surface was raised. Otherwise, the influence of a change in surface level was small.In all cases, Reynolds shear-stress varied approximately linearly with height in the lower two-thirds of the internal layer and no constant stress region was apparent near the surface, even 2 m downstream of the roughness change. When the roughness change was not accompanied by a change in level, Reynolds shear-stress values extrapolated to the surface agreed well with surface shear-stress inferred from the law of the wall.Changes in mean squares of vertical and lateral velocity fluctuations and in integral time scales, as the flow passed downstream of the roughness change, were surprisingly small.     

Velocity and Surface Shear Stress Distributions Behind a Rough-to-Smooth Surface Transition: A Simple New Model

A simple new model is proposed to predict the distribution of wind velocity and surface shear stress downwind of a rough-to-smooth surface transition. The wind velocity is estimated as a weighted average between two limiting logarithmic profiles: the first log law, which is recovered above the internal boundary-layer height, corresponds to the upwind velocity profile; the second log law is adjusted to the downwind aerodynamic roughness and local surface shear stress, and it is recovered near the surface, in the equilibrium sublayer. The proposed non-linear form of the weighting factor is equal to ln(z/z ₀₁)/ln(δ _i /z ₀₁), where z, δ _i and z ₀₁ are the elevation of the prediction location, the internal boundary-layer height at that downwind distance, and the upwind surface roughness, respectively. Unlike other simple analytical models, the new model does not rely on the assumption of a constant or linear distribution for the turbulent shear stress within the internal boundary layer. The performance of the new model is tested with wind-tunnel measurements and also with the field data of Bradley. Compared with other existing analytical models, the proposed model shows improved predictions of both surface shear stress and velocity distributions at different positions downwind of the transition.