Sensitivity of Dynamic Stall Models to Dynamic Excitation on Large-Flexible Wind Turbine Blades in Edgewise Vibrations

1. Introduction

Wind turbine loads are heavily influenced by the operating conditions and environments [1]. Therefore, it is defined by the standards [2,3] that blade designers/manufacturers need to run thousands of load cases that simulate wind turbines under various conditions in the design and certification processes. The test cases vary from normal power production operations (typically aimed at fatigue loads analysis) to the extreme conditions when the turbine is shut-down, either in parked or idling modes. The parked/idling conditions are typically applied for high wind speed cases above the cut-out wind speed of the turbine. In such conditions, wind turbine controller is normally no longer active when no backup power supply is available. Hence, the blade is prone to operate at large angle of attack depending on the yaw angle, and it remains locked at this position. This could lead to strong unsteady characteristics and the aeroelastic behavior of the blade can be dangerous for its own structural integrity. In these conditions, the unsteady aerodynamic effects on the blade can manifest as dynamic stall when the blade sections operate at high angles of attack.

The consideration of unsteady aerodynamic effects becomes even more important in the last few decades since the blade size has increased considerably, making the blade more slender and flexible. The aeroelastic characteristics of such long blades are not well understood where instability and vibrations often become huge concerns and limiting factors for designs. Despite this clear environmental challenge, wind turbine blades are required to be able to survive even in harsh conditions as defined in the set of loads calculation within the IEC [2] and DNV [3] standards.

Wind turbine aeroelastic response is strongly affected by the amount of damping on the blade itself, that can generally be divided into the structural and aerodynamic damping. When the total damping of the blade becomes negative, instability of the blade is enhanced. It is well understood that the edgewise direction is prone to negative damping conditions. The structural damping represents the structural design and the material used to manufacture the blade itself, and it is relatively easy to estimate with high confidence level. However, predicting the aerodynamic damping is not a trivial task because it can significantly change depending on the operating conditions of the turbine. This type of damping is affected by the acting forces on the blade, which alter the angle of attack due to elastic response.

The aeroelastic phenomenon in these extreme conditions is highly unsteady, and as a direct consequence, the induction-based blade element momentum theory commonly adopted in engineering tools needs to be corrected using a dynamic stall model. Examples of such models are the Beddoes-Leishman (BL) model [4], Larsen model [5], Øye model [6], Snel model [7], Risø model [8], IAG model [9,10], ONERA model [11], Boeing-Vertol model [12], etc.

Understanding how these dynamic stall models behave in challenging conditions will be an important factor to properly estimate the impacts on aeroelastic characteristics. Several studies have attempted to evaluate the dynamic stall characteristics of wind turbines in standstill conditions (parked/idling) [13,14,15,16,17,18,19,20], and it is commonly believed that dynamic stall models tend to underestimate the true aerodynamic damping due to inaccurate hysteresis predictions [21]. Hansen [14] provided a proposal for an improved modal analysis to improve the understanding of wind turbine stability. Sarkar and Bijl [22] evaluated non-linear aeroelastic behavior of an oscillating airfoil during stall-induced vibration in two-dimensional environments. Wang et al. [15], for instance, focused their studies on the linear model predictions when using the BL model and the ONERA model. Bir and Jonkman [23] investigated the aeroelastic instability on a 5 floating MW wind turbine, and pointed out the distinction between aeroelastic instability and resonance, where the latter requires external excitation at a distinct frequency. Tibaldi et al. [24] performed a numerical investigation on wind turbine resonant vibrations, where it was highlighted that depending on the external excitation, different aeroelastic modes can be excited. The responses of different dynamic stall models in wind turbine under standstill conditions are further evaluated for example by Wang et al. [15], Faber [25], Bangga et. al. [13], Andreou [26], Bangga and Yu [17]. In Ref. [27], commonly adopted engineering approaches were reviewed and discussed. Horcas et al. [28] provided an improved understanding of standstill flow conditions by adopting high fidelity computational fluid dynamics (CFD) approaches, but a direct adoption of the studies to practical engineering methods are still yet to be explored.

