Abstract

The paper presents a nonlinear approach to wind turbine (WT) using two-mass model. The main aim of the controller in the WT is to maximize the energy output at varying wind speed. In this work, a combination of linear and nonlinear controllers is adapted to variable speed variable pitch wind turbines (VSVPWT) system. The major operating regions of the WT are below (region 2) and above rated (region 3) wind speed. In these regions, generator torque control (region 2) and pitch control (region 3) are used. The controllers in WT are tested for below and above rated wind speed for step and vertical wind speed profile. The performances of the controllers are analyzed with nonlinear FAST (Fatigue, Aerodynamics, Structures, and Turbulence) WT dynamic simulation. In this paper, two nonlinear controllers, that is, sliding mode control (SMC) and integral sliding mode control (ISMC), have been applied for region 2, whereas for pitch control in region 3 conventional PI control is used. In ISMC, the sliding manifold makes use of an integral action to show effective qualities of control in terms of the control level reduction and sliding mode switching control minimization.

1. Introduction

In recent years, wind energy is one of the major renewable energy sources because of environmental, social, and economic benefits. The major classification of wind turbines (WT) is fixed speed wind turbine (FSWT) and variable speed wind turbine (VSWT). Compared with FSWT, VSWT have many advantages such as improved energy capture, reduction in transient load, and better power conditioning [1]. For any kind of WT, control strategies play a major role on WT characteristics and transient load to the network [2]. In VSWT the operating regions are classified in to two major categories, that is, below and above rated wind speed. At below rated wind speed the main objective of the controller (i.e., torque control) is to optimize the wind energy capture by avoiding the transients in the turbine components especially in the drive train. At above rated wind speed the major objective of the controller (i.e., pitch control) is to maintain the rated power of the WT. For extracting the maximum power at below rated wind speed the WT rotor speed should operate at reference rotor speed which is derived from effective wind speed. In general wind speed is measured by the anemometer which is varying along the rotor shift area so it is difficult to get the mean wind speed by this measurement. Several literatures are addressing this issue by estimating the effective wind speed. But some other literatures assume that the effective wind speed is directly available [3, 4]. In [5] a combination of proportional integral (PI) and SMC is used to adjust the turbine rotor speed for extracting maximum power without estimating the wind speed. In [6] a PI based torque control is used to control the WT, where optimal gains are achieved by particle swarm optimization and fuzzy logic theory, without estimating the wind speed. Estimation of effective wind speed by an inversion of static aerodynamic model with known pitch angle is discussed in [7]. Nonlinear static and dynamic state feedback linearization control is addressed in [8, 9] where both the single- and two-mass model are taken into consideration and the wind speed is estimated by using Newton Raphson (NR). To accommodate the parameter uncertainty and robustness a higher order sliding mode controller is proposed in [10], which ensures the stability of the controller in both the regions that is, below and above rated speed. An extended Kalman Filter based wind speed estimator is discussed in [11], which uses the conventional ATF and ISC for controlling the WT. In [12] fuzzy logic control has been proposed and the effective wind speed is estimated using sequential Monte Carlo simulator. Discrete time optimal LQR/LTR based on optimal quadratic function is studied in [13]. In [14] the effective wind speed is estimated by the frequency domain data fusion and the observer is formulated based on the mixed sensitivity problem with linear matrix Inequalities.

For above rated wind speed authors in [15] discussed the multivariable control strategy by combining the nonlinear state feedback control for region 2 with linear control for region 3. Finally the results are compared with the existing control strategy such as PID and LQG. In [16] the WT controller is the combination of linear control for blade pitch angle with nonlinear torque control. Nowadays fuzzy and neural networks (NN) are the powerful soft computing methods for controlling nonlinear systems. Authors in [17] discussed a fuzzy and neurofuzzy controller for controlling the WT at above rated wind speed. PSO is used to train the adaptive neurofuzzy controller. WT control using adaptive radial basic NN used for both pitch and torque controller is addressed in [18]. Active disturbance rejection based pitch control for variable speed WT is presented in [19], where the extended state observer is used to estimate the state variable and disturbance. In [20] the controller has two parts, that is, nonlinear feed-forward path and linear feedback path. The feed-forward path uses the information about the desired power output, wind velocity, and the turbine speed to determine the pitch angle required. In [21] NN based pitch angle control is implemented by using multilayer perceptrons with back propagation learning algorithm and radial basis function network for tuning the pitch controller.

