The Inﬂuence of Winglet Pitching on the Performance of a Model Wind Turbine: Aerodynamic Loads, Rotating Speed, and Wake Statistics

: The objective of this study is to investigate the inﬂuence of winglet pitching as an aero-brake on the performance of a model wind turbine by wind tunnel experiments. Time-resolved particle image velocimetry, force sensor, and datalogger were used to characterize the coupling between wake statistics, aerodynamic loads, and rotation speed. Results highlighted that, for a winglet with 4% of the rotor diameter length, the increase of its pitching angle can signiﬁcantly reduce the turbine rotation speed up to ∼ 28% and thrust coefﬁcient of ∼ 20%. The winglet pitching induced minor inﬂuence on the velocity deﬁcit in the very near wake regions, while its inﬂuence on accelerating the wake recovery become clear around three diameters downstream the turbine rotor. The turbulence kinetic energy exhibited a distinctive increase under large pitching angles in the near wake region at the turbine hub height due to the strong vertical ﬂow ﬂuctuations. Further investigation on the spectra of wake velocities revealed that the pitching of winglet can suppress the high-pass ﬁltering effects of turbines on wake ﬂuctuations; such large-scale turbulence facilitated the ﬂow mixing and accelerated the wake transport.


Introduction
Wind energy is one of the most abundantly available and fastest growing renewable energy sources in the US, which has increased more than six times during the last decade [1]. While wind turbines across various scales have been widely integrated in power grids, many engineering and scientific challenges still exist where the power generation cost is higher compared to traditional energy resources [2].
One of the most distinctive challenges for wind energy is the unsteadiness of incoming wind flows and therefore fluctuating power output and aerodynamic loads for wind turbines. Extensive studies [3][4][5][6][7][8] have pointed out that the atmospheric boundary layer varies continuously with season, time of the day, temperature, and other geophysical parameters. With wind tunnel and field measurement, recent works [9][10][11][12] have systematically investigated the influence of incoming flow fluctuations on the blade rotating speed and wake transport over both single turbine and wind farms. Specifically, under very high incoming velocities, the turbine rotating speed may exceed the rated value accompanied with dramatic aerodynamic loads. In such scenario, pitch control of turbine blades has been one of the mostly explored approaches for protecting the turbines under extreme wind conditions and optimize the wind turbine performance. By altering the effective angle of attacks and therefore the aerodynamic loads on blades, the pitching control allows for keeping constant turbine rotating speed and therefore power output at high wind speeds [13]. Nagai et al. [14] designed and constructed a prototype 3kW wind turbine generator with a variable pitch of the turbine blades. They found that it is not easy to maintain a constant rotational speed with wind fluctuation due to small moment of inertia of the turbine. They reported that the aerodynamic loads on the turbine can be significantly reduced by feathering the blades, which allows the structure to survive even under typhoon-like conditions. Muljadi and Butterfield [15] studied blade pitching in a variable speed wind turbine. They found that the capability of quick blade pitching is the key to maintaining constant rotor speed. To improve the pitch performance due to delays caused by hydraulic pressure driven units, Gao and Gao [16] developed a novel approach by integrating the delay-perturbation estimation and tested on turbines from kW to MW scales. In the last few decades, more attention has focused on the 'smart' blade pitching, where multiple control devices are equipped on a blade to adjust the local aerodynamic forces. Compared to the pitching of entire blade, 'smart' blade pitching allows faster responses to the incoming flow conditions and better suppresses the unsteady wind loads [17][18][19].
During recent decades, adding winglets on the blade tip has proven to be an effective approach to improve the turbine efficiency. With a small extension of the blade tip, winglets can effectively restrict the spanwise velocity, reduce the induced drag in the tip region, diminish wing tip vortices, and therefore lead to more favorable aerodynamics. Extensive efforts have focused on the role of winglets in improving the blade designs and increasing the power output. According to Whitcomb [20], the addition of a small winglet can reduce the induced drag loss more than double that reduced by increasing the wingspan with the same length. Johansen and Sørensen [21] studied different designs of winglets with the height of winglet being approximately 1.5% of rotor radius using computational fluid dynamics. They illustrated that, compared to upstream-facing winglets, the downstream-facing counterparts performed better with an increase of 1.71% in power coefficient and 1.81% in thrust coefficient at 10 m/s wind speed. They also reported that the power and thrust coefficients are highly dependent on tip-speed ratio and, for some designs, the winglets may reduce the power output at low wind speeds. Imamura et al. [22] numerically analyzed the horizontal axis wind turbine with winglets and blade tip extension. They showed that, compared to blade tip extension, rotors with winglets were more effective at decreasing the wake turbulence and increasing the power output. Gauna and Johansen [23] performed a theoretical study using an actuator cap with cylinders at the tips representing winglets. They demonstrated that the Lanchester-Betz-Joukowsky limit still holds for a wind turbine with winglets, and the increase in power production is due to the reduction in tip losses. Experimental work by Saravanan et al. [24,25] explored the effect of winglet at the tip of the blade on a wind turbine. They found that the pressure difference on the wind turbine blade increased due to the presence of the winglet. They also stated that pressure difference increases with longer winglet or smaller curvature radius of the winglet. Belferhat et al. [26] studied the flow over an isolated wing equipped with a small winglet. Their study involved different winglet cant angles of 0 ○ , 55 ○ , 65 ○ , and 75 ○ and reported that winglet cant angle of 55 ○ performed the best for increasing the power production. Wind tunnel experiments by Tobin et al. [27] found that, with a winglet of 6.7% rotor radius length, the power coefficient and thrust of the turbine can increase up to 8.2% and 15%. Zhu et al. [28] numerically investigated both upwind and downwind winglets on a horizontal axis wind turbine. They illustrated that, under tip-speed-ratio of 7 and blade pitch angle at 15 ○ , the upwind winglet produced higher power than the downwind counterpart. With both numerical and experimental efforts, Khaled et al. [29] reported that at a cant angle of 48.3 ○ and winglet length 6.32% of the rotor radius, the turbine showed the maximum power and thrust increase of 8.787%. Mühle et al. [30] did experimental investigation on a two bladed rotor with two interchangeable wing tips. The first wing tip was like a regular wind turbine with straight-cut wingtip while the second wing tip was a downstream-facing winglet. They reported that an optimized winglet design can improve the energy extraction as well as accelerate the recovery of turbulent kinetic energy in the wake by quicker tip vortex interaction. Mourad et al. [31] computationally explored the effects of winglet geometry on horizontal axis wind turbines. They found that downwind winglets with toe angles from 10 ○ to 30 ○ can reduce power outputs compared to rotors with no winglets.
Despite these progresses, the optimal control and designs of wind turbines to achieve the best performance under various incoming flow conditions remain a challenge. Specifically, there has been little investigations on the active/passive pitching of winglets, which may provide novel approaches for altering the aerodynamic characteristics of blades and protect turbines in extreme wind conditions. In this work, we explored the possibility of adjusting the performance of a wind turbine with winglets by changing its pitching angles. Wind tunnel experiments were designed to investigate how the variation of winglet pitching angles influences the aerodynamic loads, rotating speed, and wake characteristics of a modeled turbine. Results of this work will shed lights on reducing the thrust force and turbine rotating speed under strong wind, and estimating its influence on the performance of downstream turbines. The paper is organized as follows. The experimental setup is illustrated in Section 2; main results and analysis are detailed in Section 3, and final conclusions are summarized in Section 4.

