Experimental investigation of wind turbine wake and load dynamics during yaw maneuvers
- 1Univ. Orléans, INSA-CVL, PRISME EA4229, 45072 Orléans, France
- 2Ecole Centrale de Nantes, LHEEA, 1 rue de la Noë, 44321 Nantes, France
- 3Centre de Recherche de l’École de l’Air, B.A. 701, 13661 Salon-de-Provence, France
- 4Engie Digital, Le Monolithe, 59 rue Denuzière, 69285 Lyon, France
Correspondence: Sandrine Aubrun (email@example.com)
This article investigates the far wake response of a yawing upstream wind turbine and its impact on the global load variation in a downstream wind turbine. In order to represent misalignment and realignment scenarios, the upstream wind turbine was subjected to positive and negative yaw maneuvers. Yaw maneuvers could be used to voluntarily misalign wind turbines when wake steering control is targeted. The aim of this wind farm control strategy is to optimize the overall production of the wind farm and possibly its lifetime, by mitigating wake interactions. While wake flow and wind turbine load modifications during yaw maneuvers are usually described by quasi-static approaches, the present study aims at quantifying the main transient characteristics of these phenomena. Wind tunnel experiments were conducted in three different configurations, varying both scaling and flow conditions, in which the yaw maneuver was reproduced in a homogeneous turbulent flow at two different scales and in a more realistic flow such as a modeled atmospheric boundary layer. The effects of yaw control on the wake deviation were investigated by the use of stereo particle imaging velocimetry while the load variation on a downstream wind turbine was measured through an unsteady aerodynamic load balance. While overall results show a nondependence of the wake and load dynamics on the flow conditions and Reynolds scales, they highlight an influence of the yaw maneuver direction on their temporal dynamics.
The rising market demand for wind energy, together with the need to reduce costs and maximize power yield, has led to an increase in wind farm density (Pao and Johnson, 2009) with a concomitant increase in wake interactions. Wake interactions are the main cause of power losses and increase in fatigue loads in wind farms (Sanderse, 2009). In order to attenuate these negative effects of wind farm densification, different wake control strategies started to be envisaged. A wind farm control strategy consists in controlling each wind turbine individually in order to reduce its wake effects on the nearest downstream turbines and maximize the total power production. The most common solutions investigated are twofold: induction control and yaw control. The former is based on a power curtailment strategy, as reducing the power extraction of an upstream wind turbine leaves more kinetic energy available for a downstream one. The latter consists in wake steering: a wind turbine (henceforth WT) is voluntarily misaligned with respect to the wind direction in order to deviate its wake from its nominal position and hence reduce wake effects on a downstream wind turbine. If properly applied, wind farm control can improve the overall production. Its potential has already been investigated and confirmed by both simulations (e.g., Bossanyi and Jorge, 2016; Gebraad et al., 2017) and full-scale tests (e.g., Machielse et al., 2008; Wagenaar et al., 2012; Fleming et al., 2017, 2019, 2020), although the uncertainty on the net production gain still remains significant. Some studies on the effects of yaw misalignment on wind turbine wakes, mainly based on quasi-static approaches, have already been carried out describing the effect of WT yaw on the wake position, in wind tunnel conditions (e.g., Bastankhah and Porté-Agel, 2016; Grant et al., 1997; Howland et al., 2016; Schottler et al., 2018) and at full scale (e.g., Howland et al., 2020). However, analyzing yaw maneuver dynamics, by studying the transient process between the non-yawed and yawed conditions, affords new insights into wake interactions. This provides the opportunity to estimate the delay between the maneuver on an upstream WT and its effect on a downstream WT. Analysis of the time lag between the start of yaw motion and the beginning of wake deviation effects at the downstream wind turbine position can inform us about the wake advection time. Moreover, tracking wake response during the yaw maneuver of a wind turbine makes it possible to compare the maneuver duration with the time needed for wake stabilization in the deviated condition. The same comparisons can be made regarding the effects of the maneuver on the global load measured on a downstream wind turbine. Those analyses are very important to implement the proper strategy in the WT controllers. To this end, this study deals with a typical scenario of wake interaction between two wind turbines which is reproduced to investigate the effects of the yaw control strategy, paying particular attention to the WT wake and load dynamics. This work continues and completes the investigations already performed in a previous study (Macrì et al., 2018) and aims to characterize the wake and load response to the yaw variation. It is based on an experimental approach conducted in low-speed wind tunnels. Both WTs are then modeled with porous discs. Given the low geometric scales and the focus on far-wake characteristics, it is common practice to model wind turbine rotors with porous discs (e.g., Muller et al., 2015; Yu et al., 2017; Bastankhah and Porté-Agel, 2014; Van Gent et al., 2017) to isolate the unsteady wake behavior from the wind turbine geometric and kinetic parameters. All other sources of unsteadiness (blade wakes and rotation, tip vortices, controllers) are then avoided, to focus on the wake aerodynamics only. The use of the porous disc is justified by the fact that, in addition to its representativeness of the wind turbine far wake (Aubrun et al., 2013; Lignarolo et al., 2016), it is the most simplified wind turbine modeling approach and also the most used both for wake engineering models and for numerical simulations of single or multiple wind turbine wakes (see Porté-Agel et al., 2020). Modeling a wind turbine via a porous disc implies representing a fixed operational point of a wind turbine in terms of thrust coefficient and consequently velocity deficit within the wake, avoiding all the aforementioned additional sources of unsteadiness. Moreover, this approach is coherent with the general approach used for wind farm production optimization tools (i.e., FLORIS). This kind of representation permits a good reproducibility of the results and the possibility of reproducing the far wake of a wind turbine at a low geometrical scale with a simplified model. Indeed, as experimented in Macrì et al. (2020), it is very complex to achieve satisfactory statistical reliability and reproducibility of the results obtained through the use of full-scale experiments. The WT wake behavior and the load variation on a downstream WT in interaction are investigated with particular focus on the transition between the no-yaw condition (rotor facing the wind) and the yawed condition (in this study, 30∘ between the wind direction and the normal to the rotor area), in both maneuver directions: yaw increase (0–30∘) and yaw decrease (30–0∘). The same protocol is followed for each experiment, i.e., wake deviation tracking and downstream WT load variation monitoring. The wake behavior is analyzed by the use of particle image velocimetry measurements (stereo-PIV 2D-3C) while the load variations are analyzed using an unsteady balance of 6 degrees of freedom (DoF). PIV fields aim at measuring the profile for the far wake velocity distribution and permit the deduction of the wake deflection. The latter can be described by a wake skew angle that depends on the yaw angle and the downwind distance in particular (Bastankhah and Porté-Agel, 2016). In this study, a wake deviation angle, computed from the estimation of the wake center displacement will be considered. In addition, a wake deviation duration will be introduced to analyze the transient aspects of this wake deviation process during yaw maneuver. First the static wake deviation angle due to yaw misalignment and the associated global load modification on a downstream wind turbine are quantified, and then the dynamical properties of the wake steering process and of the global load variation during yaw maneuvers are studied. For this purpose, a realistic yaw maneuver is re-scaled in a wind tunnel to properly compare its duration with the associated wake response and load variation. In order to achieve a characterization of this phenomenon as complete as possible, the wind tunnel experiments are conducted in three different set-ups, varying the flow conditions and the geometric scales. Two experimental campaigns are carried out in a homogeneous and isotropic turbulent (HIT) incoming flow, varying the model scale and testing different flow velocities. A third experimental campaign is performed in a more realistic flow such as a modeled atmospheric boundary layer (ABL). The paper is organized as follows. In Sect. 2 the experimental set-ups and the data pre-processing are described. Section 3 presents the different indicators used to quantify the effects of a yawed upstream WT on the flow farther downstream, namely the upstream wake deviation, the available wind power density for a downstream WT model and the thrust applied to a downstream WT model. Static and dynamic influences of the yaw modification to these indicators are detailed in Sects. 4 and 5, respectively. Section 6 provides a summary of the outcomes and some conclusions.
2.1 Wind tunnels and flow conditions
In this section the different experimental set-ups are described, and then the methodology and the metrics applied to the study of the wake deviation, the available wind power variation and the disc load variation will be presented. Experiments were carried out in two wind tunnel facilities of the PRISME Laboratory at the University of Orléans. The first campaign in HIT conditions (HIT1) was performed in the Eiffel-type wind tunnel, while the other two campaigns were carried out in the closed-loop wind tunnel “Lucien Malavard” in the main test section for the HIT conditions (HIT2) and in the return test section for the modeled ABL conditions (ABL). The wake behavior during yaw variation was studied in different flow conditions and scales, and wind turbines were modeled with porous discs (Muller et al., 2015) and scaled at 1 : 800 for the first campaign and 1 : 320 for the others, giving a range of Reynolds numbers based on the disc diameter of around 104. Experiments were carried out with discs of two different porosity levels (henceforth P1 and P2) representing two different axial induction factors , where Uref is the reference wind speed at hub height and Uwake the velocity measured within the wake. As the disc porosity is homogeneous, the velocity deficit within the near wake is characterized by a constant value that was measured, out of the induction zone, to determine Uwake. The consequences of the yaw maneuver on wake steering and on the available wind power for a downstream WT were investigated at a fixed downstream distance (Δx). The fixed downstream distance was chosen according to the data set from a working wind farm studied in Garcia et al. (2019) and Macrì et al. (2020). The consequences on the aerodynamic loads applied to the downstream WT located at the same downstream distance were then performed. The wind turbine models and the spacing for the three experimental set-ups are summarized in Table 1.
In order to ensure a relevant comparison between the results of the three campaigns, the wind turbine models and the instrumentation were installed in an equivalent configuration. The set-ups basically differ for the flow condition and/or the geometric scale, as summarized in Table 1. The upstream wind turbine model was installed in the test section at a fixed position and its yaw motion was controlled by a Kollmorgen AKM24D-ANBNC-00 rotational servomotor and measured by a Kubler-type 8.5872.3832.G141 circular encoder. The disc center of this upstream wind turbine is the origin of the reference framework . For wake steering tracking and available wind power assessment, the velocity field was measured by stereo-PIV on a crosswise plane located at a downstream distance Δx (Fig. 1a). For the WT load variation measurements, a second but identical model was installed on a 6 DoF unsteady aerodynamic balance at the same downstream distance Δx as the laser light sheet and aligned with the upstream model (Fig. 1b). The motion system and the acquisition board were controlled via LabView, providing a continuous measurement of the angular WT position and the load variations, when measured. More details about the measurement systems will be given in the following sections.
Regarding the flow characteristics, in HIT conditions, turbulence grids were installed at the entrance of the main test section of the two wind tunnels in order to generate the desired turbulence characteristics as in Macrì et al. (2018). The reference wind speeds Uref, as well as the turbulence intensities in the empty field (no disc) at the disc location () and at the stereo-PIV measurement plane (), are summarized in Table 2. In ABL conditions, experiments were carried out in the return test section as was done in Muller et al. (2015), reproducing a neutrally stratified atmospheric boundary layer (Fig. 2). The modeled ABL represents a flow above a moderately rough terrain with a roughness length of m (full scale z0=0.02 m), a power-law exponent of α=0.14 and a friction velocity of m s−1. The reference wind speed at hub height (disc center), as well as the turbulence intensities at the disc location and at the stereo-PIV measurement plane, without the presence of the disc, is also summarized in Table 2. Particular attention was paid to the yaw variation scaling in order to properly represent realistic yaw maneuver dynamics in a wind tunnel. Taking as full-scale references a wind speed at hub height of 12 m s−1, a rotor diameter of 80 m and a nominal speed for yaw motion of 0.5∘ s−1, it is possible to retrieve the duration for a 0–30∘ and a 30–0∘ rotation. The yaw variation therefore takes 10τ0, where τ0 is an aerodynamic timescale based on the inflow velocity at hub height Uref and disc diameter D . For the wind tunnel condition, the yaw motion was scaled in order to have this 10τ0 duration, thus respecting the Strouhal similarity based on the wind turbine rotor dimension between the reduced and full-scale conditions. This similarity law is considered the most relevant when one studies unsteady phenomena in the wake of a bluff or porous body (Cannon et al., 1993).
2.2 PIV measurement system
The PIV system consisted of an Nd:YAG laser Evergreen 200 (2 × 200 mJ) emitting pulse with a wavelength 532 nm and a 2.5 Hz emission rate for non-synchronized acquisition. The light sheet was oriented in order to cross the test section transversely. Seeding particles were micro-sized olive oil droplets sprayed by a PIVTEC seeding system. Images are acquired with two LaVision Imager LX cameras (4032 px × 2688 px) with a 105 mm lens equipped with a 532 nm wavelength filter. The time delay between the images was set according to flow speeds, at 105 µs for the HIT1 set-up; at 126, 84 and 42 µs for the HIT2 set-up (for Uref of 5, 10 and 15 m s−1, respectively); and 126 µs for the ABL set-up. The images were processed with a multi-pass decreasing size (64 px × 64 px, 32 px × 32 px) interrogation window with an overlap of 50 %. An ensemble average of the collected velocity fields is then performed. For dynamic conditions, the acquisition was triggered according to the yaw motion progress. The collection of hundreds of image pairs was triggered at integer multiples of the timescale τ0, and a conditional averaging approach was then applied. Indeed, due to the PIV system characteristics, during one yaw cycle from 0 to 30∘ and vice versa, it was possible to acquire only one image pair at a chosen time delay. The yaw cycle was then reproduced either 300 times or 1000 times depending on the set-up (see Table 2), and a conditional averaging was applied in order to achieve the statistical convergence of the results. By varying the time delay from 3 to 20τ0, it is possible to reconstruct the phase-averaged velocity field evolution due to the yaw motion. The maximum statistical uncertainty of the mean wind speed and its standard deviation was defined according to Benedict and Gould (1996), with Z=1.96 (confidence interval of 95 %), Nb the number of independent samples, and the maximum turbulence intensity measured in the wake region. These parameters, together with the model size diameter D, the spacing between the model and the stereo-PIV measurement plane location Δx, and the dimensionless vector resolution RV, defined as the distance between two vectors in the PIV velocity field normalized by the disc diameters (), are summarized for HIT and ABL conditions in Table 2. It has to be mentioned that the present vector resolution was obtained after a bilinear interpolation on a mesh 2 times finer than original PIV resolution.
2.3 Load measurement system
The load measurements were performed by an unsteady six-component aerodynamic balance ATI™ model mini 40. The balance was mounted on a rigid structure located underneath the test section floor, and the downstream wind turbine model was installed on a specifically designed support. The balance has six analogic channels that, as for the encoder measuring the yaw variation, are acquired at a frequency of 2 kHz by a National Instrument card. The whole system is controlled by LabView. It should be noted that the aerodynamic balance can detect the asymmetric loading on the downstream WT. Nevertheless the analysis will be focused on the thrust measurement because that is the only loading that can be related to a notion of WT performance. In the case of static load measurements (no yaw variation in the upstream model), the sampling time was 2 min. In the case of dynamic configurations, acquisitions were designed in order to perform cycle-averaging of load for a minimum of 500 consecutive 0–30∘ backward and forward yaw displacement cycles, separated by a pause of at least the same duration as the yaw motions. A pre-processing of the acquired balance time series is needed to filter out signal fluctuations due to the natural frequency of the balance. A zero-phase digital low-pass-filtering strategy was applied to the balance signal in order to filter out the balance resonance without rejecting the first harmonics of the cyclic yawing frequency. The filtering tuning was driven by the duration metrics assessment (Sect. 5.1), in order to limit its effect on the ramp slope, despite the residual overshoot at the edges. An example of the load measurement data processing for a dynamic yaw variation is given for the HIT2 set-up in Fig. 3. Figure 3a shows the power spectral density of the yaw variation and of the low-pass-filtered and non-filtered model drag fluctuations. After subtracting the offset voltages measured before starting the experiment, the six low-pass-filtered voltage time series were converted into the six load components via a conversion matrix provided by the balance constructor. The 6 degrees of freedom load tensor is retrieved, but the single load of interest for the present study is the WT model drag force Fx. Due to the balance resolution and sensitivity, it was chosen to retrieve the drag force by the measurement of the moment along the transverse axis TY divided by the lever arm b (). By analogy with a real wind turbine, drag force will be called “thrust” T in the following.
For dynamic conditions, the last part of the data processing consists in the proper phase averaging on the rising and descending yaw ramp. Indeed, thanks to the synchronization of the encoder and the balance acquisitions, it was possible to isolate each ramp of yaw variation (henceforth positive yaw variation for 0–30∘ and negative yaw variation for 30–0∘) relying on the yaw measurement. In this way it was possible to properly and repeatedly detect the start and the end of each ramp in order to make a representative phase average of all the ramps. Figure 3b shows, as an example, the averaging results for a “cycle averaging” approach. The results highlight the post-processing protocol efficacy of the averaging effect since the shape of the time evolution of in the absence of a filtering strategy shows that the phase average itself is not sufficient to filter out the balance resonance. The shape of the shows the smoothing effect of the filtering strategy. Moreover, due to the general superposition of the two curves, it is assumed that the residual oscillations of the non-filtered cycle-averaged curve were related to the balance resonance and that filtering out its resonance frequency does not distort the main pattern. According to the balance specifications, the torque measurement uncertainty (95 % confidence level) is 1.75 % of the full-scale load (4 Nm). This is too close to the torque values measured for both porosity levels in the HIT1 configuration and for the higher porosity level in ABL conditions. Consequently, these three configurations will not be considered in the further load analysis.
3.1 Upstream WT wake center position
The wake deviation was determined by estimating the displacement of the wake center position Yc from Yc=0 when the WT model is normal to the free-stream flow (yaw angle equal to zero) to Yc≠0 when the WT model is misaligned (yaw angle non-null). The wake deviation angle θ can then be easily retrieved by simple trigonometric considerations (Fig. 4). These parameters are linked by the relation
where χ is the skew angle of the wake and γ the yaw angle of the WT model. Different methods to determine the wake center were tested (Muller et al., 2015; Howland et al., 2016; Parkin et al., 2001; Vollmer et al., 2016; Schottler et al., 2017), and the most robust for the present database was the one based on the estimation of the available wind power density of a potential downstream wind turbine, as in Vollmer et al. (2016).
The available wind power density (Pav, W m3 kg−1) corresponds to the kinetic energy density crossing a virtual downstream wind turbine rotor-swept area and normalized by the air density as in Eq. (2), where u is the streamwise velocity component over a plane normal to the flow direction. The available wind power density is calculated for all the possible rotor position in the ranges of and . The wake center position is so determined by the values of y0 and z0 that minimize the available power: and .
With this method, it is possible to analyze both crosswise coordinates of the wake center. Due to the negligible variation detected over the vertical coordinate Zc, only the horizontal displacement of the wake will be considered. A detailed comparison between most of the aforementioned approaches applied to the current paper data set has been done in Macrì (2020). This comparison led to the choice of the method based on the available power density (Vollmer et al., 2016). Indeed, this, due to the integration domain definition (see above), reduces the potential sources of bias due to the PIV possible lower quality at the boarder as detailed in Macrì (2020). As concerns measurement uncertainties, taking into account the uncertainties on the wind speed, the PIV vector resolution (Table 2) and the methods used to estimate the wake center, it was possible to estimate the maximal measurement error for both set-ups by applying usual resolution-based error estimation methods. The estimated measurement uncertainty is and m for HIT1 conditions and and m for HIT2 and ABL conditions.
3.2 Available wind power density at the downstream WT position
To complement the local indicator of wake deviation, which is the wake center position, the more integrated indicator provided by the above available wind power density at the downstream WT model location was also used. This metric is interesting to analyze because it is a good intermediate between the wake deviation of the upstream WT model and the actual load response of the downstream WT model, when it is located at the same position as the PIV measurement plane.
3.3 Downstream WT thrust coefficient
The thrust coefficient CT definition is given (Eq. 3) by
where T is the thrust force (equivalent to the drag force in the present modeling approach) measured by the aerodynamic balance, ρ the air density, Uref the free-stream speed and AD the disc area. This parameter is a direct and reliable indicator of the dimensionless load that the model is subjected to. Moreover, since the principal aim of the study is to evaluate load variation, it was chosen to non-dimension CT with the thrust coefficient at γ=0∘ or at the start of the yaw maneuver in the case of dynamic measurements.
In this section the results concerning the influence of the static yaw angle applied to an upstream WT on its wake deviation, on the available wind power for a virtual downstream WT and on the actual thrust applied to a downstream WT, will be provided. They will serve as the baseline for the remainder of the study, in which dynamics will be added to the yaw motion system.
4.1 Wake center deviation
The measurements were performed for yaw angles from γ between 0 and 30∘ by steps of 10∘. The static wake deviation was studied for all the experimental configurations listed in Table 2 and for both disc porosity levels. Figure 5a shows a summary of the wake deviation angle versus the WT yaw angle for static conditions. Several general comments can be made. Indeed, as already shown by previous studies (Parkin et al., 2001; Howland et al., 2016; Espana, 2009), the relationship between the wake deviation angle and the yaw angle is a nonlinear monotonically increasing function and the skew angle is 1 order of magnitude lower than the yaw angle. Theoretically, due to the absence of rotational entrainment in the wake of a porous disc and the absence of Coriolis force at such a reduced scale of observation, the absolute value of the wake deviation angle is identical for negative or positive yaw angles. A different behavior had been observed for a rotating wind turbine model, with a light dependence of the wake deviation angle on the direction of misalignment (Bastankhah and Porté-Agel, 2016; Bartl et al., 2018). It is assumed in the present work that this asymmetry does not play a major role in the wake dynamics and therefore is not studied. As concerns the measurement quality, the results show some scatter inherent to the propagation of cumulative errors (measurement, statistical, processing) and illustrate the difficulty of accurately determining such a small deviation angle (maximum measurement uncertainty ). Regarding the impact of the flow conditions on the wake deviation, it is not possible to detect any dependence of the results on flow conditions. Indeed, a scatter in values of the wake deviation angle or the thrust is visible, but it cannot be linked to any flow conditions. The porosity level (and hence the equivalent induction factor) affects the wake deviation: cases with a lower porosity level (higher induction factor) present a higher wake deviation. In Macrì et al. (2018), the present authors made some comparisons between similar experimental results and empirical wake deviation models. The trends were similar but the models systematically overestimated the wake deviation compared to experimental values. The results are also dependent on the wake center definition used (not shown here but mentioned in Coudou et al., 2018). Indeed, there are several wake center tracking methods, and their extrapolation to skewed wakes is still under discussion. Looking at the presented results more in detail, some specific discrepancies can be stressed. Indeed, the HIT1 P2 configuration presents a higher deviation than the other ones with the same porosity level P2, without any straightforward physical ground. This can be mostly attributed to the sensitivity of the center tracking method as discussed in Macrì (2020). Moreover, the ABL P2 configuration presents a discrepancy between its trend (especially at ) and the other results at the same porosity level. This is because the flow inhomogeneity together with the higher level of ambient turbulence make the velocity deficit generated by the higher porosity disc rather small and unsuitable to properly track the wake center. For these reasons, the ABL P2 configuration will not be discussed further. It is important to stress that the discrepancies that can be seen between the cases HIT2b P1 and HIT2c P1 and the rest of the cases are not unusual in experimental measurements (Aubrun et al., 2019) and may be due to a minor variation in the performances of one of the experimental measurement systems while changing the set-up.
4.2 Downstream WT thrust coefficient variation
Looking at the results it is possible to state that the relationship between the downstream WT thrust coefficient and the yaw angle is a nonlinear monotonically increasing function. This indicator presents less scatter than the wake deviation angle, and the results do not show any unexplained outliers. Its integrative nature (global load applied to the porous disc) plays a smoothing role, and the range of variations is also more significant. Additionally, there is a clear thrust difference depending on the porosity level. For the P1 porosity cases, the thrust gain of the downstream WT for a 30∘ yaw angle variation in the upstream WT compared to the 0∘ case is around 13 %, irrespective of the flow conditions, while for the lower porosity level, it is about 3 %. These values suggest a remarkable influence of the porosity (and hence the axial induction factor) on the consequences of a yaw modification of an upstream WT on the load applied to a downstream one.
5.1 Metrics for dynamics
To analyze the yaw maneuver effects on the upstream WT wake and downstream WT thrust dynamics, some metrics to determine the properties of the transient phenomena need to be defined. It is necessary to establish a common protocol (equally reliable) to assess the transient duration, start and end for the three phenomena: yaw maneuver, wake deviation and thrust variation. First, the cycle-averaged yaw angle, wake deviation and thrust coefficient values were normalized in order to have transient curves between 0 and 1:
with γstart and γend the γ values before and after the transient. It has to be noted that for positive yaw maneuver, γstart=0, and for negative yaw maneuver, γend=0.
with θstart and θend the θ values before and after the transient.
with and the CT values before and after the transient.
In order to facilitate the determination of the start and end of transient phenomena, fitting laws were then applied to the transient curves, depending on the yaw maneuver sign. For the positive yaw maneuver, the cycle-averaged results were fitted to Eq. (7) while for the negative yaw variation they were fitted to Eq. (8).
For both equations, c and τlag are the fitting coefficients, τ the time triggered with the yaw maneuver start, and τ0 the aerodynamic timescale. c represents the variation rate of the function and can be linked to the transient duration, whereas τlag represents the beginning of the function evolution and can be linked to the transient start. The fitting coefficients were determined by a classical nonlinear least-square fitting method. Regarding the yaw angle, a fitting procedure was not necessary because of the high time resolution and precision of the acquisition. Finally, the transient start τstart corresponds to the time when the curve crosses the 5 % threshold of the total variation and the transient end τend to the time when the curve crosses the 95 % threshold. Dimensionless values are obtained by dividing all times by the aerodynamic timescale τ0: and . The transient dimensionless duration is obtained by retrieving the time difference between these two dimensionless values: .
Other metrics defined for this study are and . The first one represents the ratio between the wake deviation duration and the maneuver duration . The second represents the ratio between the advection velocity and the reference wind speed, where the advection velocity is defined as the ratio between the streamwise spacing Δx between the upstream WT model and the downstream location where the wake deviation is investigated (or where the downstream wind turbine model is located) and the delay between the start of the maneuver and the start of the wake deviation . In order to obtain intermediate information between the upstream WT wake deviation and the response in terms of thrust modification on the downstream WT, the same metrics were also calculated for the available wind power density at the downstream WT position (see Sect. 3.2).
All the metrics applied to the dynamic wake deviation of the upstream WT, and to the thrust variation in the downstream WT, are summarized in Tables 3 and 4, respectively. It is important to mention that for cases 10 and 11, the yaw maneuver speed was doubled in order to evaluate the influence of a speed-up.
5.2 Transient durations
The first parameter to analyze in order to characterize the dynamical consequences of the yaw maneuver is , since this parameter is useful to compare the duration of the yaw maneuver with the duration of the induced wake flow and load modifications. Figure 8 shows an intuitive visual way of analyzing ; by plotting and against the maneuver duration Δτm, values far from the diagonal line indicate a dynamic behavior of the wake deviation or of the thrust variation different than the yaw maneuver. If values are below the diagonal, the response to the maneuver is faster than the yaw maneuver duration; if values are above the diagonal, the response to the maneuver is slower than the maneuver duration. The results are classified according to the disc porosity level and the yaw maneuver direction. In general, wake deviation durations are all around the line, except for one outlier (case 11) where the difference can be explained by a higher RMSE for the fitting procedure (due to a more scattered wake center evolution). No significant difference between the wake dynamics and the yaw maneuver can be noticed, regardless of the disc porosity levels and the flow conditions. In contrast, thrust variation durations are in general lower than the yaw maneuver duration. The influence of the porosity level can be noticed, and in general, cases of negative yaw variation evolve faster than cases of positive yaw variation, especially for porosity level P1. Concerning the cases of higher yaw maneuver speed, no clear trend is visible on the effect of porosity and yaw maneuver. This may be due to the smaller number of cases analyzed, but generally a higher yaw maneuver speed does not seem to significantly impact the behavior of the downstream WT thrust variation.
In Tables 3 and 4, the fitting coefficient c that represents the transient variation rate is given. There is a direct relationship between this parameter and the dimensionless transient durations Δθ*(τ) and deduced from thresholds (Sect. 5.1). This fitting parameter is not further exploited in the present study, but its robustness suggests that it could be used to set up dynamic models of the wake deviation or of the thrust variation.
5.3 Transient starts and ends
In order to obtain more detailed knowledge of the timing characteristics and the history of modification, it is important to check and for wake deviation, available wind power density and thrust variations, since values give information about the time delays before the upstream WT wake starts deflecting, before the available wind power density for the downstream WT is modified and then before the thrust applied to the downstream WT starts increasing. values inform about the delays necessary for the stabilization of all the variations in the final state.
Figure 9 shows a summary of the time parameters, through a timeline representation for the wake deviation, the available wind power density and the thrust variations. The transient duration values are also shown in a parallel plot to facilitate reading. From the analysis of the aforementioned parameters, several conclusions on the wake dynamics can be drawn. In general, the wake deviation transient (green timelines) has a slightly shorter duration and starts later for the negative yaw variation than for the positive one. As explained in Sect. 5.1, a time delay between the start of the yaw maneuver and the wake deviation is expected since the wake deviation is observed at a certain downstream distance. Air mass needs time to travel from the upstream WT model to the downstream location of interest. This wake transport velocity, or advection velocity, is generally assessed as equal to the free-stream velocity Uref (Trujillo et al., 2011), but some previous studies also proposed a lower value due to the velocity deficit within the wake of 0.8Uref (Machefaux et al., 2015) or an average between the free-stream velocity and the wake speed (Bossanyi, 2018; Keck et al., 2014). This information is checked through the ratio reported in Table 3. The ratio is presently close to 1 for a positive yaw maneuver and around 0.8 for a negative yaw maneuver. A detailed analysis of timing parameters does not provide any evidence of a major effect of the porosity on the wake deviation behavior. Concerning the cases of higher yaw maneuver speed, no particular impact of the yaw maneuver speed on the wake response is visible. Theoretically, the fitting coefficient τ*lag lag (Tables 3 and 4) can be interpreted as a time delay before the transient starts. The relationship between this parameter and will be investigated in Sect. 5.4 together with the relation between c and the transient duration.
Figure 9 also presents the transient parameters for the available wind power density deduced from the wake measurement by PIV at the downstream location where the second WT model will be installed (black timelines). Systematic differences in the timing parameters of this available wind power density compared to the wake center deviation are visible, with shorter starting time delays for the positive yaw maneuver and longer ones for the negative yaw maneuver. As this value is obtained by a space integration of the wake velocity, it illustrates that the modification of the overall wake velocity field can be different from the modification of the wake center position, which is a more local indicator and for which the range of deviation to be captured is very small. This available wind power density dynamics also illustrates how the incoming flow will dynamically impact the downstream WT model and modify its corresponding thrust. It was therefore expected that the thrust transient starts and ends (red timelines) will be systematically later than the available wind power ones. The additional time delay between both would represent the WT model response to the modification of the incoming flow. This trend is not verified in the present results since both situations can be observed. The reasoning was based on the hypothesis that the induction zone of the downstream WT model does not play a role in its dynamic process, which is clearly not true.
The thrust variation starts later for a positive yaw maneuver than for a negative one. This trend is in opposition with the wake deviation behavior and contradicts the assumption that one can use the same time delay due to advection for the upstream WT wake deviation and for the downstream WT thrust variation.
5.4 Interpretation of fitting law coefficients
In this section, an evaluation of the fitting coefficients of the exponential law is done with respect to the timing parameters used for the transient analysis. In Fig. 10a, the wake deviation duration is plotted against the c coefficient of the fitting law for both positive and negative maneuver and porosity levels, while in Fig. 10b, the same analysis is done for the thrust variation duration . Generally c has a very clear correlation with the duration for both wake deviation and thrust variation cases. Indeed, it is possible to retrieve, by linear fitting imposing the fit law to pass through the origin, a linear relationship that links c to both wake deviation and thrust variation duration. This linear fitting was done separately for the wake deviation and thrust variation but led to the same slope of 0.93:
where Δτ* is or the wake deviation duration either the thrust variation duration . These results confirm the robustness of the fit coefficient c, although the reason that c is constantly lower than the phenomenon duration has still to be investigated. In Fig. 11a, the start of the wake deviation is plotted against the coefficient of the exponential fitting laws for both positive and negative maneuver and porosity levels. In Fig. 11b, the same analysis is done for the start of the thrust variation. Several considerations on the coefficient can be made. Indeed, although a higher scatter than for the c coefficient can be observed, also shows a detectable trend. The thrust data seem to have lower scatter than the wake deviation data, and this could be related to the higher time resolution of the load measurements. There is generally a bias between the values that oscillate between 1 and 2 τ0. This could be partially attributed to the threshold value used for the determination.
In conclusion, the fitting laws implemented for the dynamic yaw conditions seem to be exploitable for a generalization of a law describing the transition phenomenon. In particular, the c coefficient has proven to be quite robust to be used in an empirical law to model the transient. Regarding the coefficient, although its physical signification is clear, its implementation in an empirical law demands more caution. Indeed, the measured values do not follow a trend as clear as the one observed for the duration coefficient. This is probably due to the specificity of this parameter, which is directly linked to the start of the phenomenon, and consequently it is more sensible to dispersion. Nevertheless, , being representative of advection, could be adjusted by making some simple advection hypothesis. The higher robustness of the c parameter, being representative of the transient duration, can probably be due to the fact that it better damps the scatter of and , which are generally concordant (both either overestimated or underestimated). At this stage, further exploitation of these data and studies have to be done in view of the implementation of an empirical model.
The aerodynamic characterization of the wake of a wind turbine model and its effect on the load of a similar downstream wind turbine model consequent to a dynamic positive (misalignment scenario) or negative (realignment scenario) yaw variation were experimentally studied for different incoming flow conditions, Reynolds scales and induction factors by using porous discs. Wake deviation, available wind power or thrust variations were the main metrics estimated from PIV and aerodynamic load measurements. First, the wake deviation angle and the thrust coefficient of the downstream wind turbine model were analyzed as a function of different yaw angles in static conditions. Then, the duration, start and end of the temporal response of the metrics to a dynamic yaw variation were estimated and compared. The main results are summarized below. Overall results do not show any noticeable influence of the flow conditions (homogeneous isotropic flows or atmospheric boundary layer flow), or of the Reynolds scales on the static and dynamical properties of the different metrics. The influence of the degree of physical modeling of the wind turbine (porous disc versus rotating wind turbine model) on the results had not been studied since it was assumed that this feature does not play a major role in the yawed far-wake dynamics, but this question needs to be further investigated. Concerning the characterization of the magnitude of the wake deviation in static yaw conditions, results show that the wake deviation is 1 order of magnitude lower than the yaw increment and that the relationship between the wake deviation and the yaw angles is nonlinear as previously found in the literature. In the case of a higher porosity level, the thrust increment is much lower than in the case of lower porosity. However, while a relevant influence of the porosity on the wake deviation angle and thrust magnitude of the downstream wind turbine model is found in static conditions, no significant influence of porosity is observed in dynamic yaw variation conditions. On the other hand, in these conditions, the analysis of the three metrics reveals different temporal characteristics depending on whether yaw variation is positive or negative. In general, the wake deviation transient has slightly shorter duration and starts later for the negative yaw variation than for the positive one. A proper knowledge of such effects could help to determine the time necessary for the wake steering to totally develop. Conversely, the thrust variation starts later for the positive yaw maneuver than for the negative one. The influence of the yaw maneuver speed was tested, and doubling the yaw maneuver speed does not seem to influence the wake or load dynamics. Finally, the study shows that the same dynamical properties of the wake of the upstream wind turbine and load variation in the downstream wind turbine cannot be generalized to any yaw variation configurations and that the advection velocity should be assumed to be different according to the yaw maneuver direction. A first step in proposing fitting coefficients was made in order to support a dynamic model of the wake deviation or the thrust for positive or negative yaw variations. However, further development will be performed to confirm the reliability of the proposed models.
Post-processed data can be made available upon request.
SA, AL and NG co-supervised the PhD work of SM. SM contributed to the design of the experimental set-ups, conducted all experiments, and processed and analyzed all data. SM, SA, AL and NG contributed to the interpretation of the results. SM, SA and AL contributed to the manuscript writing.
The authors declare that they have no conflict of interest.
The present work is part of a PhD co-funded by Labex Caprysses and Engie.
This research has been supported by Labex Caprysses (ANR-11-LABX-0006-01) and Engie SA.
This paper was edited by Raúl Bayoán Cal and reviewed by Nicholas Hamilton and one anonymous referee.
Aubrun, S., Loyer, S., Hancock, P., and Hayden, P.: Wind turbine wake properties: Comparison between a non-rotating simplified wind turbine model and a rotating model, J. Wind Eng. Ind. Aerod., 120, 1–8, 2013. a
Aubrun, S., Bastankhah, M., Cal, R. B., Conan, B., Hearst, R. J., Hoek, D., Hölling, M., Huang, M., Hur, C., Karlsen, B., Neunaber, I., Obligado, M., Peinke, J., Percin, M., Saetran, L., Schito, P., Schliffke, B., Sims-Williams, D., Uzol, O., Vinnes, M. K., and Zasso A.: Round-robin tests of porous disc models, J. Phys.-Conf Ser., 1256, 012004, https://doi.org/10.1088/1742-6596/1256/1/012004, 2019. a
Bartl, J., Mühle, F., Schottler, J., Sætran, L., Peinke, J., Adaramola, M., and Hölling, M.: Wind tunnel experiments on wind turbine wakes in yaw: effects of inflow turbulence and shear, Wind Energ. Sci., 3, 329–343, https://doi.org/10.5194/wes-3-329-2018, 2018. a
Bastankhah, M. and Porté-Agel, F.: A new analytical model for wind-turbine wakes, Renew. Energ., 70, 116–123, 2014. a
Benedict, L. and Gould, R.: Towards better uncertainty estimates for turbulence statistics, Exp. Fluids, 22, 129–136, 1996. a
Bossanyi, E.: Combining induction control and wake steering for wind farm energy and fatigue loads optimisation, in: J. Phys.-Conf. Ser, 1037, 032011, https://doi.org/10.1088/1742-6596/1037/3/032011, 2018. a
Cannon, S., Champagne, F., and Glezer, A.: Observations of large-scale structures in wakes behind axisymmetric bodies, Exp. Fluids, 14, 447–450, 1993. a
Coudou, N., Moens, M., Marichal, Y., Van Beeck, J., Bricteux, L., and Chatelain, P.: Development of wake meandering detection algorithms and their application to large eddy simulations of an isolated wind turbine and a wind farm, J. Phys.-Conf. Ser., 1037, 072024, https://doi.org/10.1088/1742-6596/1037/7/072024, 2018. a
Espana, G.: Étude expérimentale du sillage lointain des éoliennes à axe horizontal au moyen d'une modélisation simplifiée en couche limite atmosphérique, PhD thesis, Université d'Orléans, 2009. a
Fleming, P., Annoni, J., Shah, J. J., Wang, L., Ananthan, S., Zhang, Z., Hutchings, K., Wang, P., Chen, W., and Chen, L.: Field test of wake steering at an offshore wind farm, Wind Energ. Sci., 2, 229–239, https://doi.org/10.5194/wes-2-229-2017, 2017. a
Fleming, P., King, J., Dykes, K., Simley, E., Roadman, J., Scholbrock, A., Murphy, P., Lundquist, J. K., Moriarty, P., Fleming, K., van Dam, J., Bay, C., Mudafort, R., Lopez, H., Skopek, J., Scott, M., Ryan, B., Guernsey, C., and Brake, D.: Initial results from a field campaign of wake steering applied at a commercial wind farm – Part 1, Wind Energ. Sci., 4, 273–285, https://doi.org/10.5194/wes-4-273-2019, 2019. a
Fleming, P., King, J., Simley, E., Roadman, J., Scholbrock, A., Murphy, P., Lundquist, J. K., Moriarty, P., Fleming, K., van Dam, J., Bay, C., Mudafort, R., Jager, D., Skopek, J., Scott, M., Ryan, B., Guernsey, C., and Brake, D.: Continued results from a field campaign of wake steering applied at a commercial wind farm – Part 2, Wind Energ. Sci., 5, 945–958, https://doi.org/10.5194/wes-5-945-2020, 2020. a
Garcia, E. T., Aubrun, S., Coupiac, O., Girard, N., and Boquet, M.: Statistical characteristics of interacting wind turbine wakes from a 7-month LiDAR measurement campaign, Renew. Energ., 130, 1–11, 2019. a
Gebraad, P., Thomas, J. J., Ning, A., Fleming, P., and Dykes, K.: Maximization of the annual energy production of wind power plants by optimization of layout and yaw-based wake control, Wind Energy, 20, 97–107, 2017. a
Grant, I., Parkin, P., and Wang, X.: Optical vortex tracking studies of a horizontal axis wind turbine in yaw using laser-sheet, flow visualisation, Exp. Fluids, 23, 513–519, 1997. a
Howland, M. F., Bossuyt, J., Martínez-Tossas, L. A., Meyers, J., and Meneveau, C.: Wake structure in actuator disk models of wind turbines in yaw under uniform inflow conditions, J. Renew. Sustain. Energ., 8, 043301, https://doi.org/10.1063/1.4955091, 2016. a, b, c
Howland, M. F., González, C. M., Martínez, J. J. P., Quesada, J. B., Larranaga, F. P., Yadav, N. K., Chawla, J. S., and Dabiri, J. O.: Influence of atmospheric conditions on the power production of utility-scale wind turbines in yaw misalignment, J. Renew. Sustain. Energ., 12, 063307, https://doi.org/10.1063/5.00237462020. a
Keck, R., Mikkelsen, R., Troldborg, N., de Maré, M., and Hansen, K.: Synthetic atmospheric turbulence and wind shear in large eddy simulations of wind turbine wakes, Wind Energy, 17, 1247–1267, https://doi.org/10.1002/we.1631, 2014. a
Lignarolo, L., Ragni, D., Ferreira, C., and Van Bussel, G.: Experimental comparison of a wind-turbine and of an actuator-disc near wake, J. Renew. Sustain. Energ., 8, 023301, https://doi.org/10.1063/1.4941926, 2016. a
Machefaux, E., Larsen, G. C., Troldborg, N., Gaunaa, M., and Rettenmeier, A.: Empirical modeling of single-wake advection and expansion using full-scale pulsed lidar-based measurements, Wind Energy, 18, 2085–2103, 2015. a
Machielse, L., Barth, S., Bot, E., Hendriks, H., and Schepers, J.: Evaluation of “Heat and Flux” Farm Control, Final Report, Tech. rep., ECN-E–07-105, 2008. a
Macrì, S., Coupiac, O., Girard, N., Leroy, A., and Aubrun, S.: Experimental analysis of the wake dynamics of a modelled wind turbine during yaw manoeuvres, J. Phys.-Conf. Ser., 1037, 072035, https://doi.org/10.1088/1742-6596/1037/7/072035, 2018. a, b, c
Macrì, S., Duc, T., Leroy, A., Girard, N., and Aubrun, S.: Experimental analysis of time delays in wind turbine wake interactions, J. Phys.-Conf. Ser., 1618, 062058, https://doi.org/10.1088/1742-6596/1618/6/062058, 2020. a, b
Pao, L. Y. and Johnson, K. E.: A tutorial on the dynamics and control of wind turbines and wind farms, in: 2009 American Control Conference, 2076–2089, 2009. a
Parkin, P., Holm, R., and Medici, D.: The application of PIV to the wake of a wind turbine in yaw, in: Conference of Particle Image Velocimetry, Gottingen, Germany, 17–19 September 2001, 155–162, 2001. a, b
Porté-Agel, F., Bastankhah, M., and Shamsoddin, S.: Wind-turbine and wind-farm flows: a review, Bound.-Lay. Meteorol., 174, 1–59, 2020. a
Sanderse, B.: Aerodynamics of wind turbine wakes-literature review, Technical report E-09-016, Energy research Centre of the Netherlands, 2009. a
Schottler, J., Mühle, F., Bartl, J., Peinke, J., Adaramola, M. S., Sætran, L., and Hölling, M.: Comparative study on the wake deflection behind yawed wind turbine models, J. Phys.-Conf. Ser., 854, 012032, doi::10.1088/1742-6596/854/1/012032, 2017. a
Schottler, J., Bartl, J., Mühle, F., Sætran, L., Peinke, J., and Hölling, M.: Wind tunnel experiments on wind turbine wakes in yaw: redefining the wake width, Wind Energ. Sci., 3, 257–273, https://doi.org/10.5194/wes-3-257-2018, 2018. a
Trujillo, J., Bingöl, F., Larsen, G., Mann, J., and Kühn, M.: Light detection and ranging measurements of wake dynamics. Part II: two-dimensional scanning, Wind Energy, 14, 61–75, https://doi.org/10.1002/we.402, 2011. a
Van Gent, P., Michaelis, D., Van Oudheusden, B., Weiss, P.-É., de Kat, R., Laskari, A., Jeon, Y. J., David, L., Schanz, D., Huhn, F., Gesemann, S., Novara, M., McPhaden, C., Neeteson, N. J., Rival, D. E., Schneiders, J. F. G., and Schrijer, F. F. J.: Comparative assessment of pressure field reconstructions from particle image velocimetry measurements and Lagrangian particle tracking, Exp. Fluids, 58, 33, https://doi.org/10.1007/s00348-017-2324-z, 2017. a
Vollmer, L., Steinfeld, G., Heinemann, D., and Kühn, M.: Estimating the wake deflection downstream of a wind turbine in different atmospheric stabilities: an LES study, Wind Energ. Sci., 1, 129–141, https://doi.org/10.5194/wes-1-129-2016, 2016. a, b, c
Wagenaar, J., Machielse, L., and Schepers, J.: Controlling wind in ECN’s scaled wind farm, Proc. Europe Premier Wind Energy Event, Copenhagen, Denmark, 16-19 April 2012, 685–694, 2012. a
Yu, W., Hong, V., Ferreira, C., and van Kuik, G.: Experimental analysis on the dynamic wake of an actuator disc undergoing transient loads, Exp. Fluids, 58, 149, https://doi.org/10.1007/s00348-017-2432-9, 2017. a