Abstract
A set of systematical optimum operational parameters for wind turbines under various wind directions is derived by using combined momentum-energy and blade-element-energy concepts. The derivations are solved numerically by fixing some parameters at practical values. Then, the interactions between the produced power and the influential factors of it are generated in the figures. It is shown that the maximum power produced is strongly affected by the wind direction, the tip speed, the pitch angle of the rotor, and the drag coefficient, which are specifically indicated by figures. It also turns out that the maximum power can take place at two different optimum tip speeds in some cases. The equations derived herein can also be used in the modeling of tethered wind turbines which can keep aloft and deliver energy.
1. Introduction
Wind turbines have been around for centuries, from the simple water pump in Arabia to the current widely used HAWT (horizontal axis wind turbine). With more and more public attention drawn to wind energy for its ease of development and environmental friendliness, wind farms that consist of hundreds of wind turbines at the same location are now being constructed rapidly all around the world [1]. Some of the wind farms are even built offshore by using large floating platforms due to the stronger and steadier wind at sea. However, the potential impact of wind farms on the biological, physical, and human environments is still under inspection, and the economics of offshore wind farms is still being demonstrated [2, 3]. Recently, much interest has arisen in new types of wind turbines like tethered wind turbines which can fly at high altitudes while delivering energy from wind [4–6]. Since the captured wind energy is proportional to the cube of the wind speed, a tethered wind turbine at high altitude may have 500 times as much available power as a wind turbine on the surface.
Although wind turbines vary in configuration, most of them incorporate the same mechanism of extracting wind energy by using rotor blades. To determine the optimum conditions when the maximum power is generated by the individual rotor, CFD (computational fluid dynamics) is a widely used method due to its predictive accuracy. Other methods like FVM (free vortex models) have also been developed. However, these approaches become computationally intensive for the study on wind farms that involve hundreds of individual wind turbines [7], and they are of limited use in the modeling of new types of wind turbines such as tethered wind turbines. Therefore, BEM (blade element momentum) is still widely used because of its simplicity and overall accuracy.
BEM theory was formulated by Glauert and gives both the force in the axial direction and the tangential direction [7]. It analyzes the rotor blade in sections and adds up the forces on all sections to get the resultant force of the rotor [8]. However, wind direction is usually not considered in the BEM, and the wind is generally assumed to be perpendicular to the rotor. Therefore, it is still unknown what the true available power and the optimum operational condition would be for wind turbines in yawed flows. This gap has created obstacles in the consideration of wind turbines in yawed flows.
The purpose of the work described herein is to use momentum and blade-element considerations to determine the optimum operating conditions for a yawed wind turbine and to determine the maximum power that the turbine can generate under various wind directions. The work culminates in the design figures which can be used to design wind turbines. The derivation described herein also paves the way for our further study on tethered wind turbines.
2. Physical System
Figure 1 shows the physical system that will be analyzed in this paper. It should be noted that the yawed wind turbine herein is considered to be a rotor, the axis of which is fixed to some support. The wind is assumed to be horizontal, and it impinges on the rotor at a wind direction angle (between the wind direction and the axis of rotor disk), denoted as .
It is further assumed that the forces on the blades are determined solely by the lift and drag characteristics of the airfoil shape of the blades. Since momentum theory holds, the induced flow will be normal to the rotor disk. Swirl velocity is neglected in this study. No in-plane net loading is considered. Our desire is to determine the optimum operating conditions of the wind turbine under various wind directions in order to generate the most power possible. The results are shown in the figures at the end of the paper.
3. Case of a Horizontal Axis Wind Turbine
The first set of equations to be derived and solved are for the case of rotor with horizontal axis, . This case provides a good baseline and insight for the more complicated yawed case that will be derived later.
By breaking a blade down into several small sections, we derive the forces on each of these small blade elements by analyzing the airfoil of the blade (Figure 2). These forces are then integrated along the entire blade. Thus, the resultant force on the rotor is obtained from the forces on all blades [8]. This is achieved by the following.
All the forces can be normalized on : where is the air density, is the radius of the rotor, is the relative wind speed, is the slope of the lift curve, is the angle of attack between the chord line and the relative wind, is the profile drag parallel to the wind direction, is the drag coefficient, is the lift force normal to the wind, and is the lift coefficient.
Since we are interested only in the force normal to and tangential to the plane of rotation, the lift and drag forces are projected in these directions where is the pitch angle between the chord line and the plane of rotation, is the angle between the relative wind and the plane of rotation, is the thrust of each blade normal to the plane of rotation, is the drag of each blade tangential to the circle swept by the blade, is the component of the relative wind speed normal to the plane of rotation, and is the component of the relative wind speed tangential to the circle swept by the blade.
In this theoretical development, it is assumed that we are finding the optimum loading condition for a turbine in uniform air flow. For that reason, it is assumed that the airfoil sections (i.e., the chord and airfoils) have been chosen to give a reasonable life-to-drag ratio. Thus, this paper assumes that all airfoils are in the unstalled region such that the lift coefficient is a constant lift-curve slope multiplied by the angle of attack. Naturally, for a final wind turbine design, the optimum condition would have to be checked to make sure that the airfoils are not stalled and that the tip speeds are in a reasonable Mach range. From this and (1)–(6), it is possible to determine the thrust of each blade element normal to the plane of rotation, denoted as , and the drag of each blade element tangential to the circle swept by the blade, denoted as , as follows (in the appendix): where is the chord length and is the radial location of a blade element.
From the flow geometry, where is the induced flow and is the angular velocity of the blade.
Thus, the thrust force on the entire rotor, denoted as , can be obtained from integrated forces on all blades and the extracted power, denoted as , can be obtained from the work done by the drag force: where b is the number of blades.
Substituting (7) and (8) into (9), we obtain
The above two equations are derived from blade-element theory. It is obvious that the thrust and the power are both relative to the blade number and the angular velocity of the rotor.
We can also obtain the thrust from momentum theory:
Combining (10) and (12) and normalizing all velocities on , we obtain: where is the rotor property parameter, is the induced flow coefficient, which is and should be less than 0.5, is the normalized tip speed, which is .
For (11), we normalize power on and all velocities on : where is the normalized power coefficient extracted from wind.
Equations (13) and (15) are used to determine optimum operational conditions for the horizontal axis rotor. It is useful to firstly assign practical values to a and σ, both of which are only determined by the shape of blades: where and .
Equations (16) and (17) consist of five variables: , , , and . It is difficult to find the direct relationship among five variables with only two equations. To find the optimum operating conditions, we fix one of the two variables ( or ) at a certain value. By varying from 0.01 to 0.09 discretely, we use (16) and (17) to find the relationships among the three remaining variables. Thus, the discussion can be divided into two parts: fixed and fixed.
3.1. Numerical Analysis of -Fixed Case
Figures 4(a), 4(b), 4(c), 5(a), 5(b), and 5(c) show how the produced power interacts with the tip speed and the induced flow when the pitch angle is fixed at practical values. For each value of and , we use (17) to find each corresponding with varied from 0 to 0.5. Then, (16) is used to find each corresponding . It is common in these figures that each curve has one peak which indicates optimum condition, except for the one case with in Figure 4(b), which has two clear peaks. That special curve will be specifically discussed later in our yawed wind turbine. As becomes relatively larger, the phenomenon of the curve twisting is seen in Figures 5(b) and 5(c), which indicate the possibility of one induced flow value corresponding to two different powers produced.
Figures 6, 7, and 8 show the optimum operating conditions when the greatest power is produced. Each of these plots is produced by starting with a value of , stepping through values of , and getting each corresponding and . Then, it is possible to determine the maximum , the corresponding optimum and . It is obvious that the maximum power, the optimum tip speed, and the optimum induced flow decrease as the drag coefficient gets larger. This can be verified by comparing peak positions of each curve in Figures 4(a), 4(b), 4(c), 5(a), 5(b), and 5(c). Although the optimum induced flow coefficient is 1/3 for an ideal rotor in axial flows, the optimum values in Figure 8 are not always 1/3.
3.2. Numerical Analysis of -Fixed Case
Figures 9(a), 9(b), 9(c), 9(d), 10(a), 10(b), 10(c), and 10(d) show how the power interacts with the pitch angle and the induced flow when the tip speed is fixed at practical values. The plotting method is the same numerical process as that for Figures 4(a), 4(b), 4(c), 5(a), 5(b), and 5(c). As increases from Figures 9(a) to 9(d), the curves in each figure move apart from each other and the peak of each curve gradually comes out and moves rightward. However, the movement of the peak is unclear in Figures 10(a), 10(b), 10(c), and 10(d). These optimum conditions can be seen in Figure 13. It should be noticed that the two curves with being 0.08 and 0.09 in Figure 10(d) totally lie below zero and the turbine is unable to produce power in these cases.
Figures 11, 12, and 13 show the optimum operating conditions when the greatest power is produced. The plotting method is the same as that of Figures 6, 7, and 8. It is clear from Figure 11 that the maximum power decreases linearly as the drag coefficient gets larger with the tip speed drastically affecting the velocity of decreasing. In Figures 12 and 13, however, the drag coefficient does not affect the optimum pitch angle and the induced flow, both of which are only determined by the tip speed.
4. Case of a Yawed Wind Turbine
When the wind direction angle is an arbitrary value (Figure 3), we can divide the initial wind speed into components normal and parallel to the rotor axis and obtain the thrust and the power based on momentum and blade-element derivations. We generalize the momentum equations in the following: where is the initial wind speed, is the component of normal to plane of rotation, is the component of parallel to the plane of rotation, is the induced flow, is the mass flow rate at the rotor, and is the thrust on the rotor.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Substituting (18)–(20) into (21), we obtain
We then use the blade-element derivations to obtain the thrust and the drag on the element of the th blade: where is the angle between and the th blade, is the component of the relative wind speed which is tangential to the circle swept by the th blade, is the component of the relative wind speed which is normal to the plane of rotation, is the thrust on the element of the th blade which is normal to the plane of rotation, and is the drag on the element of the th blade which is tangential to the circle swept by the th blade.
Substituting (23) and (24) into (25) and (26), we find
Again, the thrust force on the entire rotor T can be obtained from integrated forces on all blades, and the extracted power P can be obtained from the work done by the drag forces: where is the number of blades.
Substituting (28) into (29), we obtain
To simplify the equation, we consider the most widely used rotor which has three blades:
Thus,
Then, we combine (22) and (32) and normalize all power on , and all velocities on : where is the induced flow coefficient, is the normalized tip speed, and is the rotor property parameter, which is .
Again, we assign a and σ the practical values: where and .
Figures 14(a), 14(b), 15(a), 15(b), 16(a), 16(b), 17(a), 17(b), 18(a), 18(b), 19(a), and 19(b) show how the produced power interacts with the tip speed and the induced flow when the pitch angle and the wind direction angle are fixed at practical values. All plots are produced by starting with a value of , stepping through values of , and getting each corresponding and based on (35). Most figures in the yawed rotor are similar to the figures in the horizontal axis rotor: as the tip speed or the induced flow gets larger, the power increases to a global maximum followed by monotonous decreases; the phenomenon of twisting comes out in versus at the relatively large pitch angle and wind direction angle.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
However, the curve with in Figure 15(a) has two global maxima, which is the same as those in Figure 4(b). For being 2° and being 0°, the maximum is 0.59 when is 8.13 and 20.52; for being 2° and being 15°, the maximum is 0.58 when is 10.43 and 16.54. It is somewhat surprising to get two optimum values of for only one maximum . To study this further, we also fixed at values around 2° and computed the maximum and its corresponding when is 0. The numerical results are shown in Table 1.
It can be concluded from the table that, at some pitch angle value around 1.8°, 1.9°, 2.0°, and 2.1°, two optimum tip speeds may coexist provided that the drag coefficient equals 0 and the wind direction angle is small.
Figures 20(a), 20(b), 20(c), 21(a), 21(b), 21(c), 22(a), 22(b), and 22(c) show the optimum operating conditions when the greatest power is produced. All plots are generated by starting at a value of , stepping through values of , and getting each corresponding and . Then, it is possible to determine the maximum , the corresponding optimum and . It is clear that the maximum power goes down monotonously as the wind direction angle gets larger. For the optimum tip speed and the optimum induced flow, most of them also decrease monotonously as the wind direction angle increases except for the ones with . It should be noted that the cases are numerically sensitive and so behave differently from the cases. As a verification of the coexisting optimum tip speeds mentioned above, the optimum jumps upside and down when in Figure 21(b).
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
5. Conclusion
Momentum-energy and blade-element-energy equations have been derived to find optimum conditions for both yawed wind turbines and new types of wind turbines using rotor blades. They are solved numerically by fixing some parameters at practical values. Then, the interactions between the produced power and the influential factors of it are generated in the figures.
It is a general principle that, under various wind direction angles, the maximum possible power produced is strongly affected by the tip speed, the pitch angle of the rotor, and the drag coefficient from the wind, which are specifically indicated by figures. For each of optimum operating conditions, most maximum powers correspond to unique optimum tip speed except for some special cases in which two optimum tip speeds coexist. As the wind direction angle increases to right angle, the maximum power decreases monotonously to zero.
Appendix
Derivation for the Differential Thrust and Drag
Substituting (1), (4), (5), and (6) into differential forms of (2) and (3), we obtain
Since is far less than , the pitch angle is very little, and the drag coefficient is far less than the lift curve slope in practical situation, we find approximations
Thus,
Nomenclature
: | Slope of lift curve |
: | Number of blades |
: | Chord length (m) |
: | Drag coefficient |
: | Lift coefficient |
: | Normalized power coefficient |
: | Thrust on blade element, which is normal to the plane of rotation (N) |
: | Drag on blade element, which is tangential to the circle swept by rotor (N) |
: | Thrust on element of the th blade, which is normal to the plane of rotation (N) |
: | Drag on element of the th blade, which is tangential to the circle swept by the th blade (N) |
: | Profile drag parallel to the wind direction (N) |
: | Thrust on blade, which is normal to the plane of rotation (N) |
: | Drag on blade, which is tangential to the circle swept by the rotor (N) |
: | Normalized tip speed |
: | Lift force on the rotor (N) |
: | Mass flow rate (kg/s) |
: | Power extracted from wind (Watts) |
: | Radius of blade (m) |
: | Thrust on the rotor (N) |
: | Induced flow (m/s) |
: | Induced flow coefficient (v/W) |
: | Wind speed (m/s) |
: | Component of the initial wind speed which is normal to the plane of rotation (m/s) |
: | Component of the initial wind speed which is parallel to the plane of rotation (m/s) |
: | Component of the relative wind speed which is normal to the plane of rotation (m/s) |
: | component of the relative wind speed which is tangential to the circle swept by the rotor (m/s) |
: | Component of the relative wind speed which is normal to the plane of rotation (m/s) |
: | Component of the relative wind speed which is tangential to the circle swept by the th blade (m/s) |
x: | Radial location of a blade element (m) |
: | Density of air (kg/m3) |
: | Angle of attack (rad) |
: | Wind direction angle, which is between the wind direction and the rotor axis (rad) |
: | Pitch angle (rad) |
: | Angle between the relative wind and the plane of rotation (rad) |
: | Blade angular velocity (rad/s) |
: | Angle between and the th blade (rad) |
: | Rotor property parameter . |