The present paper is specifically aimed at addressing the critical issue on how sensitive dynamic stall models actually are when they are exposed to dynamic excitation in edgewise direction. To the best of the author’s knowledge, comprehensive studies in this area are lacking for extremely large wind turbine blade structure. This is especially crucial when different dynamic stall models and excitation frequencies will be considered that lead to uncertainties in wind turbine blade design and certification processes. Considering the current state-of-the-art, the present paper investigates the aeroelastic sensitivity of a large wind turbine blade using three different dynamic stall models commonly applied in the designs, namely (1) the Øye model [6], (2) the incompressible version of the BL model (Risø model [8]) and the newly developed IAG model [10]. The studies start by comparing different dynamic stall models against a pitching airfoil case, where experimental data is available. Then the investigations are extended by exciting the rotor blade by means of an oscillating wind at various frequencies, ranging from below to above the first edgewise modal frequency of the blade. The growth of the angle of attack response is monitored at different wind speeds to quantify the direct influence of the unsteady aerodynamic models on the aeroelastic responses.

The present paper is organized as follows. method provides a description of the employed code, approaches and the test cases of interest. The results are discussed in resultsOSUValidation to resultstwist starting by the comparison of the pitching airfoil data, then continued with the aeroelastic response at the blade level. Finally, the paper will be concluded with suggestions for follow up studies in conclusion.

2. Methodology

2.1. Wind Turbine Design Code Bladed

Bladed is an engineering tool used in numerous wind turbine industries for designs and certifications as well as in research environment. The code solves aeroservohydroelastic problems of wind turbines by utilizing blade element momentum (BEM) theory coupled with structural and control algorithms. The rotor aerodynamic solutions are obtained using time-marching approaches of the induced velocities at each blade section. Bladed uses a full local BEM approach implementation [8,29]. The solutions are coupled with the structural dynamic models that adopt a multibody dynamics approach. A multi-part blade method is adopted to allow the structural dynamic equations to model larger deflections using a linear beam theory. Efficient solutions of these complex nonlinear problems are made possible via the fixed-step implicit Newmark-$\beta$ integrator. Several publications documented the verification studies of the Bladed code performance and accuracy [10,30,31]. Bladed version 4.16.0 is utilized in the present studies.

In the present studies, three dynamic stall models in Bladed are adopted, namely ye [6], Beddoes-Leishman (BL) [8,29] and a recently developed IAG model [10]. The incompressible version of the BL model [8,29] has been widely adopted in wind turbine designs and can be considered as the state-of-the-art model. Although the model works sufficiently well for most conditions, the accuracy of the model is constantly challenged when it operates at complex inflow conditions at large yaw misalignment. This implies that the angle of attack of the blade sections can be significantly large [10]. The IAG model has recently been developed specifically for wind turbine applications and has been tested against experimental data at various conditions [9,10]. The robustness and applicability of the IAG dynamic stall model for wind energy field were further improved in the second generation model by enabling a formulation in the state-space representation, as well as by removing the compressibility effects and by decreasing the number of required constants [10]. In Ref. [10], it was demonstrated that the IAG model performs well against experimental data and it outperforms other models when the airfoil operates at large angles of attack beyond stall, and the tests were carried out at various reduced frequencies and amplitudes. The IAG model was further compared against other dynamic stall models for various load cases in Refs. [13,17].

2.2. Blade Specification and Airfoil Data Treatment

The present studies mainly adopt a blade constructed by the S809 airfoil. This airfoil is chosen because there is experimental data available both for static and dynamic conditions [32]. Thus, validation of the results for pitching airfoil case will be possible. This airfoil has a relative thickness of $t/c$ = 21% and was measured at a Reynolds number of a million. For the present studies, results based on the airfoil equipped with a leading edge grit (turbulator) will be adopted. This is done to simulate the "soiled" effects on a wind turbine blade in field.

In the present studies, the blade geometry of the IEA 15 MW reference wind turbine version 1.1 is adopted [33]. The blade is selected since it is representative for the current blade design that is associated with all complexity and highly flexible structures. The blade has a length of about 117 m. Starting from 91.33 m of the blade length, the relative thickness of the airfoil section (relative to chord) is about 21%. In the present studies, the IEA 15 MW reference wind turbine blade is modified by replacing the applied airfoils at all blade sections. First, the airfoil data from 91.33 m of the blade length is replaced with the S809 airfoil. The remaining parts of the blade up to the root area are then generated based on interpolation (according to relative thickness). The root of the blade adopts a cylinder with a constant drag coefficient of about 0.35. This process is illustrated in Figure 1 and Figure 2. It is noted that most investigations will be focused on the dynamic response of the blade section at 93.68 m except stated otherwise. A more realistic test case can be designed by using the actual blade airfoils without modification. However, this will involve complex behavior of different airfoils with different stall characteristics. This makes it difficult to draw conclusions and to analyze further. Thus, the simplification done in the present paper is used as a starting point.

