Simulating the aerodynamic performance and wake dynamics of a vertical‐axis wind turbine

The accurate prediction of the aerodynamics and performance of vertical‐axis wind turbines is essential if their design is to be improved but poses a significant challenge to numerical simulation tools. The cyclic motion of the blades induces large variations in the angle of attack of the blades that can manifest as dynamic stall. In addition, predicting the interaction between the blades and the wake developed by the rotor requires a high‐fidelity representation of the vortical structures within the flow field in which the turbine operates. The aerodynamic performance and wake dynamics of a Darrieus‐type vertical‐axis wind turbine consisting of two straight blades is simulated using Brown's Vorticity Transport Model. The predicted variation with azimuth of the normal and tangential force on the turbine blades compares well with experimental measurements. The interaction between the blades and the vortices that are shed and trailed in previous revolutions of the turbine is shown to have a significant effect on the distribution of aerodynamic loading on the blades. Furthermore, it is suggested that the disagreement between experimental and numerical data that has been presented in previous studies arises because the blade–vortex interactions on the rotor were not modelled with sufficient fidelity. Copyright © 2010 John Wiley & Sons, Ltd.     

Sound propagation through the wake flow of a hilltop wind turbine—A numerical study

The sound propagation from a wind turbine situated on the top of a hill into the downwind domain is studied by numerical simulations for 13 cases with varying hill geometry and inflow conditions. The influence of the hill on the atmospheric flow and the wake due to the rotor are simulated by precursory large‐eddy simulations. In addition to the combined consideration of hill and turbine wake effects, these effects are also separately evaluated. The results show that placing the turbine on top of a hill leads to slightly lower sound levels on the downwind plane, although the wake alone supports downward refraction and tends to increase the sound impact near the ground at greater distance. Variations of the hill geometry and the inflow conditions do not have significant effects on the near‐ground sound levels in the downwind domain.     

A vortex panel model for the simulation of the wake flow past a vertical axis wind turbine in dynamic stall

A 2D vortex panel model with a viscous boundary layer formulation has been developed for the numerical simulation of a vertical axis wind turbine (VAWT), including the operation in dynamic stall. The model uses the ‘double wake’ concept to reproduce the main features of the unsteady separated flow, including the formation and shedding of strong vortical structures and the wake–blade interaction. The potential flow equations are solved together with the integral boundary layer equations by using a semi‐inverse iterative algorithm. A new criterion for the reattachment of the boundary layer during the downstroke of a dynamically stalled aerofoil is implemented. The model has been validated against experimental data of steady aerofoils and pitching aerofoils in dynamic stall at high and low Reynolds numbers (Re = 1.5 × 10⁶ and Re = 5 × 10⁴). For the low Reynolds number case, time‐resolved 2D particle image velocimetry (PIV) measurements have been performed on a pitching NACA 0012 aerofoil in dynamic stall. The PIV vorticity fields past the oscillating aerofoil are used to test the model capability of capturing the formation, growth and release of the strong leading edge vortex that characterizes the dynamic stall. Furthermore, the forces extracted from the PIV velocity fields are compared with the predicted ones for a quantitative validation of the model. Finally, the model is applied to the computation of the wake flow past a VAWT in dynamic stall; the predicted vorticity fields and forces are in good agreement with phase‐locked PIV data and CFD‐DES available in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of inflow conditions on turbine loading and wake structures predicted by large eddy simulations using exact geometry

A large eddy simulation was performed on an National Renewable Energy Laboratory (NREL) phase VI wind turbine (10 m diameter), using the exact blade geometry, to determine the influence of different inflow conditions on the aerodynamic loadings and the near wake characteristics. The effects of the three inflow conditions, uniform inflow, linear wind shear and linear wind shear with turbulence, are investigated. Wind shear causes periodic variations in power and aerodynamic loading with an additional force component exerted along the lateral direction. Significant separation occurs in the high wind region on the suction side of the blades, resulting in unstable loading in off‐design inflow conditions. Because of the shear effect between the near‐blade tip vortex and ambient flow, the strong vortex core in the helical structure dissipates and transforms into a continuous vorticity sheet when x/D > 1.5. The combination of inflow turbulence and wind shear enhances the turbulence generation mechanism in the near wake, where energy is withdrawn from large wake structures and converted into energy of small‐scale structures. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of a LES technique to characterize the wake deflection of a wind turbine in yaw

