The Circulation Distribution on the Lifting Line for a Given Extracted Power

Abstract

Presently, there exist few numerical methods which treat the inverse problem for the determination of the geometry of wind turbine blades. In this work, authors intend to solve the inverse optimum project for horizontal axis wind turbine in which the selection of the circulation distribution is obtained by resolving two variational problems: the first consists in sorting the circulation distribution on the lifting line, which, for a given power extracted by the wind turbine, minimizes the loses due to the induced velocity. In the second, the optimal circulation distribution is selected such that the kinetic energy of the wind downstream of the rotor disc is minimum, when the energy extracted by the wind turbine for one rotating period is imposed. A code has been developed which incorporates the real pitch of the helicoidal vortex wake. Very promising results have been obtained: the circulation distribution for a given extracted power and the chord lengths distribution law along the blade span.

1. Introduction

The blade element momentum (BEM) model [1] is the most common model used in aerodynamic and aeroelastic codes for the prediction of wind turbine characteristics. This is due to its simplicity, which makes it possible to implement the model as a fast numerical algorithm. The accuracy of the BEM model is in general reasonable for a wind turbine in normal operating conditions. A verification study [2] comprising the most common aeroelastic codes in Europe showed a typical difference of 5%–10% between measured and simulated mean blade loads and 5%–20% difference in dynamic blade loads. However, BEM simulations depend heavily on the input aerofoil sectional data, as was seen in the NREL blind test comparison of aerodynamic codes with data from the NASA-Ames wind tunnel experiment [1, 3]. Usually a correction of the aerofoil data for three-dimensional flow effects is necessary in order to improve the correlation of numerical results with experimental data.

There are some fundamental shortcomings of the BEM model which can be ascribed to the assumptions introduced in its derivation. The equation relating the induction at the disc with the axial thrust is derived for a stream tube enclosing the whole rotor disc, but in the BEM model the same relation is used for differential stream tubes at different radial positions (strip theory). Comparison of the BEM model with more accurate induction models, such as the numerical actuator disc model with the flow solved with a CFD model, indicates that this assumption is not accurate in regions of the blade with strong radial variation of the loading, as, for example, towards the tip and at the root [1, 4].

Owing to the assumptions in the BEM model, more physically based aerodynamic models such as vortex line and panel methods have been developed for some problems (rotor shapes or flow conditions). These models reflect detailed information about the induced flow field at the rotor disc and have been successfully applied, for example, to studies of dynamic inflow situations [1, 5].

Though numerous studies have been effectuated on the resolution of the inverse problem on rotating machines, there exist few numerical methods which treat the inverse problem for the determination of the geometry of a wind turbine blade.

If the wake of the blade and its slipstream were to be known exactly, then the field of induced velocities and other relevant quantities could be computed more readily. This can be done, in fact, by appropriate manipulation of Biot-Savart equation for the inviscid vortex line. The classical expression of a vortex-induced velocity is the Biot-Savart law that is the fundamental relationship between a vortex, its shape, and the velocity that it induces. This method proves extremely valuable, because it can track the main pattern of the vortex system.

The main idea of the present article consists in determining the circulation distribution which corresponds to a given extracted power. It is evident that the power coefficient should be chosen under Betz limit moreover under the limit curve for 3-blade Betz rotor [6]. We try to develop the inverse “OPHWT” model [7, 8]. This model can provide us the optimal aerodynamic and geometric characteristics for a given extracted power. So after validating the “OPHWT” in [7], we try to exploit it really by using some of existing power coefficient for the machine.

The selection of the circulation distribution for the inverse problem is obtained by resolving two variational problems: the first consists in sorting the circulation distribution on the lifting line, which, for a given power extracted by the wind turbine, minimizes the loses due to the induced velocity “OPHWT” [7, 8]. In the second, the optimal circulation distribution is selected such that the kinetic energy of the wind downstream of the rotor disc is minimum, when the energy extracted by the wind turbine for one rotating period is imposed.

It is shown that the optimal circulation distribution corresponding to the best energetic efficiency of the wind turbine induces a velocity which satisfies the slip condition. This actually occurs in the case where the helicoidally vortex sheet becomes rigid and rotates with the same speed as the wind turbine, but in the opposite sense.

The method used in calculating the induced velocities is publicized in [7].