Since the S809 airfoil was measured only to a limited range of the angle of attack ($\alpha$), it is a common practice to extrapolate the data to cover $\alpha$ = -180^∘ to $\alpha$ = +180^∘ regime. In this paper, a modified Viterna approach [34] is adopted and the results are plotted in figPolarExtrapolated. The description of the method is given in Appendix A. The usage of the extrapolated data will affect the accuracy of wind turbine calculations. The reconstruction accuracy deteriorates the shorter the data one has. A 360^∘ extrapolation can be used in early stage of wind turbine designs to give a rough idea about the turbine loads before any further detailed assessments are made. It is highlighted in Ref [35] how the dynamic calculations are influenced by the steady data reconstruction.

The extrapolated airfoil is then interpolated based on the relative thickness distribution shown in Figure 3. Toward the root area, cylinder is applied and the interpolation results at various relative thicknesses are shown in Figure 2. The chord distribution of the IEA 15 MW reference wind turbine is applied without change, but twist distribution is modified to zero (untwisted blade) to simplify the problem of interest. However, it is understandable that this is not realistic for real wind turbine blades. Therefore, the influence of twist distribution will be given in Section 6. The structural properties like mass, bending stiffness, shear stiffness, torsional stiffness and axial stiffness distributions are adopted without any change. This includes the shear center distribution, neutral axis, principal axis orientation and mass axis orientation, but these parameters are not shown in this paper for conciseness.

2.3. Description of the Test Cases

In the present studies, the test case is limited only to cover the blade characteristics to simplify the physical interactions and also to enable isolation of the blade excitation effects. Therefore, tower and all other structures are set to rigid and only one blade is modeled. The blade is parked at an azimuth angle of zero (blade pointing up) with the blade being pitched to 71.45^∘, which effectively creates a mean angle of attack of 18.55^∘ with the wind coming from upstream. This mean angle of attack is chosen to match the experimental data for a pitching airfoil case defined in [32]. The rotational speed of the rotor is set to zero to simulate the turbine condition in parked operation. The calculation for the momentum theory is no longer relevant for this case, therefore the induction calculations and the associated engineering sub-models like tip loss correction and dynamic wake are switched off. However, dynamic stall effects are included since they are highly relevant for the case. In the present studies, three dynamic stall models are tested: ye, BL and IAG models. The blade is exposed to sinusoidal variation of the wind velocity in the following form:

$$\begin{array}{cc}\hfill u\left(t\right)& =U_{\infty }cos\alpha \left(t\right)\hfill \\ \hfill \end{array}$$

(1)

$$\begin{array}{cc}\hfill v\left(t\right)& =U_{\infty }sin\alpha \left(t\right)\hfill \end{array}$$

(2)

where

$$\alpha \left(t\right)=\Delta \alpha sin\left(2\pi ft\right)$$

(3)

Here, $\Delta \alpha$ represents the excitation amplitude and is set to be constant at 10.45^∘, which again corresponds to the experimental data for a pitching airfoil case defined in [32]. The values of the incoming wind speed $U_{\infty }$ and excitation frequency f are varied from 10 m/s - 50 m/s and 0.086 Hz - 1.03 Hz, respectively. Depending on the test cases being investigated, it is often more convenient to represent the frequency in non-dimensional form. In the studies, this is done either in form of the reduced frequency:

$$k={\displaystyle \frac{2\pi fc}{2U_{\infty }}}$$

(4)

or as a ratio to the natural frequency. All simulations are performed for 100 s with a timestep size of 0.02 s.

Since the studies are focused on the edgewise resonance, it is a logical choice to use the natural frequency of the blade edgewise mode ($f_0$). For this purpose, a modal analysis was performed in Bladed which yields the modes listed in Table 1. It can be seen that the 1st edgewise modal frequency is about 0.687 Hz. This means that the frequencies being evaluated (0.086 Hz - 1.03 Hz) will be within the range of 0.125$f_0$ to 1.5$f_0$.