When a wind turbine works in yaw, the wake intensity and the power production of the turbine become slightly smaller and a deflection of the wake is induced. Therefore, a good understanding of this effect would allow an active control of the yaw angle of upstream turbines to steer the wake away from downstream machines, reducing its effect on them. In wind farms where interaction between turbines is significant, it is of interest to maximize the power output from the wind farm as a whole and to reduce fatigue loads on downstream turbines due to the increase of turbulence intensity in wakes. A large eddy simulation model with particular wind boundary conditions has been used recently to simulate and characterize the turbulence generated by the presence of a wind turbine and its evolution downstream the machine. The simplified turbine is placed within an environment in which relevant flow properties like wind speed profile, turbulence intensity and the anisotropy of turbulence are found to be similar to the ones of the neutral atmosphere. In this work, the model is used to characterize the wake deflection for a range of yaw angles and thrust coefficients of the turbine. The results are compared with experimental data obtained by other authors with a particle image velocimetry technique from wind tunnel experiments. Also, a comparison with simple analytical correlations is carried out. Copyright © 2009 John Wiley & Sons, Ltd.     

Impact of different forest densities on atmospheric boundary‐layer development and wind‐turbine wake

The aim of this work is to investigate the atmospheric boundary‐layer (ABL) flow and the wind turbine wake over forests with varying leaf area densities (LAD). The forest LAD profile used in this study is based on a real forest site, Ryningsnäs, located in Sweden. The reference turbine used to model the wake is a well‐documented 5‐MW turbine, which is implemented in the simulations using an actuator line model (ALM). All simulations are carried out with openFOAM using the Reynolds averaged Navier‐Stokes (RANS) approach. Twelve forest cases with leaf area index (LAI) ranging from 0.42 to 8.5 are considered. Results show that the mean velocity decreases with increasing LAI within the forest canopy, but increases with LAI above the hub height. Meanwhile, the turbulent kinetic energy (TKE) varies nonmonotonically with forest density. The TKE increases with forest density and reaches to its maximum at an average LAI of 1.70, afterwards, it decreases gradually as the density increases. It is also observed that the forest density has a clear role in the wake development and recovery. Comparisons between no‐forest and forest cases show that the forest characteristics help in damping the added turbulence from the turbine. As a consequence, the forest with the highest upstream turbulence has the shortest wake downstream of the turbine.     

Quantification of power losses due to wind turbine wake interactions through SCADA,meteorological and wind LiDAR data

Power production of an onshore wind farm is investigated through supervisory control and data acquisition data, while the wind field is monitored through scanning light detection and ranging measurements and meteorological data acquired from a met‐tower located in proximity to the turbine array. The power production of each turbine is analysed as functions of the operating region of the power curve, wind direction and atmospheric stability. Five different methods are used to estimate the potential wind power as a function of time, enabling an estimation of power losses connected with wake interactions. The most robust method from a statistical standpoint is that based on the evaluation of a reference wind velocity at hub height and experimental mean power curves calculated for each turbine and different atmospheric stability regimes. The synergistic analysis of these various datasets shows that power losses are significant for wind velocities higher than cut‐in wind speed and lower than rated wind speed of the turbines. Furthermore, power losses are larger under stable atmospheric conditions than for convective regimes, which is a consequence of the stability‐driven variability in wake evolution. Light detection and ranging measurements confirm that wind turbine wakes recover faster under convective regimes, thus alleviating detrimental effects due to wake interactions. For the wind farm under examination, power loss due to wake shadowing effects is estimated to be about 4% and 2% of the total power production when operating under stable and convective conditions, respectively. However, cases with power losses about 60‐80% of the potential power are systematically observed for specific wind turbines and wind directions. Copyright © 2017 John Wiley & Sons, Ltd.     

An experimental study on the effects of winglets on the tip vortex interaction in the near wake of a model wind turbine

An experimental study of the near wake up to four rotor diameters behind a model wind turbine rotor with two different wing tip configurations is performed. A straight‐cut wing tip and a downstream‐facing winglet shape are compared on the same two‐bladed rotor operated at its design tip speed ratio. Phase‐averaged measurements of the velocity vector are synchronized with the rotor position, visualizing the downstream location of tip vortex interaction for the two blade tip configurations. The mean streamwise velocity is found not to be strongly affected by the presence of winglet tip extensions, suggesting an insignificant effect of winglets on the time‐averaged inflow conditions of a possible downstream wind turbine. An analysis of the phase‐averaged vorticity, however, reveals a significantly earlier tip vortex interaction and breakup for the wingletted rotor. In contradistinction, the tip vortices formed behind the reference configuration are assessed to be more stable and start merging into larger turbulent structures significantly further downstream. These results indicate that an optimized winglet design can not only contribute to a higher energy extraction in a rotor's tip region but also can positively affect the wake's mean kinetic energy recovery by stimulating a faster tip vortex interaction.     