The resolution is achieved by the singularities method based on discrete distributions of circulation: it consists in dividing the blade span-length 𝑅 in to a number of small spanwise discrete segments of elementary lengths Δ𝑟, such that the divisions are progressively more concentrated towards the root and the tip and less concentrated towards the middle of the blade [7].

2. Inverse OPHWT Model

In this inverse “OPHWT” model, we intend in the reality to find the efficacy rotor. In fact, for a given wind site characterized by the wind speed and for a given number of the blades, we aim to use the machine with the minimum losses of power. These wishes consist in imposing some law for the induced velocities on the lifting line, and we determinate the optimum circulation distribution which induce these velocities.

The selection of the circulation distribution is based on the hypothesis of limiting the span length of the blade. An optimal circulation distribution is sorted by minimizing the power losses due to the effect of the induced velocity downstream of the turbine disc. It is a variational problem which consists in finding the circulation distribution Γ(𝜉) for which a power output 𝑊 is as high as possible and which leads to a minimal losses 𝑃 due to the effect of the induced velocity [7, 8].

It should be noted that the wake shed along the blades is assumed to form a helicoidal vortex sheet originating from the trailing edge. The system of free vortices emanating from the trailing edge is specified to persist without disintegration up to a relatively far distance downstream of the rotor disc, by following a local stream line [9–12].

The method used in calculating the induced velocities by the vortex sheet is derived directly from the general lifting lines and lifting surfaces theories without any simplifications made. It was elaborated by Goldstein [9]. He was of the opinion that the rotation of the wake in the case of a finite number of blades constituted a loss of energy that would have to be provided by the propeller’s motor and incorporated the loss in his theory. In order to correct for the losses, Goldstein used a complicated method based on the development of the circulation in terms of Bessel’s functions and trigonometric series whose convergence is not assured especially towards the blades tips and the values of the circulation obtained are only poor approximations [7]. A recent development of these calculations was applied to marine screw propellers where the lifting line hypothesis was replaced by that of the lifting surface, because of the larger blades by Luu. The calculations of Morgan and Wrench [13] shall be adopted in the present paper. The vortex theory is applicable in the case of small attack angles where the lift curve of the blade sections is a linear function of the attack angle. The vortex sheet is as a result stationary with respect to the blades [7].

The integral equation for determining the optimum spanwise distribution of circulation Γ(𝜉) on the blade is given by [7]^𝑅_𝑟𝑚Γ_𝑀⋅_𝑁𝑃^+∞_−∞𝑔()𝑟𝑄,𝑀^𝑀𝑟_𝑄𝑉²₀Ω²+𝑟²_𝑄𝑑𝛼_𝑄𝑑𝑟_𝑀+^𝑅_𝑟𝑚Γ_𝑀⋅_𝑁𝑃^+∞_−∞𝑔(𝑀,𝑄)𝑉²₀Ω²+𝑟²_𝑀𝑑𝛼_𝑀𝑑𝑟_𝑀+𝐶𝑉₀=0. (1)

A semi-infinite helicoidal vortex is made up of a vortex filament which is moving from the rotor disc plane to downstream infinity while carrying the vortex intensity (Γ/Ω𝑅²)_𝑗 (see Figure 1). The component (𝑤/Ω𝑅) of the velocity induced according to Morgan and Wrench [13] is given by𝑤Ω𝑅_𝑖=𝑁_𝑏4𝜋𝜆^∗1+𝐺₁ΓΩ𝑅²_𝑗if𝜉_𝑖≺𝜉_𝑗,𝑤Ω𝑅_𝑖𝑁=−^𝑏4𝜋𝜆^∗𝐺₂ΓΩ𝑅²_𝑗if𝜉_𝑖≻𝜉_𝑗, (2) where 𝐺₁ and 𝐺₂ are given by𝐺₁=𝑧_𝑗𝑧_𝑖^1/2𝜇,𝐺1−𝜇−𝐴log1−𝜇₂=𝑧_𝑗𝑧_𝑖^1/21𝜇−1+𝐴log1−𝜇⁻¹. (3)(4) With 𝑦=𝜉/𝜆^∗, 𝑧=1+𝑦² and 𝑧𝜇=^𝑖−1𝑦_𝑖𝑦_𝑖𝑧_𝑗−1^𝑁𝑝𝑒^{𝑁𝑝(𝑧𝑖−𝑧𝑗)}. (5)