This paper uses a PI control for pitch control and two nonlinear controllers such as SMC and ISMC for torque control. The paper is organized as follows. Section 2 discusses the modeling of two-mass model. Problem formulation and control objectives are discussed in Section 3. The proposed controllers for all the regions are discussed in Section 4. Section 5 discusses the validation of the results using FAST simulator. Finally in Section 6 a conclusion is drawn from the obtained results, which concludes that the proposed controllers are working fine for controlling the WT at below and above rated wind speed.

2. WT Model

WT is a device which converts the kinetic energy of the wind in to electric energy. Simulation complexity of the WT purely depends on the type of control objectives. In case of WT modelling complex simulators are required to verify the dynamic response of multiple components and aerodynamic loading. Generally dynamic loads and interaction of large components are verified by the aeroelastic simulator. For designing a WT controller, instead of going with complex simulator, the design objective can be achieved by using simplified mathematical model. In this work WT model is described by the set of nonlinear ordinary differential equation with limited degree of freedom. This paper describes the control law for a simplified mathematical model with the objective of optimal power capture at below and above rated wind speed. The proposed controllers are validated for FAST WT model in three different cases of wind speed profile such as below rated wind speed (region 2), above rated wind speed (region 3) and a smooth transition between these two wind speeds (region 2.5). The parameter of the two-mass model is given in Appendix A. Generally VSWT system consists of the following components, that is, aerodynamics, drive trains, and generator, shown in Figure 1.

Equation (1) gives the nonlinear expression for aerodynamic power capture by the rotor From (1) it is clear that the aerodynamic power () is directly proportional to the cube of the wind speed. The power coefficient is the function of blade pitch angle () and tip speed ratio (). The tip speed ratio is defined as ratio between linear tip speed and wind speed:

Generally wind speed is stochastic nature with respect to time. Because of this, tip speed ratio gets affected, which leads to variation in power coefficient. The relationship between aerodynamic torque () and the aerodynamic power is given in (3) where is the torque coefficient given as Substituting (5) in (4) we get In above equation the nonlinear term is which can be approximated by the 5th order polynomial given in (7) where to are the WT power coefficient.

The values of approximated coefficients are given in Table 1. Figure 2 shows the versus curve.

Figure 3 shows the two-mass model of the WT. Equation (8) represents dynamics of the rotor speed with rotor inertia driven by the aerodynamic torque (): Breaking torque acting on the rotor is low speed shaft torque () which can be derived by using stiffness and damping factor of the low speed shaft given in (9): Equation (10) represents dynamics of the generator speed with generator inertia driven by the high speed shaft torque () and braking electromagnetic torque (): Gearbox ratio is defined as From (10) the high speed shaft torque can be expressed as Putting the values of from (13) in (12) we get

3. Problem Formulation

Figure 4 shows the various operating region in VSVPWT.

Region 1 represents the wind speed below the cut in wind speed. Region 2 represents the wind speed between cut in and rated wind speed. In this region the main objective is to maximize the energy capture from the wind with reduced oscillation on the drive train. Region 2.5 represents the wind speed nearer to the rated wind speed, that is, transient period. Region 3 describes the wind speed above the rated wind speed. In this region pitch controller is used to maintain the WT at its rated power.

3.1. Control Structure

The time response of WT electrical system is much faster than the other parts of the WT. This makes it possible to decouple the generator and the aeroturbine control designs and thus define a cascaded control structure around two control loops.(1)The inner control loop consists of electrical generator with power converters.(2)The outer loop has the aeroturbine control which gives the reference to the inner loop.

Figure 5 shows the WT control levels. In this paper we made an assumption that the inner loop is well controlled.

3.2. Control Objective

In this work we consider two objectives. First objective is torque control for optimal power extraction with reduced oscillation on the drive train at below rated wind speed. To achieve the above objective (Region 2) the blade pitch angle () and tip speed ratio () are set to be its optimal value. In order to achieve the optimal tip speed ratio the rotor speed must be adjusted to the reference/optimal rotor speed () by adjusting the control input that is, generator torque (). Equation (15) defines the reference/optimal rotor speed: Second objective is the pitch control for above rated wind speed, where torque control input is consider as maximum value. In region 3 only the pitch angle is varying according to the difference between the nominal speed and output generator speed. This work addresses both controller designs to achieve the control objectives for below and above wind speed.