Experimental Setup
Model wind turbines with winglets under various pitching angles were placed on the bottom wall of the Boundary Layer and Subsonic Tunnel (BLAST) in the University of Texas at Dallas. The test section of BLAST is 30 m long, 2.1 m high, and 2.8 m wide; during the experiment, incoming air flow temperature was controlled by a cooling system and maintained within 22.5 ± 0.3 ○ C to ensure constant incoming air properties. In order to well develop the turbulent boundary layer (TBL) for turbines, cubic blocks with 2.5 cm height located 0.2 m between each other were placed on the bottom surface of the test section ( Figure 1a). The rotation frequency of the fan of BLAST was fixed at 120 Hz; this leads to the mean velocity of incoming flow at hub height of the turbine as U hub = 6.43 m s −1 where the turbulence intensity is I u = σ u U hub = 12%, which is similar to those observed with full-scale turbines [32]. Here, σ u denotes the standard deviation of the streamwise velocity. The TBL reaches the freestream velocity at z ≈ 55 cm and therefore boundary layer thickness of δ z hub ≈ 2.75. Details for the non-dimensional profiles of the incoming velocity U U hub , streamwise turbulence intensity σ u U hub and kinematic shear stress −ū ′ v ′ U 2 hub within the TBL are presented in Figure 2. The horizontal-axis wind turbines were designed on a base model from Sandia National Labs [33,34] and manufactured from the Stratasys F370 3D printer in University of Texas at Dallas. All modeled turbines share the same rotor diameter of d T = 200 mm and the hub height z hub = 200 mm, leading to a Reynolds number of Re = U hub d T ν = 1.3 × 10 5 ; here, ν is the kinematic viscosity of air. The turbine tower was made from the M5 threaded rod. Details of the turbine blade geometry across various sections are provided in Tobin et al. [9,27] and not duplicated here briefly. The blockage ratio of the turbine was 0.53% based on the rotor sweeping area and wind tunnel cross section, which led to neglectable blockage effects. In this work, special attention is focused on the influence of winglet pitching on the wake and turbine performance. The winglets were designed by extending the blade's tip with l w = 8 mm (4% of rotor diameter) in a plane perpendicular to that of the rotor. Here, the geometry of the winglet is a rectangular plate with its width b w the same as the blade tip (Figure 1c,d) and resembles a wing-tip flap. It is worth pointing out that the winglet length tested here is slightly longer than those applied in other works [27,29]; this is limited by the resolution of force/flow measurements where longer winglet facilitates to highlight the influence of its pitching in altering the turbine performance and wakes. The pitching angle of the winglets α is defined by the angle between blade tip edge and the surface of winglets (Figure 1b). In this work, four winglet pitching angles with α = 0 ○ , 30 ○ , 60 ○ and 90 ○ were investigated. A 24 mm-diameter DC motor from RS PRO with a constant resistance through the experiments was used as the loading to control the rotating speed [11,27,35]. This led to the tip-speed-ratio (TSR) of λ = 2πωR U hub ≈ 4.2 for the base case (α = 0 ○ ), where R is the turbine radius. The measured power coefficient for the base case is C p ≈ 0.127. Both TSR and C p from current work are close to previous studies with similar rotor design [9,11,27]. The instantaneous turbine rotating speed was inferred from the output voltage of the DC motor measured directly via a USB-6210 datalogger from National Instruments (Austin, TX, USA). For each experiment, the voltage was sampled at a frequency of 10 kHz for a period of 100 s. Complementary measurements of time-averaged aerodynamic loads were performed with an ATI high-resolution force sensor connected to the turbine tower. The sensor was embedded within the bottom roughness to avoid any disturbs on the wake development. During each experiment, the instantaneous thrust forces were captured at a frequency of 1 kHz for periods of 60 s. The uncertainty of the sensor was less than 1.2% according to the minimum measured thrust. The incoming TBL and wake statistics along central axis of each turbine were characterized by a time-resolved particle image velocimetry (PIV) system from TSI. Two high-speed Phantom VEO440 cameras with 4 MP (2560 × 1600 pixels) resolutions were horizontally aligned to create a combined field of view (FOV 1) of 500 mm × 340 mm which includes ∼20 mm gap to investigate the mean wake flows. The FOV was located at the turbine wake region within streamwise distances x d T ∈ [1.85, 3.05], [3.15, 4.35] and heights z z hub ∈ [0. 3,2], where the coordinate system was defined with the origin coincident with the rotor plane at the bottom wall. For each experiment, 2500 image pairs were collected at a frequency of 30 Hz. The characterization of unsteady flow dynamics focuses on the wake region near the turbine hub height. Here, FOV 2 with the same streamwise regions but heights of z z hub ∈ [0.8, 1.2] were explored, where the sampling frequency was set at 500 Hz with a period of 50 s (i.e., 25,000 image pairs) to capture the instantaneous wake fluctuations. This sampling frequency is constrained by the storage of the PIV camera; as will be shown in the following analysis, this frequency allows for well characterizing both energy containing and inertial sub-range of turbulent flows. For both mean and unsteady wake measurements, the FOVs were illuminated by a 1 mm thick laser sheet generated from a 30 mJ/pulse laser, where air flow was seeded by 15 µm-diameter soap bubbles from the TSI bubble generator. The soap bubbles were well mixed by the recirculating wind tunnel before each experiment. The image pairs were processed using an Insight4G software package (Version: 11.1.1.0 TSI, Shoreview, MN, USA) from TSI with a multipass scheme. The final interrogation window size was 32 × 32 pixels with 50% overlap, resulting in a final vector grid spacing ∆x = ∆y = 2.4 mm. The overall uncertainty of the identified seeding particle locations was ∼0.1 pixel [36]; this led to the uncertainty of flow velocity measurement of ∼1.4% given the bulk particle displacement of 7 pixels between two successive images.