Figure 1 

Helicoidal vortex sheet originating from the lifting line.

𝐴 is given by Morgan and Wrench:1𝐴=24𝑁_𝑝9𝑦²_𝑗+2𝑧³_𝑗+3𝑦²_𝑖−2𝑧³_𝑖. (6)

The power extracted by the wind turbine is written in the following form:𝑊=^𝑅_𝑟𝑚𝜌Ω𝑟_𝑀𝑉₀+𝑤Γ_𝑀𝑑𝑟_𝑀=𝑊−𝑃, (7) where 𝑊=^𝑅_𝑟𝑚𝜌Ω𝑟_𝑀𝑉₀Γ_𝑀𝑑𝑟_𝑀,−𝑃=^𝑅_𝑟𝑚𝜌Ω𝑟_𝑀𝑤Γ_𝑀𝑑𝑟_𝑀, (8)𝑊 represents the kinetic energy of wind (power extracted without any losses) and 𝑃 represents the losses resulting from the induced velocities.

The power coefficient 𝐶_𝑃 is given by𝐶_𝑃=∫^𝑅_𝑟𝑚𝜌Ω𝑟_𝑀𝑉₀+𝑤Γ_𝑀𝑑𝑟_𝑀(1/2)𝜋𝜌𝑅²𝑉³₀. (9) Details of this model are presented in [7].

3. Minimization of the Kinetic Energy Model

The aim of this method is to minimize the kinetic energy of the wind down stream of the rotor disc. This leads to increasing the induced velocity downstream in order to attain the optimum efficiency which is under the limit curve for 3-blade Betz rotor [6]. The optimization of wind turbines using helicoidal vortex model with minimization of the kinetic energy was made by Chattot in 2003 [14]. His solution differs from the present work because he used pure numerical simulations with a discretization of the helical vortex lines instead of the Morgan and Wrench’s approximation. In this paper, numerical model is achieved by the singularities method based on discrete distributions of circulation: it consists in dividing the blade distance 𝑟_𝑚𝑅 in to a number small spanwise segments, such that the divisions are progressively more concentrated towards the root and the tip and less concentrated towards the middle of the blade [6, 7].

A variational calculation is used to show that the optimal distribution of the circulation Γ(𝜉) corresponding to a high-energy output by the wind turbine will induce a normal velocity 𝑣_𝑛, whose expression is similar to that obtained by applying the slip condition, in the case where the helicoidal vortex sheet becomes stiff and turning with a rotational speed Ω𝑢⃗_𝑧. The optimal condition is written as follows [9, 13]:𝑣_𝑛=𝐶Ω𝑢⃗_𝑧×𝑟𝑢⃗_𝑟⋅𝑛⃗. (10)

Indeed, if we designate by Σ₀ a portion of the helicoidal vortex sheet corresponding to one complete functioning period 𝑇, the kinetic energy 𝐸 imparted to the fluid downstream of the rotor disc during the period is given by𝐸(Γ)𝜌=−2_Σ0Γ𝜙^′_𝑛𝑑Σ. (11)

The power extracted by the wind turbine for one functioning period is given by𝑈(Γ)=^𝑅_𝑟𝑚^𝑇₀𝜌Ω𝑢⃗_𝑧×𝑟𝑢⃗_𝑟⋅Γ𝑑𝑟𝑢⃗_𝑟×𝑊𝑑𝑡𝑢⃗_𝑠𝑑Σ=𝜌Ω_Σ0𝑟𝑢⃗_𝜃⋅𝑢⃗_𝑛Γ𝑑𝑟. (12)

The target now is to look for the condition to be satisfied so that the kinetic energy 𝐸 is minimal when the value of 𝑈 in the preceding (5) is imposed. Let us suppose that there is a variation of the circulation distribution Γ(𝑟); by designating this variation by 𝛿Γ(𝑟), the variation of 𝑈 which results from (5) is given by𝛿𝑈(𝛿Γ)=𝜌_Σ0Ω𝑟𝑢⃗_𝜃⋅𝑢⃗_𝑛𝛿Γ𝛿Σ. (13)

The corresponding variation of 𝐸 is given by ()𝜌𝛿𝐸𝛿Γ=−2_Σ0𝛿Γ𝜙^′_𝑛+Γ𝛿𝜙^′_𝑛𝛿Σ. (14)

