Wind Turbine Power Maximization Using Log-Power Proportional-Integral Extremum Seeking

: This paper proposes a Log-Power Proportional-Integral Extremum Seeking Control (LP-PIESC) framework for maximizing the power capture of a wind turbine operating at below-rated wind speeds, i.e., the so-called region-2 of a turbine’s power curve. Extremum seeking control (ESC) has emerged as a viable algorithm to maximize energy capture for a wind turbine operating in region-2. Despite the encouraging results of early ESC strategies, the basic algorithm suffers from slow and inconsistent convergence behavior under changing wind speed within region-2. It has been shown that replacing the power signal with its logarithm results in an algorithm that is robust and predictable even when the mean wind speed varies. In addition, new studies have suggested that replacing conventional ESC with proportional plus integral ESC (PIESC) results in faster convergence to optimal conditions. In the current paper, the idea of log-power feedback is merged with the PIESC scheme and is applied to tune the parameters of the region-2 torque controller for the NREL 5-MW turbine reference model. The results of this new algorithm are compared with the ESC with log-of-power feedback using NREL OpenFAST simulations. The log-power feedback PIESC is also implemented for the blade pitch set-point angle. Energy capture over the course of the simulations and damage equivalent loads calculated with MLife are used to assess the results. The simulations performed under different turbulent intensity cases demonstrate the rapid convergence of the log-power feedback PIESC.


Introduction
Wind turbines extract power in two different wind regimes: above-rated wind speed when the machine is regulated to produce its rated power while maintaining the rated rotor angular speed despite wind variations, and below-rated wind speed when the angular rotor speed is actively controlled to maximize the power extracted from the wind. These two modes of operation are loosely referred to as region-3 (above-rated wind condition) and region-2 (below-rated wind conditions).
The rotor power in region-2 is proportional to the available wind power, and the proportionality parameter, i.e., the power coefficient, or aerodynamic efficiency, is optimized to maximize the output power. In variable-speed variable-pitch turbines, the power coefficient (C p ) is a unimodal function of the tip-speed ratio (λ) and the blade-pitch angle (β) [1], which implies there is an optimal tip-speed ratio and blade-pitch angle for achieving the optimal C p and thus maximum output power. In practice, below-rated (region-2) operation entails keeping the blade-pitch angle constant at its ideal value β opt while adjusting the rotor speed to maintain the optimal tip-speed ratio λ opt despite wind changes [1,2]. Optimal tip-speed ratios and blade-pitch angles can differ from the manufacturer's specification due to site-specific complex aerodynamic effects and other factors. Furthermore, the optimal settings for these turbine parameters may change over time because of the blade wear and tear, variations in environmental conditions, etc.
Extremum Seeking Control (ESC) has emerged as an effective real-time optimization strategy to maximize power capture in region-2 despite a lack of knowledge of optimal turbine parameters. The ESC is a gradient-based steepest ascent search algorithm to find (local) optima in real-time. ESC-based region-2 control for tuning the torque gain, blade pitch angle, and yaw angle was investigated by Creaby et al. [3]. They used rotor power as the feedback. Xiao et al. [4] provided field test results in region-2 using the NREL CART3 turbine. Despite the encouraging outcomes of these ESC investigations, one disadvantage of these techniques is inconsistent convergence rates to the optimal control settings as wind speed changes within region-2. Rotea [5] analyzed the root cause of inconsistent convergence; he proposed a "log-power feedback strategy" for ESC-based wind turbine control. The improved and consistent convergence of the log-power ESC was demonstrated in [5] using a simplified reduced-order model. Afterwards, it was shown through highaccuracy large-eddy simulations that the log-power feedback ESC calibrated at a given wind speed exhibits consistent and robust performance across all wind speeds within region-2 under realistic turbulent wind conditions [6].
Guay et al. [7] recently proposed an ESC algorithm that estimates the gradient using a variant of the recursive least-squares algorithm for parameter estimation. The technique in [7] introduces a proportional-integral approach to improve the transient performance of the extremum-seeking controller. The algorithm in [7] shall be referred to as PIESC. This algorithm is faster than the conventional ESC algorithm.
In this paper, we propose a log-of-power feedback PIESC scheme for region-2 power maximization. The algorithm adjust the torque gain parameter k g in the "classical τ g = k g ω 2 g " (τ g is the reference generator torque, ω g is a filtered version of the measured generator speed) control law [1] and the blade pitch angle β for region-2 power maximization individually. Simulation results demonstrate that the combination of the log-of-power transformation with the PIESC scheme yields a fast algorithm for optimal tuning of the torque gain with a less invasive dither signal than currently possible with the conventional log-of-power ESC. We also show results for the optimal tuning of the blade pitch angle with the log-of-power PIESC scheme. Such tuning is useful when changes in blade aerodynamics or the environment result in a different optimal blade pitch angle than the one provided by the manufacturer. This paper is an extension of the work in [8], which contained initial comparative results for the LP-PIESC and LP-ESC (i.e., a conventional ESC driven by the log-of-power signal) for torque-gain k g tuning only. Here, we extend these initial results to provide a detailed analysis of energy capture and loads when the LP-ESC and LP-PIESC are used to determine the optimal torque gain. This detailed comparative study demonstrates the faster convergence of the LP-PIESC over the LP-ESC. In addition to faster convergence, we show that the LP-PIESC requires a smaller dither amplitude than the LP-PIESC. We also introduce the use of the LP-PIESC to determine the optimal blade pitch angle in below-rated wind conditions and demonstrate its rapid convergence to the optimal angle. All simulations are done using the NREL 5-MW reference turbine [2] for a range of wind speeds and turbulence intensities. In order to focus on the evaluations in belowrated wind speeds, we do not implement the baseline NREL 5-MW region switching strategy; instead, we allow the turbine to operate in region-2 for the wind speeds considered. The torque control always follows the τ g = k g ω 2 g control law for all wind speeds and turbulence intensities.
The organization of the paper is as follows. Background information for the PIESC algorithm and log-of-power feedback is given in Section 2. Section 3 describes the implementation, design parameters, and simulation results of the LP-ESC and LP-PIESC for torque gain control. The LP-PIESC implementation, design parameters, and simulation results for the case of blade pitch angle control are shown in Section 4. Final remarks are given in Section 5.

