Various parameters used in the axisymmetric boundary layer.
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 (
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 (
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,
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
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
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 . In such a weak interaction flow, the basic assumptions made are:
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.
where the non-dimensional velocity could be obtained by a viscous-inviscid interaction procedure for flow past a 2-D airfoil (Drela, 1989). Since the primary concern of the present work is to investigate the rotational 3-D effects on prestalled and stalled blades by means of the boundary layer method, the inviscid pressure distribution is simply considered as
is the local pressure coefficient,
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
Figures 5 and 6 illustrate such chordwise inviscid velocity distributions and the corresponding variations of the peripheral skin-friction coefficient
Fig 6. Variation of a) the peripheral skin-friction coefficient; b) the boundary-layer shape parameter.
Momentum, θ component:
An order of magnitude analysis shows that and , and thereby these terms can be neglected as a first approximation. The widely used Pohlhausen (Schlichting, 1979) and Mager (Mager, 1951) velocity profiles and the associated closure relations are introduced to solve the laminar integral boundary layer equations (Dumitrescu et al., 2007). For the turbulent boundary layer a power law-type of velocity profile is assumed for the mean stream velocity profile
The Mager cross flow profile is assumed (Mager, 1951)
where is the limiting streamline parameter (tan ).
Equations (11) and (12) can now be written in terms of the parameters
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 is an empiric constant (a value of 0.52 is used), is the Reynolds number, based on the streamwise velocity at the edge of the boundary layer and the streamwise momentum thickness , and is the distance between the leading edge and the transition point along the streamwise direction (the laminar separation point is used).
The pressure gradients cause large changes in velocity profiles and, consequently, in the shape parameter
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 is a function of the factor . represents the volume flow rate per unit area through the surface and is the rate of entrainment of inviscid external flow into the boundary layer.
The entrainment function for 3-D flow is not yet available and hence, the correlation for 2-D flow is used (Head, 1958):
Also, the similarity solutions have shown that
is a function of the streamwise boundary-layer shape parameter
Then, it is assumed that the variation of the entrainment rate with follows the same relationship for three-dimensional flows.
Equations (15), (16) and (22) are to be solved for
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
are calculated on the prestall (
) and stall (
) conditions, with a Reynolds number of 106. 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 (
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
3. Flow at high wind speeds
At high wind speeds conditions ( ), though the stall incidences are exceeded on the whole blade, the flow is partially separated on the inboard half-span of the blade and massively separated on the outboard half-span, Fig. 9 .
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 for much of their operational time. The flow is neither two-dimensional nor steady, and is affected by rotation. As the root is approached the flow causes high values of tangential velocity and it behaves rather like the rotating flow over a stationary disk, Fig. 10 a), b). This is the result of a strong nonlinear interaction between the high wind speed and the rotational flow induced by a constant-speed rotor, which can be written as
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
The angular velocity of the disk Ω
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.
The boundary conditions at the surface of the disk and the edge of boundary layer are
The outer flow radial momentum equation becomes then and the pressure in the outer flow is
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
is assumed to be some given function, Eq. (24), while the radial velocity vanishes. The variables are the boundary-layer thickness
The Blasius shear laws for three-dimensional flow (Schlichting, 1979) cover both the laminar and turbulent case
is equal to the maximum value of
This is the final result of the analytical formulation and the equations (32) and (33) must be integrated numerically for a particular circulation function Γ
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 (
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 decrease rapidly with the radius towards the edge of the disk ( ) and have maximum values near the hub (here ).
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 ) the wake-blade pressure interaction for the increased wind speed ( ) is a pressure redistribution process of pressures ( ) along the chord at constant circulation (lift). The rule of pressure redistribution is sketched as is shown in Fig. 14.
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 , chordwise pressure redistribution at the constant value of circulation at the start regime ( ), for , followed by circulation decay for increasing radius. The pressure redistribution process is companied by the move of the pressure center to the midchord.
The sucking suction effect is affected by two key non-dimensional parameters . Figures 15a and 15b illustrate the influence of these parameters on the chordwise pressureFig 14. The rule of the suction pressure with superimposed rotation effect in a normalized plan:
distribution for the previous hypothetic blade, with constant post-stall incidence and chord along the span
. The figures also indicate the boundary layer state on the upper surface:
4. Model for delayed stall regime
Stall occurs at an airfoil when flow separates at high angles of attack, typically>150, 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
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 of each blade. Then the blade vortices from span will each be a line vortex running axially to the center of the rotor. The direction of rotation of all span vortices will be the same, forming a core (or root) vortex of total strength . The root vortex is primarily responsible for inducing the rotational 3-D effects on blades. The acquisition of the tangential component of air velocity is compensated for by a fall in the static pressure (suction) in the wake.
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 . In this working mode the centrifugal force has a strong stabilizing effect on the boundary layer through the favorable radial pressure gradient which occurs. In contrast with the pumping - working mode at high tip speed ratios, now the centrifugal forces produce a centerwise suction effect, resulting mainly in boundary layer stabilizing against the separation. On the other hand, in the presence of a solid surface (blade), this balance of forces is disrupted by the friction retardation of the tangential flow close to the boundary. Under such a condition, the radial pressure gradient proceeds to drive the retarded boundary flow along the surface towards the center, resulting in large radial velocities residing close to the surface. Above the boundary layer where the cyclostrophic balance is still intact, the radial velocity inevitably falls to zero. This gives rise to an inflexion point in the radial velocity profile. At some location close to the axis, this inward flowing air moves away from the solid boundary and effuses into a vortex core. The radial pressure gradient has a very strong stabilizing effect on the boundary layer and acts to revert it to its reattaching, and even more to its laminar state upstream of the effusing core.
Then, as the blade section moves through the air a circulation develops around it. In order to comply with Kelvin’s theorem,
a starting vortex must exist such that the total circulation around a line that surrounds both the airfoil and the wake remains unchanged. Since the circulation of the sectional airfoil, where the phenomenon of stall was initiated at , is practically the blade circulation due to rotation, (at the sectional airfoil is stalled for ), then .
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 , inviscid (pressure) load redistribution along the airfoil chord at constant circulation or of the start regime ( ) for increased suction pressure at , followed by spanwise circulation or lift decay involving the stretching of separation bubble surface all the way to the trailing edge for increasing radius; beyond that blade section ( ) the flow is separated over the whole airfoil and the leading-edge stall occurs.
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
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 ( ) is more stable, while a thinner boundary layer with more concentrated vorticity ( ) is less stable. Since the lift decay law depending on the stability of the boundary layer, is different from that suggested by Fig. 13, tending to be slower for higher λ and steeper for lower λ, an average decay corresponding to is assumed. Therefore, the decay is assumed to be terminated for all regimes at midspan, i.e. the end of the flow regime for . This conclusion is consistent with the result found in (Dumitrescu & Cardoş, 2010), which indicates the full separation at .
|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, followed by suction pressure redistribution and movement of pressure center to midchord; no correction to drag|
|3||Thick skewed 3-D boundary layer with stretching bubble||Gradual decay of lift for; no correction to drag|
|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 , where the suction pressure equalizes to the kinetic energy of wake (Fig. 16) resulting mainly in boundary layer reattaching and the rise of a leading edge separation bubble, a stable focus as topological entity.
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, CNS, 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.
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
is imposed and the hub is at approximately 0.20
The comparison of predicted and available measured
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.
Banks W. Gadd G. 1963Delayed effect of rotation on laminar separation, 1 4
Bjorck A. 1995, FFA TN- 1995 31.
Breton S. Ph Coton F. N. Moe G. 2008A study of rotational effects and different stall delay models using a prescribed wake vortex scheme and NREL phase VI experiment data, , 11 459 482.
Burton T. Sharpe D. Jenkins N. Bossanyi E. 2001, Wiley, Chichester.
Cebeci T. Cousteix J. 1999, Springer-Verlag, Berlin.
Chaviaropoulos P. K. Hansen M. O. L. 2000Investigating three-dimensional and rotational effects on wind turbines blades by means of quasi-3-D Navier-Stokes solver, , 122 330 336.
Corrigan J. J. Schlichting J. J. 1994Empirical model for stall delay due to rotation, , San Francisco, CA, 8.4-(1-15).
Corten G. P. 2001Flow separation on wind turbine blades, , University of Utrecht.
Drela M. 1989XFOIL: An analysis and design systems for low Reynolds number airfoils, , University of Notre Dame, June 1989.
Du Z. Schling M. S. 1998A 3-D stall-delay model for horizontal axis wind turbines performance prediction, -98-0021.
Du Z. Selling M. S. 2000The effect of rotation on the boundary layer of a wind turbine blade, 20 167 181.
Duque E. P. N. Burkland M. D. Johnson W. 2003Navier-Stokes and comprehensive analysis performance predictions of the NREL phase VI experiment, , 125 4 457 467.
Dumitrescu H. Cardoş V. 2004Rotational effects on the boundary layer flow in wind turbines, , 42 2 408 411.
Dumitrescu H. Cardoş V. Dumitrache A. 2007Modelling of inboard stall delay due to rotation, The Science of Making Torque from Wind, , 75IOP Publication doi:10.1088/1742-6569/1/012022.
Dumitrescu H. Cardoş V. 2009Inboard boundary layer state on wind turbine blades, , 89 3 163 173.
Dumitrescu H. Cardoş V. 2010Analysis of leading-edge separation bubbles on rotating blades, , 47 5 1815 1819.
Fletcher T. M. Brown R. E. Kim D. H. Kvon O. J. 2009Predicting wind turbine blade loads using vorticity transport and RANS methodologies, , Marseille, March.
Fogarty L. E. Sears W. R. 1950Potential flow around a rotating, advancing cylindrical blade, , 17 9 599
Fogarty L. E. 1951The laminar boundary layer on a rotating blade, , 18 4 247 252.
Glauert H. 1963Windmills and fans, in , 14edited by Durand W.F., Dover ed.
Hansen A. C. Butterfield C. P. 1993Aerodynamics of horizontal-axis wind turbines, , 25 115 149.
Head R. M. 1958Entrainment in turbulent boundary layers, , Sept. 1958.
Himmelskamp H. 1947Profile investigations on a rotating airscrew, (PhD dissertation, Götingen, 1945), 832
Lakshminarayana B. Govindan T. R. 1981Analysis of turbulent boundary layer on cascade and rotor blades of turbomachinery, 19 10 1333 1341.
Lock R. C. Williams B. R. 1987Viscous-inviscid interactions in external aerodynamics, , 24 2 51 171.
Ludwieg H. Tillman W. 1949Investigation of the wall shearing stress in turbulent boundary layers, 1285.
Mager A. 1951Generalization of boundary layer momentum-integral equations to three-dimensional flows including these of rotating systems, 1067.
Mc Croskey W. J. Dwyer H. A. 1969Methods of analyzing propeller and rotor boundary layers with crossflow. NASA SP-228, 473 514.
Narramore J. C. Vermeland R. 1992Navier-Stokes calculations of inboard stall delay due to rotation, , 29 1
Ronsten G. 1992Static pressure measurements on a rotating and non-rotating 2.375 m wind turbine blade. Comparison with 2-D calculations, 39 1-3, 105 118.
Rott N. Lewellen W. S. 1966Boundary layers and their interactions in rotating flows, in , 7 111 144.
Schepers J. G. Brand A. J. Bruining A. Graham J. M. R. Hand M. M. Infield D. G. Paynter R. J. H. Simms D. A. 1997Final report of IEA Annex XIV: field rotor aerodynamics, ECN-C- 97 027.
Schlichting H. 1979, McGraw-Hill, New York.
Sears W. R. 1950Potential flow around a rotating cylindrical blade, , 17 9 183 184.
Shen W. Z. Sørensen J. N. 1999Quasi-3-D Navier-Stokes model for a rotating airfoil, , 150 518 548.
Snel H. Houwink R. van Bussel G. J. W. Bruining A. 1993Sectional prediction of 3-D effects for stalled flow on rotating blades and comparison with measurements, , Lübeck-Travemünde, Germany, 395 399.
Sørensen J. N. 1986Prediction of the three-dimensional stall on wind turbine blade using three-level, viscous-inviscid interaction model, , 429 435.
Sørensen N. N. Michelsen J. A. Schreck S. 2002Navier-Stokes predictions of NREL phase VI rotor in the NASA Ames 80ft×120ft wind tunnel, 5 2-3, 151 169.
Tangler J. L. Kocurek J. D. 1993Wind turbine post stall airfoil performance characteristics guidelines for blade-element momentum methods, (591), 1-10.
Tangler J. L. 2002The nebulous art of using wind tunnel airfoil data for predicting rotor performance, , 5 245 257.
Tangler J. L. 2004Insight into wind turbine stall and post-stall aerodynamics, 7 247 260.
Vatistas G. H. Kozel V. Mih W. C. 1991A simpler model for concentrated vortices, , 11 73 76.