A low‐order model for analysing effects of blade fatigue load control

A new low‐order mathematical model is introduced to analyse blade dynamics and blade load‐reducing control strategies for wind turbines. The model consists of a typical wing section model combined with a rotor speed model, leading to four structural degrees of freedom (flapwise, edgewise and torsional blade oscillations and rotor speed). The aerodynamics is described by an unsteady aerodynamic model. The equations of motion are derived in non‐linear and linear form. The linear equations of motion are used for stability analysis and control design. The non‐linear equations of motion are used for time simulations to evaluate control performance. The stability analysis shows that the model is capable of predicting classical flutter and stall‐induced vibrations. The results from the stability analysis are compared with known results, showing good agreement. The model is used to compare the performance of one proportional–integral–derivative controller and two full‐state feedback controllers. Copyright © 2006 John Wiley & Sons, Ltd.     

An improved BEM model for the power curve prediction of stall‐regulated wind turbines

Blade element momentum (BEM) theory is the standard computational technique for the prediction of power curves of wind turbines; it is based on the two‐dimensional aerodynamic properties of aerofoil blade elements and some corrections accounting for three‐dimensional wing aerodynamics. Although most BEM models yield acceptable results for low‐wind and pitch‐controlled regimes where the local angles of attack are small, no generally accepted model exists up to date that consistently predicts the power curve in the stall regime for a variety of blade properties and operating conditions. In this article we present a modified BEM model which satisfactorily reproduces the power curves of four experimental wind turbines reported in the literature, using no free fit parameters. Since these four experimental cases comprehend a great variety of conditions (wind tunnel vs field experiments, different air densities) and blade parameters (no twist and no taper, no taper but twist, both twist and taper, different aerofoil families), it is believed that our model represents a useful working tool for the aerodynamic design of stall‐regulated wind turbines. Copyright © 2004 John Wiley & Sons, Ltd.     

Development of an aeroelastic code based on three‐dimensional viscous–inviscid method for wind turbine computations

Aerodynamic and structural dynamic performance analysis of modern wind turbines are routinely estimated in the wind energy field using computational tools known as aeroelastic codes. Most aeroelastic codes use the blade element momentum (BEM) technique to model the rotor aerodynamics and a modal, multi‐body or the finite‐element approach to model the turbine structural dynamics. The present work describes the development of a novel aeroelastic code that combines a three‐dimensional viscous–inviscid interactive method, method for interactive rotor aerodynamic simulations (MIRAS), with the structural dynamics model used in the aeroelastic code FLEX5. The new code, called MIRAS‐FLEX, is an improvement on standard aeroelastic codes because it uses a more advanced aerodynamic model than BEM. With the new aeroelastic code, more physical aerodynamic predictions than BEM can be obtained as BEM uses empirical relations, such as tip loss corrections, to determine the flow around a rotor. Although more costly than BEM, a small cluster is sufficient to run MIRAS‐FLEX in a fast and easy way. MIRAS‐FLEX is compared against the widely used FLEX5 and FAST, as well as the participant codes from the Offshore Code Comparison Collaboration Project. Simulation tests consist of steady wind inflow conditions with different combinations of yaw error, wind shear, tower shadow and turbine‐elastic modeling. Turbulent inflow created by using a Mann box is also considered. MIRAS‐FLEX results, such as blade tip deflections and root‐bending moments, are generally in good agreement with the other codes. Copyright © 2017 John Wiley & Sons, Ltd.     

New vortex‐lift and tangential‐force models for HAWT aerodynamic load prediction

Horizontal axis wind turbines (HAWTs) experience three‐dimensional rotational and unsteady aerodynamic phenomena at the rotor blades sections. These highly unsteady three‐dimensional effects have a dramatic impact on the aerodynamic load distributions on the blades, in particular, when they occur at high angles of attack due to stall delay and dynamic stall. Unfortunately, there is no complete understanding of the flow physics yet at these unsteady 3D flow conditions, and hence, the existing published theoretical models are often incapable of modelling the impact on the turbine response realistically. The purpose of this paper is to provide an insight on the combined influence of the stall delay and dynamic stall on the blade load history of wind turbines in controlled and uncontrolled conditions. New dynamic stall vortex and nonlinear tangential force coefficient modules, which integrally take into account the three dimensional rotational effect, are also proposed in this paper. This module along with the unsteady influence of turbulent wind speed and tower shadow is implemented in a blade element momentum (BEM) model to estimate the aerodynamic loads on a rotating blade more accurately. This work presents an important step to help modelling the combined influence of the stall delay and dynamic stall on the load history of the rotating wind turbine blades which is vital to have lighter turbine blades and improved wind turbine design systems.     