Log-Power Maximization Using Gradient-Based Algorithms
Extremum Seeking Control (ESC) and PIESC are gradient-based algorithms that tune parameters in real-time in order to maximize a system's performance index. These algorithms do not require comprehensive physical models to work. The gradient of the performance index with respect to the parameters is estimated by perturbing the tunable parameters with a zero-mean periodic dither signal. By demodulating a measurement of the perturbed performance index (ESC) or by using variants of recursive least squares (PIESC), a gradient estimate can be obtained. Then, the values of the tunable parameters are updated according to the estimated gradient. The optimum is reached asymptotically, starting from a suboptimal guess of the parameter values.
Rotea [5] conducted a study for power maximization with gradient-based algorithms that demonstrates the advantages of using the logarithm of the power signal rather than the power itself as the performance index. The intuition behind the result in [5] is relatively simple to explain, and it is summarized as follows.
The power available in the wind is given by where R is rotor radius in m, ρ is air density in kg/m 3 , and V is the incoming wind speed perpendicular to the rotor plane in m/s. The rotor power produced can then be given by where C p denotes the non-dimensional power coefficient of the wind turbine, and u is a control parameter to adjust the power coefficient. In practice, u can be the torque gain of the generator torque control signal, the blade pitch angle, or some other controllable parameter influencing C p (u). The goal is to maximize C p by proper selection of an optimal u-value using measurements of the power P only. When the rotor power is the performance index, the gradient of the performance index (rotor power) with respect to u is which is proportional to the power available in the wind P wind . However, when the logarithm of the power signal is utilized as the performance index, i.e., J = ln P, the gradient becomes ∂J ∂u which does not depend on the available wind power P wind , and it is only a function of the power coefficient. This result, however simple, has remarkable implications. Note from (1) and (3) that the convergence properties of gradient-based algorithms using a direct power measurement is a function of V 3 . Thus, without re-tuning, such algorithms will exhibit slow convergence at low wind speeds and may become unstable at high wind speeds [5,6]. On the other hand, if a gradient-based algorithm is used to maximize the log-power performance index, then it follows from (4) that convergence and stability are independent of the wind speed V. This results in an easier-to-tune optimization algorithm which performs consistently throughout region-2. This idea from in [5] is patentpending [9]. The consistency and robustness of ESC with the log-power performance index are demonstrated in [6] using detailed large-eddy simulations.