Results and Discussion
In this section, we discuss in detail the distinctive influences of blade winglet pitch angles on the statistics of turbine rotation speed, thrust coefficients, and wake development of the model wind turbines.

Rotation Speed and Thrust Force
Specific insight on the effect of winglet pitch angles on the turbine rotation speed can be obtained from the mean and associated fluctuation statistics in each case. Figure 3 shows a comparison between time-averaged tip-speed-ratio λ and thrust coefficient C T under all different α normalized by that of the base case (α = 0 ○ ), λ o = 4.2, and C T o = 0.585. Here, the thrust coefficient is defined as where T is the mean thrust force and ρ a is the air density. Generally, with the growth of winglet pitching angle, λ presented distinctive decrease up to ∼28% under α = 60 ○ and 90 ○ . The evolution of C T shows a similar trend where the maximum decrease reached ∼20%. The power coefficients of turbines with α = 30 ○ , 60 ○ and 90 ○ reduced to 0.088, 0.067, and 0.066, respectively. As a dominant impact, with the growth of winglet pitching, its projected area impinged by the induced velocity due to turbine rotation increased monotonously; this led to higher local drag on the blade tip region and slowed down the turbine rotations. Interestingly, continuous growth of α from 60 ○ to 90 ○ induced a minor influence on both λ and C T . This can be explained from two aspects. First, as pointed out in Modi et al. [37], the difference of drag coefficient for a flat plat inclined at 60 ○ and 90 ○ is small. In addition, the winglet pitching can alter the local flow characteristics with its schematics shown in Figure 3b. Here, with the growth of α, the winglet can no longer block the tip vortices (red arrows) well, which reduces the blade efficiency. At the same time, due to turbine rotation, additional pressure difference is produced between the windward side (surface A) and leeward side (surface B) of the winglet; such pressure difference can trigger secondary flow motions (blue arrows) which enhances the tip vortices and further slows down the turbine rotation. Note that, with α → 90 ○ , the secondary flow is nearly perpendicular to the blade edges and therefore no longer facilitates the tip vortices. It is worth pointing out that the mass of the winglet portion is only ∼1.8% of the entire blade; therefore, compared to traditional pitching of the whole blade, the winglet pitching provides potential with a much faster response of turbines under instant variations of wind velocities. Specifically, with the same generator resistance load, our work demonstrated that the turbine rotation speed can remain constant under a ∼40% increase of incoming wind velocity by adjusting the winglet pitching angle from 0 ○ to 60 ○ . This provides new approaches for efficiently protecting the turbines under extreme wind conditions. The influence of winglet pitching angles on turbine rotation speed fluctuations is illustrated in Figure 4; here, the power spectra (Φ ω ) of instantaneous rotation angular velocity ω are analyzed for selected cases, and the spectrum of incoming turbulence (Φ u ) at turbine hub height is added in the subplot for discussion. In general, the spectra under 0 ○ and 90 ○ presented a similar trend including two distinctive regions. In the first region with f d T U hub < 0.06 (R 1 ), the turbulence is within the range of energy containment, where the energy cascade decay of Φ ω is only related to the aerodynamic characteristics of the rotor [9,38]. Under higher frequency regions (R 2 ), the turbulence is within the inertial sub-range following Φ u ∝ f −5 3 which induced a faster decay rate of Φ ω . For both winglet pitching angles, Φ ω plateaus at f d T U hub ∼ 1, indicating that the turbine rotation was no longer affected by the small scale turbulence. It is worth pointing out that both spectra presented minor differences in R 2 , while, in general , Φ ω of α = 90 ○ showed lower energy under R 1 . This reveals that the growth of winglet pitching angle reduces low-frequency turbine rotation speed fluctuations.
To further highlight the influence of winglet pitching angle on the rotation speed fluctuations in low-frequency regions, the integrated energy was calculated as shown in Figure 5. Overall, small winglet pitching angles (α = 0 ○ and 30 ○ ) led to larger σ 2 ω across all frequency ranges, corresponding to the higher Φ ω found in Figure 4. Note that, for all investigated α, the curves of σ 2 ω plateaus at f d T U hub ≈ 0.06, i.e., the beginning of R 2 , indicating that turbine rotation fluctuations are mostly modified by low-frequency components.