Self‐induced vibrations of a DU96‐W‐180 airfoil in stall

This work presents an analysis of two‐dimensional (2D) and three‐dimensional (3D) non‐moving, prescribed motion and elastically mounted airfoil computational fluid dynamics (CFD) computations. The elastically mounted airfoil computations were performed by means of a 2D structural model with two degrees of freedom. The computations aimed at investigating the mechanisms of both vortex‐induced and stall‐induced vibrations related to a wind turbine blade at standstill conditions. In this work, a DU96‐W‐180 airfoil was used in the angle‐of‐attack region potentially corresponding to stall‐induced vibrations. The analysis showed significant differences between the aerodynamic stability limits predicted by 2D and 3D CFD computations. A general agreement was reached between the prescribed motion and elastically mounted airfoil computations. 3D computations indicated that vortex‐induced vibrations are likely to occur at modern wind turbine blades at standstill. In contrast, the predicted cut‐in wind speed necessary for the onset of stall‐induced vibrations appeared high enough for such vibrations to be unlikely. Copyright © 2013 John Wiley & Sons, Ltd.     

Experimental and computational analysis of stall cells on rectangular wings

The present paper is the second part of a combined (experimental and computational) study on stall cells (SCs) on a rectangular wing. In the first part, tuft data were used in order to geometrically characterize a stabilized SC resulting from a localized spanwise disturbance introduced by a zigzag tape. Here, pressure measurements on the model and in the wake and aerodynamic polars at midspan are reported. The wing model had an aspect ratio value of 2, the Reynolds number was 10⁶ and the range of angles of attack (α) was from ?6^° to 16^°. Experimental results confirm previous findings. Furthermore, two‐dimensional and three‐dimensional Reynolds Averaged Navier‐Stokes RANS simulations are used in order to better understand the structure of SCs. 3D simulations reproduce the experimental data with a 3° delay in α and permit a qualitative analysis. It is found that the SC vortices start normal to the wing surface and extend downstream in the wake; the evolution of the SC vortices in the wake is in strong interaction with the separation line vortex and the trailing edge line vortex; as the SC vortex develops downstream in the wake, its centreline is contracted towards the SC centre; the wing wake is pushed upstream at the centre of the SC and downstream at the sides by the SC vortices; spanwise lift and drag distributions always attain their minimum at the SC centre. Copyright © 2013 John Wiley & Sons, Ltd.     

Combined analytical/FEA-based coupled aero structure simulation of a wind turbine with bend–twist adaptive blades

The simulation of wind turbines with bend–twist adaptive blades is a coupled aero-structure (CAS) procedure. The blade twist due to elastic coupling is a required parameter for wind turbine performance evaluation and can be predicted through a finite element (FE) structural analyser. FEA-based codes are far too slow to be useful in the aerodynamic design/optimisation of a blade. This paper presents a combined analytical/FEA-based method for CAS simulation of wind turbines utilising bend–twist adaptive blades. This method of simulation employs the induced twist distribution and the flap bending at the hub of the blade predicted through a FEA-based CAS simulation at a reference wind turbine run condition to determine the wind turbine performance at other wind turbine run conditions. This reduces the computational time significantly and makes the aerodynamic design/optimisation of bend–twist adaptive blades practical. Comparison of the results of a case study which applies both combined analytical/FEA-based and FEA-based CAS simulation shows that when using the combined method the required computational time for generating a power curve reduces to less than 5%, while the relative difference between the predicted powers by two methods is only about 1%.     

Inverse design of horizontal axis wind turbine blades using a vortex line method

