Open access peer-reviewed chapter

Synthesis and Optimization of Wind Energy Conversion Devices

By Janis Viba, Vitaly Beresnevich and Martins Irbe

Submitted: February 12th 2019Reviewed: December 10th 2019Published: April 15th 2020

DOI: 10.5772/intechopen.90819

Downloaded: 42

Abstract

An approximate method for the analysis of interaction between wind flow and rigid flat blades is considered. The method allows synthesis and optimization of wind energy conversion systems without using space-time-programming procedures. By this method, the action of wind flow on the blade is subdivided on frontal pressure and vacuum (depression) on leeward side. The method was tested by computer simulation and experiments in wind tunnel. Examples of optimization tasks are solved in application to blades with simple shape. New wind energetic device with controlled orientation of flat blades to air flow is developed. Theoretical and experimental analysis of blade’s interaction with airflow is performed. Aerodynamic coefficients for blade’s drag and lifting forces are determined experimentally in wind tunnel. Optimization of system parameters is made. To increase the efficiency of energy transformation, it is proposed to change, by special law, the orientation of blade’s working surface relative to airflow during rotation of the rotor. It is shown that the optimal angular rotation frequency ratio between rotor and blade is equal to two. Serviceability and main advantages of the proposed method are confirmed by experiments with physical model of airflow device.

Keywords

  • wind energy conversion
  • rigid blades
  • airflow
  • experimental aerodynamics
  • simulation
  • optimization

1. Introduction

Extracting energy from wind flow is one of the most important renewable energy sources. Operation principle of wind devices is based on air flow action on blades mounted on special wheel and further transformation of air flow kinetic energy into the mechanical energy of wheel rotation [1, 2]. In order to increase the efficiency of wind energy transformation, it is proposed to change orientation of blades to air flow during rotation of main rotor. Position of blade relative to air flow is considered as an optimal, if resulting aerodynamic force gives maximal torsional moment about longitudinal axis of the rotor. This chapter considers solution of the problem stated by the theoretical and experimental analysis of blade’s interaction with air flow in different aerodynamic conditions.

Solid body interaction with air flow is studied in many works. Lift and drag coefficients for flat plate interacting with air flow in the range of Reynolds number from 104 to 106 were analyzed by Haibo et al. [3]. It was shown that at small angles of attack lift–drag relation is parabolic, but at higher angles of attack, it is a circle with radius of unity. Stationary regimes with constant vectors of relative velocities are well researched and mathematically described (for example, see [4]).

Nonstationary regimes of interaction between continuous air flow and rigid body, as shown in [5], mathematically are more complicated. Such kind of interaction can be analyzed by numerical simulation with space time programming methods [6].

In this chapter, a simplified method for the analysis of interaction between air flow and rigid body is proposed. Method is based on separate consideration of flow-body interaction on pressure and suction (vacuum) sides.

2. Simplified model of air flow interaction with flat plate

To simplify the analysis, optimization and synthesis of wind energetic devices, it is proposed to use approximated model of flow-plate interaction without considering viscous effects of air medium. For this purpose, air flow interaction with the plate is subdivided on two phases: one interaction is on the pressure side and the other one-on the vacuum side (Figure 1).

Figure 1.

Air flow interaction with rectangular flat plate.

To analyze air flow-plate interaction, the theorem on change of momentum of mechanical system in differential form, is used [7]. In accordance with this theorem, on the side of pressure the following equation can be written:

dmVcosβ=dNdt,E1

where dm is an elemental mass of air flow; V is a velocity of air flow; β is the angle between air flow direction and normal to the plate surface L1 (Figure 1); dN is the impact force in the normal direction to the plate surface L2; and t is time.

Elemental mass dm of air flow can be mathematically described by the equation:

dm=ρVBcosβdtdLE2

where dL is elemental length of plate’s surface; B is width of the plate (in the case of two-dimensional task B = const); and ρ is the air density.

By the integration of Eq. (1), extra pressure and force along x-axis can be determined. Mathematically, the pressure distribution can be described by the expressions:

pL1=V2ρcosβ2E3
pL2=V2ρsinβ2E4

Vacuum side of the plate is loaded by constant pressure, which can be determined by formula:

p2=+V2ρCE5

where C is constant to be found experimentally or by computer simulation.

By this way, it is possible to find total force applied to the body in the air flow. For example, projection on x-axis of total force F 1 applied to the rectangular plate is described by the following equation:

F1x=HBV2ρC+cosβ3+d·sinβ3cosβ+d·sinβ,E6

