Various parameters used in the axisymmetric boundary layer.

## 1. Introduction

In the design process of improved rotor blades the need for accurate aerodynamic predictions is very important. During the last years a large effort has gone into developing CFD tools for prediction of wind turbine flows (Duque et al., 2003; Fletcher et al., 2009; Sørensen et al.,2002). However, there are still some unclear aspects for engineers regarding the practical application of CFD, such as computational domain size, reference system for different computational blocks, mesh quality and mesh number, turbulence, etc. Thus, in the design process and in the power curve prediction of wind turbines, the aerodynamic forces are calculated with some form of the blade element method (BEM) and its extensions to the three-dimensional wing aerodynamics. The results obtained by the standard methods are reasonably accurate in the proximity of the design point, but in stalled condition the BEM is known to underpredict the forces acting on the blades (Himmelskamp, 1947). The major disadvantage of these methods is that the airflow is reduced to axial and circumferential flow components (Glauert, 1963). Disregarding radial flow components present in the bottom of separated boundary layers of rotating wings leads to alteration of lift and drag characteristics of the individual blade sections with respect to the 2-D airfoils (Bjorck, 1995).

Airfoil characteristics of lift (_{L}) and drag (_{D}) coefficients are normally derived from two-dimensional (2-D) wind tunnel tests. However, after stall the flow over the inboard half of the rotor is strongly influenced by poorly understood 3-D effects (Banks & Gadd, 1963; Tangler, 2002). The 3-D effects yield delayed stall with _{L}higher than 2.0 near the blade root location and with correspondingly high _{D}* .*Now the design of constant speed, stall-regulated wind turbines lacks adequate theory for predicting their peak and post-peak power and loads.

During the development of stall-regulated wind turbines, there were several attempts to predict 3-D post-stall airfoil characteristics (Corrigan & Schlichting, 1994; Du & Schling, 1998; Snel et al., 1993), but these methods predicted insufficient delayed stall in the root region and tended to extend the delayed stall region too far out on the blade.

The present work aims at giving a conceptualization of the complex 3-D flow field on a rotor blade, where stall begins and how it progresses, driven by the needs to formulate a reasonably simple model that complements the 2-D airfoil characteristics used to predict rotor performance.

Understanding wind turbine aerodynamics (Hansen & Butterfield, 1993) in all working states is one of the key factors in making improved predictions of their performance. The flow field associated with wind turbines is highly three-dimensional and the transition to two-dimensional outboard separated flow is yet not well understood. A continued effort is necessary to improve the delayed stall modeling and bring it to a point where the prediction becomes acceptable. In the sequel, based on previous computed and measured results, a comprehensible model is devised to explain in physical terms the different phenomena that play a role and to clarify what can be modeled quantitatively in a scientific way and what is possible in an engineering environment.

In 1945 Himmelskamp (Himmelskamp, 1947) first described through measurements the 3-D and rotational effects on the boundary layer of a rotating propeller, finding lift coefficients much higher moving towards the rotation axis. Further experimental studies confirmed these early results, indicating in stall-delay and post-stalled higher lift coefficient values the main effects of rotation on wings. Measurements on wind turbine blades were performed by Ronstend (Ronsten, 1992), showing the differences between rotating and non-rotating pressure coefficients and aerodynamic loads, and by Tangler and Kocurek (Tangler & Kocurek, 1993), who combined results from measurements with the classical BEM method to properly compute lift and drag coefficients and the rotor power in stalled conditions.

The theoretical foundations for the analysis of the rotational effects on rotating blades come at the late 40’s with Sears (Sears, 1950), who derived a set of equations for the potential flow field around a cylindrical blade of infinite span in pure rotation. He stated that the spanwise component of velocity is dependent only upon the potential flow and it is independent of the span (the so-called independence principle). Then, Fogarty and Sears (Fogarty & Sears, 1950) extended the former study to the potential flow around a rotating and advancing blade. They confirmed that, for a cylindrical blade advancing like a propeller, the tangential and axial velocity components are the same as in the 2-D motion at the local relative speed and incidence. A more comprehensive work was made once more by Fogarty (Fogarty, 1951), consisting of numerical computations on the laminar boundary layer of a rotating plate and blade with thickness. Here he showed that the separation line is unaffected by rotation and that the spanwise velocities in the boundary layer appeared small compared to the chordwise, and no large effects of rotation were observed in contrast to (Himmelskamp, 1947). A theoretical analysis done by Banks and Gadd (Banks &Gadd, 1963), focused on demonstrating how rotation delays laminar separation. They found that the separation point is postponed due to rotation, and for extreme inboard stations the boundary layer is completely stabilized against separation.