A new inverse design process for horizontal axis wind turbine blades is developed to account for three‐dimensional blade features such as non‐planar wing tip. The multidimensional Newton iteration method combined with a vortex line method is used to provide blade geometry parameters given desired aerodynamic behaviors such as lift coefficient and axial induction. The Jacobian matrix is visualized to show the effect of the change of the blade twist and chord on the change of the aerodynamic behaviors. The method is validated for a canonical straight blade with uniform lift coefficient and axial induction distributions. The results show an excellent agreement with those obtained by PROPID, which is a blade element momentum theory‐based inverse design code. The National Renewable Energy Laboratory Phase VI blade is used to validate the method for a straight blade with non‐uniform distributions of the lift coefficient and axial induction. The method is also applied successfully to a non‐straight blade design with a non‐planar wing tip. A noticeable change in the twist and chord for this non‐straight blade is seen compared with a straight blade. Finally, the inverse design code is used to make a large rotor blade, and the power output generated by this blade is computed. Copyright © 2014 John Wiley & Sons, Ltd.     

Fluid structure interaction of a morphed wind turbine blade

One serious challenge of energy systems design, wind turbines in particular, is the need to match the system operation to the variable load. This is so because system efficiency drops at off‐design load. One strategy to address this challenge for wind turbine blades and obtain a more consistent efficiency over a wide load range, is varying the blade geometry. Predictable morphing of wind turbine blade in reaction to wind load conditions has been introduced recently. The concept, derived from fish locomotion, also has similarities to spoilers and ailerons, known to reduce flow separation and improve performance using passive changes in blade geometry. In this work, we employ a fully coupled technique on CFD and FEM models to introduce continuous morphing to desired and predetermined blade design geometry, the NACA 4412 profile, which is commonly used in wind turbine applications. Then, we assess the aerodynamic behavior of a morphing wind turbine airfoil using a two‐dimensional computation. The work is focused on assessing aerodynamic forces based on trailing edge deflection, wind speed, and material elasticity, that is, Young's modulus. The computational results suggest that the morphing blade has superior part‐load efficiency over the rigid NACA blade. Copyright © 2013 John Wiley & Sons, Ltd.     

压气机分层优化气动设计体系研究

基于压气机分层次气动优化设计思想,结合了优化算法、CFD(计算流体力学)技术与压气机气动设计程序,利用商业软件,建立了压气机气动优化设计平台,将气动优化设计思想融入到了压气机气动设计多个阶段,以实现压气机气动设计使性能最优化的目的.应用商业通用优化平台进行了轴流压气机一维、S2反问题气动优化设计,并采用NUMECA提供的全三维优化设计平台对压气机进行了全三维气动优化设计.计算结果表明:分层气动优化设计是提高压气机气动性能的有效手段,先进的优化算法在气动设计的各阶段都能够比传统设计手段更大限度地实现压气机性能的最优化.     

Laminar‐turbulent flow simulation for wind turbine profiles using the γ– transition model

The accurate prediction of the laminar‐turbulence transition process is fundamental in predicting the aerodynamic performance of wind turbine profiles. Fully turbulent flow simulations have been shown to over‐predict the aerodynamic performance and thereby negatively impacting the design of airfoils in flow regimes where the possible presence of laminar flow could be exploited to improve the performance of wind turbine rotors. Correlation‐based transition modelling offers a fully computational fluid dynamics compatible approach, where the model integrates completely with the existing turbulence model, allows for the prediction of various transition mechanisms, is applicable to three‐dimensional flows and compatible to adjoint‐based design optimization frameworks. The present paper addresses several modifications necessary for a robust transition model and investigates the accuracy of the model for a wide range of angles of attack and Reynolds numbers, which are necessary for a thorough validation of the correlation‐based transition model for wind turbine profiles. The transition model was employed to predict the transition locations; and an assessment of the various transition mechanisms, Reynolds number effects, sectional characteristics and aerodynamic performance for the NLF(1)‐0416 and S809 airfoils is presented with comparisons to experimental data and numerical solutions. Copyright © 2013 John Wiley & Sons, Ltd.     

Structural investigation of composite wind turbine blade considering various load cases and fatigue life

