## 5 Results from BEM Computations
At the NTM7.5 wind speed realization, the largest increment average power output was reached by the case with the MTs implemented in the station 8 of the blade, with an increase of the power in compared with the clean wind turbine of 4.425%. At the NTM10 wind speed realization the largest average power value is reached by the case with the MTs mounted on st8 with an increase with respect to the clean case of 3.282%. The largest increases for the three wind speed realization used in the present work have been performed by the case with the MTs in the blade station 8, corresponding to the case DU91W(2)250 MT9520 st8. Figure 14 illustrates in a bar plot the increments in the power output for every case in comparison with the clean case. The effect of mounting MTs on the blade stations studied in the present work is more significant at low wind speeds than at wind speed close to the nominal power. At those low wind speeds, the MTs can help to increase notably the wind turbine power output performance. Note that even thought the case DU91W(2)250 MT9520 st8 st9 has more MTs mounted along the blade, its power performance in terms of average power output is quite similar in comparison with the other cases with MTs on stations 8 and 9. After applying the BEM algorithm to all control volumes, the tangential and normal load distribution is known and global parameters such as thrust and bending moment at the root of the blade can be computed. In the current work, the mean values of both thrust and bending moment have been calculated for each wind speed realization. The thrust has been calculated by Equation (5) taking into account the trust distribution along the blade and the bending moment has been determined by Equation (6). \[T=\int_{0}^{R}4\cdot\!\tau\cdot\!r_{i}\!\cdot\!\rho\!\cdot\!(V_{0})^{2}\cdot\! \omega\!\cdot\!a_{i^{\prime}}(1-a_{i})\cdot\!F\cdot\!dr_{i} \tag{5}\]  
\[M=\int_{0}^{R}4\cdot\!\tau\cdot\!r_{i^{\prime}}^{3}\cdot\!\rho\!\cdot\!V_{0} \cdot\!\omega\!\cdot\!(1-a_{i})\cdot\!a_{i^{\prime}}^{4}\cdot\!F\cdot\!dr_{i} \tag{6}\]  
Figure 14: Average wind turbine power output increases for all wind speed realizations of the microtabed wind turbine in comparison with the clean one. The Prandtl's tip loss correction factor \(F\) has been determined by the Equation (7), for both thrush and bending moment calculations:  
\[F=\frac{2}{\pi}\cos^{-1}\left(e^{-f}\right)\text{where}f=\frac{B}{2}\frac{R-r_{ i}}{r_{i}\sin\Phi} \tag{7}\]  
The variables used for thrust and bending moment estimations and the corresponding dimensions are shown in Table 5. All calculations are based on the 5 MW reference wind turbine described in [22]. Table 6 represents the mean values of thrust calculations for all MT cases according to the matrix described in Table 2. The thrust was calculated by integrating the thrust distribution along the blade by the Equation (5) for each wind speed realization. Afterwards, the average thrust was determined according to the wind realization duration of Figure 11. The increments in the average thrust of the microtabel blade cases have been illustrated in Table 6 in comparison with the clean case. The average thrust has been experienced an increase in every case with MTs implemented and once again, the case with the MTs mounted on the blade station 8 presents larger increments than the other cases, which is in concordance with the results of power output presented in Figure 14.  
The bending moment in the root of the blade has been determined by Equation (6) taking into account the bending moment distribution along the blade. The bending moment in the root of the blade was determined for each wind speed realization and computed along blade according to Equation (6). Afterwards, the mean bending moment was determined according to the wind realization duration of Figure 13:0 The increments in the average bending moment of the microtabel blade cases have been illustrated in Table 7 in comparison with the clean case. No extraordinary increment in the mean