3. Comparison of Dynamic Stall Model against Pitching Airfoil Data

To verify the accuracy of dynamic stall models, Bladed calculations are compared with experimental data of a pitching S809 airfoil. This is done specifically at the blade station located at 93.68 m. Since the pitching airfoil data obtained from Ref. [32] does not include flexibility (measured on rigid airfoil in a wind tunnel), flexibility setting for the blade is turned off in Bladed. Equation (1) and Equation (2) are applied through a Wind DLL module in Bladed with an angle of attack variation in Equation (3) being set such that the reduced frequency in Equation (4) at 93.68 m matches the measurement campaign. This effectively simulates a pitching airfoil case at k = 0.079.

The results for all three force components - lift ($C_L$), drag ($C_D$) and pitching moment about a quarter chord ($C_M$) - are presented in Figure 4. Both static and dynamic data sets obtained from experimental campaign are presented for comparison. It can be seen that the hysteresis effects of $C_L$, $C_D$ and $C_M$ are massively underestimated by the BL and Øye models. The IAG model clearly shows a better match with the experimental data at this deep stall situation. The improved prediction accuracy can be attributed to a more accurate estimation of the leading edge vortex effects in the IAG model [9,10].

4. Dynamic Responses of the Aerodynamic Coefficients

In Section 4.1 and Section 4.2, the effects of blade flexibility are taken into account. Here the dynamic characteristics of lift, drag and pitching moments are evaluated against two different conditions. First, the same reduced frequency as applied in Section 3 will be adopted in the discussion given in Section 4.1, but now with the blade flexibility being activated. Second, resultspolarresonance discusses the dynamic responses of the aerodynamic forces when the excitation frequency f is set to match the 1st edgewise modal frequency of the blade, see tabmodalfrequencies.

4.1. Exposed to a Constant Reduced Frequency

The simulations at a constant reduced frequency of k = 0.079 (defined at 93.68 m blade station) are performed within $U_{\infty }$ = [10,20,30,40,50] m/s. Here the magnitude of the wind speed becomes important on top of the excitation frequency because this drives the blade deflection. For the sake of conciseness, only the results for $U_{\infty }$ = 10 m/s, 30 m/s and 50 m/s for lift are presented in Figure 5, Figure 6 and Figure 7, respectively. Results for different dynamic stall models for both rigid and flexible blades are presented. The complete plot for lift, drag and pitching moment are available in Appendix B.1.

For the smallest wind speed of 10 m/s, it can be seen in Figure 5 that blade flexibility has practically almost negligible influences. The difference between rigid and flexible settings remains small. Now, by increasing the wind speed to 30 m/s in Figure 6, it can be seen that the results with flexible blade start to deviate from the rigid blade setting. This is particularly seen at upstroke regime (e.g., the higher lift part of the hysteresis). This behavior stems from the fact that the angle of attack increases within this part of the dynamic motion and dynamic solutions are more sensitive to engineering models being applied. The IAG model shows increased vibration characteristics compared to other models. Looking back at Figure 4, it seems that this is how the blade will likely react in reality (note that the same reduced frequency is simulated here).

As for the highest evaluated wind speed of 50 m/s, results start to become much more radical as seen in Figure 7. It can be seen that Øye model, BL model and IAG model behave fairly differently. The Øye model shows some increased variations of $C_L$ in terms of the force magnitudes and the angle of attack range. In contrast to the Øye model, BL model shows somewhat reduced hysteresis curves. A similar behavior is also observed in the IAG model to some degree. However, since the vortex effects in the IAG model are more pronounced than the BL model, it can be seen that force fluctuations are actually increased locally at high angles of attack. The average dynamic polar gradient of the IAG model seems to increase, which might be attributed to the impulsive effects.

4.2. Exposed to Edgewise Resonance Frequency

In this section, the blade is exposed to a constant frequency that matches the 1st edgewise modal frequency of the blade, which excites resonance. The simulations were done at varying wind speeds from 10 m/s - 50 m/s. This means that the reduced frequency itself decreases with the wind speed. Similar to the discussion in Section 4.1, only the results for $U_{\infty }$ = 10 m/s, 30 m/s and 50 m/s for lift are presented for the sake of conciseness. The complete plot for lift, drag and pitching moment are available in Appendix B.2.