This study proposes a structural design for developing a medium scale composite wind turbine blade made of E-glass/epoxy for a 750 kW class horizontal axis wind turbine system. The design loads were determined from various load cases specified at the IEC61400-1 international specification and GL regulations for the wind energy conversion system. A specific composite structure configuration, which can effectively endure various loads such as aerodynamic loads and loads due to accumulation of ice, hygro-thermal and mechanical loads, was proposed. To evaluate the proposed composite wind turbine blade, structural analysis was performed by using the finite element method. Parametric studies were carried out to determine an acceptable blade structural design, and the most dominant design parameters were confirmed. In this study, the proposed blade structure was confirmed to be safe and stable under various load conditions, including the extreme load conditions. Moreover, the blade adapted a new blade root joint with insert bolts, and its safety was verified at design loads including fatigue loads. The fatigue life of a blade that has to endure for more than 20 years was estimated by using the well-known S–N linear damage theory, the service load spectrum, and the Spera's empirical equations. With the results obtained from all the structural design and analysis, prototype composite blades were manufactured. A specific construction process including the lay-up molding method was applied to manufacturing blades. Full-scale static structural test was performed with the simulated aerodynamic loads. From the experimental results, it was found that the designed blade had structural integrity. In addition, the measured results of deflections, strains, mass, and radial center of gravity agreed well with the analytical results. The prototype blade was successfully certified by an international certification institute, GL (Germanisher Lloyd) in Germany.     

A CFD study of coupled aerodynamic‐hydrodynamic loads on a semisubmersible floating offshore wind turbine

The prediction of dynamic characteristics for a floating offshore wind turbine (FOWT) is challenging because of the complex load coupling of aerodynamics, hydrodynamics, and structural dynamics. These loads should be accurately calculated to yield reliable analysis results in the design phase of a FOWT. In this study, a high‐fidelity fluid‐structure interaction simulation that simultaneously considers the influence of aero‐hydrodynamic coupling due to the dynamic motion of a FOWT has been conducted using computational fluid dynamics based on an overset grid technique. The DeepCwind semisubmersible floating platform with the NREL 5‐MW baseline wind turbine model is considered for objective numerical verification with the NREL FAST code. A state‐of‐the‐art computational model based on the coupled computational fluid dynamics and dynamic structure analysis is constructed and analyzed to solve multiphase flow, 6 degrees of freedom motions of OC4 semisubmersible FOWT. A quasi‐static mooring solver is also applied to resolve the constraint motion of floater because of a 3‐line mooring system. The influence of tower shadow on the unsteady aerodynamic performance and loads is also demonstrated. Finally, complex unsteady flow fields considering blade and tower interference effects among blade‐tip vortices, shedding vortices, and turbulent wakes are numerically visualized and investigated in detail.     

Innovative design approaches for large wind turbine blades

A preliminary design study of an advanced 50 m blade for utility wind turbines is presented and discussed. The effort was part of the Department of Energy WindPACT Blade System Design Study with the goal to investigate and evaluate design and manufacturing issues for wind turbine blades in the 1–10 MW size range. Two different blade designs are considered and compared in this article. The first is a fibreglass design, while the second design selectively incorporates carbon fibre in the main structural elements. The addition of carbon results in modest cost increases and provides significant benefits, particularly with respect to blade deflection. The structural efficiency of both designs was maximized by tailoring the thickness of the blade cross‐sections to simplify the construction of the internal members. Inboard the blades incorporate thick blunt trailing edge aerofoils (flatback aerofoils), while outboard more conventional sharp trailing edge high‐lift aerofoils are used. The outboard section chord lengths were adjusted to yield the least complex and costly internal blade structure. A significant portion of blade weight is related to the root buildup and metal hardware for typical root attachment designs. The results show that increasing the number of studs has a positive effect on total weight, because it reduces the required root laminate thickness. The aerodynamic performance of the blade aerofoils was predicted using computational techniques that properly simulate blunt trailing edge flows. The performance of the rotor was predicted assuming both clean and soiled blade surface conditions. The rotor is shown to provide excellent performance at a weight significantly lower than that of current rotors of this size. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of seismic and aerodynamic load interaction on structural dynamic response of multi-megawatt utility scale horizontal axis wind turbines