Evaluation of the effects of turbulence model enhancements on wind turbine wake predictions

The modelling of wind turbine wakes is investigated in this paper using a Navier–Stokes solver employing the k–ω turbulence model appropriately modified for atmospheric flows. It is common knowledge that even single‐wind turbine wake predictions with computational fluid dynamic methods underestimate the near wake deficit, directly contributing to the overestimation of the power of the downstream turbines. For a single‐wind turbine, alternative modelling enhancements under neutral and stable atmospheric conditions are tested in this paper to account for and eventually correct the turbulence overestimation that is responsible for the faster flow recovery that appears in the numerical predictions. Their effect on the power predictions is evaluated with comparison with existing wake measurements. A second issue addressed in this paper concerns multi‐wake predictions in wind farms, where the estimation of the reference wind speed that is required for the thrust calculation of a turbine located in the wake(s) of other turbines is not obvious. This is overcome by utilizing an induction factor‐based concept: According to it, the definition of the induction factor and its relationship with the thrust coefficient are employed to provide an average wind speed value across the rotor disk for the estimation of the axial force. Application is made on the case of five wind turbines in a row. Copyright © 2010 John Wiley & Sons, Ltd.     

Near wake Reynolds‐averaged Navier–Stokes predictions of the wake behind the MEXICO rotor in axial and yawed flow conditions

In the present paper, Reynolds‐averaged Navier–Stokes predictions of the flow field around the MEXICO rotor in yawed conditions are compared with measurements. The paper illustrates the high degree of qualitative and quantitative agreement that can be obtained for this highly unsteady flow situation, by comparing measured and computed velocity profiles for all three Cartesian velocity components along four axial transects and several radial transects. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of the near‐wake between actuator‐line simulations and a simplified vortex model of a horizontal‐axis wind turbine

The flow around an isolated horizontal‐axis wind turbine is estimated by means of a new vortex code based on the Biot–Savart law with constant circulation along the blades. The results have been compared with numerical simulations where the wind turbine blades are replaced with actuator lines. Two different wind turbines have been simulated: one with constant circulation along the blades, to replicate the vortex method approximations, and the other with a realistic circulation distribution, to compare the outcomes of the vortex model with real operative wind‐turbine conditions (Tjæreborg wind turbine). The vortex model matched the numerical simulation of the turbine with constant blade circulation in terms of the near‐wake structure and local forces along the blade. The results from the Tjæreborg turbine case showed some discrepancies between the two approaches, but overall, the agreement is qualitatively good, validating the analytical method for more general conditions. The present results show that a simple vortex code is able to provide an estimation of the flow around the wind turbine similar to the actuator‐line approach but with a negligible computational effort. Copyright © 2015 John Wiley & Sons, Ltd.     

雾霾之后的反思

新年过后的第二个周末,浓重的雾霾已在全国多个城市肆虐,这让人们的心情变得糟糕。数据显示,截至1月13日零时,全国有33个城市的部分监测点PM2.5浓度超过300微克/立方米,个别城市出现PM2.5"爆表",比如北京的PM2.5浓度最高达到950微克/立方米。环保专家称,如此严重的空气质量污染,可以说已近人类所能承受的极限。于是,人们看到了政府有关方面发布的紧急预案,比如通知市民减少户外活动,要求学校停止户外体育锻炼,这体现了政府的责任意识。但我们是满足于制定完美的灾情应对预案,还是谋求从根本上消除灾难?     

Investigation of terrain and wake effects on the performance of wind farms in complex terrain using numerical and experimental data