Now the velocity 𝜙^′_𝑛 induced by a doublet distribution of intensity Γ into the points 𝑃 and 𝑀 on the sheet Σ₀ is given by𝜙^′_𝑛𝑀=_Σ0Γ(𝑃)𝜕𝜕𝑛_𝑀𝜕𝜕𝑛_𝑃𝑔(𝑃,𝑀)𝑑Σ_𝑃=_Σ0Γ(𝑃)𝑘(𝑃,𝑀)𝑑Σ_𝑃, (15) where 𝑔(𝑃,𝑀)=1/(4𝜋𝑟_𝑃𝑀) and the function 𝑘(𝑃,𝑀) is defined by:𝑘()=𝜕𝑃,𝑀𝜕𝑛_𝑀𝜕𝜕𝑛_𝑃𝑔(𝑃,𝑀), (16)

𝑘(𝑃,𝑀) is a symmetric function with respect to the variables 𝑃 and 𝑀; as a result we can write()𝜌𝛿𝐸𝛿Γ=−2_Σ0(𝛿ΓΓ𝑘(𝑃,𝑀)+Γ𝛿Γ𝑘(𝑃,𝑀))𝛿Σ_𝑃,𝛿𝐸(𝛿Γ)=−𝜌_Σ0𝜕Γ𝜙^′_𝑛𝑑Σ_𝑃. (17)

Let 𝜕Γ and 𝜕Γ^′ be arbitrary variations of the initial distribution of Γ(𝑃); the corresponding variations of 𝑈 are 𝜕𝑈(𝜕Γ) and 𝜕𝑈(𝜕Γ^′). The two variations of Γ can be superposed in such a way that the resulting variation of 𝑈 is zero. This implies that 𝑈 remains unchanged as expressed in the following equation:𝜕𝑈(𝜕Γ)−𝑘𝜕𝑈𝜕Γ^′=0. (18)

The resulting variation of Γ(𝑃) is 𝜕Γ−𝑘𝜕Γ^′.

The corresponding variation of 𝐸 is equally given by 𝜕𝐸(𝜕Γ)−𝑘𝜕𝐸(𝜕Γ^′)

If 𝐸(Γ) corresponding to the initial distribution of Γ is extremum; that is, it is either minimal or maximal, then the corresponding variation of 𝐸 has to be zero. This last condition enables us to obtain an expression for 𝑘 given by 𝑘=𝜕𝐸(𝜕Γ)/𝜕𝐸(𝜕Γ^′).

By substituting 𝑘 in (10), we obtain𝜕𝑈(𝜕Γ)−𝜕𝑈𝜕Γ^′𝜕𝐸𝜕Γ^′𝜕𝐸(𝜕Γ)=0. (19)

By introducing the Lagrange constant: 𝐶=[𝜕𝑈(𝜕Γ^′)]/[𝜕𝐸(𝜕Γ^′)], the preceding equation becomes: 𝜕𝑈(𝜕Γ)−𝐶𝜕𝐸(𝜕Γ)=0.

By considering the expressions which give 𝜕𝑈(𝜕Γ) and 𝜕𝐸(𝜕Γ), that is, (6) and (9), the optimal condition is written as𝜌_Σ0Ω𝑟𝑢⃗_𝜃⋅𝑢⃗_𝑛−𝐶𝜙^′_𝑛𝜕Γ𝑑Σ_𝑃=0. (20)

We notice that this is a homogenous equation. If 𝜕Γ must be different from zero, then the content of the bracket must be zero; this leads to𝐶𝜙^′_𝑛=Ω𝑟𝑢⃗_𝜃⋅𝑢⃗_𝑛, (21) in which (𝜙^′_𝑛) is recognized as a speed 𝑣_𝑛 resembling the slip condition in the case of a stiff helicoidal vortex sheet rotating at a suitable angular velocity. The optimal condition in (3) is written now as𝑣_𝑛=𝐶Ω𝑟𝑢⃗_𝜃⋅𝑛.⃗ (22)

But the component in the direction of 𝑢⃗_𝜃 of the vector normal to the vortex sheet is given by𝑢⃗_𝜃⋅𝑛⃗=𝑛_𝜃=𝜆^∗𝜆^∗2+𝜉^21/2, (23) where 𝜆^∗ is the parameter characterising the real advance pitch of the vortex system.