LP-ESC and LP-PIESC
In this paper, we use conventional extremum seeking control (ESC) in continuous time with the logarithm of the rotor power as the performance index to be maximized following the design methodology in [4]. We shall refer to this algorithm as the LP-ESC. A typical LP-ESC is shown in Figure 1, where F in (s) and F out (s) are unity-gain linear timeinvariant (LTI) approximation of the input and output dynamics, respectively, and f (u) represents a static map of the performance index (rotor power). A periodic dither signal d 1 (t) = a sin(ωt) and a demodulation signal d 2 (t) = sin(ωt + θ) are applied to estimate the gradient ∂ f /∂u, where a is the dither amplitude, ω is the dither frequency, and θ represents the phase compensation angle. The high-pass filter F HP (s) is applied to the logarithm of the performance index ln(y) to eliminate the steady-state term of the objective function. Following the multiplication by d 2 (t), a low-pass filter F LP (s) extracts the gradient information. The LP-ESC loop is closed by an integrator with gain k to drive ∂ f /∂u to zero asymptotically. Our goal is to develop an improved version of the LP-ESC using the proportional plus integral ESC (PIESC) developed by Guay et al. [7]. We shall refer to this gradient-based algorithm as the LP-PIESC.
The LP-PIESC holds the promise of faster convergence than the LP-ESC due to two main modifications: the use of a recursive least squares-like method to estimate the gradient and the addition of a proportional term to accelerate convergence to the optimum. In the proposed proportional-integral extremum seeking control, the gradient estimation is done using parameter estimation instead of the conventional dither-demodulation technique. The development of the LP-PIESC is based on a time-varying estimate of the gradient of the cost and a PI control law to drive the system to its optimal operating point.
The proportional-integral extremum seeking controller proposed in [7] is given by (5) and can be represented as in Figure 2: where k p is the proportional gain, τ I is the time constant, andθ 1 is proportional to the gradient ∂ f /∂u. The parameter to be estimated as presented in [7] is θ = [θ 0 , θ 1 ] T , where the superscript T denotes transpose. In our application, θ is a real-valued vector of dimension two. Only the scalar component θ 1 contains the gradient information we need. The parameter θ 0 is a bias component needed for the parameter estimation to work properly. The parameter estimation update law is given bẏθ where Proj(·) is a Lipschitz projection operator designed to ensure that the estimates are bounded within the constraint set and guarantee stability. This projection algorithm was implemented as discussed in Appendix E in [10] and the constraint set adaptation was adopted as per Adetola et al. [11]. The forcing term in the differential Equation (6) is the prediction error e given by e = y −ŷ (7) where y is the logarithm of the performance index to be optimized andŷ represents its prediction. The equation for the predicted outputŷ is given by the following dynamics:ŷ where φ = [1, (u −û)] T and K is a positive constant to be determined. The time-varying parameter c in (8) is the solution of the differential equatioṅ An auxiliary variable η is introduced and its filtered estimate is given bẏη The matrix Σ is the solution of the matrix differential equatioṅ with Σ(0) = αI and α > 0, σ > 0, k T > 0, i.e., they are positive constants. The inverse of Σ is given by the solution of the matrix differential equatioṅ with Σ −1 (0) = 1 α I. The estimation algorithm propagates Σ −1 , which is used in the parameter dynamics (6).

Implementation of LP-ESC and LP-PIESC for Torque Control in Region-2
The 5MW NREL reference turbine model, simulated with OpenFAST, is used in this study [12]. Table 1 lists the major parameters of this turbine model. The generator torque (τ g ) is generally used as the control variable to maximize power in region-2 and is computed as a tabulated function of the filtered generator speed in NREL's baseline control strategy. The main control law to maximize power below-rated wind speeds is given by where ω g denotes the generator speed and k g is the generator torque gain used as the control parameter for the extremum seeking algorithms. The theoretical optimal value of k g is obtained from [1] k opt = 1 2 where A is the cross-sectional area swept by the rotor and ρ is the air density given by 1.22 kg/m 3 . For this turbine, the theoretical optimal tip-speed ratio is λ opt = 7.55 and the corresponding maximum power coefficient is C p,max = 0.48 [2]. Then, the optimal generator torque gain takes the value k opt = 2.34 Nm/(rad/s) 2 .
The control parameter u is the non-dimensional normalized torque gain introduced in [5]. The optimal normalized torque gain value is A high-level block diagram of the region-2 torque gain LP-ESC/LP-PIESC implementation is depicted in Figure 3. Input to the LP-ESC/LP-PIESC is the logarithm of the rotor power P normalized with respect to the rated power P r = 5 MW. A rate limit of 15 kNm/s is applied on the torque reference as per the definition of the 5-MW turbine to avoid sudden change in the torque command [2].