Wake Characteristics
The influence of winglet pitching on the wake development is first assessed by comparing the normalized mean streamwise velocity distributions U U hub shown in Figure 6. In general, under higher α, the wake velocity deficit is reduced especially at the hub height region. Such phenomena are more clearly illustrated with the selected velocity deficit profiles ∆U = U inc (z) − U(z) in Figure 7, where U inc is the time-averaged incoming flow speed. It proves that, in the very near wake x d T ≈ 2, the wake recovering is not yet affected by the winglet pitching angles, while distinctive differences of wake velocities occur at x d T ≈ 3 and gradually diminish in the far wake. It is worth pointing out that the induced velocity deficit is highly coupled with the thrust coefficient of the turbines; to better analyze the relation between aerodynamic loads and wake recovering, the Bastankhah and Porté-Agel's [39] analytical wake model is included to compare with our measurements. This model applied mass and momentum conservation to a control volume around a turbine with an assumption of self-similarity for the wake velocity profiles as: where y and z are spanwise and vertical coordinates, and A represents the wake development in streamwise direction as: where k * is the wake growth rate, and β is a function of C T , given as: In this work, specific attention is focused on the wake development in the turbine hub height where the influence of winglet pitching is most distinctive. Therefore, Equation (3) can be simplified with z = z hub and y = 0:  The comparison between measured (solid line) and modeled (dashed line) hub-height wake velocities are summarized in Figure 8. Key to determine the wake recovery rate is estimating the value of k * . Previous works have shown that k * is a complex function of incoming turbulence and turbine designs [40,41]. In this work, k * is determined via complementary experiments of a turbine without winglet, where the wake flows and thrust forces were measured; this information is presented in Figure 8b, which shows that k * = 0.0324 is the best fit for the wake development. Note that C T = 0.5744 for the turbine without winglet and is slightly lower than the base case, which follows the same trend as reported in Tobin et al. [27]. In addition, Bastankhah and Porté-Agel [39] pointed out that the velocity deficit profiles present a self-similar shape after a certain downstream distance (x d T ≳ 3); therefore, the model is compared with experimental results only at the second FOV (x d T ≥ 3.15). In general, the model well presented the wake recovery trend under the influence of winglet pitching. Based on the model prediction, the winglet pitching from α = 0 ○ to 60 ○ produces ∼6% higher hub-height wake velocity at x d T = 5 and ∼3.5% at x d T = 7; with a first order approximation of P ∼ U 3 , this benefits the power production of downstream turbines with 19% and 11% located at x d T = 5 and 7, respectively. It is worth pointing out that, for all investigated configurations, especially under high α, this model led to a somewhat slower wake recovery rate dU dx prediction compared to the experimental results (Table 1). This indicates that, compared to traditional turbines sharing the same thrust coefficient, the pitching of winglet brings additional influence on wake recovery and induced stronger wake mixing to accelerate flow recovery [42,43].  The quantification of wake fluctuations provides insights for estimating the unsteadiness of power outputs of downstream turbines; this information is characterized by the turbulence kinetic energy TKE = ⟨u ′2 + w ′2 ⟩ 2U 2 hub shown in Figure 9; here, u ′2 and w ′2 are the velocity fluctuations in streamwise and spanwise directions and ⟨.⟩ denotes the time-average. Due to the limitation of two-dimensional PIV measurement, the velocity fluctuations in the y-axis component is not included [35,44]. This quantity illustrates the strong fluctuations of wake velocity due to the influence of α. Overall, distinctive wake fluctuations occur near the top-tip regions of the turbine where the mean velocity shear is significant; such phenomenon is similar to those observed in traditional turbine wakes from both numerical simulations and experiments [45,46]. It is worth pointing out that high α induces a distinctive increase of TKE in the hub-height region, especially within the near wake; this is more clearly shown by the TKE distribution profiles in Figure 10, where the difference of TKE across various α reaches the maximum at x d T = 2, z z hub ≈ 1.2 and gradually diminishes along the streamwise direction. To further investigate the production of strong TKE under high α, wake fluctuations normalized by hub-height incoming velocity at streamwise (I u ) and vertical (I v ) directions are analyzed as shown in Figures 11 and 12 for selected cases. While the variation of α induced little influence on the distribution of I u (Figure 11), higher α produced significantly stronger vertical flow fluctuations especially near the hub-height region within x d T ≲ 3 ( Figure 12); this contributed to the overall stronger TKE shown in Figure 9.  Insights on the influence of α for the structure of turbulence may be characterized by the difference between pre-multiplied streamwise velocity spectra in the wake and incoming flow at hub height along the stream-wise direction shown in Figure 13, which is given as where f is the frequency and Φ is the corresponding energy spectra of flow velocity. Overall, across all α, the turbines increased turbulent motions within f d T U ≳ 10 −1 and suppressed wake fluctuations with larger scales; similar 'high-pass filtering' effects of turbines on the wake flows were observed in previous experimental work [47]. Interestingly, for the very near wake flow (x d T ≲ 3), such suppression is most distinctive under α = 0 ○ but gradually diminishes with the growth of α. As discussed in Figure 3, the winglet pitching significantly reduces the tip-speed-ratio and therefore the blades became less 'efficient' in 'high-pass filtering' the wake flows, especially in the very near wake regions.  Finally, we discuss the wake turbulence transport under different winglet pitching; this helps to analyze the correlation of power output fluctuations between upstream and downstream turbines [11]. This quantity is analyzed by the cross-correlations on the time-series of streamwise velocities at two locations: Equation (8) describes the coherent motions between flow velocities at x = 2d T and downstream locations x = 2d T + δx along the hub-height. Selected cases for the evolution of η are shown in Figure 14a. It is clear that, with a specific time delay (∆t), the value of η reaches the local maximum as highlighted by the red dashed lines, where the footprint of upwind flow fluctuations is transported to the downstream locations. As expected, with longer distance (i.e., higher δx), the corresponding transportation time increased where the strength of correlation (i.e., maximum value of η) decreased.
To further investigate the influence of α on the wake transport speed, the distribution of ∆t normalized by δx and U hub (i.e., the ratio between averaged wake transport speed within the distance δx and the incoming velocity) across various α and δx is summarized in Figure 14b. Interestingly, the wake transport speed is in general higher than the time-averaged wake velocity at the turbine hub height (U U hub ) from Figure 8), especially under high α cases. This further highlighted that the pitching of winglet facilitates the wake mixing, where the wake transport is not barely determined by the local wind velocity but also influenced by the mixing with the surrounding background TBL, which is always higher than the local wake flow velocities [48].