Horizontal axis wind turbines can experience significant time varying aerodynamic loads that has the potential to cause adverse effects on structural, mechanical, and power production. The progress in the wind industry has caused the construction of wind farms in areas prone to high seismic activity. With the advances in computational tools, a more realistic representation of the behavior of wind turbines should be performed. One of the simulation platforms was developed using the 5 MW NREL utility scale reference turbine model. The performed simulations will be used to evaluate the effects of aerodynamic and seismic load coupling on the power generation and structural dynamics behavior of this structure. Different turbine operational scenarios such as (i) normal operational condition with no earthquake, (ii) idling condition with the presence of seismic loads, (iii) normal operational condition with earthquake, and (iv) earthquake-induced emergency shutdown will be simulated with various loading conditions to show the differences in generated power and dynamic response. The results of this paper provide formulations for calculating generated power and design deriving parameters by considering different intensity measures. Moreover, the effects of aerodynamic damping and pitch control system are presented to shows reduction in the resulting design demand loads.     

Numerical validation of a finite element thin‐walled beam model of a composite wind turbine blade

This paper presents a numerical validation of a thin‐walled beam (TWB) finite element (FE) model of a realistic wind turbine rotor blade. Based on the theory originally developed by Librescu et al. and later extended to suit FE modelling by Phuong, Lee and others, this computationally efficient yet accurate numerical model is capable of capturing most of the features found in large blades including thin‐walled hollow cross section with variable thickness along the section's contour, inner reinforcements, arbitrary material layup and non‐linear anisotropic fibre‐reinforced composites; the present application is, for the time being, restricted to linearity. This one‐dimensional (1D) FE model allows retaining information of different regions of the blade's shell and therefore approximates the behaviour of more complex three‐dimensional (3D) shell or solid FE models more accurately than typical 1D FE beam models. A 9.2 m rotor blade, previously reported in specialized literature, was chosen as a case study to validate the static and dynamic behaviour predicted by a TWB model against an industry‐standard 3D shell model built in a commercial software tool. Given the geometric and material complexities involved, an excellent agreement was found for static deformation curves, as well as a good prediction of the lowest frequency modes in terms of resonance frequencies, mode shapes and frequency response functions; the highest (sixth) frequency mode shows only a fair agreement as expected for an FE model. It is concluded that despite its simplicity, a TWB FE model is sufficiently accurate to serve as a design tool for the recursive analyses required during design and optimization stages of wind turbines using only readily available computational tools. Copyright © 2011 John Wiley & Sons, Ltd.     

某1 MW水平轴风力机叶片气动设计及载荷计算

基于片条理论,考虑了叶尖损失、叶根损失、叶柵影响和重载荷下对片条理论参数修正的情况下,完成了某1 MW水平轴风力机叶片的气动设计,并对其气动性能进行了评估;最后根据IEC规范对叶片在不同风况状态下进行载荷计算,所得结果可为同类风力机气动设计和结构设计提供参考。     

Demonstrating that power and instantaneous loads are decoupled in a vertical‐axis wind turbine

The flow in the meridian plane of a high aspect ratio vertical‐axis wind turbine (VAWT) can be described as two dimensional. The wake that is generated by the VAWT in a two‐dimensional flow consists of shed vorticity and is a result of the temporal variation of bound circulation on the blades, following Kelvin's theorem. The strength and location of the vorticity that is produced by the VAWT in a two‐dimensional flow are thus independent of the average bound circulation on the blade. Two independent computational models—a potential flow panel model and a method that is based on the vorticity–velocity formulation of the Navier–Stokes equations—have been used to show that the VAWT can produce the same power for different azimuthal distributions of the blade aerodynamic loading. It is thus demonstrated that the instantaneous blade aerodynamic loading and the power conversion of a VAWT are decoupled. This observation has, potentially, significant impact on the design of the VAWT and reopens the research on asymmetric blade shapes in order to optimize the performance of this turbine configuration. Copyright © 2013 John Wiley & Sons, Ltd.