Parameter Design
LP-ESC parameters are designed following the guidelines in [13,14]. The dither frequency is chosen within the bandwidth of the plant dynamics. The rotor inertia and the actuator dynamics yield the input dynamics, while the output dynamics is due to the sensor dynamics and/or signal conditioning. The input and output dynamics are merged at the plant input and approximated using open-loop step responses under constant wind input to simplify the algorithm design. The response of the rotor speed (ω r ) under staircase step changes in normalized torque-gain indicates a first-order dynamics, as shown in Figure 4. The test was performed for the hub-height mean wind speed of 8 m/s with no turbulence. The estimated time constant from the normalized torque-gain to the rotor speed ranges from 7 to 8 s. The largest time constant was adopted for ESC design, i.e., 8 s. The corresponding bandwidth of the combined input-output dynamic for the torque-gain ESC is 0.125 rad/s. As dither frequency should be selected within the estimated bandwidth, it was selected as 0.02 rad/s. Figure 5 shows the Bode plots for the estimated plant dynamic and the filters in the ESC design. The low-pass F LP and high-pass F HP filters (see Figure 1) are first-order filters. The dither amplitude was selected using trial and error. The dither amplitude is about 33% of the optimal normalized torque gain in (15). The integrator gain k was designed to achieve a settling time T s = 15 min, approximately for the 8 m/s wind speed with 10% TI case and then the same integrator gain was used for all other cases. The settling time T s is defined as the time it takes to reach optimum torque gain value for the first time after turning on the dithering (control on). The LP-ESC parameters used are listed in Table 2.   Table 2. LP-ESC Parameters of NREL 5MW Region 2 controller (see Figure 1). The LP-PIESC scheme was designed to have the same dither frequency as the LP-ESC. After trial and error, the dither amplitude could be reduced significantly to 25% of the dither amplitude value of LP-ESC. This result is important because the external dither is a perturbation to the overall turbine. One plausible explanation for the smaller dither amplitude is that the method for gradient estimation in LP-PIESC, Equations (6)- (12), is more sophisticated than the one used in the conventional LP-ESC. The parameters for the LP-PIESC scheme are listed in Table 3. Table 3. LP-PIESC Parameters of NREL 5 MW Region 2 controller. See Equations (6)- (12).

Simulation Results
The LP-ESC and the LP-PIESC controllers with the algorithm parameters shown in Tables 2 and 3 are evaluated with OpenFAST simulations for different hub-height mean wind speeds of 4 m/s, 8 m/s, and 12 m/s, no wind shear, and under turbulence intensities (TI) of 10% and 15%, respectively.
The turbulent wind profiles were generated using NREL TurbSim [15]. TurbSim uses IEC 61400-1 [16] to generate the wind profiles. In order to generate the wind input files for our simulations, we specify the following five parameters in the meteorological boundary conditions section of the TurbSim input file (defined in Appendix A of [15]): (1) Turbulence model is set as Kaimal (IECKAI), (2) IEC turbulence type is chosen to be the Normal Turbulence Model (NTM), (3) hub-height of 90 m for NREL 5MW reference turbine model, (4) mean wind speed, and (5) percentage turbulence intensity. Other parameters were left to their default values in TurbSim. Mathematical expression for the Kaimal spectrum can be found in Appendix C.3 of [17]. Figure 6 shows the hub-height wind time series for the three mean wind speeds (i.e., 4 m/s, 8 m/s, and 12 m/s) used in the simulation study with 10% TI. These wind profiles are then used to evaluate the performance of LP-ESC and LP-PIESC in Figures 7 and 8, respectively. The top plot shows the normalized torque gain (u), the middle plot shows the tip speed ratio (λ), and the estimated power coefficient (C p ) is shown in the bottom plot. Tip speed ratio (λ) and the estimated power coefficient (C p ) plots are obtained by applying a 100 s moving average filter on the OpenFAST outputs for these variables. The vertical dashed lines in Figure 7 show the time LP-ESC is turned on and the time for practical convergence to the optimum power coefficient C p . Both LP-ESC and LP-PIESC are turned on at 150 s.
It can be observed that the LP-ESC takes approximately 850 s (i.e., ≈15 min as designed) to converge to the optimum torque gain for 8 m/s and 12 m/s mean wind speeds. In stark contrast, the LP-PIESC converges to the optimum almost instantaneously for all cases.
Simulations with increased turbulence intensity (15% TI) are also used to evaluate the performance of LP-PIESC. The simulation results are shown in Figure 9. It can be observed that increasing the turbulent intensity does not affect the transient or steady-state performance of the LP-PIESC.  Table 2 and hub-height wind from  Table 3 and hub-height wind from Figure 6. Normalized torque gain u (top), tip speed ratio λ (middle), estimated power coefficient C p (bottom). The dashed horizontal lines indicate optimal parameters u opt , λ opt and C P,max . The LP-PIESC is turned on at 150 s.  Table 3 and wind input with 15% TI. Normalized torque gain u (top), tip speed ratio λ (middle), and estimated power coefficient C p (bottom). The dashed horizontal lines indicate optimal parameters u opt , λ opt and C P,max . The LP-PIESC is turned on at 150 s.
The results in this section provide evidence that LP-PIESC is faster to converge to the optimum than LP-ESC and also robust to the variations in mean wind speed and turbulence intensities. In addition, the LP-PIESC requires a dither amplitude significantly lower than that for LP-ESC.