Conclusions
In this study, the influence of winglet pitching as an aero-brake on the performance and wake characteristics of a modeled turbine was investigated with wind tunnel experiments. Despite a small part of the entire blade (4% of turbine diameter length), the pitching of winglet can significantly alter the turbine rotation speed, thrust coefficients, wake recovery rate, and turbulence transport. Key concluding points in this work are summarized as follows: (i) The pitching of winglets can effectively decrease the turbine rotation speed and thrust coefficient up to ∼28% and ∼20% as the pitching angle is increased to 60 ○ . The unsteady turbine rotations are dominated by the low-frequency fluctuations, which can be suppressed under large winglet pitching angles.
(ii) The pitching of winglet leads to smaller velocity deficit which is induced by both the variation of thrust coefficient and blade aerodynamic features. The growth of winglet pitching angle significantly increases the vertical velocity fluctuations; this is clearly reflected by the distinctively higher turbulence kinetic energy at hub height of the near wake regions.
(iii) Large winglet pitching angles reduce the high-pass filtering effects of wake fluctuations in the near wake. While the mean wake velocity presents distinctive variation at the hub-height along the streamwise direction, the speed of wake fluctuation transport under large winglet pitching angles demonstrates much smaller changes.
To sum up, this work provided evidence that, for turbines with winglets, altering their pitching angles can effectively adjust the turbine performance, influencing the wake flows and therefore downstream turbine operations. This concept is similar to the control of flaps which changes the local aerodynamic characteristics of the blades. Future work will explore in detail the influence of winglet lengths and shapes, as well as the possibility of applying different winglet pitching angles on each blade to balance the induced torque acting on the hub due to nonuniform incoming flows.