In this work, the combination effects on wind turbine performances of wakes and terrain‐driven flow are investigated. The test case is a subcluster of four turbines from a wind farm sited in southern Italy in a very complex terrain. The layout, the inter‐turbine distance and the wind rose result in a challenging performance scenery. The subcluster is analyzed, when the wind blows from the west, through computational fluid dynamics numerical simulations and experimental supervisory control and data acquisition data mining. Two wind intensity regimes and several simulation setups are employed. It is shown that the main effect of the terrain is the northward distortion of the wake of the upstream turbine. This explains the non‐trivial yawing patterns of the cluster and the fact that the wake line affects the overall performances of the subcluster less than it would do in flat terrain. It is further shown that the presence of the rest of the subcluster in operation southward deviates the wake line of the upstream turbine. The dependency on wind intensity of these directional distortions allows to estimate the relative importance of wakes and terrain‐driven flow. A bijective feedback between models and data is established and a convincing framework is constructed, for separating and assessing the effect of the terrain and of the single and multiple wake. Copyright © 2017 John Wiley & Sons, Ltd.     

Wake modelling of a VAWT in dynamic stall: impact on the prediction of flow and induction fields

Dynamic stall is a relevant phenomenon in the design and operation of a vertical axis wind turbine (VAWT) as it impacts loading, control and wake dynamics. Although streamtube models and single‐wake vortex models are commonly used for VAWT simulation, they either do not explicitly simulate the distribution of vorticity in the wake (streamtube models) or simplify it into a single‐wake release point (single‐wake vortex models). This can lead to inaccurate predictions of the vorticity distribution and wake dynamics, and therefore of the induction field, rotor loading and wake development, including wake mixing and re‐energizing. In this work, we use a double‐wake panel model developed for the simulation of dynamic stall in a VAWT to analyse (i) what is the flow field in dynamic stall, including the induction field, (ii) what is the error due to assuming a simplified wake, in both vorticity distribution and induction and (iii) how an incorrect simulation of the vorticity distribution can affect the prediction of the dynamics of the near and far wake. The results demonstrate that for mild separation (tip speed ratio λ≥3), single‐wake models can produce acceptable results. However, for lower tip speed ratios (λ < 3), the inaccuracy in the prediction of loads, induction field and vorticity distribution becomes significant because of an inadequate representation of the wake dynamics. These results imply that using lower order models can lead to inaccurate estimations of loads, performance and power control requirements at low tip speed ratios. Copyright © 2014 John Wiley & Sons, Ltd.     

Analytical and numerical fluid–structure interaction study of a microscale piezoelectric wind energy harvester

In this paper, an analytical approach and two numerical models have been developed to study an energy‐harvesting device for micropower generation. This device uses wind energy to oscillate a cantilevered beam attached to a piezoelectric layer for generating electric energy output. The analytical approach and the first numerical model consider the fluid–structure interaction phenomenon in the harvester performance. The equations governing beam oscillations and airflow have been coupled to a set of four differential equations in the analytical approach. This set of equations has been solved to determine the beam deflection and the air pressure variation with time. The numerical methods have been conducted by employing a commercial software. The results of the analytical method and the first numerical model have been compared in different working conditions, and their credibility has been discussed. In the second numerical model, the electromechanical performance of the piezoelectric material has also been incorporated in the harvester device analysis. This model has been verified against experimental data for the output voltage and power of the device available in the literature. Finally, the effect of different geometrical parameters has been studied on the harvester performance, and suggestions have been made to improve the harvester efficiency.     

Predictions of turbulent flow and heat transfer in gas–droplets flow downstream of a sudden pipe expansion

Droplets-laden turbulent flow downstream of a sudden pipe expansion has been investigated by using Euler/Euler two-fluid model for the gaseous and dispersed phases. Significant increase of heat transfer in separated flow at the adding of evaporating droplets has been demonstrated (more than 2 times compare with one-phase flow at the value of mass concentration of droplets M_L1  0.05). Addition of dispersed phase to the turbulent gas flow results in insignificant increase of the reattachment length. Low-inertia droplets (d₁  50 μm) are well entrained into the circulation flow and present over the whole pipe section. Large particles (d₁ ≈ 100 μm) go through the shear layer not getting into the detached area. Comparison with experimental data on separated gas–droplets flows behind the plane backward-facing step has been carried out.     

Characterization of the unsteady flow in the nacelle region of a modern wind turbine