In the NASA report done by McCroskey and Dwyer (McCroskey & Dwyer, 1969), the so-called secondary effects in the laminar incompressible boundary layer of propeller and helicopter rotor blades are widely studied, by means of a combined, numerical and analytical approach. They showed that approaching the rotational axis, the Coriolis force in the cross flow direction becomes more important. On the other hand the centrifugal pumping effect is much weaker than was generally thought before, but its contribution increases particularly in the region of the separated flow. The last two decades have known the rising of computational fluid dynamics and the study of the boundary layer on rotating blade has often been carried on through a numerical approach. Sørensen (Sørensen, 1986) numerically solved the 3-D equations of the boundary layer on a rotating surface, using a viscous-inviscid interaction model. In his results the position of the separation line still appears the same as for 2-D predictions, but where separations are more pronounced a larger difference between the lift coefficient calculated for the 2-D and 3-D case is noticed. A quasi 3-D approach, based on the viscous-inviscid interaction method, was introduced by Snel et al (Snel et al., 1993) and results were compared with measurements. They proposed a semi-empirical law for the correction of the 2-D lift curve, identifying the local chord to radii (* c*/

*) ratio of the blade section as the main parameter of influence. This result has been confirmed by Shen and Sørensen (Shen & Sørensen, 1999), and by Chaviaropoulos and Hansen (Chaviaropoulos & Hansen, 2000), who performed airfoil computations applying a quasi 3-D Navier-Stokes model, based on the streamfunction-vorticity formulation and respectively a primitive variables form. Du and Seling (Du & Selling, 2000), and Dumitrescu and Cardoş (Dumitrescu & Cardoş, 2010), investigated the effects of rotation on blade boundary layers by solving the 3-D integral boundary-layer equations with the assumed velocity profiles and a closure model. Dumitrescu and Cardoş studying the boundary layer behavior very close to the rotation center (*r

*/*r

*< 1) stated that the stall delay depends strongly on the leading edge separation bubbles formed on inboard blade segment due to a suction effect at the root area of the blade (Dumitrescu & Cardoş, 2009).*c

Recently, with the advent of the supercomputer, computational fluid dynamics (CFD) tools have been employed to investigate the stall-delay for wind turbines (Duque et al., 2003; Fletcher et al., 2009; Sørensen et al., 2002). These calculations in addition to Narramore and Vermeland’s results (Narramore & Vermeland, 1992) showed that the 3-D rotational effect was particularly pronounced for the inboard sections, but the genesis of this phenomenon is a problem still open. The present concern aims at giving explanations in physics terms of the different features widely-observed experimentally and computationally in wind turbine flow.

## 2. Flow at low wind speeds

At low wind speed conditions, i.e. at high speed ratios (TSR>3), the visualization of the computed flow indicates that the flow is well-behaved and attached over much of the rotor, Fig. 1 . Figure 1 shows the separated area and radial flow on the suction side of a commercial blade with 40 m length at the design tip speed ratio (TSR=5); the secondary flow is strongest at approximately 0.17* R*and reaches up to 0.31

*, where*R

*is the rotor radius. The local air velocity relative to a rotor blade consists of free-wind velocity*R

V

_{w}defined as the wind speed if there were no rotor present, that due to the blade motion Ω

_{ b}

*and the wake induced velocities; at high TSR, a weak wake (Glauert, 1963) occurs and its rotational induction velocity can be neglected. Wind turbine blade sections can operate in two main flow regimes depending on the size of the rotation parameter defined as the ratio of wind velocity to the local tangential velocity*r

Therefore, at low wind speed, the main rotational effect is due to the Coriolis force which delays the occurrence of separation to a point further downstream towards the trailing edge, and by this the suction pressures move towards higher levels as * r*/

*decreases. The pumping effect is much weaker than was generally thought before.*c

The pressure field created by the presence of the turbine is related to the incoming flow field around the blade, taken as being composed of the free wind velocity and the so-called induction velocity due to the rotor and its wake. Thus, the incoming field results from a weak interaction between two different flows: one axial and the other rotational

the radial independence principle is applied to flow effects, i.e. induction velocities used at a certain radial station depend only on the local aerodynamic forces at that same station;

the mathematical description of the air flow over the blades is based on the 2-D flow potential independent of the span, and on corrections for viscosity and 3-D rotational effects.

These assumptions suitable to BEM methods reduce the complexity of the problem by an order of magnitude yielding reliable results for the local forces and the overall torque in the proximity of the design point, at high tip speed ratios. In order to estimate the 3-D rotational effects, usually neglected in the traditional BEM model, the flow around a hypothetic blade with prestall/stall incidence and chord constant along the whole span is considered in the sequel.

### 2.1. Representation of flow elements

The set of equations including a simplified form of the inviscid flow and the full three-dimensional boundary-layer equations are used to identify the influence of the3-D rotational effects at low wind speeds.

* A. Inviscid flow*. In order to find the velocity at the airfoil surface in absence of viscous effects, the reference velocity at a point on a rotating wind turbine blade is

where _{w}is the wind speed, * R*is the radius of the turbine and

*is the axial induced velocity interference, a function of the speed ratio λ (Burton et al., 2001). Starting from the idea of Fogarty and Sears (Fogarty & Sears, 1950), an inviscid edge velocity can be calculated as*a

where * U*,

*and*V

*represent the velocity components in the cylindrical coordinate system*W

_{ b}(Fig. 3).

*can be derived from the local 2-D velocity potential. However, this spanwise component is very small and thus neglected in the present work. The potential edge velocity components can be approached as*V

where the non-dimensional velocity

where * x*/

*is the airfoil abscissa by airfoil chord and subscripts*c

*and*U

*indicate properties on the upper surface and lower surface respectively. The scale parameter*L

C

_{p}and is only dependent on the local aerodynamic field: the airfoil shape, incidence and Reynolds number. This one-parameter pressure distribution family ranges a large gamut from the prestall distributions, values of

Because a wind turbine can operate a long time at low tip speed ratios, the inboard regions of the blade are stalled, and there leading-edge separation bubbles can occur (Dumitrescu & Cardoş, 2010). Behind the mild separation it is assumed that the flow relaxes with the vanishing skin friction. Therefore, the following flow is applied to solve

where * H*are, respectively, the momentum thickness and the boundary-layer shape parameter in the streamwise direction.

Figures 5 and 6 illustrate such chordwise inviscid velocity distributions and the corresponding variations of the peripheral skin-friction coefficient _{fx}and the boundary layer shape parameter * H*, calculated for various values of

*/*r

*= ∞ (2-D). Laminar separation takes place at*c

x

_{t}is the point that corresponds to the minimum skin-friction, which sometimes can be the point of laminar separation.

Fig 6. Variation of a) the peripheral skin-friction coefficient; b) the boundary-layer shape parameter.

* Momentum integral equations.*The flow in the boundary layer on a rotating blade in attached as well as in stalled conditions is represented using the integral formulation developed in (Dumitrescu & Cardoş, 2010) for analyzing separated and reattaching turbulent flows involving leading-edge separation bubbles. The 3-D incompressible steady boundary layer equations are written in the cylindrical coordinate system

*the axial and*z

*the radial (spanwise) direction. The equations used early (Cebeci & Cousteix, 1999) are shown in their conservative laminar form:*r

Continuity

Momentum, θ component:

Momentum, * r*component:

where * u*,

*and*v

*stand for velocity components in θ,*w

*and*r

*directions, respectively,*z

*for the local radius measured from the centre of rotation, ρ for the fluid density, ν for the kinematic viscosity, and*r

*is a pressure like term including the centrifugal effect*P

with * p*denoting the static pressure.

The equations (7)-(9) are integrated with respect to z (normal to blade) from 0 to δ (boundary layer thickness) and the integral forms of equations are obtained as

where (_{e}, _{e}=0) are the inviscid freestream velocity components and

An order of magnitude analysis shows that

where * H*is the local boundary layer shape factor.

The Mager cross flow profile is assumed (Mager, 1951)

where

Equations (11) and (12) can now be written in terms of the parameters* H*and

The skin-friction relation for flows with pressure gradients and rotation effects is based on the experimental data for a turbulent boundary layer in a rotating channel (Lakshminarayana & Govindan, 1981),

This correlation is a modified version of the correlation developed by Ludwieg and Tillmann (Ludwieg & Tillman, 1949), which includes the effect of rotation. In this relation

The pressure gradients cause large changes in velocity profiles and, consequently, in the shape parameter * H*. The variation of

*cannot be neglected and an additional equation is required. Out of the available auxiliary equations, only the energy integral equation and the entrainment equation have been suitable for the turbulent boundary layers.*H

The entrainment equation is chosen to model the rotor boundary layer, which in the coordinate system used here can be written as

where the entrainment coefficient

The entrainment function

Also, the similarity solutions have shown that * H*. This relationship, which results from a best fit to experimental data (Lock & Williams, 1987) is

*< 4.(21)*H

Then, it is assumed that the variation of the entrainment rate with

Equation (19) is written in a form similar to Eqs. (15) and (16)

Equations (15), (16) and (22) are to be solved for* H*(

*Eq. (21)) with the prescribed boundary conditions. At the leading edge,*H,

*.*H

### 2.2. Application

A hypothetic blade, with constant incidence and chord along the span is chosen to analyze the rotational effects on the boundary layer of a wind turbine blade. In this example, an external flow with the distribution of velocity on the surface of blade given by Eqs. (5) and (6) is imposed.

The curves plotted in Fig. 7 are calculated on the prestall (^{6}. Relationships between the separation points and the blade radius are shown in Fig. 7, in which two apparent trends can be seen. Firstly, as compared to the 2-D prestall condition (* r*/

*= ∞), the separation point is postponed because of the effect of Coriolis force, which acts as a favorable pressure gradient tending to delay the separation to a point further downstream towards the trailing edge. It is shown that, as the radial location*c