Energy Capture
The energy capture for both the LP-ESC and the LP-PIESC controllers is evaluated and compared with the baseline controller (i.e., optimal torque gain value at all times). The controllers are evaluated for the hub-height mean wind speeds of 4 m/s, 8 m/s, and 12 m/s, no wind shear and turbulence intensities (TI) of 10% and 15%, respectively. TurbSim [15] is used to generate wind profiles from six distinct seeds for each turbulent intensity case. The average energy capture over those six wind profiles will be presented here. All the calculations are done using data from the time the extremum seeking controllers are turned on (150 s) till the end of the simulation (3000 s).
First, the LP-ESC is compared with the baseline controller. The result for average energy output can be seen in Figure 10. Using LP-ESC results in an average energy loss of approximately 13% (TI = 10%) and 16% (TI = 15%), respectively.  The energy capture with LP-PIESC was also compared with the baseline controller. The average change in energy capture is shown in Figure 11.  The results suggest that the LP-PIESC gives an energy output very close to the baseline controller, where the optimal tuning parameter (torque gain for our study) value is known. The maximum energy loss was 0.7% (reduction of around 1 kWh) which was seen for the low wind speed of 4 m/s. Averaging results across all wind conditions for both LP-ESC and the LP-PIESC, we obtain Table 4.  Table 4 suggests that due to the faster convergence of the torque gain to its optimal value, the LP-PIESC could produce nearly optimal energy output. That is, the LP-PIESC could reach very close to the maximum energy output without knowing the optimal setting. This optimal setting in the baseline controller could change over time due to the wear and tear of blade surfaces and variations in atmospheric and environmental conditions such as icing, insects, etc. In this case, the baseline controller can result in a loss of energy. The LP-PIESC could aid in finding the new optimal torque gain value without any significant loss in energy capture.

Structural Loads Analysis
This section shows the change in loads (tower, drive train, blades) due to the change in torque gain imparted by the LP-PIESC. We investigate the impact on transient loads caused by the sudden change in torque gain as well as the impact of the dither signal on the wind turbine loads. The loads considered are shown in Table 5, along with the natural frequencies of the relevant structural modes obtained from [2,18]. Table 5. Blade, drivetrain, and tower loads considered for analysis with respective natural frequencies.

Load First Mode Natural Frequency (Hz)
Blade-root flapwise bending moment (BRFW) 0.5-1 Blade-root edgewise bending moment (BREW)) 0.9-1. TurbSim [15] is used to generate wind profiles from six distinct seeds for each turbulent intensity case. The time series of the loads are analyzed with MLife to calculate aggregate damage equivalent loads (DELs) [19]. The DELs obtained from the time series of the LP-PIESC are then compared against those obtained from the baseline controller. All the evaluations are done on the data from 150 s till the simulation end time, i.e., 3000 s.
As shown in Figures 8 and 9, when the LP-PIESC is turned on, there is a sudden change in torque gain from the initial setting to the optimal value. This sudden change will propagate to a step-change in the generator torque because the rotor speed cannot change instantaneously. Figure 12 illustrates this situation for the wind speeds shown in Figure 6. The response of the baseline constant torque gain is shown as well. As expected, the generator torque and the rotor angular speed converges to the optimal time series, eventually.

