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

Authors: Ali Helali; Badreddine Kamoun; David Afungchui; Mohamed Jomaa Safi
Journal: International Journal of Rotating Machinery
 (2012)

---

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


---

## Body

## 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-lengthR in to a number of small spanwise discrete segments of elementary lengths Δr, 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 W¯ is as high as possible and which leads to a minimal losses P 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](1)∫rmRΓM⋅∑NP∫-∞+∞g(Q,M)rMrQV02Ω2+rQ2dαQdrM+∫rmRΓM⋅∑NP∫-∞+∞g(M,Q)V02Ω2+rM2dαMdrM+CV0=0.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(Γ/ΩR2)j (see Figure 1). The component (w/ΩR) of the velocity induced according to Morgan and Wrench [13] is given by(2)(wΩR)i=Nb4πλ*(1+G1)(ΓΩR2)jifξi≺ξj-,(wΩR)i=-Nb4πλ*G2(ΓΩR2)jifξi≻ξj-,
where G1 and G2 are given by(3)G1=(zjzi)1/2[μ1-μ-Alog⁡(1-μ)],(4)G2=(zjzi)1/2[1μ-1+Alog⁡(1-μ-1)].
With y=ξ/λ*, z=1+y2 and
(5)μ=(zi-1yiyizj-1)NpeNp(zi-zj).Figure 1
Helicoidal vortex sheet originating from the lifting line.A is given by Morgan and Wrench:(6)A=124Np[9yj2+2zj3+3yi2-2zi3].The power extracted by the wind turbine is written in the following form:(7)W̅=∫rmRρΩrM(V0+w)ΓMdrM=W-P,
where
(8)W=∫rmRρΩrMV0ΓMdrM,-P=∫rmRρΩrMwΓMdrM,W represents the kinetic energy of wind (power extracted without any losses) and P represents the losses resulting from the induced velocities.The power coefficientCP is given by(9)CP=∫rmRρΩrM(V0+w)ΓMdrM(1/2)πρR2V03.
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 rmR 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 vn, 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 Ωu⃗z. The optimal condition is written as follows [9, 13]:(10)vn=(CΩu⃗z×ru⃗r)⋅n⃗.Indeed, if we designate byΣ0 a portion of the helicoidal vortex sheet corresponding to one complete functioning period T, the kinetic energy E imparted to the fluid downstream of the rotor disc during the period is given by(11)E(Γ)=-ρ2∫Σ0Γϕn′dΣ.The power extracted by the wind turbine for one functioning period is given by(12)U(Γ)=∫rmR∫0Tρ(Ωu⃗z×ru⃗r)⋅(Γdru⃗r×Wdtu⃗s)dΣ=ρΩ∫Σ0ru⃗θ⋅u⃗nΓdr.The target now is to look for the condition to be satisfied so that the kinetic energyE is minimal when the value of U in the preceding (5) is imposed. Let us suppose that there is a variation of the circulation distribution Γ(r); by designating this variation by δΓ(r), the variation of U which results from (5) is given by(13)δU(δΓ)=ρ∫Σ0Ωru⃗θ⋅u⃗nδΓδΣ.The corresponding variation ofE is given by
(14)δE(δΓ)=-ρ2∫Σ0(δΓϕn′+Γδϕn′)δΣ.Now the velocityϕn′ induced by a doublet distribution of intensity Γ into the points P and M on the sheet Σ0 is given by(15)(ϕn′)M=∫Σ0Γ(P)∂∂nM∂∂nPg(P,M)dΣP=∫Σ0Γ(P)k(P,M)dΣP,
where g(P,M)=1/(4πrPM) and the function k(P,M) is defined by:(16)k(P,M)=∂∂nM∂∂nPg(P,M),k(P,M) is a symmetric function with respect to the variables P and M; as a result we can write(17)δE(δΓ)=-ρ2∫Σ0(δΓΓk(P,M)+ΓδΓk(P,M))δΣP,δE(δΓ)=-ρ∫Σ0∂Γϕn′dΣP.Let∂Γ and ∂Γ′ be arbitrary variations of the initial distribution of Γ(P); the corresponding variations of U are ∂U(∂Γ) and ∂U(∂Γ′). The two variations of Γ can be superposed in such a way that the resulting variation of U is zero. This implies that U remains unchanged as expressed in the following equation:(18)∂U(∂Γ)-k∂U(∂Γ′)=0.The resulting variation ofΓ(P) is ∂Γ-k∂Γ′.The corresponding variation ofE is equally given by ∂E(∂Γ)-k∂E(∂Γ′)IfE(Γ) corresponding to the initial distribution of Γ is extremum; that is, it is either minimal or maximal, then the corresponding variation of E has to be zero. This last condition enables us to obtain an expression for k given by k=∂E(∂Γ)/∂E(∂Γ′).By substitutingk in (10), we obtain(19)∂U(∂Γ)-∂U(∂Γ′)∂E(∂Γ′)∂E(∂Γ)=0.By introducing the Lagrange constant:C=[∂U(∂Γ′)]/[∂E(∂Γ′)], the preceding equation becomes: ∂U(∂Γ)-C∂E(∂Γ)=0.By considering the expressions which give∂U(∂Γ) and ∂E(∂Γ), that is, (6) and (9), the optimal condition is written as(20)ρ∫Σ0(Ωru⃗θ⋅u⃗n-Cϕn′)∂ΓdΣP=0.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)Cϕn′=Ωru⃗θ⋅u⃗n,
in which (ϕn′) is recognized as a speed vn 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)vn=CΩru⃗θ⋅n⃗.But the component in the direction ofu⃗θ of the vector normal to the vortex sheet is given by(23)u⃗θ⋅n⃗=nθ=λ*(λ*2+ξ2)1/2,
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(24)vn=Cλ*Ωr(λ*2+ξ2)1/2.Now, the normal component of the induced velocity is linked to the axial componentvz, along u⃗z by(25)vn=1ξ(λ*2+ξ2)1/2vz.It comes out that the optimal condition concerning the componentvz of the induced velocity is given by(26)vzΩR2=Cλ*Ωξ2λ*2+ξ2.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]:(27)∑j=1j=11(VZIJ)i,j(Γ*ΩR2)j=ξ2λ*2+ξ2,
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 X, 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 Figures2 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: λ0=V0/ΩR) and for a given number of the blades (Nb=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 forNb=3 for the inverse OPHWT model.Figure 3
Circulation Coefficient on the lifting line forNb=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 Figure4. 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 CP=0.37.For the same coefficient powerCp=0.37, the difference in the value notified in Figure 4 reaches 30% in the reference chord spanwise as defined at 0.75 R. 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.Figures5 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 Cp. We make some variations on the power coefficient with the functioning parameter lambda (λ0=V0/ΩR) 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 forNb=3 and α=8°.Figure 6
Span wise variation of the chord in the kinetic energy model forNb=3 and α=8°.Figure6 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λ0 (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. Figure7 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−1 and for a coefficient power Cp=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 Cp=0.48.Figure 8
Spanwise variation of the chord in the inverse OPHWT model and the Near wake model forNb=3 and α=8°.As can be observed in Figure7, 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−1 corresponding to a tip speed of 80 m s−1 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 Γ=k∞C(r)W⃗(r)α(r). 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λ0, a fixed power coefficient CP 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 CP.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.



---

*Source: 476908-2012-01-24.xml*