where d is the ratio of edges L2/L1; H is section height of the plate in the direction perpendicular to the flow (Figure 1, H=L1sinβ+L2cosβ).

As follows from the results obtained, an approximate analytical method can be applied to solve air flow and blade interaction problems. Validation of the results of analytical calculations can be performed using computer simulation and experiments. Examples will be discussed in the following sections.

3. Theory and simulation of plane motion of flat plate

To determine coefficient C, numerical simulations in Solid Works for flat plate interaction with air flow were performed [8, 9]. Modeling results for a rectangular flat plate were compared with analytical formula (6). It was shown that approximate value of coefficient C is about 0.5. The estimation of accuracy of formula (6) for C = 0.5 in comparison with simulation results is presented in graphical form in Figure 2.

Figure 2.

The accuracy of formula (6) for C = 0.5 as function of angle β.

From the analysis of Figure 2, it can be concluded that interaction force fluctuates within certain limits, and value of fluctuation is dependent on angle β. At smaller angles, the value of fluctuations is also small, but with the increasing of β fluctuation of interaction force becomes greater. The relative difference is not very large, and its mean value is about 3.6%.

Therefore, engineering analysis by the proposed method does not require a step by step space-time calculations to find interaction coefficient C. A new nonstationary flow interaction formula was obtained, and this formula includes the object shape, state, and flow rate direction parameters. Coefficient C is taken as C = 0.5. The efficiency of the proposed method is illustrated by the analysis of flat plate interaction with air flow.

A two-dimensional (2D) model of translation motion of thin flat plate (thickness d ∼ 0) in co-ordinate plane xy is shown in Figure 3. The model includes linear elastic element with stiffness coefficient c and linear viscous damper with damping coefficient b.

Figure 3.

Model of thin flat plate to obtain renewable energy from air flow.

In accordance with methods of classical mechanics [7], relative interaction velocity Vr must be taken into account:

Vr=V+v,E7

where V is velocity of air flow; and vis velocity of flat plate along x-axis.

In this case, differential equation of plate motion along x-axis can be written in the following form:

mx¨=cxbẋ0.5+cosβ2cosβ(V+ẋ)2sign(V+ẋ),E8

where A = LB is a surface area of the plate; ρ is the air density; β is plate angle against air flow; and m is mass of the plate.

The renewable energy is generated due to the action of damping force (bẋ).Therefore, momentary power can be determined by formula

P=b(ẋ)2.E9

The average power Pa during time t is determined by integration of Eq. (9):

Pa=0tb(ẋ)2dtt.E10

By the analysis of Eq. (8), it can be concluded that five parameters can be used to control the efficiency of this system. These parameters are as follows: c, b, A, β, and V.

Mathematical simulation of Eq. (8) was performed with program MathCad assuming the following values of main system’s parameters: A = 0.04 m2 (length L = 0.2 m and width B = 0.2 m); V 0 = 10 m/s; ρ = 1.25 kg/m3 (at temperature 10°C).

Results of simulation for control action by angle β=π2.5sin7tare presented in Figures 4 and 5.

Figure 4.

Displacement x as function of time t.

Figure 5.

Average power Pa of generator force b x ̇ during short transient process (t = 4 s).

Average power Pa (Figure 5) is presented as percentage of maximal power, which can be achieved under the plate’s velocity less than one third of flow velocity.

Results of simulation for the case of control action on the system by harmonic variation of angle β (Figures 4 and 5) are as follows:

  • stationary regime of motion occurs very quickly, practically within two to three cycles;

  • it is possible to synthesize the optimal parameters of the system (for example, stiffness c, area A, frequency and amplitude of action, etc.), which would provide the maximum power within the given limitations;

  • further increasing of the efficiency can be achieved by the use of more complex control variation of angle β (biharmonic, polyharmonic, etc.).

Results of simulation of plate’s motion for the case of control action by velocity V=V020.5sin10tare presented in Figures 6 and 7 (it is assumed that V0=10m/s).

Figure 6.

Velocity v of plate’s central point as function of time t.

Figure 7.

Average power Pa of generator force b x ̇ during long transient process (t = 30 s)

As it is seen from the analysis of graphs presented (Figures 6 and 7), almost stationary oscillatory regime with maximal average power Pa can be achieved after some cycles of transient process.

Results of simulation for the case of control action on the system by harmonic variation of velocity V (Figures 6 and 7) are as follows:

  • a new opportunity to generate energy by the harmonic variation of flow rate is discovered;

  • new variants of control action on the system by the flow rate variation using more complex laws (biharmonic, polyharmonic, etc.) are opened for the further research.