Hence, the normal component of the induced velocity corresponding to the optimal condition is given by𝑣_𝑛=𝐶𝜆^∗Ω𝑟𝜆^∗2+𝜉^21/2. (24)

Now, the normal component of the induced velocity is linked to the axial component 𝑣_𝑧, along 𝑢⃗_𝑧 by𝑣_𝑛=1𝜉𝜆^∗2+𝜉^21/2𝑣_𝑧. (25)

It comes out that the optimal condition concerning the component 𝑣_𝑧 of the induced velocity is given by𝑣_𝑧Ω𝑅²=𝐶𝜆^∗Ω𝜉²𝜆^∗2+𝜉². (26)

The optimal distribution of Γ^∗(𝜉) is obtained to a multiplicative constant by resolving the following linear system of equations deduced from a combination of the Biot Savart’s law and the calculations of Morgan and Wrench [13]:^𝑗=11_𝑗=1(VZIJ)_𝑖,𝑗Γ^∗Ω𝑅²_𝑗=𝜉²𝜆^∗2+𝜉², (27) where Γ^∗(𝜉)=Γ(𝜉)/𝜆^∗.

The resolution of (20) is achieved by the singularities method based on discrete distributions of circulation. The initial values of the circulation distributions and the components of the induced velocities obtained can be further improved iteratively by introducing a multiplicative factor 𝑋, calculated by maximizing the torque coefficient. See [7, appendix  3]. Details of these calculations are presented in [7].

4. Results and Discussions

For a given extracted power, we present in Figures 2 and 3 the optimal circulation coefficient for the two models. In fact, for a given wind site characterized by the wind speed to tip speed ratio at the extremity of the blade (lambda: 𝜆₀=𝑉₀/Ω𝑅) and for a given number of the blades (𝑁_𝑏=3), we determine the optimum circulation distribution which corresponds to the maximum power extractable by the wind turbine.

Figure 2 

Circulation Coefficient on the lifting line for 𝑁𝑏=3 for the inverse OPHWT model.

Figure 3 

Circulation Coefficient on the lifting line for 𝑁𝑏=3 obtained by the kinetic energy model.

This result is very important in designing the blade of a wind turbine. In effect, by using this code we can determine the optimum circulation for every given power coefficient; later on we can determine the chord lengths distribution along the blade span leading to the optimum shape of the blade. This is the main the advantage of the OPHWT model.

A comparison of the next two figures reveals some discrepancy between the two models. This is presented in Figure 4. The circulation coefficient distribution of the OPHWT model appears to be quasiconstant on the lifting line. But for the kinetic energy model this distribution gives the impression to be an elliptic function. It is further seen that the circulation for the two models goes to zero at both ends of the blade, which is well known from finite aspect ratio wing theory.

Figure 4 

Comparison of the circulation coefficient distribution on the lifting line between the OPHWT model and the kinetic energy model for 𝜆0=0.5 and 𝐶𝑃=0.37.

For the same coefficient power 𝐶_𝑝=0.37, the difference in the value notified in Figure 4 reaches 30% in the reference chord spanwise as defined at 0.75 𝑅. This can be attributed to the higher induced velocity in the case of the model where the kinetic energy of the wind downstream of the rotor disc is minimized, and for the reason that in the OPHWT model we reduce the effect of the normal component of the induced velocity.

Figures 5 and 6 present the variation of the chord lengths of the blade sections along the span for an incident angle of 8° for the two models for a fixed power 𝐶_𝑝. We make some variations on the power coefficient with the functioning parameter lambda (𝜆₀=𝑉₀/Ω𝑅) for a wind turbine of 3 blades. This point is perfectly intended by the constructers who suggest to manufacture the specific wind turbine for a given wind site.

Figure 5 

Spanwise variation of the chord in the inverse OPHWT model for 𝑁𝑏=3 and 𝛼=8∘.

Figure 6 

Span wise variation of the chord in the kinetic energy model for 𝑁𝑏=3 and 𝛼=8∘.

Figure 6 reveals clearly the most resistant form of a blade characterized by larger chords towards the hub which progressively becomes less large towards the tip in order to maintain the quasiconstant circulation along the blade (see Figure 2) and resulting in the possibility of the blades enhancing higher loads. It is the well-known shape of the blade of the wind turbine.