3.3. Classical Control Techniques

In region 2 in order to compare the results of proposed and existing conventional controllers a brief description of the well-known control techniques, that is, ISC (indirect speed control) and ATF (aerodynamic torque feed-forward) are discussed in this section. In ISC, it is assumed that the WT is stable around its optimal aerodynamic efficiency curve. The two-mass model control signal is given in (16) where where is the low speed shaft damping coefficient brought up to the high speed shaft.

In ATF, proportional control law is used to control the WT. The rotor speed and the aerodynamic torque () are estimated using Kalman filter, which is used to control the WT [22]. The control law is given in (19) The optimal value of proportional gain is found to be . The above existing control techniques have three major drawbacks, that is, the ATF control having more steady state error, so an accurate value of is needed, and in ISC the WT has to operate at its optimal efficiency curve which introduces more power loss for high varying wind speed. Both the controllers are not robust with respect to disturbances. To avoid the above drawbacks two nonlinear controllers, that is, ISMC and SMC, are proposed in region 2. For region 3 control conventional proportional Integral (PI) control is considered.

3.4. Wind Speed Estimation

The estimation of effective wind speed is related to aerodynamic torque and rotor speed provided the pitch angle is at optimal value: The aerodynamic power coefficient is approximated with 5th order polynomial as given in (7): The estimated wind speed can be obtained by solving (24) using MNR. The above equation has unique solution at below rated region. With known the optimal rotor speed is calculated by using (15). The above wind speed estimation is only considered for below rated wind speed condition. Appendix B gives the steps for MNR algorithm.

4. Control Techniques for Below and Above Rated Wind Speed

4.1. Sliding Mode Control (SMC) for Optimal Power Extraction

The proposed control strategy combines MNR based estimator with second order sliding mode controller. SMC is one of the effective nonlinear robust approaches with respect to system dynamics and invariant to uncertainties. Lyapunov stability approach is used in SMC to keep the nonlinear system under control. The objective of the controller is to optimize the energy capture from the wind by tracking the reference rotor speed. In general designing SMC has two steps. First to find the sliding surface second is to develop the control signal .

Time varying sliding surface is defined as where is defined as the difference between rotor speed and reference rotor speed, is the order of the system, and is positive constant. Consider Finally the sliding surface is defined as

Generally three types of reaching law are proposed [23]. Direct switching function method is applied in this work having the condition:

The control should be chosen in such a way that the following candidate Lyapunov function satisfies Lyapunov stability criteria. Lyapunov candidate function is defined as A sufficient condition is that the system output should stay on the sliding surface, that is, . By taking the derivate of the Lyapunov function, only when . From this it is clear that the defined control law is able to track the time varying reference signal, that is, rotor speed reference. If then , from this so The convergence condition is given by the Lyapunov equation which makes the sliding surface attractive and invariant. At steady state the rotor should track the optimal rotor speed asymptotically; that is, as . Consider meets the following condition: The values of the control variable are to be set in such a way that the system will be stable. These control variables are given as follows: In order to derive the control input that is, , the following conversion has been made. By substituting the in (31) we get

By using the relationship given in (14) finally we will get

The control structure is defined as

Generally the SMC have two parts, that is, equivalent control and switching control . The switching is used to avoid the parameter uncertainty and disturbances and the equivalent control is used to control the overall system behaviour, that is, tracking control: The switching control is defined in two ways:

Generally the SMC with signum function introduces the chattering phenomenon in the system. This chattering introduces high frequency dynamics in the WT system, that is, control action to the system in undesirable. To neglect this chattering a smooth control discontinuity is introduced. As the signum function varies between −1 and +1 discontinuously, it is replaced by a tangent hyperbolic function (tanh). Finally the torque control structure is given in (39) where “” is the sliding gain which is chosen based on the empirical results from the simulation. In summary the controller performance depends on the sliding gain and the boundary layer thickness. By using trial and error method the sliding gain is found to be 0.2 and boundary layer thickness is 1 for simulation. If the gain value increases more than 0.2, it introduces more oscillation on the low speed shaft. The boundary layer is chosen arbitrarily that depends on the system performance.

4.2. Integral Sliding Mode Control (ISMC) for Optimal Power Extraction

To improve the sliding surface and overcome the steady state error the integral action is included in the sliding surface.

