ALBERT

Description: We show that the upscaling of wind turbines from rotor diameters of 15–20 m to presently large rotors of 150–200 m has changed the requirements for the aerodynamic blade element momentum (BEM) models in the aeroelastic codes. This is because the typical scales in the inflow turbulence are now comparable with the rotor diameter of the large turbines. Therefore, the spectrum of the incoming turbulence relative to the rotating blade has increased energy content on 1P, 2P, …, nP, and the annular mean induction approach in a classical BEM implementation might no longer be a good approximation for large rotors. We present a complete BEM implementation on a polar grid that models the induction response to the considerable 1P, 2P, …, nP inflow variations, including models for yawed inflow, dynamic inflow and radial induction. At each time step, in an aeroelastic simulation, the induction derived from a local BEM approach is updated at all the stationary grid points covering the swept area so the model can be characterized as an engineering actuator disk (AD) solution. The induction at each grid point varies slowly in time due to the dynamic inflow filter but the rotating blade now samples the induction field; as a result, the induction seen from the blade is highly unsteady and has a spectrum with distinct 1P, 2P, …, nP peaks. The load impact mechanism from this unsteady induction is analyzed and it is found that the load impact strongly depends on the turbine design and operating conditions. For operation at low to medium thrust coefficients (conventional turbines at above rated wind speed or low induction turbines in the whole operating range), it is found that the grid BEM gives typically 8 %–10 % lower 1 Hz blade root flapwise fatigue loads than the classical annular mean BEM approach. At high thrust coefficients that can occur at low wind speeds, the grid BEM can give slightly increased fatigue loads. In the paper, the implementation of the grid-based BEM is described in detail, and finally several validation cases are presented. Comparisons with blade loads from full rotor CFD, wind tunnel experiments and a field experiment show that the model can predict the aerodynamic forces in half-wake, yawed flow, dynamic inflow and turbulent inflow conditions.

Description: Laminar-turbulent transition behavior of a wind turbine blade section is investigated in this study by means of field experiments and 3-D computational fluid dynamics (CFD) rotor simulations. The power spectral density (PSD) integrals of the pressure fluctuations obtained from the high-frequency microphones mounted on a blade section are analyzed to detect laminar-turbulent transition locations from the experiments. The atmospheric boundary layer (ABL) velocities and the turbulence intensities (T.I.) measured from the field experiments are used to create several inflow scenarios for the CFD simulations. Results from the natural and the bypass transition models of the in-house CFD EllipSys code are compared with the experiments. It is seen that the bypass transition model results fit well with experiments at the azimuthal positions where the turbine is under wake and high turbulence, while the results from other cases show agreement with the natural transition model. Furthermore, the influence of inflow turbulence, wake of an upstream turbine, and angle of attack (AOA) on the transition behavior is investigated through the field experiments. On the pressure side of the blade section, at high AOA values and wake conditions, variation in the transition location covers up to 44 % of the chord during one revolution, while for the no-wake cases and lower AOA values, variation occurs along a region that covers only 5 % of the chord. The effect of the inflow turbulence on the effective angle of attack as well as its direct effect on transition is observed. Transition locations for the wind tunnel conditions and field experiments are compared together with 2-D and 3-D CFD simulations. In contrast to the suction side, significant difference in the transition locations is observed between wind tunnel and field experiments on the pressure side for the same airfoil geometry. It is seen that the natural and bypass transition models of EllipSys3D can be used for transition prediction of a wind turbine blade section for high-Reynolds-number flows by applying various inflow scenarios separately to cover the whole range of atmospheric occurrences.