A three‐dimensional Navier–Stokes solver has been used to investigate the flow in the nacelle region of a wind turbine where anemometers are typically placed to measure the flow speed and the turbine yaw angle. A 500 kW turbine was modelled with rotor and nacelle geometry in order to capture the complex separated flow in the blade root region of the rotor. A number of steady state and unsteady simulations were carried out for wind speeds ranging from 6 m s^?1 to 16 m s^?1 as well as two yaw and tilt angles. The flow in the nacelle region was found to be highly unsteady, dominated by unsteady vortex shedding from the cylindrical part of the blades, which interacted with the root vortices from each blade, generating high‐velocity gradients. As a consequence, the nacelle wind speed and the nacelle flow angle were found to vary significantly with the height above the nacelle surface. The nacelle anemometry showed significant dependence on both yaw and tilt angles with yaw errors of up to 10 degrees when operating in a tilted inflow. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of rotor–tower interaction,tilt angle,and yaw misalignment on the aeroelasticity of a large horizontal axis wind turbine with composite blades

This paper presents a detailed analysis of the rotor–tower interaction and the effects of the rotor's tilt angle and yaw misalignment on a large horizontal axis wind turbine. A high‐fidelity aeroelastic model is employed, coupling computational fluid dynamics (CFD) and structural mechanics (CSM). The wind velocity stratification induced by the atmospheric boundary layer (ABL) is modeled. On the CSM side, the complex composite structure of each blade is accurately modeled using shell elements. The rotor–tower interaction is analyzed by comparing results of a rotor‐only simulation and a full‐machine simulation, observing a sudden drop in loads, deformations, and power production of each blade, when passing in front of the tower. Subsequently, a tilt angle is introduced on the rotor, and its effect on blade displacements, loads, and performance is studied, representing a novelty with respect to the available literature. The tilt angle leads to a different contribution of gravity to the blade deformations, sensibly affecting the stresses in the composite material. Lastly, a yaw misalignment is introduced with respect to the incoming wind, and the resulting changes in the blade solicitations are analyzed. In particular, a reduction of the blade axial displacement amplitude during each revolution is observed.     

Validation of four LES and a vortex model against stereo-PIV measurements in the near wake of an actuator disc and a wind turbine

In this paper we report the results of a workshop organised by the Delft University of Technology in 2014, aiming at the comparison between different state-of-the-art numerical models for the simulation of wind turbine wakes. The chosen benchmark case is a wind tunnel measurement, where stereoscopic Particle Image Velocimetry was employed to obtain the velocity field and turbulence statistics in the near wake of a two-bladed wind turbine model and of a porous disc, which mimics the numerical actuator used in the simulations. Researchers have been invited to simulate the experimental case based on the disc drag coefficient and the inflow characteristics. Four large eddy simulation (LES) codes from different institutions and a vortex model are part of the comparison. The purpose of this benchmark is to validate the numerical predictions of the flow field statistics in the near wake of an actuator disc, a case that is highly relevant for full wind farm applications. The comparison has shown that, despite its extreme simplicity, the vortex model is capable of reproducing the wake expansion and the centreline velocity with very high accuracy. Also all tested LES models are able to predict the velocity deficit in the very near wake well, contrary to what was expected from previous literature. However, the resolved velocity fluctuations in the LES are below the experimentally measured values.     

Tidal farm analysis using an analytical model for the flow velocity prediction in the wake of a tidal turbine with small diameter to depth ratio

A tidal turbine is a device converting hydrodynamic power into electrical power. Lately, more and more projects have been developed in order to optimize the productivity of this kind of energy. In such research with industrial interest, under the impact of the wake effect on the output power, the analysis of a tidal farm layout is regarded as the first priority. Simple approaches such as those developed for wind farms could be used in tidal turbine arrangement optimization. These methodologies can be improved by taking into account the turbulence in tidal farms and tidal turbines' mechanical characteristics. The goal of this work is to propose a predictive analytical model to estimate the tidal speed in the far wake of tidal turbines with small diameter to depth ratio (20% here). It is a first step prior to integrate the wake model in a tidal farm layout optimization algorithm. The wake model development is achieved reanalyzing the far wake's equations used in wind farm applications. A turbine represented by an Actuator Disc (AD) in conjunction with a Computational Fluid Dynamics (CFD) numerical model is used as a reference for this purpose. The CFD-AD model has been validated with experimental results from literature. The novelty of the present work consists in expressing the far wake's radius expansion as a function of the ambient turbulence and the thrust coefficient. The proposed equation is used in conjunction with the Jensen's model in a manner that the velocities downstream a tidal turbine can be estimated. The velocity distribution in the far wake of a single turbine obtained by the proposed model is in good agreement with the CFD numerical model. As a matter of fact, the model provides satisfactory accuracy in the cases of two parks: one with five aligned turbines and one with ten staggered turbines.