4. Parametric optimization of flat plate translational motion

The model of air flow interaction with flat plate is shown in Figure 8. In the given optimization problem, the power P generated due to action on the plate of air flow with velocity V is of interest. Parametric optimization aims to achieve maximum power P (criterion) providing plate’s translational motion with constant speed u and in condition of braking with force Q (parameter).

Figure 8.

Optimization model: V–air flow velocity; u–velocity of flat plate translation motion; and n–normal axis relative to plate.

In accordance with Eq. (6), the motion of a thin plate (d = 0) is described by the following differential equation:

mx¨=LBρ1+CVsinγαẋsinα2·signVsinγαẋsinαsinαQ,E11

where m is mass of the plate; L and B are dimensions of the plate; α and γ are angles of flow and plate orientation against x axis.

The optimization task is solved under the condition of plate motion with a constant velocity ẋ= u. In this case, using Eq. (11), the braking power P of force Q can be expressed in the following form:

P=LBρ1+CVsinγαusinα2·signVsinγαusinαusinα.E12

Parametric optimization problem was solved on the base of analysis of Eq. (12) with the aid of program MathCad. The aim of the optimization: to look for the parameters u, α, and γ that give the maximum power P at specified limits. The example of response surface for the criterion P as a function of plate angle α and velocity u is shown in Figure 9. In addition, Figure 10 presents the section of this response surface corresponding to the constant velocity u = 8 m/s.

Figure 9.

Response surface for the optimization criterion P as a function of plate angle α and velocity u (for the case of L = 0.5 m; B = 0.5 m; ρ = 1.25 kg/m3; C = 0.5; V = 10 m/s; γ = π/2 rad).

Figure 10.

Section of the response surface P(α, u) corresponding to the constant velocity u = 8 m/s.

The analysis of Eq. (12) shows that there is an extreme change of criterion P under the parameter u corresponding to the following condition:

u=Vsinγα3sinαE13

Under the velocity u corresponding to the condition (13) the maximal power P is achieved:

Pmax=4C0V3sinγα327,E14

where C0=LBρ1+C.

For the other initial data, the optimal parameters must be sought numerically.

5. Wind energy conversion device with controlled orientation of blades to air flow

Efficiency of energy conversion at great extent is dependent on blade’s orientation relative to air flow. In majority of existing air flow devices, there is no possibility for special change of blade’s orientation relative to air flow during rotation of main wheel, and therefore, position of blade can be optimal only in specific time instants [2].

Special variation of flat blade’s turning angle during rotation of main wheel is realized in the device described in [10]. But interrelation between turning angles of blade and main wheel are not optimal in this device. Due to this, potential possibilities to increase efficiency of energy extraction are not fully realized.

Theoretical analysis of air flow interaction with rotating flat blade in different aerodynamic conditions is considered in [11, 12]. It is shown that position of blade is optimal, if resulting aerodynamic force gives maximal torsion moment about longitudinal axis of main wheel.

This section deals with the determination of optimal control law for variation of blade’s turning angle during rotation of central wheel.

5.1 Analysis of air flow interaction with flat blade

The considered model of air flow device is shown in Figure 11. Flat blades 2 are hinged to the rotor 1, besides longitudinal axes O1 and O2 of rotor and blades are mutually parallel. Position of blade 2 relative to air flow is given by angle α, but rotation of rotor 1 is evaluated by angle φ. In order to increase the efficiency of wind energy transformation, it is necessary to find optimal relations between angles φ and α during operation of the system.

Figure 11.

Principle model of wind device: 1–central wheel; 2–flat blade.

In accordance with aerodynamics theory [13], flat blade 2 placed in air flow is subjected to action of aerodynamic force R (Figure 11). Force R can be resolved into two components: drag force Fx (acts along flow direction) and lift force Fy (acts in direction perpendicular to air flow). The following formulae are used to calculate these forces [14]:

Fx=0.5CxV2;Fy=0.5CyV2E15

where Cx and Cy –dimensionless drag and lift aerodynamic coefficients; A–area of blade’s working surface; ρ–density of air medium; and V–velocity of air flow.

Aerodynamic coefficients Cx and Cy are dependent on blade’s geometry, its orientation relative to air flow and dimensionless Reynolds number [13]. Coefficients Cx and Cy for the considered flat blade were determined experimentally in wind tunnel ARMFIELD. Principle diagram of experimental setup is shown in Figure 12.

Figure 12.

Diagram of experimental setup.