A sliding surface is defined as where is the integral gain.

The order of the system then the sliding surface modified as

By taking the same Lyapunov function as mentioned in SMC with the same condition meets the following condition: In order to derive the control input, that is, , the following conversion has been made. By substituting the in (41) we get

Finally the torque control structure is given in The necessary condition of ISMC is to eliminate the steady state error in conventional SMC. The linear sliding surface is By taking the Laplace transform of above equation The steady state error is calculated by applying the final value theorem (FVT) to conventional sliding surface: Equation (48) represents the conventional SMC having steady state error which is proportional to boundary layer thickness () and inversely proportional to the coefficient and sliding gain .

The steady state error is calculated by applying the final value theorem (FVT) to integral sliding surface: By taking the Laplace transform of above equation

4.3. Pitch Controller

For above rated wind speed the torque control output is fixed to its rated value that is, maximum value of the control input. In this condition a linear pitch control is introduced. The gain scheduled pitch control is allowed to maintain the generator speed around its nominal value [24]. To achieve the above objective, proportional + integral (PI) control is introduced: where is the nominal value of the generator speed, that is, 1600 rpm, and is the generator speed.

5. Validation Results

The numerical simulations are performed on the CART (controls advanced research turbine) WT model and the characteristics are given in Table 2. The CART is located at the NREL (National Renewable Energy Laboratory) National Wind center near Golden, Colorado. The CART 3 is a three-blade variable speed, variable pitch WT with nominal power rating of 600 kW. It has mainly three parts, that is, rotor, tower, and nacelle. The rotor includes blade and their attachment points are called hub which is maintained by the low speed shaft. The tower is a cantilever beam which supports a yaw bearing and nacelle. The yaw bearing allows the turbine to rotate in the wind speed direction. The nacelle houses the complete drive train assembly. It contains the gearbox, generator, and low speed shaft. The gearbox is directly connected to the squirrel cage induction generator through high speed shaft. The generator is connected to the grid through power electronics that can directly control generator torque [25]. The power electronics consists of three-phase PWM (pulse width modulation) converters with a constant dc link voltage. The main objective of the grid side converter is maintaining the dc link voltage constant [26, 27].

5.1. Description for Simulator

FAST is an aeroelastic WT simulator which is developed by NREL. It can able to model both two- and three-blade horizontal axis wind turbine (HAWT). This FAST code can predict both extreme and fatigue loads. The tower and flexible blade are modeled by using the assumed mode method. Other components are modeled as rigid bodies. An advanced certified code is used in FAST to model the aerodynamic behavior of the WT. WT loads are calculated by using BEM (blade element momentum) and multiple component of wind speed profile [28]. FAST code is approved by the Germanischer Lloyd (GL) Windenergie GmbH for calculating onshore WT loads for design and certification [29]. Due to the above advantages and exact nonlinear modeling of the WT, the proposed controllers are validated by using FAST. In general three blade turbine has 24 DOF (degree of freedom) to represent the wind turbine dynamics. In this work 3 DOF is considered for WT, that is, variable generator, rotor speed, and blade teeter. FAST codes are interface with -function and implemented with Simulink model. FAST uses an AeroDyn file as an input for aerodynamic part. AeroDyn file contains aerodynamic analysis routine and it requires status of a WT from the dynamic analysis routine and returns the aerodynamic loads for each blade element to the dynamic routine [30]. Wind profile acts as the input file for AeroDyn. The wind input file is generated by using TurbSim which is developed by the NREL.

In order to analysis the performance of the WT three cases are chosen.

5.2. Below Rated Wind Speed (Region 2)

The test wind profile with full field turbulence is generated by using TurbSim developed by NREL. Figure 6 shows the hub height wind speed profile. In general any wind speed consists of two components, that is, mean wind speed and turbulence component. The test wind speed consists of 10 min dataset that was generated using Class A Kaimal turbulence spectra. It has the mean value of 7 m/s at the hub height, turbulence intensity of 25%, and normal IEC (International Electrotechnical Commission) turbulence type. The above wind speed is used as the excitation of WT.

