Load Mitigation and Optimal Power Capture for Variable Speed Wind Turbine in Region 2
This paper proposes the two nonlinear controllers for variable speed wind turbine (VSWT) operating at below rated wind speed. The objective of the controller is to maximize the energy capture from the wind with reduced oscillation on the drive train. The conventional controllers such as aerodynamic torque feedforward (ATF) and indirect speed control (ISC) are adapted initially, which introduce more power loss, and the dynamic aspects of WT are not considered. In order to overcome the above drawbacks, modified nonlinear static state with feedback estimator (MNSSFE) and terminal sliding mode controller (TSMC) based on Modified Newton Raphson (MNR) wind speed estimator are proposed. The proposed controllers are simulated with nonlinear FAST (fatigue, aerodynamics, structures, and turbulence) WT dynamic simulation for different mean wind speeds at below rated wind speed. The frequency analysis of the drive train torque is done by taking the power spectral density (PSD) of low speed shaft torque. From the result, it is found that a trade-off is to be maintained between the transient load on the drive train and maximum power capture.
In recent years, wind energy is one of the major renewable energy sources because of environmental, social, and economic benefits. The major classifications of wind turbines (WT) are fixed speed wind turbine (FSWT) and VSWT. Compared with FSWT, VSWT has many advantages such as improved energy capture, reduction in transient load, and better power conditioning . For any kind of WT, control strategies play a major role on WT characteristics and transient load to the network . In VSWT, the operating regions are classified into 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. In , the maximum power for VSWT is achieved by PI (proportional integral) controller, which is based on the fuzzy system. Error is taken as the input to the controller, that is, difference between the actual and optimal rotor speed, and the output of the controller is generator torque. Fuzzy logic systems (FLS) are used for tuning the PI controller gains for various wind speed. PI gains are optimized for different wind speed by particle swarm optimization (PSO). In , radial bias function neural network (RBFNN) and torque observer based control algorithm are used to control the WT for optimal energy capture. RBFNN is trained online by using MPSO (modified particle swarm optimization) training algorithm. In order to achieve the maximum power, the difference between the actual and optimal rotor speed is to be minimized. In , a new maximum adaptive algorithm for extraction of optimal power is proposed for small WT. Perturb and absorb scheme is adapted for different wind speed to obtain optimum relationship for regulating the maximum power point. In , two control strategies are developed for optimal power extraction with reduced mechanical stress. The first one is tracking controller with wind speed estimator which ensures the optimal angular speed of the rotor. In the second one, a robust power tracking is developed by nonhomogenous quasicontinuous high order sliding mode controller without considering wind velocity. Maximum power extraction from VSWT is achieved by a Takagi-Sugeno-Kang (TSK) fuzzy model which is based on data driven model . In TSK model, a combination of fuzzy clustering method and genetic algorithm (GA) is used for portioning the input-output space and least square (LS) algorithm is used for parameter estimation. Nonlinear static and dynamic state feedback linearization control are addressed in [8, 9], where both the single and two-mass model are taken into consideration and the wind speed is estimated by Newton Rapshon (NR) method. To accommodate the parameter uncertainty and robustness, a higher order sliding mode controller is proposed in , which ensures the stability of the controller in both regions, that is, below and above rated speed. Feedback torque control is applied for mathematical model FSWT for maximum power extraction . In order to achieve the maximum power point in the WT, FLC tuned by GA is discussed in . The width of the membership function in FLC is adjusted by GA. In , sliding mode controller (SMC) and integral sliding mode controller (ISMC) are designed for all the regions of variable speed variable pitch wind turbine (VSVPWT) with FAST simulator.
The objective of this paper is to prove the efficacy of nonlinear controllers which considers the dynamic aspect of the wind and aero turbine, without the wind speed measurement. Finally, the objective is to track the reference rotor speed asymptotically. This paper is organized as follows. The objective of the work is discussed in Section 2. Section 3 discusses the modeling of the two-mass model. The conventional and proposed controllers are discussed in Sections 4 and 5. In Section 6, FAST model results are analyzed. Finally, a conclusion is drawn from the obtained results in Section 7, which shows the proposed method is having better performance compared to other existing controllers.
2. Problem Formulation
Generally, WT is classified into two types, that is, fixed and variable speed WT. Variable speed WT has more advanced and flexible operation than fixed speed WT. Operating regions in variable speed WT are divided into three types. Figure 1 shows the various operating region in variable speed WT.
Region 1 represents the wind speed below the cut-in wind speed. Region 2 represents the wind speed between cut-in and cut-out. In this region, the main objective is to maximize the energy capture from the wind with reduced oscillation on the drive train. Region 3 describes the wind speed above the cut-out speed. In this region, pitch controller is used to maintain the WT at its rated power.
Figure 2 shows the WT control scheme. 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 (1) defines the reference/optimal rotor speed:
3. WT Model
A WT is a device which converts the kinetic energy of the wind into 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 equations with limited degree of freedom. In this paper, the control law is designed based on simplified mathematical model with the objective of optimal power capture at below rated wind speed and reduced oscillation of the drive train. The proposed controllers are tested with different wind profiles. Finally, the controllers are validated for FAST WT model. The parameters of the two-mass model are given in . Generally, VSWT system consists of the following components; that is, aerodynamics, drive trains, and generator are shown in Figure 3.
Equation (2) gives the nonlinear expression for aerodynamic power capture by the rotor:From (2), 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 (4):where is the torque coefficient given asSubstituting (6) in (5), we getIn the above equation, the nonlinear term is which can be approximated by the 5th-order polynomial given in the following:where to are the WT power coefficients.
The values of approximated coefficients are given in Table 1. Figure 4 shows the versus curve.
Figure 5 shows the two-mass model of the WT. Equation (9) 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 the following: Equation (11) 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 (11), the high speed shaft torque can be expressed as Putting the values of from (14) in (13), we get
4. Conventional Controllers
In order to compare the results of proposed and existing conventional controllers, a brief description of the well-known control techniques, that is, ISC and ATF, is 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 the following:wherewhere 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 . The control law is given in (18): 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; 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, MNSSFE and TSMC, are proposed.
4.1. 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 (8): The estimated wind speed can be obtained by solving (23) using MNR. The above equation has unique solution at below rated region. With known , the optimal rotor speed is calculated by using (1).
5. Proposed Nonlinear Controllers
5.1. Terminal Sliding Mode Control (TSMC) for Optimal Power Capture
Let us consider the linear sliding surface:The first-order derivative of the above equation can be obtained asA nonsingular terminal sliding mode manifold  is first designed aswhere is a design constant and and are the positive integer, which satisfy the following condition:The linear sliding mode is combined with nonsingular terminal manifold . After reaches zero in finite time, both and will also reach finite time; then, the tracking errors and can asymptotically converge to zero. Once reaches zero, it will stay on zero by using the control law. Then, the sliding surface will converge to zero in finite time. The total time from to can be calculated by using the equations , from which time taken from to is obtained as By taking the derivative of (26),The stability of the system is investigated by choosing the following Lyapunov function:By taking the derivative of the above equation, we will getTo make the controller more adaptive to uncertainty and disturbances, we introduce the parameter called , where contains both modelling error and external disturbances: From the two-mass model system equation, By substituting (33) in (31),According to the above equation, the terminal sliding mode control law is defined as By substituting (35) into (34), the following equation can be obtained. Since the lumped parameter uncertainty and disturbance are unknown in practical application and the upper bound is very difficult to determine, the adaptive control law is adapted for lumped uncertainty .
The Lyapunov candidate function is chosen as where and is a positive constant. By taking the derivative of (36),The adaptive control rule can be selected asAccording to the above equation, the terminal sliding mode control law is defined asIn order to make the chattering free control law, the signum function can be replaced by : Substituting (38) into (37), the following can be obtained:
5.2. Modified Nonlinear Static State Feedback Linearization with Estimator (MNSSFE) for Maximum Power Extraction
In , the authors have explained nonlinear static state feedback with estimator, where the second derivative of the rotor speed and the first derivative of the low speed shaft are considered to express the control law. In order to avoid higher derivatives and complex control law, a modified nonlinear static state feedback with estimator is proposed in this section. Rearranging the terms in (9), we will getBy using the relationship given in (14) and (15), we will get By separating the control input , finally, the control torque can be expressed aswhere is approximated by the new input :The first-order error dynamics can be written asFrom (45), the new input is defined asBy substituting in (46), we get the final control law for the WT two-mass model:
6. Validation Results
CARTs (Control Advanced Research Turbines) are located in the center of the national wind NREL (National Renewable Energy Laboratory), near Golden, Colorado. The CART3 is a three-bladed variable speed and variable pitch wind turbine and has a rating of 600 kW. It mainly consists of three parts, namely, the rotor, the tower, and the nacelle. The generator is connected to the grid through power electronics that can directly control generator torque . The power electronics consist of three-phase PWM (Pulse Width Modulation) converters with a constant dc link voltage. The main objective of the grid side converter is to maintain the dc link voltage constant [17, 18].
6.1. Simulation Using FAST Model
FAST was developed by the NREL; it is used for WT aeroelastic simulator. The modelling of two- and three-blade horizontal axis wind turbines (HAWT) is obtained by FAST. This FAST code can be able to predict extreme and fatigue loads. Tower and flexible blade server are modelled by “assumed mode method.” WT loads are calculated by using BEM (Blade Element Momentum) and multiple component of wind speed profile . FAST code is approved by the Germanischer Lioyd (GL) WindEnergie GmbH for calculating onshore WT loads for design and certification . Due to the above advantages and exact nonlinear modeling of the WT, the proposed controllers are validated by using FAST. In general, three-blade turbines have 24 DOF (degrees of freedom) to represent the wind turbine dynamics. In this work, 3 DOF are considered for WT, that is, generator, rotor speed, and blade teeter. FAST codes are interfaced with S-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 . Wind profile acts as the input file for AeroDyn. The wind input file is generated by using TurbSim which is developed by the NREL. 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. 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 the following: 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 (50).(2)The reduced oscillation on the drive train and control torque smoothness are measured by the STD (standard deviation) and maximum value.The abovementioned values for all the controllers are given in Table 2. The rotor speed comparisons for FAST simulator are shown in Figures 7 and 8. 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 allows reducing the steady state error. In ISC, during fast transient wind speed, it introduces more power loss. Moreover, these controllers are not robust with respect to high turbulence wind speed profile. To overcome the above drawbacks, TSMC and MNSSFE are proposed. Figure 8 shows the rotor speed comparisons for MNSSFE and TSMC. From this figure, it is clear that at initial wind condition 0–30 sec both the controllers are not able to track the optimal rotor speed due to the initial setting in the AeroDyn input file. At high wind speed variations 220, 350, and 550 sec, the TSMC is almost tracking the reference rotor speed compared to MNSSFE. Except MNSSFE and TSMC, all the other controllers are having more tracking error in rotor speed. To achieve more power capture, the rotor speed should closely track the optimal rotor speed.
Table 2 gives the performance analysis of all the conventional and proposed controllers. From Table 2, it is clear that the STD of and is the lowest for TSMC and the highest for ATF controller. This ensures that the smoothness of the control input in TSMC is better compared to other controllers. ISC has very less STD of and ; at the same time, the efficiency is very low compared to all the controllers. Also the ISC control only depends on the generator speed which is not accurate and for higher variation in wind speed and introduces significant power loss. So a trade-off should be made between the efficiency and the fatigue load on drive train. Except TSMC, ATF and MNSSFE are having more standard deviation which ensures more drive train transient load. Compared to ISC and ATF, the aerodynamic and electrical efficiency of the proposed TSMC and MNSSFE are better.
To analyze the controller performances in a more detailed fashion, Figures 9 and 10 show the box plot for low speed shaft torque and generator torque with the mean, median, quartiles (notch boundaries), quartiles (box ends), bounds, and the outliers. From the size of the boxes shown, it is clear that the ISC experiences minimum variation compared to others. It ensures that ISC has the minimum transient load on the drive train; at the same time, we can find from Table 2 that the efficiency of ISC is not comparable with other controllers. Comparing the box plot of TSMC and MNSSFE, with TSMC having less variation in low speed shaft torque and generator torque, this indicates smoothness of the controller and reduction in transient load.
Figure 11 shows the box plot for rotor speed for FAST simulator. From this figure, it is clear that TSMC has almost the same variation as the variation in the reference rotor speed. It is observed that, apart from TSMC and MNSSFE, other controllers such as ATF and ISC are having more variations with respect to reference speed. This indicates that, for ATF and ISC, the obtained rotor speed is not able to track the reference rotor speed.
The frequency analysis is carried out by using the PSD on the low speed shaft torque which is shown in Figure 12. As the MNSSFE plot is completely above the TSMC plot, it is clear that low speed shaft torque variation is more for MNSSFE than TSMC. This ensures that TSMC gives minimum excitation to the drive train.
As shown in Figure 13, the MNSSFE controller has improved power capture by 2.53% compared to TSMC. Fast tracking introduces more variation in control input and drive train. So an intermediate tracking has to be chosen and a compromise has been made between efficiency and load mitigation. From this analysis, even though MNSSFE gives a little better efficiency than TSMC, by considering transient load on drive train and smooth control input, TSMC is found to 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 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. Therefore, a compromise should be made between the power capture and transient load reduction. The better optimal speed tracking leads to better power capture for TSMC controller. The simulations are performed with different wind speed profiles with the mean wind speed at below rated wind speed. The results are given in Tables 3 and 4. From these tables, it is observed that, with an increase in mean wind speed, the maximum value of the control input () also increases. In all the cases, both TSMC and MNSSFE controllers are having almost the same efficiency but the transient load reduction is better for TSMC. As the mean wind speed increases, the standard deviation also increases for MNSSFE compared with TSMC. It is observed that when the wind speed undergoes high variation the TSMC can be able to produce better power capture with reduced transient load on the drive train.
Figure 14 shows the electrical power comparison for industrial baseline controller, with MNSSFE and TSMC. From this figure, it is observed that industrial baseline controller has more oscillation compared to MNSSFE and TSMC. Both MNSSFE and TSMC are almost having the same power. From Table 5, it is found that industrial baseline has more oscillation in control torque compared with other controllers.
Figure 15 shows the rotor speed comparison for MNSSFE and TSMC with constant additive disturbance of 1 kNm/. From this figure, it is clear that both controllers are robust with respect to disturbance. With reference to Table 6, the STD is less for TSMC as compared to MNSSFE; at the same time, the efficiency of these controllers is similar.
This paper deals with the problem of controlling the maximum power generation at below rated wind speed of VSWT. The objective is to design a robust controller that maximizes the energy extraction from the wind while reducing the transient loads. For the above purpose, two nonlinear controllers, that is, TSMC and MNSSFE, which have the ability to reject disturbance and accommodate parameter uncertainty, are proposed in this study. Finally, it is concluded that a trade-off is to be maintained between the efficiency and mechanical stress on the drive train. The performances of these controllers are compared with the conventional ATF and ISC using FAST aeroelastic simulator. The proposed controllers are found to produce satisfactory results in achieving the control objectives.
|:||Low speed shaft stiffness (Nmrad−1)|
|:||Generator inertia (kg·m2)|
|:||Rotor inertia (kg·m2)|
|:||Generator external damping (N·m·rad−1·s−1)|
|:||Low speed shaft damping (N·m·rad−1·s−1)|
|:||Rotor external damping (N·m·rad−1·s−1)|
|:||Aerodynamic power (W)|
|:||Electrical power (W)|
|:||Rotor radius (m)|
|:||Aerodynamic torque (N·m)|
|:||Generator (electromagnetic) torque (N·m)|
|:||High speed shaft torque (N·m)|
|:||Low speed shaft torque (N·m).|
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
T. Burton, D. Sharpe, N. Jenkins, and E. Bossanyi, Wind Energy Handbook, Wiley Publications, New York, NY, USA, 2001.
F. D. Bianchi, F. D. Battista, and R. J. Mantz, Wind Turbine Control Systems: Principles, Modelling and Gain Scheduling Design, Springer, Buenos Aires, Argentina, 2nd edition, 2006.
M. Sheikhan, R. Shahnazi, and A. N. Yousefi, “An optimal fuzzy PI controller to capture the maximum power for variable-speed wind turbines,” Neural Computing and Applications, vol. 23, no. 5, pp. 1359–1368, 2013.View at: Publisher Site | Google Scholar
C.-M. Hong, C.-H. Chen, and C.-S. Tu, “Maximum power point tracking-based control algorithm for PMSG wind generation system without mechanical sensors,” Energy Conversion and Management, vol. 69, pp. 58–67, 2013.View at: Publisher Site | Google Scholar
I. Kortabarria, J. Andreu, I. M. de Alegría, J. Jiménez, J. I. Gárate, and E. Robles, “A novel adaptative maximum power point tracking algorithm for small wind turbines,” Renewable Energy, vol. 63, pp. 785–796, 2014.View at: Publisher Site | Google Scholar
J. Mérida, L. T. Aguilar, and J. Dávila, “Analysis and synthesis of sliding mode control for large scale variable speed wind turbine for power optimization,” Renewable Energy, vol. 71, pp. 715–728, 2014.View at: Publisher Site | Google Scholar
V. Calderaro, V. Galdi, A. Piccolo, and P. Siano, “A fuzzy controller for maximum energy extraction from variable speed wind power generation systems,” Electric Power Systems Research, vol. 78, no. 6, pp. 1109–1118, 2008.View at: Publisher Site | Google Scholar
B. Boukhezzar, H. Siguerdidjane, and M. Maureen Hand, “Nonlinear control of variable-speed wind turbines for generator torque limiting and power optimization,” Journal of Solar Energy Engineering, vol. 128, no. 4, pp. 516–530, 2007.View at: Publisher Site | Google Scholar
B. Boukhezzar and H. Siguerdidjane, “Nonlinear control of a variable-speed wind turbine using a two-mass model,” IEEE Transactions on Energy Conversion, vol. 26, no. 1, pp. 149–162, 2011.View at: Publisher Site | Google Scholar
B. Beltran, T. Ahmed-Ali, and M. E. H. Benbouzid, “High-order sliding-mode control of variable-speed wind turbines,” IEEE Transactions on Industrial Electronics, vol. 56, no. 9, pp. 3314–3321, 2009.View at: Publisher Site | Google Scholar
M. Liao, L. Dong, L. Jin, and S. Wang, “Study on rotational speed feedback torque control for wind turbine generator system,” in Proceedings of the International Conference on Energy and Environment Technology (ICEET '09), vol. 1, pp. 853–856, Guilin, China, October 2009.View at: Publisher Site | Google Scholar
H. M. Amine, H. Abdelaziz, and E. Najib, “Wind turbine maximum power point tracking using FLC tuned with GA,” Energy Procedia, vol. 62, pp. 364–373, 2014.View at: Publisher Site | Google Scholar
S. Rajendran and D. Jena, “Control of variable speed variable pitch wind turbine at above and below rated wind speed,” Journal of Wind Energy, vol. 2014, Article ID 709128, 14 pages, 2014.View at: Publisher Site | Google Scholar
H. Vihriala, R. Perela, P. Makila, and L. Soderlund, “A gearless wind power drive: part 2: performance of control system,” in Proceedings of the Wind Energy for the New Millennium European Conference (EWCE '01), pp. 1090–1093, Copenhagen, Denmark, July 2001.View at: Google Scholar
S. Mondal and C. Mahanta, “Adaptive second order terminal sliding mode controller for robotic manipulators,” Journal of the Franklin Institute, vol. 351, no. 4, pp. 2356–2377, 2014.View at: Publisher Site | Google Scholar | MathSciNet
L. J. Fingersh and K. Johnson, “Controls advanced research turbine (CART) commissioning and baseline data collection,” Tech. Rep., National Renewable Energy Laboratory (NREL), 2002.View at: Google Scholar
R. Ottersten, On control of back-to-back converters and sensorless induction machine drives [Ph.D. thesis], Chalmers University of Technology, Gothenburg, Sweden, 2003.
R. Peña, R. Cardenas, R. Blasco, G. Asher, and J. Clare, “A cage induction generator using back to back PWM converters for variable speed grid connected wind energy system,” in Proceedings of the 27th Annual Conference of the IEEE Industrial Electronics Society (IECON '01), vol. 2, pp. 1376–1381, IEEE, Denver, Colo, USA, December 2001.View at: Publisher Site | Google Scholar
M. O. L. Hansen, J. N. Sørensen, S. Voutsinas, N. Sørensen, and H. A. Madsen, “State of the art in wind turbine aerodynamics and aero elasticity,” Progress in Aerospace Sciences, vol. 42, no. 4, pp. 285–330, 2006.View at: Publisher Site | Google Scholar
A. Manjock, “Design codes FAST and ADAMS for load calculations of onshore wind turbines,” Tech. Rep., National Renewable Energy Laboratory (NREL), Golden, Colo, USA, 2005.View at: Google Scholar
D. J. Laino and A. C. Hansen, “User's guide to the wind turbine aerodynamics computer software aerodyn,” Tech. Rep., National Wind Technology Center, 2003.View at: Google Scholar