During experiments, angle α between air flow and blade’s flat surface was varied from 0° (working surface of the blade is parallel to air flow) to 90° (working surface of the blade is perpendicular to air flow), while the air velocity varies from 5 to 20 m/s. Aerodynamic forces Fx and Fy for different angles α and for the given flow velocity V were measured with the aid of weight mechanism (Figure 12). After that, dimensionless aerodynamic coefficients Cx and Cy were calculated using Eq. (15).

On the base of mathematical processing of experimental results, coefficients Cx and Cy are plotted as functions of angle α (Figures 13 and 14, blue lines). To simplify the application of experimental data in engineering calculations, approximation of curves Cx (α) and Cy (α) is made (Figures 13 and 14, red lines) using program EXCEL.

Figure 13.

Drag aerodynamic coefficient Cx as function of the incidence angle α.

Figure 14.

Lift aerodynamic coefficient Cy as function of the incidence angle α.

Approximation functions Cx (α) and Cy (α) mathematically can be expressed with the following equations:

Cxα=1.31sin0.0137α0.01E16
Cyα=0.5094sin0.00315α0.28318α0.6143E17

Coefficient C for resulting aerodynamic force R can be determined by the formula

C=Cx2+Cy2E18

For each value of rotor’s turning angle φ, it is necessary to find the optimal incidence angle α. Angle α is considered as an optimal, if resulting aerodynamic force R gives maximal torsional moment M about longitudinal axis O1 of the rotor (Figure 15). But moment M is determined by projection of force R on tangent line τ. Therefore, angle α is optimal, if projection becomes maximal. To satisfy this condition, it is necessary to maximize a projection of resulting aerodynamic coefficient C.

Figure 15.

Decomposition of aerodynamic force R on components.

Projection of resulting coefficient C on tangent line can be determined by formula (Figure 15)

Crα=Ccosβφ,E19

where β=arctgCxα/Cyαthe angles between vectors Rand Fy.

Optimal values of angle α were calculated with computer program MathCad using Eqs. (16)(19). Calculations were made varying a turning angle φ of rotor from 0° till 360° with the step Δφ = 5°. The results of calculations are presented in Table 1.

φ[deg]α[deg]φ[deg]α[deg]φ[deg]α[deg]φ[deg]α[deg]
0359090180145270180
5379591185148275180
103910092190150280181
154110593195152285182
204311099200154290183
2546115105205156295185
3049120110210158300188
3552125115215160305190
4055130119220163310193
4558135122225165315195
5061140126230167320197
5565145129235170325200
6070150132240172330202
6575155134245175335204
7081160137250177340206
7587165139255178345208
8088170141260179350210
8589175143265180355212
9090180145270180360215

Table 1.

Interconnection between rotor’s turning angle φ and incidence angle α.

Using data of Table 1, a graph α = f(φ) is plotted (Figure 16). Curve α = f(φ) can be used in designing of wind devices to determine optimal positions of blade (angle α) for different possible values of rotor’s turning angle φ.

Figure 16.

The relationship between the optimal incidence angle α and rotor’s turning angle φ.

5.2 Prototype model of air flow device

To achieve maximal efficiency of wind energy transformation, it is necessary to take into account earlier determined optimal relations between turning angles φ and α (Figure 16). But exact curve α = f(φ) is a nonlinear one, therefore its practical realization in wind devices is very difficult.

To simplify practical realization of optimal angles α in real mechanisms, it is proposed to use linear approximation of exact (experimental) curve α = f(φ). Exact and approximated curves α = f(φ) are compared graphically in Figure 17.

Figure 17.

Relationship between angles φ and α (experimental exact curve and its linear approximation).

Mathematically the linear relation between angles φ and α can be described by the following equation:

α=45°+0.=π/4+0.E20

As follows from Eq. (20), one revolution of blade takes place after two whole revolutions of rotor. Therefore, during one cycle both opposite working surfaces of each blade take up a running air flow in turn.

Calculation error on application of approximate linear function (20) has been evaluated. For this purpose, values of angle α determined by formula (20) were inserted into the earlier used mathematical model (16)(19). It was shown that in the most adverse case a relative error in determination of force does not exceed 10%. Herewith, as follows from Figure 17, approximate value of angle α during some stages of rotor motion is slightly smaller than optimal one, but in some other stages–a little more of it. Therefore, on average in one revolution, the difference between approximate and optimal values of angle α (or between generated powers) will be very small.

Besides, application of linear approximation (20) makes it possible to simplify practical realization of optimal operation regime in wind device. In accordance with Eq. (20), ratio of angular rotation frequencies between rotor and blades must be constant and equal to 2. Such constant frequency ratio can be realized using simple transmission (gear, belt or chain). Taking account of these considerations, it has been found expedient to use a linear approximation (20) in designing of wind device.