The proposed and conventional controllers are implemented using FAST interface with MATLAB Simulink. Table 3 shows the comparisons of proposed and conventional controllers. The main objectives of the controllers are to maximize the energy capture with reduced stress on the drive train. The efficiency of the controllers is compared by using the following terms, that is, aerodynamic () and electrical () efficiency given in (53) where is the optimal aerodynamic power for the wind speed profile. The following objectives are used to measure the performance of the controllers.(1)Maximization of the power capture is evaluated by the aerodynamic and electrical efficiency which is defined in (53).(2)The reduced oscillation on the drive train and control torque smoothness are measured by the STD (Standard Deviation) and maximum value.

Figure 7 shows the estimation of effective wind speed by using the Modified Newton Raphson (MNR) estimator. The MNR estimator gives the correct reference which ensures the dynamic aspect of wind.

The rotor speed comparisons for FAST simulator are shown in Figures 8 and 9. The conventional controllers such as ATF and ISC are not able to track the optimal reference speed. ATF has only single tuning parameter, that is, , which cannot minimize the steady state error. In ISC, during fast transient wind speed, it introduces more power loss. More over these controllers are not robust with respect to high turbulence wind speed profile. To overcome the above drawbacks SMC and ISMC are proposed. Figure 9 shows the rotor speed comparisons for ISMC and SMC. It is found that ISMC and SMC track the reference rotor speed better than ISC and ATF.

Table 3 gives the performance analysis of all the conventional and proposed controllers. From Table 3, it is clear that the STD of and is less for ISMC compared to SMC and ATF. This ensures the smoothness of the control input and low speed shaft torque in ISMC compared SMC and ATF. Table data shows that ISC having lowest STD of and , but its efficiency is very low (69.73%) compared to all other controllers. So a trade-off should be made between the efficiency and the fatigue load on drive train.

To analyze the controller performances in a more detailed fashion, Figures 10 and 11 show the box plot for low speed shaft torque and generator torque with the mean, median, ±25% quartiles (notch boundaries), ±75% quartiles (box ends), ±95% bounds, and the outliers. From the size of the boxes shown, it is clear that the ISC experiences minimum variation than others. It ensures that ISC having the minimum transient load on the drive train; at the same time, we have seen from Table 3 that the efficiency of ISC is not comparable with other controllers. Comparing the box plot of ISMC and SMC, ISMC has less variation in low speed shaft torque and generator torque; this indicates smoothness of the controller and reduction in transient load.

Figure 12 shows the boxplot for rotor speed for FAST simulator. From this figure it is observed that ISMC and SMC have almost same variation in the reference rotor speed. Apart from ISMC and SMC, other controllers such as ATF and ISC are having more variations in reference speed.

The frequency analysis is carried out by using the PSD on the low speed shaft torque which is shown in Figure 13. As the SMC plot is completely above the ISMC plot, it can be concluded that low speed shaft torque variation is more for SMC than ISMC. This indicates that ISMC gives minimum excitation to the drive train compared to SMC.

Figure 14 shows comparison of low speed shaft torque and control torque for ATF, SMC, and ISMC by considering ISC as baseline control. All the low speed shaft torques (LSSTq) and control torques (ControlTq) are having higher values compared to baseline controller.

As shown in Figure 15 the SMC controller has improved power capture by 0.75% compared to ISMC. An intermediate tracking has been chosen and a compromise has been made between efficiency and load mitigation. From the above analysis and results given in Table 3 it is observed that even though SMC gives slightly more efficiency than ISMC, by considering transient load on drive train and smooth control input, ISMC found be optimal.

In order to avoid the torsional resonance mode by choosing the proper tracking dynamics a trade-off is made between power capture optimization with smooth control and reduced transient load on low speed shaft torque. A good dynamic tracking; that is, similar to WT fast dynamics gives better power capture but it requires more turbulence in control torque. Conversely slow tracking gives smooth control action with less power capture. The simulations are performed with different mean wind speed at region 2. The results are given in Tables 4 and 5. From these tables it observed that with an increase in mean wind speed the maximum value of the control input () also increases. In all the cases both SMC and ISMC controllers are having almost same efficiency but the transient load reduction is better for ISMC. As the mean wind speed increases the rate of increase of STD is more for SMC than ISMC. It observed that when the wind speed undergoes high variation the ISMC can produce better power capture with reduced transient load on the drive train.

5.3. Transition between Below and Above Rated Wind Speed (Region 2.5)