Drivetrain Loads
The first turbine component affected by sudden changes in generator torque is expected to be the drivetrain. DEL result obtained from the time series data of drivetrain torsion (DTT) from ESC start, i.e., 150 s to simulation end time, i.e., 3000 s is as shown in Figure 13. DTT does not show much change for 8 m/s and 12 m/s wind speeds at 10% and 15% TI but increases for 4 m/s wind speed by 70% at 10% TI and 28% at 15% TI, respectively. The high percentage change is largely due to the smaller values at a low wind speed of 4 m/s. The actual difference in DEL values is 5.6 kNm for 10% TI and 3.4 kNm for 15% TI wind.
To better understand the large percentage change in the DTT DEL at 4 m/s, we examine the time series for the wind profile in Figure 6, which corresponds to one of the input wind files used to calculate damage equivalent loads. Figure 14 shows the time series of the drivetrain torsional moment. Note the large increase in the amplitude of the drivetrain torsion signal at 150 s (LP-PIESC turned on). These larger oscillations have a frequency of approximately 1.7 Hz, which corresponds to the drivetrain natural frequency shown in Table 5. This result suggests that the percentage change in DELs at low wind speed, when the dimensional moments are small, can be attributed to the natural transient response of the drive train to a sudden change in one of its inputs, i.e., the generator torque. To increase confidence in this explanation, the DELs are also computed using the time series from 300 s to the end of the simulation (3000 s). The results are shown in Figure 15. It can be observed that once the transients die out ( Figure 14) due to the sudden change in torque (Figure 12), the DELs obtained from the LP-PIESC are almost identical to that of the baseline controller, which supports our hypothesis. In addition, from Figure 15 we note the continuous presence of the dither signal does not affect the drive train DELs, that is, the DELs with the LP-PIESC are almost identical to the ones with the constant torque gain.

Tower Loads
Next, we studied the two fundamental mechanical moments that must be considered, one is related to the forward and backwards motion of the tower (tower-base fore-aft bending moment (TBFA)) and the other is related to the side-side tower motion (towerbase side-to-side bending moment (TBSS)). DEL results for the tower loads are shown in Figures 16 and 17.  For TBSS (Figure 17), there is no substantial change in DEL values at 8 m/s and 12 m/s. On the other hand, at 4 m/s DEL values increase 22% for TI = 10% and 9% for TI = 15%. Again, this percentage value is higher because of the smaller numbers at 4 m/s wind speed. That is, the actual difference in DEL values is 32 kNm (10% TI) and 19 kNm (15% TI). The TBSS experiences two effects: (i) as the LP-PIESC starts with a highly sub-optimal initial torque gain value, the initial transient is high, and it takes more than 150 s to converge, and (ii) the effect of switching on the LP-PIESC at 150 s. Therefore, examining the TBSS DEL after these two effects have died out would be meaningful as it is difficult to separate their effects. Therefore, we also evaluated the TBSS DELs using the time series from 300 s to the end of the simulation (3000 s). The results are shown in Figure 18. We can observe that the DEL values for both the LP-PIESC and the baseline are almost identical. The continuous presence of the dither signal in the LP-PIESC does not affect the TBSS DELs.

Implementation of LP-PIESC for Blade Pitch Tuning in Region-2
Given the promising results LP-PIESC showed for the torque gain control in the previous section, in this section, we design and evaluate the LP-PIESC algorithm for the control parameter of blade pitch angle, keeping the normalized torque gain at its optimal value as given in (15). As in the prior section, the design and analysis are done using the NREL 5 MW turbine model in OpenFAST. A high-level block diagram of the blade pitch tuning for region-2 is shown in Figure 22, where u is the blade pitch angle reference commanded by the LP-PIESC. A rate limit of 8 • /s is applied on the blade pitch angle reference as per the definition of the 5-MW reference turbine model. The optimal blade pitch angle β for this turbine is 0 • .

