Abstract
Wind farm control has demonstrated power production improvements using yaw-based wake steering compared to individual turbine optimization. However, slower yaw actuation rates in response to rapid inflow changes lend to impracticality of yaw-based steering, as it causes time-varying downstream rotor–wake overlap, power production fluctuation, and consequent reduction. Therefore, closed-loop wake control is required to mitigate wake deflection uncertainty. To respond to rapid inflow variations, rotor speed actuated is investigated here. Furthermore, wake position information is required as feedback for closed-loop control function. For field-installed turbines, nacelle-based Light detection and ranging (LIDAR) is expected to provide this information. So far, LIDAR-derived wake position has been determined through model-based field reconstruction of scattered LIDAR data. However, this requires sophisticated, economically prohibitive LIDARs. To incorporate inexpensive, two-beam LIDAR for wake detection, a tip vortex-based approach was developed and is also presented here. These contributions can be considered as intermediate steps toward realization of a novel closed-loop wake control.
Introduction
To make wind power competitive with other conventional energy sources, it is important to minimize the cost of wind energy, thus maximize the power production. Grouping wind turbines in large multirow facilities referred to as wind farms has contributed in achieving this goal given the benefits of economies of scale. However, wind turbines within a wind farm are often operating in the wake of another turbine, leading to a reduction in power production and an increase in structural loads. To minimize wake effects, wind turbines are typically spaced between 7 and 10 rotor diameters D. However, full-scale wake measurements have revealed that the wake influence can extend in excess of 20–30D (Hirth and Schroeder, 2013; Hirth et al., 2012). On the contrary, the current wind farm control industry practice relies on the greedy control approach to optimize the power production of individual wind turbines (e.g. Abdullah et al., 2012; Koutroulis and Kalaitzakis, 2006; Kumar and Chatterjee, 2016) without considering the impact of wake effects on neighbor turbines. This approach is suboptimal for the wind farm performance as wind farm wakes account for up to 20% power loss (compared to the case where all the turbines are exposed to freestream inflow
Furthermore, yaw-based wake steering optimization procedures have been proposed with the goal of finding a set of yaw angles among multiple wind turbines that produce a certain wake configuration within the wind farm leading to maximum power production (Gebraad et al., 2016). It is important to note that wake models used to describe the flow for wind farm optimization are static in nature (i.e. assuming steady-state conditions), which overestimate the expected power production enhancement. Although dynamic model considerations are available in recent literature (Munters et al., 2016), simplified models are preferred by the industry, since in their experience, the cost-to-benefit ratio of higher-fidelity models is yet to be realized (Ulker, 2018). In real operating conditions, the wake propagation is dynamically disturbed by the continuously changing wind direction. These wake deflection dynamics might lead to continuous overlapping with downstream rotors, increasing structural loads, causing fluctuation in overall power production, and a consequent reduction in average power production compared to the predicted power output from the static models. In this regard, closed-loop wake steering control can be an effective approach to guide the wake to a steady position in response to changing wind conditions, as illustrated in Figure 1.

Wake trajectory stabilization. Solid lines: desired wake trajectory; dotted lines: instantaneous disturbed wake trajectory.
One aspect of the closed-loop wake steering control is the control action to actuate the wake. Currently, the wake deflection strategies are driven by active turbine yaw. However, the efficacy of yaw-based wake steering to control the wake position is reduced at small timescales, as the wind variations are too rapid to follow by the low-acting yaw actuator (Vollmer et al., 2016). Studies in the field (Bromm et al., 2018) and wind tunnel (Castillo et al., 2017) have demonstrated that yaw-based wake steering is more prominent at larger downstream distances (>7D). In addition, parameterizations, such as Guntur et al. (2012), show that wake deflection significantly increases with downstream distances. A small change in mean wind direction can produce a significant variation in wake conditions of downstream turbines, which leads to large changes in overall wind farm power output (Porté-Agel et al., 2013). This emphasizes consideration of large-scale mean-flow transients (mesoscale, 10–100 km or 10–100 min, weather changes) to accurately predict wind farm power for typical meteorological conditions (Munters et al., 2016). Furthermore, it is critical that the wake steering control action has a time response that is less than the timescales at which the mesoscale transients are observed to be effective. In addition, the yaw mechanism actuation, although relatively quick (0.5°/s), must rest at the new position for a certain time period. References state a resting period of 10 min (Mikkelsen et al., 2013; Thorson, 2013) and consequent yaw activation only upon meeting predetermined yaw error threshold (Kim and Dalhoff, 2014). This slow movement and frequent pause are required to minimize the effect of destabilizing gyroscopic forces (Kim and Dalhoff, 2014). Therefore, although inter-turbine wake or flow propagation timescale is within the yaw system’s time response, a minor yaw-misalignment local error will likely lead to a high downstream wake deflection. This will reduce the efficacy of the yaw-based wake steering, since the local wind direction is not accurately predicted by upstream measurements. Therefore, to dynamically steer the wake in response to wind variations, the effect of a much faster control action such as rotor speed on the wake deflection is discussed in this article.
In addition, the ability to sense the wake position that could be used to provide the feedback required to adjust the wake trajectory is considered here. Wake position has been characterized in terms of the location of the wake center in several wind tunnel studies using particle image velocimetry (PIV) (Howland et al., 2016). For field-installed turbines, Raach et al. (2016b) proposed local tracking of the wake center trajectory using nacelle-mounted Light detection and ranging (LIDAR) looking downstream. Since the wake center is a function of the geometry of the wake profile (Howland et al., 2016; Vollmer et al., 2016), it cannot be derived directly from single LIDAR line-of-sight measurements. In the wake center tracking approach described in the study by Raach et al. (2016a), the wind velocity field and hence the location of the wake center were estimated through the model-based wind field reconstruction approach. This wake detection approach was tested with high-fidelity simulation SOWFA (Simulator for Off/Onshore Wind Farm Applications), but validation at full scale has not been reported yet. Full-scale wake measurements, where the wake characteristics were derived through the field reconstruction approach, have been reported by Fleming et al. (2016) and Herges et al. (2017) using a ground-based modified LIDAR mounted on the nacelle measuring downstream. However, such LIDAR requires the ability to focus at different ranges and perform complex scanning patterns, which represent a cost issue limiting its commercial application. Therefore, to enable wake detection with simple, low-cost, off-the-shelf nacelle-mounted LIDAR, a wake detection strategy based on the location of the tip vortices is presented in this article. However, in this wind tunnel-based study, explained in section “Methodology,” downwind measuring hot-wire anemometers that emulate nacelle-mounted LIDAR function are utilized for the LIDAR-based wake detection construct validation.
Methodology
Given the inherent unsteady nature of the turbulent inflow, dynamic changes in rotor speed, yaw, and pitch angle should be considered (Westergaard, 2012). High-fidelity simulation tools have been used to study open-loop wake control, requiring substantial computational resources and time (Annoni et al., 2014; Munters et al., 2016). However, it will require significant efforts to accommodate and study unsteady inflow combined with required dynamic controls, substantially increasing computing resources and time. On the contrary, exploring innovative strategies in full-scale turbines is impractical and expensive and may not be available for iterative development of wake control strategies. Therefore, the wind tunnel-based “Hyper Accelerated Wind Farm Kinematic-Control Simulator” (HAWKS), deployed in the study by Castillo et al. (2016, 2017) to characterize wake flow and test wake control strategies, was adapted to test the wake actuation and the wake detection strategies separately.
Wake deflection actuation
Several studies (Guntur et al., 2012; Jiménez et al., 2010; Park et al., 2013) have pointed out that the degree of wake deflection is proportional to the thrust coefficient CT changes. When a wind turbine operates in yawed conditions, the thrust force T exerted by the wind on the turbine rotor can be divided in two components

Simplified schematic of the forces exerted on the wake of a yawed turbine.
Since the thrust force component
Wake detection
LIDAR technology proposed in the studies by Raach et al. (2016a, 2016b) has been adapted to focus at different ranges and perform complex scanning patterns. However, for practical implementation, wake detection approach development based on a simple, inexpensive fixed focus, two-beam LIDAR is considered here. The continuous-wave, two-beam nacelle-mounted LIDAR model WindEye, developed by Windar Photonics A/S, is used as the primary reference for the development of the wake detection methodology. This LIDAR is feasible due to its lower power laser, smaller optics, and no moving mechanical parts. The LIDAR scanning pattern geometry is shown in Figure 3, where the angle between the two beams is 13° and the focus distance is 80 m.

Instantaneous vorticity contours under dynamic yaw misalignment illustrating the wake tip vortex movement relative to fixed LIDAR focus point (purple dot).
A characteristic of the near wake is the presence of the tip vortices shed by the blade tips, which rotate and propagate downstream, defining the boundary between the wake and the outer flow, in other words, defining the wake envelope. The wake deflection in the near wake area might be defined by the tip vortices as reported by Haans et al. (2007).
Therefore, as the yaw angle offset changes due to the inherent variability in the inflow, the wake interface defined by the tip vortices moves in and out of the boundary set by the LIDAR focus points as illustrated in Figure 3. The location of the tip vortices was estimated from turbulent statistics based on instantaneous velocity field measurements (e.g. Castillo et al., 2017; Lynum, 2013). These studies are wind tunnel-based using both PIV and hot-wire anemometry, respectively. On the contrary, in a wind tunnel study reported in the study by Van Dooren et al. (2017), hot-wire anemometry showed good agreement with LIDAR measurements, implying its potential to emulate LIDAR wake detection function in wind tunnel using hot-wire anemometry. Therefore, to estimate the position of the tip vortices and hence the wake envelope relative to the LIDAR focus points, hot-wire anemometry will be incorporated in the HAWKS platform to define the instantaneous wake position.
Experimental setup
The HAWKS platform was developed in the National Wind Institute (NWI) closed-loop wind tunnel at Texas Tech University. The test section has a cross section of 1.2 m × 1.8 m. The mean stream-wise incoming velocity (free stream) was kept constant at
The velocity field across a planar area of the flow field downstream of the turbine was measured with the PIV technique, and its setup is shown in Figure 4(a) to (d). Since PIV camera placement was feasible along the wind tunnel sidewall and not the bottom, to obtain wake cross-section measurements in a plane that is horizontal in the field, the turbine tower was mounted to the tunnel sidewall and not its floor. Figure 4(a) shows the 90° rotated wind turbine, where the tower appears to be a cantilever beam, and Figure 4(a) and (b) shows the laser light sheet location that represents the measurement plane at hub height. A 150 mJ/pulse Nd:YAG laser (532 nm) coupled with laser sheet forming optics was used to illuminate the flow seeded by a particle generator using propylene glycol oil droplets. The particle images were recorded with an array of four 14-bit LaVision charge-coupled device (CCD) cameras (Figure 4(c)) with a resolution of 1600 × 1200 pixels each using 50 mm focal length lens that were synchronized with the laser pulses. The time interval dt between two laser pulses was set to 300 μs. However, the HAWKS PIV system characteristics described in the study by Castillo et al. (2017) were modified in order to accommodate a larger field of view. The CCD cameras were 2.75 m away from the light sheet plane, so that the image yielded a total inspection area of about 2.24D × 7.5D (Figure 4(b)), with an image resolution of 2.52 pixels/mm. The image processing was carried out with Davis software (v8.2.2) from LaVision. The size of the interrogation window in the correlation calculation was changed to 16 × 16 pixels with an overlap of 50%. A summary of the HAWKS platform characteristics is listed in Table 1.

HAWKS PIV setup: (a) model wind turbine when observed from upstream, (b) wind tunnel test section, (c) schematic of the PIV system, and (d) four-camera array.
Summary of HAWKS platform characteristics.
HAWKS: Hyper Accelerated Wind Farm Kinematic-Control Simulator; PIV: particle image velocimetry.
To describe the wake flow turbulent behavior in detail, a Dantec Dynamics 54T42 MiniCTA (constant temperature anemometer) equipped with 55P16 hot-wire probe was utilized to emulate the LIDAR function. The data acquisition is carried out with a 16-bit NI 9215 DAQ. The hot-wire probes were mounted on the traversing system (Figure 5) to move the probe around the wake flow field.

Hot-wire traverse system located downstream the wind turbine.
Results
This section outlines the results of the tests performed in the HAWKS platform. Before a closed-loop algorithm is tested in HAWKS, its components that form an open loop were tested. The closed-loop control requires an input of wake location and an ability to respond to the input to actuate wake deflection of desired location. Each of these aspects, wake detection, and actuation were tested independently. This section discusses the experiments undertaken to test these for the HAWKS platform.
Wake deflection actuation
To illustrate the deflection of wakes of the yawed turbine with variable rotor speed, measurements were made for static yaw offset

Contours of the normalized mean stream-wise velocity
Figure 6 shows that the wake is more pronounced as the rotor speed

Comparison of normalized mean stream-wise velocity
Profiles shown in Figure 7 reveal wake shift as ω increases due to an increasing thrust force. To further illustrate this, the location of the wake center

Comparison of measured wake deflections
In the three cases, Figure 8 shows that the empirical model applied to cases CT = 0.55, 0.66, and 0.85 works reasonably well with the measured wake deflection at
Detection task
As the wake deflection changes, the fixed LIDAR focus point sweeps across the wake interface. The power spectral density (PSD) of the measured stream-wise velocity u at the focus point will vary depending on the relative wake interface location. The spectra were calculated with discrete Fourier transform (DFT). At each measurement point, the total number of samples was divided into 10 successive windows, and the DFT of each successive window was then calculated and finally averaged to find the mean spectrum Euu(f). The resulting spectrum was normalized by the free-stream velocity variance σu2. As a starting point, only zero-yawed case is analyzed here.
The tip speed ratio of the model turbine was set to 3.91, corresponding to 3000 r/min. The hot-wire probes acquired data at 20 kHz, which provides 400 samples per revolution. The sampling time of each window was 0.1024 s, giving a total of 2048 samples, corresponding to about 5 revolutions. First, the hot-wire probes were traversed across the wake interface very close to the rotor at x/D = 0.2 on one side of the wake to perform measurements at three locations. Figure 9(a) shows a comparison of the spectral characteristics at each location. The peaks in Figure 9(a) demonstrate that at every position across the wake interface, the dominant frequencies are multiples of the rotor rotational frequency, 3000 r/min or 50 Hz. The once-per revolution occurrence is denoted as 1 P. The second dominant peak is close to 100 Hz, meaning two occurrences per revolution and is denoted as 2 P. The third peak is close to 150 Hz, meaning three occurrences per revolution and is denoted as 3 P. The peak occurrences shown in the spectrum represent tip vortices.

Comparison of spectral characteristics across the wake interface at x/D = 0.2: (a) PSD and (b) 1, 2, and 3 P spectral energy.
To better interpret the peaks resulting from the PSD shown in Figure 9(a), Figure 9(b) shows the spectral energy at 1, 2, and 3 P frequencies of the different points across the wake interface at x = 0.2D. At y/D = 0.45 (the closest to the rotor center), the peaks are barely visible and hence there is no evidence of the presence of the tip vortices at this location. At y/D = 0.5, the three peaks are clearly present. At y/D = 0.55, the only peak at 3 P is an indication that the three tip vortices are clearly present at this location. When comparing the spectral energy at 3 P between y/D = 0.5 and 0.55, there is almost no change.
However, the peaks at 1 and 2 P increase the most at y/D = 0.5 compared to y/D = 0.55 and y/D = 0.44, indicating the presence of wake interface. Furthermore, the peak at 1 P is an indication that the tip vortices have merged. The peak at 2 P is also prominent, perhaps indicating an intermediate state during the vortex merging at this measurement location, which is close to the rotor plane (x = 0.2D).
The hot-wire probes were traversed across the wake interface at x/D = 3, the plane where the nacelle-mounted LIDAR would be focused, to perform measurements at six locations. Figure 10(a) shows a comparison of the spectral characteristics at each location. Figure 10(b) clearly reveals a slight difference among the spectral energy across the wake interface at both 2 and 3 P. On the contrary, the peak at 1 P, which is the most evident peak in every position, shows a significant variation across the wake interface reaching a maximum at blade tip y/D = 0.5. The fact that 1 P peak dominates at x/D = 3 at a wider range might be an indication that the tip vortices have completely merged and further diffused.

Comparison of spectral characteristics across the wake interface at x/D = 3: (a) PSD and (b) 1, 2, and 3 P spectral energy.
Figure 11 illustrates the effect of changing the rotor speed on the spectrum characteristics at x/D = 3. Figure 11 shows that the dominant peak in all the spectrums is at 1 P. It also reveals two features when the rotor speed increases: (1) the 1 P peaks occur at higher frequencies (shifted to the right) and (2) the 1 P peak spectral energy increases. The former is expected as the rotational frequency of the rotor is higher. The tip vortex circulation is proportional to the turbine blade lift force described by the Kutta–Joukowski law (Lynum, 2013). Therefore, the latter is caused by an increase in the lift force on account of the angle of attack increase due to rotor speed increase.

Effect of tip speed ratio on the spectral characteristics at x/D = 3 and y/D = 0.5.
Subsequently, the PSD magnitude for 1 P will be used to calibrate the location of the wake interface relative to the hot-wire location. This information will be used as an input to the actuation system. The rotation speed will be adjusted such that the wake is moved to a desired location. These experiments will be reported in future studies.
Conclusion and future work
In this article, experimental validation rotor speed-based wake steering and position detection in a controlled wind turbine wake simulator was presented. First, the influence of the rotor speed
Footnotes
Acknowledgements
The experiments were performed at the TTU’s National Wind Institute (NWI) closed-loop wind tunnel at Texas Tech University.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was funded by the Governor’s Office in the State of Texas Emerging Technology Funds (ETF), Global Laboratory of Energy Assets Management and Manufacturing (GLEAMM), and the Binational Industrial Research and Development (BIRD) Foundation.