In Figure 8, the results for $U_{\infty }$ = 10 m/ are presented for three different dynamic stall models. It shall be noted that this case has the highest reduced frequency of k = 0.598, which in theory shall generate the strongest vibration. In contrast, the wind speed is the smallest (10 m/s) which limits the flexibility effects. It can be observed in Figure 8 that the Øye model yields the weakest vibrations. The IAG model is characterized by the strongest hysteresis effects. This is mainly caused by the impulsive effects being dominant at high reduced frequency. On top of that, the blade flexibility, albeit at the relatively small wind speed, enhances the dynamic stall effects which yield in increased vibrations.

Increasing the wind speed to 30 m/s yields a drastic change of the vibration characteristics. This test case has a reduced frequency of about 0.199 that is still considerably high. On top of that, blade flexibility effects start to become stronger. It can be seen that the vibrations in the Øye model are amplified massively, which occurs mainly only during the vortex breakdown phase, see that the discrepancies of the rigid with the flexible case occur mainly starting from $\alpha$ = 18^∘. The BL model is severely affected by the resonance effects throughout all considered excitation range. It can be seen that $C_L$ varies violently from the lowest to the highest angle of attack. The IAG model, on the other hand, shows somewhat a more stabilized vibration behavior and the gradient of the dynamic polar is larger than the other models. At this reduced frequency, it is expected that the strength of the impulsive and circulatory effects is fairly comparable.

For the highest considered wind speed of 50 m/s, the vibrations from the ye and BL models are now stabilized to a smaller level than the 30 m/s case. This comes from the fact that the actual reduced frequency is smaller than the 30 m/s case to the 50 m/s case. The vibration characteristics for the IAG model is fairly comparable with the 30 m/s case. For this flow condition, the circulatory effects are getting stronger while the impulsive effects are becoming smaller.

Figure 8. Dynamic lift polar calculated at a constant resonance frequency of f = 0.687 Hz at 10 m/s (at k = 0.598).

Figure 8. Dynamic lift polar calculated at a constant resonance frequency of f = 0.687 Hz at 10 m/s (at k = 0.598).

Figure 9. Dynamic lift polar calculated at a constant resonance frequency of f = 0.687 Hz at 30 m/s (at k = 0.199).

Figure 9. Dynamic lift polar calculated at a constant resonance frequency of f = 0.687 Hz at 30 m/s (at k = 0.199).

Figure 10. Dynamic lift polar calculated at a constant resonance frequency of f = 0.687 Hz at 50 m/s (at k = 0.12).

Figure 10. Dynamic lift polar calculated at a constant resonance frequency of f = 0.687 Hz at 50 m/s (at k = 0.12).

5. Angle of Attack Range Variations

In this section, the range of the angle of attack change from the flexible blade simulation to the rigid blade simulation is compared. This is formulated as:

$$\Delta {\alpha }^{Flex}-\Delta {\alpha }^{Rigid}=\left({\alpha }^{Flex,MAX}-{\alpha }^{Flex,MIN}\right)-\left({\alpha }^{Rigid,MAX}-{\alpha }^{Rigid,MIN}\right),$$

(5)

which highlights how much blade vibration grows over the time. Similar with the discussion in Section 4, two cases are considered: first, at the same reduced frequency as applied in Section 3. For the second case, the excitation frequency is varied from low to high frequency levels within $0.125\le f/f_0\le 1.5$.

5.1. Exposed to a Constant Reduced Frequency

Figure 11 displays the vibration growth for the case at a constant reduced frequency of k = 0.079 (calculated at 93.68 m of the blade length). It can be seen that the range of the angle of attack is not altered too significantly for the low wind speed case up to 30 m/s. A larger difference between the rigid case and the flexible case starts to occur starting from the wind speed of 40 m/s. It can also be observed that the IAG dynamic stall model has the largest increase of the range up to 40 m/s, before it decreases at the wind speed of 50 m/s. On the other hand, the BL model consistently shows the smallest level of range change across all wind speeds. It is also notable that the case without dynamic stall model shows a consistent increase of the range with increasing wind speed.

5.2. Exposed to Varying Frequencies