Kinematic diagram of the proposed air flow device is shown in Figure 18. Flat blade 2 is hinged to the rotor 1. Besides, longitudinal axes O1 and O2 of rotor and blade are mutually parallel and kinematically connected each other with belt transmission 3. Belt pulley 4 is attached to blade 2, but pulley 5–to rotor 1. And in accordance with Eq. (20), angular rotation frequency ratio between rotor 1 and blade 2 is taken as 2.

Figure 18.

Kinematic diagram of air flow device: 1–rotor; 2–flat blade; 3–belt transmission; 4 and 5–pulleys; 6–tail.

To follow possible changes of air flow direction, the proposed device is equipped with tail 6 fastened to pulley 5. Due to the change of air flow direction, tail 6 is turned till becomes parallel to air flow. Simultaneously pulley 5 and rotor 1 are also turned, as the result blade 2 takes optimal orientation to air flow. The proposed method for energy extraction from air flow and device for its realization are patented [15].

In accordance with the proposed kinematic diagram, a prototype model of wind device is made (Figure 19).

Figure 19.

Prototype model of wind energy conversion device.

Prototype model of wind device (Figure 19) contains four identical flat blades, which are kinematically connected with rotor by a toothed belt transmission. Ratio of angular rotation frequencies between rotor and blades is taken as 2; due to this the orientation of blades relative to air flow is changed in accordance with condition (20).

Experiments with prototype device are made in wind tunnel ARMFIELD (Figure 20). Stable and effective operation of the device is demonstrated for the range of flow velocities from 5 till 20 m/s.

Figure 20.

Experiments in wind tunnel ARMFIELD.

Experimental investigations confirm the principal possibility to increase with the proposed method the angular velocity of central wheel and extracted power from air flow under the same wind velocity (in comparison with traditional air flow devices [14]).

Besides, experiments indicate on some shortcomings in operation of the prototype model shown in Figure 19. These shortcomings became especially evident with the increasing of number of blades. In such case the front blade (from the air flow side) can interfere for propagation of air flow to rearward blades. Air flow arriving to rearward blades can change its direction and intensity. Forming of air vortexes inside the space between blades has become possible. All these factors have negative influence on the efficiency of wind energy transformation.

In order to improve the efficiency of wind energy extraction, new design of air flow device is developed. Computer model of this wind device is shown in Figure 21.

Figure 21.

New design of air flow device.

Wind device contains three identical flat blades, which are fastened on axes hinged in central wheel. Blade’s axes are kinematically connected with central rotor by a system of cylindrical gear transmission. In accordance with the condition (20), ratio of angular rotation frequencies between central wheel and blades is taken as 2. Besides, it is proposed to move blades relative each other along their rotation axes. Thanks to this, it became possible to make around each blade free space in radial direction for air flow access.

During rotation blades do not interfere each other (Figure 21). Each blade has its own open channel for passing of air flow. Therefore, application of this device will make it possible to increase the efficiency of wind energy transformation. More detailed quantitative analysis of operation of this device would be critical in the future.

6. Conclusions

New method for approximate analysis of air flow interaction with rigid body is developed. This method allows to solve problems of synthesis and optimization of wind devices in a simplified way, without using intensive and laborious space–time-programming procedures.

The aerodynamic drag and lift coefficients for flat blade are determined experimentally for different positions of blade relative to air flow.

Based on theoretical and experimental analysis, new designs of wind energy conversion devices with controllable orientation of flat blades to air flow are proposed.

A prototype model of the developed air flow device is given. Experiments confirm a serviceability of the device and efficiency of wind energy transformation.

Conflict of interest

The authors confirm that this section of monograph has no conflicts of interest.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Janis Viba, Vitaly Beresnevich and Martins Irbe (April 15th 2020). Synthesis and Optimization of Wind Energy Conversion Devices, Design Optimization of Wind Energy Conversion Systems with Applications, Karam Y. Maalawi, IntechOpen, DOI: 10.5772/intechopen.90819. Available from:

chapter statistics

42total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Design Optimization of Wind Energy Conversion Systems with Applications

Edited by Karam Maalawi

Next chapter

Optimal Design and Operational Monitoring of Wind Turbine Blades

By Francis Xavier Ochieng, Craig Matthew Hancock, Gethin Wyn Roberts and Julien Le Kernec

Related Book

First chapter

Introductory Chapter: An Introduction to the Optimization of Composite Structures

By Karam Maalawi

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us