*/*r

*decreases, the separation point moves towards the trailing edge, and consequently the suction side distribution of the pressure coefficient moves towards lower levels (Chaviaropoulos & Hansen, 2000), Fig. 8. The drop of the pressure coefficient along the suction side can be related to the separated area on the blades, which reduces in comparison with the 2-D (*c

*/*r

*= ∞).*c

Secondly, in the stall incidence condition, leading edge separation bubbles are formed on the upper surface of the blade, at the root area, which delays the occurrence of massive separation. Beginning at the root of the blade, the bubble continues to stretch towards the trailing edge and at approximately * r*/

*= 4 where the bubble stretches all the way to the trailing edge, the flow separates over the whole airfoil (Dumitrescu & Cardoş, 2010) In the next section it is shown that there is a suction effect at the hub up to the mid-span which reduces the leading-edge bubble volume and produces a significant pressure drop along the suction side of the airfoils increasing, thus, the loading of the blade.*c

## 3. Flow at high wind speeds

At high wind speeds conditions (

This behavior, which is characterized by significantly increased lift coefficients as compared to the corresponding 2-D case, and by a delay of the occurrence of flow separation to higher angles of attack, is caused by the rotational augmentation at the root area of the blade. The knowledge extracted from experimental and computational visualizations of these rotational effects on the flow field along the blade can be used to develop a model for their prediction or the so-called stall-delay phenomenon.

The wind turbine blade sections, often operate with stall at low tip speed ratios and can undergo a delayed stall phenomenon on the inner part of span. As the wind velocity increases, the inboard regions of the blades are stalled, having

a form that suggest an outer rotational flow with angular velocity higher than that of the constant-speed rotor,

In contrast with the Glauert weak wake (Glauert, 1963), the strong potential vortex-like wake is responsible for all the 3-D rotational effects, affecting the boundary layer on the turbine blades once the phenomenon of stall is initiated. The enhanced rotational outer flow and the bottom boundary layer are the cause of the high-lift effect of the airfoils at high angle of attack, i.e. the stall delay.

### 3.1. Boundary layer beneath a Rankine-like vortex

The axisymmetric angular velocity (Ekman) boundary layer developing on the suction side of the rotor disk is used to illustrate the contribution of the suction effect of a strong wake on the chordwise surface pressure distribution and volume flux in the radial direction. Thus, consider an axisymmetric rotating flow of a viscous incompressible fluid over a disk of radius _{0} ( Fig. 11 ).

The angular velocity of the disk Ω_{ b}is constant and that of the fluid far away from the disk Ω_{ v}may be an arbitrary function of the radius * r*, provided it is stable, i.e.

where _{h}is the radius of forced vortex core.

This vortex model avoids the discontinuity of the velocity derivatives at the point of transition from free to forced modes and gives the possibility of the description of the flow near the hub.

For the inertial system of cylindrical coordinates

The boundary conditions at the surface of the disk and the edge of boundary layer are

The outer flow radial momentum equation becomes then

where _{ps}is the suction pressure coefficient normalized by the pressure at the vortex center, and the subscripts * c*and 0 indicate properties at the vortex center and the edge of disk, respectively.

The tangential and radial momentum-integral equations may be derived from equations (26) and (27) for radial and angular velocity profiles which satisfy the boundary conditions on the disk surface and the smoothness requirements at the edge of the boundary layer

Outside the boundary layer, the tangential velocity _{δ} and the radial volume flux

where λ_{1}, λ_{2}, λ_{3} and λ_{4} are profile form parameters defined in Table 1 (Rott & Lewellen, 1966)

Flow | ||||||||

Laminar | 6.75 | 2.5 | 1.372 | 0.467 | 0.667 | 2 | 12 | |

Turbulent | 1.69 | 4.93 | 1.63 | 0.222 | 0.250 | 0.0225 | 0.0513 |

The Blasius shear laws for three-dimensional flow (Schlichting, 1979) cover both the laminar and turbulent case

where * u*and following Blasius,

*= 0.0225 (turbulent flow).*C

This is the final result of the analytical formulation and the equations (32) and (33) must be integrated numerically for a particular circulation function Γ_{ v}(Eq. (24)) with an appropriate choice of the boundary layer form parameters.

The onset of vortical flows like tornadoes can be visible on the constant-rotational speed rotor at low speed ratios (λ < 3.0). At the root area of the blade these flows generally behave like a vacuum pump which is featured by volume flow rate (* Q*) and suction pressure (

C

_{ps}).

Figures 12 and 13 show such characteristics induced by the wake modeled as a Rankine-like vortex with its axis normal to the rotor disk for different tip speed ratios. The suction characteristics show that the effects of the disk defined by

On the other hand, large radial volume fluxes close to the hub along with the previous results (Fig. 7a) indicate a 3-D attached boundary layer with small leading-edge separation bubbles at high angles of attack.

Therefore, at low tip speed ratios and close to the hub, there is essentially a 3-D attached flow field and a theoretical description in this region can be performed by viscous-inviscid interactive flow models.

### 3.2. Blade-wake interaction

Flow at low tip speed ratios is dominated by a vortex-like wake which triggers the stall delay phenomenon on the inboard half-rotor disk. The dual working mode of the rotor blades as a stationary disk in a fast rotating flow near the hub-sucking mode, and as the rotating disk in the outboard separated flow-pumping mode, are clearly visible from their pressure distributions when the wind turbine operates at

On the upper side of the inner blade sections there are two suction pressure fields: one is azimuthally uniform generated by circulation augmentation of wake and the other is the aerodynamic non uniform pressure field with a suction spike at the leading-edge (Fig. 4). These two pressure fields interact and give rise to the radial effects which contribute to stall delay, as well as to a higher lift coefficient at the root area of the blade. There are two main effects: superimposed pressure fields of wake ( Fig. 12 ) and of blade (Fig. 4), which decrease the suction spikes at the leading edge until they cancel, and a radial flow effect induced by the volume flow merger in the boundary layers of wake and blade, resulting in velocity profile skewing and boundary-layer energizing (no separation). Assuming the conservation of circulation once the phenomenon of stall is initiated (here at

As we will show in the next section, the phenomenon of stall-delay can be described as a three step process: generation of the strong wake by the clustering of shed vortices from the inboard blade segments at

### 3.3. Application

The sucking suction effect is affected by two key non-dimensional parameters

distribution for the previous hypothetic blade, with constant post-stall incidence and chord along the span* S*-separation and

*- reattachment.*R

## 4. Model for delayed stall regime

Stall occurs at an airfoil when flow separates at high angles of attack, typically>15^{0}, beyond this angle the 2-D lift coefficient drops significantly. It has long been known (Himmelskamp, 1947; Ronsten, 1992; Snel et al., 1993), however, that the observed power curve under conditions where most parts of the blade are stalled is consistent with significantly higher lift coefficients. The lack of a conceptual model for the complex 3-D flow field on the rotor blade, where stall begins, and how it progresses, has hindered the finding of an unanimously accepted solution. This section aims at giving a better understanding of the delayed stall events.

To interest an engineer in stall-delay one might cite the fact that the inboard regions of the blades are stalled for much of their operational time. Thus predicting blade loads and the power output during stall is very important in making good predictions of wind turbine performance. With a physicist one could reason that rotational stall admits interesting transient states from 3-D inboard delayed stall in partially separated flow to 2-D outboard massively separated flow. Here, the detailed physical understanding of the viscous flow surrounding the hub area is dominated by the interaction between the vorticity and pressure fields and their dynamics. However, all the existent models have been developed more or less intuitively without a basic fluid dynamics support. Six different models mostly used to correct the airfoil characteristics for stall delay, including a wide range of different assumptions, were recently analyzed (Breton et al., 2008). The conclusion was that none of the six models studied correctly represented the flow physics, and that this was ultimately responsible for their lack of generality.

### 4.1. New physics-based model

* Flow Topology of Delayed Stall.*There are two main flow regimes pertinent to wind turbine rotors: one is a flow close to rotating disk in axial flow at low velocity, practically at rest, (λ>3, Fig. 2), and another is a rotating flow over a stationary or slowly rotating disk (Fig. 9). The last can produce delayed stall on the wind turbine blade sections operating at low tip speed ratios (λ≤3).

The key to the stall-delay phenomenon is the rotational flow surrounding the root area of the blade, i.e. the wake rotation. At low tip speed ratios (or high wind speeds), just as a vortex is shed from each tip blade, a vortex is also shed beginning from the section

There are two common characteristics which most strong potential vortices exhibit, like those produced at the blade root; the first is that the vortex velocity field above the surface boundary layer is always dominated by the tangential velocity component, while the second shows that the balance of forces between the flow outside the boundary layer and vortex core is governed by a balance between the centrifugal force and radial pressure gradient

* Modeling of Stall-Delay*. The acquisition of the enhanced rotational flow (tangential component of velocity) by the increase in the kinetic energy of wind, associated with a strong suction of the air at the hub is governed by conservation of angular momentum, like the potential vortex, and conservation of energy:

Then, as the blade section moves through the air a circulation

a starting vortex

The wake’s circulation, which is constant up to the hub induces a suction pressure along the upper side of blade (Fig. 16) stabilizing the boundary layer against the separation at post-stall angles of attack (angles of 30 degrees exist currently at the blade root), followed by forming of a leading-edge separation bubble. The bubble at the innermost sectional airfoil has almost no effect on integrated loads, because it is never more than a few percent of the chord in length, and the airfoil acts as in an ideal fluid flow.

The phenomenon of stall-delay can be described as a three step process: rise of strong wake by clustering of vortices shed from the inboard blade span at

The specific mechanism for circulation (or lift) decay is not presently known. However, a measure of the degree of radial instability of the boundary layer and implicitly the lift decay can be obtained from the sucked volume flow rate (Fig. 13) into the boundary layer beneath the Rankine vortex-like wake. Thus, the maximum sucked flow rate would correspond to the minimum volume of the separation bubble and_{L}decay can be supposed to follow the flow rate decrease,

boundary layer, namely larger values of the flow rate signify smaller bubble volumes and thinner boundary layers for reducing tip speed ratios. Further, it can be easily found that a thicker boundary layer with distributed vorticity (

Point | Flow Structure | Forces |

1 | The attached 3-D boundary layer with leading-edge separation bubble | Exceed 2-D CLmax, extrapolate linear regime |

2 | The separation bubble is sucked and its volume reduces | Maximum lift, CL,INV at |

3 | Thick skewed 3-D boundary layer with stretching bubble | Gradual decay of lift for |

4 | The bubble breaks away from the trailing edge and forms a free shear layer | Readjust to stall regime |

The various stages of the delayed stall process are summarized schematically in Fig. 17. and Table 2. Stage 1 represents the delay in the onset of separation in response to reduction in adverse pressure gradients produced by the influence of the strong vortex wake at the root. The phenomenon is initiated at

Stage 2 of the delayed stall process involves the suction pressure redistribution of the start pressure distribution with the conservation condition of its integrated load and the simultaneous move of pressure center to the midchord; as the suction/tip speed ratio increases/decreases the volume of the leading-edge separation bubble is reduced. The start pressure distribution is known from the inviscid aerodynamic loading solution over the chord; the pressure redistribution continues up to no suction peak exists at the leading-edge and even more no suction surface pressure gradient exists. On the other hand, the pressure redistribution is associated with a vorticity concentration process, which influences the stability of boundary layer.

Stage 3 is the sectional circulation (lift) decay which occurs at greater radii than the start location, where the separation bubble surface stretches all the way to the trailing edge. Beyond that zero flow rate location the flow is separated over the whole airfoil so that full stall occurs. By reason of the boundary layer stability it is inferred that the stall-delay effects are present up to a 50% of span for all the regimes.

Stage 4 is the relaxation to stall regime, which begins just as the integrated normal load of 2-D stall, C_{NS}, is reached (end of delayed stall) and continues as far as the recovery of the specific pressure distribution (suction peak at the leading edge) is achieved.

### 4.2. Application

The combined experimental rotor (Schepers et al. 1997) of the Renewable Energy Laboratory (NREL) is chosen to illustrate the stall-delay model proposed for a wind turbine blade. The rotor uses the NREL S 809 airfoil and a simpler BEM method. In this example, the one parameter pressure distribution on the upper surface of blade at the root with * R*. The curves plotted in Fig. 19 are calculated with these conditions: the onset of stall-delay at

The comparison of predicted and available measured _{N}, _{T}and * P*are shown in Fig. 19 (a, b, c), for wind speeds 13 m/s, 15 m/s and 25 m/s. Excepting predicted

C

_{T}at the blade root, the agreement between predictions and measurements is reasonably good. Most likely, at the root area of the blade there are computational and experimental uncertainness, and better more exist uncertainness for predictions. However, the present stall-delay model is the most comprehensive one able to capture much of the actual physics of flow.

## 5. Chapter review

The phenomenon of stall-delay has been shown to be an important consideration in wind turbine design because, in the presence of turbulence, it ultimately generates the loads with the highest peak-to-peak fatigue cycles for both blade and rotor shaft bending (Tangler, 2004). It has been shown that stall-delay is characterized by a favorable delay in onset of flow separation to higher angles of attack for the inboard regions of the blades. This is followed by the less favorable event of leading-edge separation bubble stretching. As long as the bubble surface stays over airfoil, it acts to decay the lift gradually. Then, the leading-edge separation occurs. The stall-delay begins at the hub and crosses varied partially separated flow regions over the inboard half-rotor disk, depending on the wind speed. Therefore, the consideration of stall-delay phenomenon represents a necessary refinement in the rotor design process which will more accurately define the performance and power control at high wind speeds.

While the prediction of the conditions for stall-delay onset and their subsequent effects clearly forms an essential part of any rotor design process, it has been shown that this is a problem not yet fully understood, nor easily predicted. For engineering analyses, the modeling of stall-delay still remains a particularly challenging problem. This is mainly because of the need to balance physical accuracy with computational efficiency and/or the need to formulate a model of stall-delay in a particular mathematical form. To this end, a number of semi-empirical models have been developed for use in wind turbine design work. A discussion of these semi-empirical methods has been recently presented, along with the conclusion that none of the six most known models correctly represented the physics of flow (Breton et al., 2008). Generally, predictions are good when measurements are available for validation or empirical refinement of the model, but their capabilities for general blade parameters (twist, taper and airfoil shape) are less certain.

A conceptualization of the complex 3-D flow field on a rotor blade, where stall begins, and how it progress, has been proposed, along with a reasonably simple model that complements the 2-D airfoil characteristics used to predict rotor performance. This concept considers that the 3-D flow field near the hub behaves like the rotational flow over a stationary disk, inducing a strong suction pressure towards the center. Here, the centrifugal force has a strong stabilizing effect on the boundary layer, in which the rotor is acting as a vacuum pump on the separated volume of bubble. The rough model has shown encouraging results, but some uncertainties remain in the prediction of the sectional circulation decay and also in the proper validation of predictions with more measured and CFD computed airloads on the rotor. The more correct physical description provided by this present model should give more coherence and consistence of future experimental and computational studies on the stall-delay phenomenon.