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

## Abstract

### Keywords

• wind energy conversion
• 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).

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 dmis an elemental mass of air flow; Vis a velocity of air flow; βis the angle between air flow direction and normal to the plate surface L1 (Figure 1); dNis the impact force in the normal direction to the plate surface L2; and tis time.

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

dm=ρVBcosβdtdLE2

where dLis elemental length of plate’s surface; Bis 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 Cis 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 F1 applied to the rectangular plate is described by the following equation:

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

where dis the ratio of edges L2/L1; His 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 Cis 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.

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 Cis 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 xyis shown in Figure 3. The model includes linear elastic element with stiffness coefficient cand linear viscous damper with damping coefficient b.

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

Vr=V+v,E7

where Vis 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¨=cxbẋ0.5+cosβ2cosβ(V+ẋ)2sign(V+ẋ),E8

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

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

P=b(ẋ)2.E9

The average power Paduring time tis determined by integration of Eq. (9):

Pa=0tb(ẋ)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); V0 = 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.

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).

As it is seen from the analysis of graphs presented (Figures 6 and 7), almost stationary oscillatory regime with maximal average power Pacan 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 Pgenerated due to action on the plate of air flow with velocity Vis 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).

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

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

where mis mass of the plate; Land Bare dimensions of the plate; α and γ are angles of flow and plate orientation against xaxis.

The optimization task is solved under the condition of plate motion with a constant velocity ẋ= u. In this case, using Eq. (11), the braking power Pof force Qcan 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 Pat specified limits. The example of response surface for the criterion Pas a function of plate angle α and velocity uis shown in Figure 9. In addition, Figure 10 presents the section of this response surface corresponding to the constant velocity u= 8 m/s.

The analysis of Eq. (12) shows that there is an extreme change of criterion Punder the parameter ucorresponding to the following condition:

u=Vsinγα3sinαE13

Under the velocity ucorresponding to the condition (13) the maximal power Pis 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.

In accordance with aerodynamics theory [13], flat blade 2 placed in air flow is subjected to action of aerodynamic force R(Figure 11). Force Rcan 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 Cxand 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 Cxand Cyare dependent on blade’s geometry, its orientation relative to air flow and dimensionless Reynolds number [13]. Coefficients Cxand Cyfor the considered flat blade were determined experimentally in wind tunnel ARMFIELD. Principle diagram of experimental setup is shown in Figure 12.

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 Fxand Fyfor different angles α and for the given flow velocity Vwere measured with the aid of weight mechanism (Figure 12). After that, dimensionless aerodynamic coefficients Cxand Cywere calculated using Eq. (15).

On the base of mathematical processing of experimental results, coefficients Cxand Cyare 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.

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 Cfor resulting aerodynamic force Rcan 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 Rgives maximal torsional moment Mabout longitudinal axis O1 of the rotor (Figure 15). But moment Mis determined by projection Rτof force Ron tangent line τ. Therefore, angle α is optimal, if projection Rτbecomes maximal. To satisfy this condition, it is necessary to maximize a projection Cτof resulting aerodynamic coefficient C.

Projection of resulting coefficient Con 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 φ.

### 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.

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 Rτ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.

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).

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.

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.

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.

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited.

## How to cite and reference

### 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:

### Related Content

#### 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

#### Optimum Composite Structures

Edited by Karam Maalawi

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.