The aeroelastic responses of the blade section at varying frequencies are discussed in the present section. To enable a clear comparison about the relative relationship between the excitation frequency (f) with the 1st edgewise natural frequency ($f_0$), the excitation frequency is reported as the normalized frequency ($f/f_0$), and this ranges from 0.125 to 1.5. To enable a comprehensive evaluations of the aeroelastic response, the results are plotted in two different ways: Figure 12 groups the plots based on the excitation frequency and the data are plotted with respect to the wind speed, meanwhile in Figure 13 the results are grouped based on the wind speed and the relationship is plotted against the excitation frequency itself.

Figure 12 clearly illustrates how the reduced frequency affects the vibration characteristics of the blade. Note that the same excitation frequency yields a different reduced frequency (k) for a different wind speed. At the excitation frequency level $f/f_0$ of 0.25 in Figure 12a, one can see that the IAG model generally shows the largest deviation compared to the rigid blade case before it is overtaken by the no dynamic stall case at 50 m/s. Here it can be seen that the BL model case generally shows the lowest deviation compared to the rigid case. By increasing the excitation frequency to half of the 1st edgewise natural frequency ($f/f_0$ = 0.25) in Figure 12b, the IAG model now shows the lowest deviation against the rigid case, followed by the BL model, Øye model and no dynamic stall model cases. This observation is mostly pronounced with increasing wind speed up to 40 m/s, before it decreases again when reaching 50 m/s case. The reason is associated with the reduction of the reduced frequency level becoming more significant than the blade flexibility effect itself.

A similar story is found at the non-dimensional excitation frequency of 0.75 in Figure 12c, where the discrepancy against the rigid case increases up to a certain wind speed. Here it is also observed that while all other cases generally show an increased level of discrepancy against the rigid blade case, the BL model is largely not affected by this particular frequency. Having observed the magnitude of the frequency in more detail, the applied frequency is actually close to the 1st flapwise natural frequency of the blade. Therefore, the flapwise mode is actually what is being excited.

Now looking closer at the actual resonance frequency of the 1st edgewise mode in Figure 12d, the IAG model shows the least deviation compared to the rigid blade case and is followed by no dynamic stall, Øye model and the BL model cases. This frequency excitation is considered the most challenging case because generally the edgewise mode is prone to negative damping and aeroelastic codes tend to underestimate the aerodynamic damping. The IAG model clearly shows an increased damping level under this condition.

Further increasing the excitation frequency in Figure 12e and Figure 12f, it can be seen that generally the results between the rigid and flexible cases do not deviate too much. This shows that the effects of resonance is not prominent and the blade remains in its stable conditions.

Figure 13 is aimed specifically to assess the characteristics of the dynamic stall model on the aeroelastic response for varying frequencies. At the smallest wind speed case of 10 m/s in Figure 13a, the excitation frequency largely has no impact on the range of angle of attack. The range only slightly increases for the no dynamic stall case and the IAG model case at the edgewise resonance frequency. By increasing the wind speed to 30 m/s, the blade flexibility plays a stronger role in the aeroelastic characteristics as shown in Figure 13b. It can be clearly observed that the excitation occurs at $f/f_0$ = 0.5 and $f/f_0$ = 1.0 due to the edgewise excitation while at $f/f_0$ = 0.75 due to flapwise excitation. For the BL model and Øye model cases, the blade section is excited strongly at the resonance frequency and this leads to the angle of attack variation to grow significantly.

Further increasing the wind speed to 50 m/s in Figure 13c, it is notable that most cases deviate from the rigid blade case for all frequency ranges, except for the BL model. However, it is clearly observed that the deviation is smaller than the one observed in Figure 13b. The explanation lies in the fact that the blade flexibility effects are prominent for this wind speed case, but the excitation reduced frequency itself is actually smaller than the the case with a wind speed of 30 m/s. This is also supported by the discussion made in Figure 12d.

6. Twist Angle Distribution Influence

It has been discussed above that dynamic stall model plays a significant role in the aeroelastic effects depending on the excitation frequency. Now, to gain deeper insights into design considerations, it is interesting to see if twist angles affect the characteristics. Therefore, the original twist angle distribution in Figure 3 is now applied. First of all, to evaluate how much twist angle changes the modal properties (keeping all other structural constraints the same), modal analysis of the blade is carried out and listed in Table 2. It can be seen that the modal frequencies are altered slightly, but not to a significant degree. The same inflow conditions are then imposed, but the blade set angle is adjusted to maintain the same angle of attack range at r = 93.68 m for the rigid blade case in comparison with the untwisted blade.