However for the same functioning parameter 𝜆₀ (lambda), those predicted by the optimum project are larger because of the hypothesis of minimization of the losses due to the induced velocity. In effect, the OPHWT model is applicable to the wind turbines while the minimization of the kinetic energy is used for the propellers of highest efficiency in which the motive power supplies the torque which is necessary to maintain the rotation. The chord in the case of the kinetic energy model has an elliptic form (Figure 6) which bears a resemblance to the parabolic shape turning on the highest velocities.

The OPHWT is confirmed very efficient to calculate the optimal performances and geometry of a wind turbine. In effect, it is applied in the case of reel functioning conditions. Figure 7 presents the circulation distribution of the OPHWT model compared to that of the near wake model of Beddoes [1, 15] for a 40 m blade rotating at a free wind speed of 7 ms⁻¹ and for a coefficient power 𝐶_𝑝=0,48. The corresponding distribution of the plan-form is shown in Figure 8.

Figure 7 

Comparison of the circulation Γ (m2/s) between the inverse OPHWT model and the near wake model for the 40 m rotating blade at a free wind speed of 7 ms−1 and for 𝐶𝑝=0.48.

Figure 8 

Spanwise variation of the chord in the inverse OPHWT model and the Near wake model for 𝑁𝑏=3 and 𝛼=8∘.

As can be observed in Figure 7, the circulation distributions follow the same trend but there are some discrepancies in the values of the circulation distributions at the tip of the blade. This divergence attains 29% at the section with abscissa 30 m on the blade. In effect, the OPHWT model is based on a variational problem aiming at reducing the energy losses due to the induced velocity by the vortex sheet. The near wake model on the other hand is in principle a lifting line model for the rotating blade, where only a quarter revolution of the wake system behind the blade is taken into account. In this model, Beddoes [15] considers that only the first part of the trailed vortex system behind the blade has an important influence on the shaping of the loading on the blade as well as on the time variation of the downwash at the blade. The circulation distribution of the OPHWT model is thus believed to be the most closest to reality because this circulation seems to be quasiconstant on the lifting line except towards the root and at the tip of the blade where it falls to zero.

The Near wake model was tested on a 40 m rotating blade with an angular velocity of 2.0 rad s⁻¹ corresponding to a tip speed of 80 m s⁻¹ and with a plan-form as shown in Figure 8. The plan-form of the Near wake model is compared to the plan-form obtained by the inverse OPHWT model. The spanwise variations of the chords follow the same trend but there are discrepancies in the values at the bottom of the blade. In effect, the chords in the present inverse problem correspond to a power coefficient of 0.48 obtained by using the circulation equation Γ=𝑘_∞𝐶(𝑟)𝑊⃗(𝑟)𝛼(𝑟). But for the near wake model, the plan-form presented in Figure 8 is an input used to determine the circulation distribution presented in Figure 7.

5. Conclusion

In this paper, an inverse problem has been resolved by two methods: the first consists in the minimization of the power losses due to the induced velocity (inverse OPHWT model) and in the second we have minimized the kinetic energy imparted to the fluid downstream of the rotor disc (prop fan).

Confrontation of the results obtained with that of Helge Aaggard and Flemming (A near wake model for Beddoes) has revealed the potential of the optimum project in predicting improved and higher rotor performances even at relatively high wind speeds.

In effect, for a given wind site characterized by a wind-rose, that is a given functioning parameter 𝜆₀, a fixed power coefficient 𝐶_𝑃 is required for the wind turbine. The inverse OPHWT code can provide us with the optimal circulation distribution as well as the most adequate plan-form for the construction of the blade corresponding to the given 𝐶_𝑃.

As a follow-up of this work we intend to construct a complete blade of a wind turbine in 3D by coupling the inverse OPHWT model with an inverse code which solves the inverse problem of the flow over a 2D airfoil [16, 17].

The aim of the project is to visualize the actual level of wind turbine technology in Tunisia (the Sidi Daoud wind site, Tunisia): essentially the blades and the aerodynamic efficiency in order to set up a local wind turbine industry. Tunisia intends to install other wind sites in order to reach an energy production of about 200 MW by 2012. This paper can help the decision makers in the choice of the blades suitable for a given wind site, which can be purchased, and also to design new forms of blades for the small installations of wind turbines for isolated sites.