Parameter Design
As done in Section 3.1, LP-PIESC parameters are designed for blade pitch angle control. The response of the rotor speed under staircase step changes in the blade pitch angle indicates a first-order dynamics, as shown in Figure 23. The test was performed for the hub-height mean wind speed of 8 m/s with no turbulence.The estimated time constant from the blade pitch angle (β) to the rotor speed (ω r ) ranges from 7 to 10 s. The largest time constant was adopted for ESC design, i.e., 10 s. The corresponding bandwidth of the combined input-output dynamic for the blade-pitch ESC is 0.1 rad/s. As the dither frequency should be selected within the estimated bandwidth, it was selected as 0.02 rad/s. The Bode plot for the estimated plant dynamic is shown in Figure 24. The dither amplitude was selected using trial and error. The parameters for the LP-PIESC scheme are listed in Table 6.

Simulation Results
As done in Section 3.2, the LP-PIESC controller with the parameters shown in Table 6 is evaluated with OpenFAST simulations for different hub-height mean wind speeds of 4 m/s, 8 m/s, and 12 m/s, no wind shear and turbulence intensities (TI) of 10% and 15%, respectively. The wind profiles in Figure 6 with 10% TI are used to determine the performance of the LP-PIESC in Figure 25. The top plot shows the blade pitch angle (β), the middle plot shows the tip speed ratio (λ), and the estimated power coefficient (C p ) is shown in the bottom plot. Tip speed ratio (λ), and the estimated power coefficient (C p ) plots are obtained by applying a 100 s moving average filter on the OpenFAST outputs for these variables. The horizontal dashed lines indicate the optimal values of β, λ, and C p . The LP-PIESC is turned on at 150 s.
It can be observed that the LP-PIESC converges to the optimum blade pitch angle almost instantaneously for all cases. The tip speed ratio (λ) and the estimated power coefficient (C p ) converges in less than 100 s.  Table 6 and hub-height wind from Figure 6. Blade pitch angle β (top), tip speed ratio λ (middle), and estimated power coefficient C p (bottom). The dashed horizontal lines indicate optimal parameters β opt , λ opt and C P,max . The LP-PIESC is turned on at 150 s.
Simulations with increased turbulence intensity (15% TI) are also used to evaluate the performance of LP-PIESC. The simulation results are shown in Figure 26. It is observed that increasing the turbulent intensity does not affect the transient or steady-state performance of the LP-PIESC.
The results in this section provide further evidence that LP-PIESC has fast convergence to the optimum and is also robust to the variations in mean wind speed and turbulence intensities.  Table 6 and wind input with 15% TI. Blade pitch angle β (top), tip speed ratio λ (middle), estimated power coefficient C p (bottom). The dashed horizontal lines indicate optimal parameters β opt , λ opt and C P,max . The LP-PIESC is turned on at 150 s.

Energy Capture
The energy capture for the LP-PIESC controller is evaluated and compared with the baseline controller (i.e., 0 • optimal blade pitch angle at all times). The controller is evaluated for hub-height mean wind speeds of 4 m/s, 8 m/s, and 12 m/s, no wind shear and, turbulence intensities (TI) of 10% and 15%, respectively. TurbSim [15] is used to generate wind profiles from six distinct seeds for each turbulent intensity case. The average energy capture is calculated over those six wind profiles. All the calculations are done using data from the time the LP-PIESC is turned on (150 s) till the end of the simulation (3000 s).
The average change in energy capture is shown in Figure 27.  The results suggest that the LP-PIESC gives an energy output very close to the baseline controller, where the optimal tuning parameter (blade pitch angle) value is known. The maximum energy loss was 0.5% (reduction of around 25 kWh) which was seen for the high wind speed of 12 m/s (15% TI). Averaging results across all wind conditions, we obtain Table 7.  Table 7 suggests that due to the faster convergence of the blade pitch angle to its optimal value, the LP-PIESC could produce nearly optimal energy output. That is, the LP-PIESC could reach very close to the maximum energy output without knowing the optimal pitch angle setting. This optimal setting in the baseline controller could change over time due to the wear and tear of the blade surface, and variations in atmospheric and environmental conditions such as icing, insects, etc. Thus, whenever the current pitch angle is no longer optimal, the LP-PIESC can find the new optimal blade pitch angle without any significant loss in energy capture.