The aeroelastic response of the blade is presented in Figure 14 for different dynamic stall modeling approaches. For all cases, it can be seen that the addition of twist reduces the growth of the angle of attack range when the blade flexibility is considered. The effects are mostly prominent for the no dynamic stall case and the case employing the Øye model. This is followed by the IAG model and the weakest effect is observed for the BL model. The twist distribution changes the acting angle of attack and pushes the operating range further away from the stall. The twist angle distribution will likely become an important aspect to enhance large wind turbine blade stability during parked or idling conditions.

7. Conclusions and Outlook

The present studies have been conducted to evaluate the sensitivity of different dynamic stall models to edgewise vibrations occurring on large-flexible wind turbine blades. The investigations were done systematically by starting from the rigid response of the blade section with respect to an oscillating wind conditions. The flexibility of the blade is then applied, causing a dynamic response of the angle of attack. From the discussions, the following aspects can be summarized:

• It is clear that dynamic stall model plays a strong role in the aeroelastic characteristics. The angle of attack range starts to grow away from the rigid blade case when the excitation frequency is between $0.5\le f/f_0\le 1.0$. It is observed that in between this range, the flapwise mode is one souce of the excitation, though the strongest cause is related to the edgewise mode at $f/f_0$ = 0.5 and $f/f_0$ = 1.0. This can be used as a rough guidance when performing wind turbine design and running load cases especially in idling or parked conditions.
• The IAG model shows the best agreement with the pitching airfoil data, which represents the rigid conditions. Under the influence of blade flexibility, the IAG model is most sensitive to frequency excitation far from the edgewise resonance. However, when the excitation is close to the edgewise resonance frequency, the IAG model shows a higher aerodynamic damping which increases the blade stability. This explains the characteristics of this model in generating lower loads for the certain idling cases as discussed in [13,17].
• The characteristics of the IAG model for different excitation frequencies are also dependent upon the impulsive contributions. A higher reduced frequency typically increases the impulsive terms.
• Twist distribution is crucial in determining the aeroelastic stability of the blade. Adding a twist distribution (assuming all structural properties the same) only slightly alters the blade modal frequencies, but affects the operating range of the angle of attack.

Future studies might be aimed at investigating the dynamic stall sensitivity at higher excitation frequency levels to gain a better picture of the higher harmonic effects. Furthermore, incorporating the torsional effects will also be important to evaluate the coupling between the edgewise and torsional terms and how the aerodynamic damping actually influences the aeroelastic characteristics.

Appendix B.1. Exposed to Constant Reduced Frequency

Figure A4. Dynamic polar calculated at a constant reduced frequency of k = 0.079 at 10 m/s.

Figure A4. Dynamic polar calculated at a constant reduced frequency of k = 0.079 at 10 m/s.

Figure A5. Dynamic polar calculated at a constant reduced frequency of k = 0.079 at 30 m/s.

Figure A5. Dynamic polar calculated at a constant reduced frequency of k = 0.079 at 30 m/s.

Figure A6. Dynamic polar calculated at a constant reduced frequency of k = 0.079 at 50 m/s.

Figure A6. Dynamic polar calculated at a constant reduced frequency of k = 0.079 at 50 m/s.

Appendix B.2. Exposed to Edgewise Resonance Frequency

Figure A7. Dynamic polar calculated at a constant resonance frequency of f = 0.687 Hz at 10 m/s (at k = 0.598).

Figure A7. Dynamic polar calculated at a constant resonance frequency of f = 0.687 Hz at 10 m/s (at k = 0.598).

Figure A8. Dynamic polar calculated at a constant resonance frequency of f = 0.687 Hz at 30 m/s (at k = 0.199).

Figure A8. Dynamic polar calculated at a constant resonance frequency of f = 0.687 Hz at 30 m/s (at k = 0.199).

Figure A9. Dynamic polar calculated at a constant resonance frequency of f = 0.687 Hz at 50 m/s (at k = 0.12).

Figure A9. Dynamic polar calculated at a constant resonance frequency of f = 0.687 Hz at 50 m/s (at k = 0.12).