Figure 16(a) shows the test wind speed consists of 10 min dataset that was generated by the binary file format. In this wind profile the wind speed changes in step of every 50 sec starts from 6 m/sec to 18 m/sec. This shows both the above and below rated wind speed are included in the wind profile. For below rated wind speed torque control comes into action with constant pitch angle and above rated wind speed pitch control comes into action with rated torque. Figure 16(b) shows the generator speed for SMC and ISMC for below and above rated wind speed. Both controllers achieve the nominal value of the generator speed at 250 sec. The corresponding wind speed is around 11.5 m/sec which can be seen from Figure 16(a). As the wind speed approaches towards the rated speed the WT generator speed reaches the nominal value, that is, 167.55 rad/sec. Figure 16(c) shows the electrical power comparison for SMC and ISMC for the transition period. At region 2.5 ISMC can able to extract the more power than SMC with almost same mechanical stress on the drive train. Figure 16(d) shows the generator torque comparison in region 2.5 for SMC and ISMC. It can be observed that ISMC produces more generated torque compared to SMC in region 2.5, that is, 250 sec to 400 sec. As generator speed remains constant from 200 sec onwards it is obvious that ISMC captures more power compared to SMC. The STD of generated torque for SMC and ISMC is found to be 1.060 kNm and 1.07 kNm respectively which is almost same. Figure 16(e) shows the pitch angle comparison for SMC and ISMC at region 2.5. Pitch variation found to be more for SMC compared to ISMC. So the pitch actuator needs more control action for SMC + PI.

Figure 17(a) shows the test wind speed consists of 10 min dataset that was generated by the binary file format. In this wind profile the wind speed changes with vertical wind profile starts from 6 m/sec to 18 m/sec. This shows both the above and below rated wind speed are included in the wind profile. Figure 17(b) shows the generator speed for SMC and ISMC for below and above rated wind speed. Both controllers achieve the nominal value of the generator speed at 280 sec. The corresponding wind speed is around 11.5 m/sec which can be seen from Figure 17(a). As the wind speed approaches towards the rated speed the WT generator speed reaches the nominal value, that is, 167.55 rad/sec. Figure 17(c) shows the electrical power comparison for SMC and ISMC for the transition period. At region 2.5 ISMC can extract more power than SMC with almost the same mechanical stress on the drive train. Figure 17(d) shows the generator torque comparison in region 2.5 for SMC and ISMC. It can be observed that ISMC produces more generated torque compared to SMC in region 2.5, that is, 280 sec to 415 sec. As generator speed remains constant from 280 sec onwards it is obvious that ISMC captures more power compared to SMC. Figure 16(e) shows the pitch angle comparison for SMC and ISMC at region 2.5. Pitch variation was found to be more for SMC compared to ISMC. So the pitch actuator needs more control action for SMC + PI.

Figure 18 shows the rotor speed comparison of SMC and ISMC controller with the rotor inertia parameter uncertainty (+30%) and disturbance of 5/ kNm. From this figure it is concluded that with the disturbance and uncertainty ISMC can track the reference rotor speed which ensures the robustness with respect to uncertainty and disturbance.

6. Conclusion

In this paper a combination of linear and nonlinear control for VSVPWT has been proposed. The proposed nonlinear controller such as SMC and ISMC is designed for wide range of below rated wind speed profiles. The main aim of the controller is to capture the maximum power with reduced oscillation on the drive train at below rated speed. The simulation of the controllers is performed by NREL CART 3, 600 kW WT. From the analysis it is concluded that the proposed ISMC (region 2) with conventional PI controller (region 3) can achieve the maximum power in all the regions of wind speed.

Appendices

A. Two-Mass Model Parameters

Rotor radius:  m,air density:  kg/m3,rotor inertia:  kg·m2,generator inertia:  kg·m2,shaft damping coefficient:  Nm/rad,shaft stiffness coefficient:  Nm/rad,rotor friction coefficient:  Nm/rad/sec,generator friction coefficient:  Nm/rad/sec,gear ratio: ,cut in wind speed: 6 m/sec,rated wind speed: 13 m/sec,cut out wind speed: 25 m/sec.

B. Modified Newton Raphson Algorithm

The estimation of effective wind speed is done by using the modified Newton Raphson. At instant the effective wind speed is obtained from rotor speed and aerodynamic torque:Step 1: ,Step 2: ,Step 3: ,Step 4: Stop if or go to Step 2,Step 5: ,where is the result of the first iterations and is the sampled rate fixed here as 1 sec.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.