Structural Loads Analysis
This section shows the change in loads (blades, tower, drive train) due to the change in blade pitch angle imparted by the LP-PIESC. We investigate the impact on transient loads caused by the sudden change in the blade pitch angle as well as the effect of the dither signal on the wind turbine loads. The loads considered are shown in Table 5.
The LP-PIESC is evaluated for the hub-height mean wind speeds of 4 m/s, 8 m/s, and 12 m/s, no wind shear, and turbulence intensities (TI) of 10% and 15%, respectively. TurbSim [15] is used to generate wind profiles from six distinct seeds for each turbulent intensity case. The time series of the loads are analyzed with MLife to calculate aggregate damage equivalent loads (DELs) [19]. The DELs obtained from the time series of the LP-PIESC are then compared against those obtained from the baseline controller, i.e., β = 0 • . All the evaluations are done on the data from 150 s till the simulation end time of 3000 s. As shown in Figures 25 and 26, when the LP-PIESC is turned on, there is a sudden change in the blade pitch angle from the initial setting to the optimal value. This sudden change will propagate to a sudden change in the blade-root flap-wise bending moment (BRFW). Figure 28 illustrates this situation for the wind speeds shown in Figure 6.  Figure 6).
The response of the baseline constant blade pitch angle is shown as well. As expected, the BRFW and the rotor angular speed converges to the optimal time series, eventually.

Blade Loads
The first turbine component affected by the sudden change in the blade pitch angle should be the blades. For each wind condition, we calculate the average DEL over all three blades. The average DEL for the blade-root flapwise bending moment (BRFW) is shown in Figure 29. Using LP-PIESC, the BRFW DELs increase for all wind conditions compared to the baseline controller. This increase can be attributed to the sudden change in the blade pitch angle (from sub-optimal to optimal), which leads to the sudden change in the BRFW as shown in Figure 28. From Figure 29, we can conclude that the percentage DEL change decreases with the increase in mean wind speed or the increase in turbulence intensity. The average increase across the three wind speeds ranges from 13% (10% TI) to 5% (15% TI). This could be because, at higher turbulence, fluctuations due to the wind could be relatively more significant than the effect of a step change in the blade pitch angle. The highest percentage change for both the 10% TI and 15% TI wind can be observed at the low wind speed case of 4 m/s which is likely due to the small actual DEL values.
We also calculated BRFW DELs for the time series from 300 s (once the transients die out) till the simulation end time (3000 s). The results are shown in Figure 30. Now, the average change in BRFW DELs is just 0.3% (10% TI) and 0.4% (15% TI). Thus, the DEL values with the LP-PIESC are essentially identical to those with the baseline controller. Furthermore, the presence of the continuous dither signal in the LP-PIESC does not affect the BRFW DELs.    For both the TBFA and TBSS, there is no substantial change in DELs, which means tower loads are not affected by the LP-PIESC.

Drivetrain Loads
We also evaluated the drivetrain torsion for both the LP-PIESC and the baseline controller. DEL results for drivetrain torsion are shown in Figure 34. Similar to the tower loads, drivetrain torsion is not affected by either the sudden change in blade pitch angle or the dither signal of the LP-PIESC for all the wind conditions. The average percentage change at both 10% TI and 15% TI wind is 0.6%.

Conclusions
This paper presented a log-power feedback PIESC (LP-PIESC) strategy for wind turbine power maximization in region-2. The proposed algorithm is compared with a conventional ESC with log-of-power feedback (LP-ESC) scheme. The performance of both strategies is evaluated using the NREL 5MW reference turbine model in OpenFAST. Simulations were performed under different below-rated wind conditions. The major improvements of the LP-PIESC, over the conventional LP-ESC, are (i) faster convergence to optimal parameters, and (ii) smaller dither amplitude to estimate the gradient of log-of-power with respect to the control parameters.
These initial results provide evidence supporting the hypothesis that the amalgamation of log-power feedback and PIESC technique could result in faster and consistent convergence across different below-rated wind speed conditions. Design of the PIESC algorithm requires tuning more parameters than the conventional ESC. Additional work is needed to establish practical guidelines for parameter design customized to control of wind turbines and wind farms. In addition, evaluation of the LP-PIESC performance under wind conditions such as ramps or gusts should be performed to evaluate potential risks of this algorithm.