1. Introduction
Wind turbine rotors extract energy from wind that passes through them. This energy extraction makes changes in the wind flow so that the wind speed is reduced and the turbulence intensity is increased downstream. This phenomenon, known as the wind turbine wake, reduces the power production and increases the fatigue load of the downstream wind turbines [
1,
2]. Especially in a wind farm composed with many wind turbines, this wake effect generally reduces the power production efficiency of the wind farm by 10 to 20%, so it is common to optimize the layout of the wind farm to minimize the wake effect [
1,
3]. However, since further optimization of the wind farm layout is not possible after constructing a wind farm, the active control of the wake using blade pitching and nacelle yawing has been investigated recently [
4,
5,
6,
7,
8,
9,
10].
The application of a wind farm control algorithm to an actual wind farm for verification is difficult to be realized, so most research on this topic is performed through simulation [
4,
5,
6,
7,
8,
9,
10]. Also, high fidelity computational fluid dynamics (CFD) codes cannot be used for this purpose because of large computational loads, and therefore, simplified wind farm simulation codes including wind turbine dynamics and control are used [
11].
The accuracy of the wake model in wind farm simulation codes is important because the performance of wind farm control is largely affected by the wind speed prediction from the wake model [
12,
13,
14]. For this reason, many researchers are making efforts to verify the accuracy of the wake models and to improve it through wind tunnel tests or high fidelity numerical simulations. Wind tunnel tests and numerical simulations are very useful methods to develop new wake models or verify their accuracy because they can have a controlled flow environment repeatedly [
15,
16].
The wake model that can be implemented to wind farm simulation codes can be divided into two, the kinematic wake model and the field wake model [
17]. The kinematic wake model, also known as the explicit model, basically predicts the wind speed deficit in the wake based on the momentum equation. The kinematic wake model also features a self-similar wind speed profile. The kinematic wake model is also widely used in research fields such as predicting the annual power production of wind farms or wind farm design because it is possible to predict the wake field very quickly through simple equations [
4,
5,
6,
7,
8,
9,
10].
The kinematic wake models have limitations in that the effect of turbulent momentum change on the wake growth and expansion cannot be considered. Therefore, the kinematic wake models fill the physically missing parts with empirical wind speed profiles and other parameters that can take into account the turbulence effects [
17]. The Jensen wake model is one of the oldest kinematic wake models and is still widely used. The Jensen wake model predicts the wake using an unrealistic wind speed profile called the ‘top hat’ and a wake decay constant which is the wake recovery rate parameter. However, a relatively recent study by Kim showed that the Jensen wake model is not less accurate than other classical kinematic wake models with more realistic wind speed profiles such as the wake models by Larsen or Frandsen [
18,
19,
20]. The wake decay constant in the Jensen model is an empirical parameter used to consider the effect of the ambient turbulence intensity due to the surface roughness on the wake [
21]. The original Jensen wake model was modified by Katic so that the wake behind a wind turbine with power regulation by stall or pitch could be calculated differently by considering the change in the thrust coefficient [
18].
Several studies were carried out to construct a wake model with a higher accuracy while retaining the advantages of the kinematic wake model by correcting the wake decay constant or the wind speed profile [
1,
22,
23,
24,
25,
26]. Wind Atlas Analysis and Application Program (WAsP), and recommends using the empirically derived wake decay constants of 0.05 and 0.075 for onshore and offshore, respectively. However, Barthelmie and Jensen found that the power production of a wind farm can be changed due to the atmospheric conditions through the measurement in Nysted offshore wind farm. They also found that the wake decay constant of 0.03 is suitable for Nysted wind farm. This value is smaller than the wake decay constant proposed by WAsP for offshore [
1]. Peña et al. proposed a new wake decay constant calculation formula that can take into account the atmospheric stability to increase the accuracy of the Jensen wake model [
27]. However, they did not do enough to fully validate the proposed model [
21]. Tian pointed out that the wake growth is greatly affected by the ambient turbulence intensity which is not considered in the Jensen model, and he proposed a new wake decay constant formula considering the effect of the atmospheric turbulence intensity. The proposed model, however, showed improved prediction compared to the original Jensen wake model only for the below rated wind speed region. The proposed wake model was found to have large errors for different thrust coefficients (above rated wind speed region).
Kinematic models have also been modified to be more accurate by combining momentum conservation with the data obtained from wind tunnel tests or numerical simulations. Bastankhah and Porté-Agel (BP) proposed a new kinematic wake model with a Gaussian wind profile using the data from Large Eddy Simulation (LES) model and verified it by LES simulations with various roughness lengths [
28]. However, their wind turbine in the LES simulation was different from a typical megawatt class wind turbine in that the tip speed ratio was only between 3 and 4 which is significantly smaller than that of modern wind turbines, which are between 7 and 10 [
29]. Also, the validation was made only for a single thrust coefficient at a fixed tip speed ratio corresponding to the below rated wind speed region.
Ishihara conducted a wind tunnel test using a 1/100 scale model of Mitsubishi’s MWT-1000 in 2004 and proposed a new wake model with a Gaussian wind profile by combining the measured wake data and the momentum conservation equation [
25]. He constructed a wake model through the wind tunnel test using a scaled wind turbine with a rotor diameter of 57 cm considering variations of the thrust coefficient and the ambient turbulence intensity. However, since the blade pitch angle of the wind turbine model was fixed, the tip speed ratio of the model had to be adjusted to change the thrust coefficient to cover the wake in the above rated wind speed region. This operating condition is different from that of general modern megawatt class wind turbines with active pitching.
Unlike kinematic wake models, field wake models are implicit models that are based on time or space averaged Navier-Stokes equations. They yield the magnitude of the velocity over the entire flow field as a solution of the implicit equation [
17]. Therefore, field wake models require more computing power than kinematic wake models. In recent years, the remarkable increase in computation speed has led many researchers to use computational fluid dynamics (CFD) to characterize the wake effect and develop wake reduction strategies [
30]. The CFD-based wake analysis is known to provide high accuracy because the wind turbine rotor is modeled and the flow field around the rotor is calculated with Navier–Stokes equations [
31]. The Energy Research Center of the Netherlands (ECN) developed a three dimensional wake model called WakeFarm based on the UPMWAKE proposed by Crespo et al. [
24]. The basic background flow of UPMWAKE is known to be modeled by the atmospheric wind profile model based on the Monin-Obukhov similarity theory [
32]. The initial WakeFarm used a hat-shaped initial wind profile similar to that of the Jensen wake model.
Schepers tried to improve the initial WakeFarm, in the Efficient Development of Offshore Wind Farms (ENDOW) project and proposed a Gaussian initial wind profile obtained from a wind tunnel test. He verified the model with the test data from the wind tunnel and from the field. Although the model was in good agreement with the measured data at the location where the initial profile was applied, the error increased with the downstream distance.
The Ainslie eddy viscosity wake model is one of the oldest field models and is a wake model based on parabolized Reynolds averaged Navier-Stokes (RANS) equations [
33]. Ainslie simplified the equation by assuming axisymmetric flow and by ignoring the pressure gradient, and constructed a model that can be solved relatively quickly and accurately. However, this caused the disadvantage that the model cannot be applied to the region close to the rotor where the pressure gradient is dominant.
Unlike WakeFarm, the Ainslie eddy viscosity model uses two empirical models such as the initial wind speed profile and the filter function to reduce the error of the wake prediction. The two empirical models are derived from the wind tunnel tests that take into account the change of the thrust coefficient and the downstream distance [
34]. However, since the wind turbine scale model used to construct the empirical wake model was not pitch regulated as the wind turbine model used by Ishihara, the tip speed ratio was adjusted to change the thrust coefficients of the wind turbine model.
In this study, a new wake model based on the Ainslie eddy viscosity wake model was investigated to improve the current weaknesses of the model. For this, a wind tunnel test was performed with a scaled wind turbine which has an active pitching mechanism similar to that of modern large capacity wind turbines.
Although research articles with a similar scope are available in the literature, the originality of this study can be summarized as follows. The first is the wind turbine scale model used in this study. It was simply a scaled wind turbine with pitch actuators. Therefore, the wake could be measured both below and above the rated wind speed regions to cover different thrust coefficients by pitching the blades. Adjusting the tip speed ratio in the below rated wind speed region only for wake measurement purposes was not implemented any more.
The second is the control of the scaled wind turbine. A Bachmann programmable logic controller (PLC) with torque and pitch control algorithms embedded was used as a controller for the scaled wind turbine. Therefore, the operating of the wind turbine while the wake was measured in the downstream of the turbine was very similar to that of modern large capacity wind turbines such as the maximum power point tracking in the below rated wind speed region and the closed loop power regulation not to exceed the rated value using a PI controller in the above rated wind speed region. The third is the modification of the initial wake profile and the filter function of the original the Ainslie eddy viscosity model to obtain a more accurate wake model that can be used for both below and above rated wind speed regions. To show the performance of the proposed wake model more clearly, the result was compared with those of the wake models proposed by other researchers having a similar scope.
4. Results and Discussion
All of the wake data presented below were measured at a constant air density. Also, the proposed wind profile is the normalized result using the mean wind speed measured by hot-wire probes 1.57 D ahead of the wind turbine.
4.1. Comparison with the Wake Model
The wind profile changes in the wake region were measured at a fixed thrust coefficient of the wind turbine.
Figure 5a shows the measured wake profile of the scaled wind turbine model operating at below rated wind speed and the comparison result of the predicted wake profiles through wake models.
Figure 5b shows measured wake profile of the above rated wind speed case.
As shown in
Figure 5a, the wind speed reduction due to the kinetic energy extraction was mainly centered at the mid-span (±0.25 D–0.3 D), which is critical for the aerodynamic performance of wind turbine blades [
42]. The wind speed reduction decreased away from the mid-span. This phenomenon resulted in a double-dip wind speed reduction, which was observed up to 3–4 D. Thus, it had a relatively dominant effect on pressure gradient. As the downstream distance increased, the wind speed reduction profile became parabolic and eventually disappeared. Since the wake profile measured in this wind tunnel test was measured by moving the hot-wire probe, the meandering of the wake and the centerline offset were not clearly confirmed in the measured data.
As shown in
Figure 5a, none of the wake models considered in this study could simulate the double-dip shape of the wake profile near the wind turbine rotor because they ignored the pressure gradient in the region close to the wind turbine. In particular, the Jensen wake model differed from the measured wind profile because it assumed the hat-shaped wind profile. However, other wake models which assumed a Gaussian wind profile such as the Larsen wake model, had profiles similar to the measured one in the far-wake region, where the pressure gradient disappeared.
The most remarkable point of this wind tunnel test is the accuracy of the wake model which uses the empirical model developed by a wind tunnel test or is developed through Large eddy simulation (LES) analysis. The Ainslie eddy viscosity wake model uses the empirical initial wake profile obtained from the wind tunnel test as mentioned above. The initial wake profile model is applied to the downstream 2 D, and then the wake calculation based on the parabolized RANS is performed. For this reason, it can be seen that the wake profile at the 2 D downstream from the turbine is very similar to the measured wake profile. However, the wake prediction accuracy of the Ainslie eddy viscosity wake model decreases gradually as the wake propagates downstream. This wake prediction error occurs because the Ainslie eddy viscosity wake model cannot sufficiently simulate the physical phenomena existing in the wake region.
The BP wake model based on the wind tunnel test and the LES analysis predicts the wind speed deficit at the wake center line very accurately. The BP wake model, however, showed a high accuracy only in the wake centerline, and overestimated the width of the wake. On the other hand, the Ishihara wake model developed through the wind tunnel test showed a very similar wake profile to the actual measured wake from 3.66 D, where the double dip phenomenon disappears. The Ishihara wake model showed very accurate prediction in the wake deficit at the wake center line as well as the wake width. On the other hand, if the thrust coefficient of the wind turbine is reduced as shown in
Figure 5b, double dip phenomenon does not occur even if it is close to the region. In addition, it can be found that the wake width expansion and the wind speed recovery are similar to the below rated case in accordance with the increase in the downstream distance.
Figure 6 shows the comparison of the mean wind speed deficit from the wake models and the wind tunnel test. The mean deficit was calculated in the wake region bounded by
U/
U0 < 0.95 which is the boundary used to calculate the mean deficit in Parkin’s study [
43].
As shown in
Figure 6a, most of the wake models assuming the Gaussian wake, profile including the recently developed Ishihara wake model, predict the average wind speed deficit and recovery in the wake region and are very similar to the actual measured wake when input wind speed is lower than the rated wind speed. On the other hand, the Jensen wake model and the BP wake model predicted excessive wind speed deficit compared to the actually measured wake, as shown
Figure 6a. The comparison results show that a wake model can improve its accuracy if it is developed or modified through full size wind turbine simulations or wind tunnel tests with scaled wind turbines.
Figure 6b shows the comparison of the measured wind speed deficit with the actual measured wake data when the input wind speed is higher than the rated wind speed. As shown in
Figure 6a, most of the wake models considered in this study showed relatively high accuracy when the input wind speed was lower than the rated wind speed. However, if the thrust coefficient of the wind turbine changed, the accuracy of the wake prediction was decreased as shown in
Figure 6b.
The Ainslie eddy viscosity wake model considered in this study uses the parabolized RANS and the empirical model at the same time to predict the wake, but it is not much different in the mean deficit errors from the other wake models. In particular, the error of the Ainslie eddy viscosity wake model was larger than that of the recently developed Ishihara wake model.
4.2. Wake Model Modification
The Ainslie eddy viscosity wake model calculates the wind speed variation in the wake region by using an initial profile, which is originally given by
where
DM and
b are calculated as follows:
where
I0 is the ambient turbulence intensity in percentage. The original Ainslie eddy viscosity wake model uses Equation (23) to calculate the initial wind speed profile which is applied to 2 D downstream. In this study, the initial profile and filter function were modified based on the measured wake profile and the Nelder Mead simplex algorithm. The Nelder Mead simplex method is one of the function minimization algorithms used to find the parameters that can minimize the value of the objective function. The detailed operation principle of the algorithm can be found in the literature [
44,
45].
The objective function of the algorithm to calculate the initial wind speed profile is shown in Equation (26). The initial profile of the wake was calculated at 3.66 D behind the wind turbine, which is the location where the double dip shape in the wake profile disappeared:
In Equation (26),
J is the objective function used in the function minimization algorithm,
UExp3.66 is the measured wake profile at 3.66 D and
UModel is the calculated wind profile.
UModel can be obtained by Equation (27)
where
P1–
P6 are the dimensionless parameters obtained from the function minimization algorithm to match the actual wake profile.
Table 2 shows the parameters obtained for the initial wake profile.
Similarly, the following objective function was used to construct the filter function:
Note that
J in Equation (28) is the objective function used in the function minimization algorithm to calculate the parameters of the modified filter function. Note that it is not the same as Equation (26).
J1 and
J2 in Equation (28) can be calculated by Equations (29) and (30).
where
δ is the wake width. The Ainslie eddy viscosity wake model does not have a separate formula for the wake width unlike the general analytical wake model, but the wake width can be modified by modifying Equation (22) which was derived from Ainslie’s wind tunnel test. Therefore, in this study, the filter function in Equation (22) was modified to be
where
P1,
P2, and
P3 were found from the function minimization algorithm to match the actual wake profile.
Table 3 presents the parameters obtained from the function minimization algorithm.
Figure 7 shows the comparison of the wakes measured by the wind tunnel test and predicted through the wake models.
Figure 7 shows the wake prediction accuracy of the modified Ainslie eddy viscosity wake model with respect to the wind turbine thrust coefficient variation.
Figure 7a,b shows the mean wind speed deficit in the wake region by distance for the below and above wind speed regions, respectively. As shown in
Figure 7a, the modified Ainslie eddy viscosity wake model showed higher prediction accuracy than other classical kinematic wake models. In particular, the modified Ainslie eddy viscosity showed high accuracy not only at 3.66 D, which is the location where the empirical initial profile was applied, but also in the subsequent locations where the wake predictions were made based on the parabolized RANS.
Figure 7b also shows that the modified Ainslie eddy viscosity wake model has high accuracy even when the thrust coefficient of the wind turbine is decreased.
As shown in
Figure 8a,b, the Ainslie eddy viscosity wake model with the modified initial wind speed profile and the filter function predicted wakes very close to the measured wake data. In particular, the modified Ainslie eddy viscosity predicts the mean wind speed deficit very close to the actual measured wake data even when the thrust coefficient of the wind turbine changes (above rated wind speed region).
Figure 8c,d shows that the modified Ainslie eddy viscosity wake model improves the accuracy of the wake width as well as the wind speed mean deficit in the wake region.
Figure 8c shows the predicted wake width for each wake model when the input wind speed is lower than the rated wind speed. As shown in the
Figure 8, most of the wake models have a wider wake region than the actual wake width. However, the Ishihara wake model predicts a wake width very close to that obtained from the wind tunnel tests in this study. On the other hand, the BP wake model based on the LES analysis showed the largest wake width prediction error among the wake models considered in this study. This error could be attributed to the fact that the considered wind turbine model when developing the BP wake model was not realistic in terms of the tip speed ratio.
Figure 8d shows the result of comparing the wake width generated when the input wind speed is higher than the rated wind speed, compared with the measured data and the wake predicted by wake models. As shown in
Figure 8c,d, the modified Ainslie eddy viscosity wake model as well as the Ishihara wake model have a relatively accurate wake width prediction even though the thrust coefficient of the wind turbine is changed. On the other hand, other wake models mostly predicted a wider wake width than the measured value. The Jensen wake model, which predicted the wake width accurately when the input wind speed was lower than the rated wind speed, showed good accuracy even when the thrust coefficient of the wind turbine was changed (above rated wind speed region).
Table 4 shows the error of the wake prediction accuracy of the wake models considered in this study, including the newly proposed modified Ainslie eddy viscosity wake model. The accuracy of the wind speed deficit is presented by the root mean square error (RMSE) and wake width error was not considered for RMSE calculation. In this study, accurate prediction accuracy of the wake shape and wind speed deficit of each wake model is confirmed through RMSE comparison.
Ainslie eddy viscosity wake model modified by this study starts wake calculation from 3.66 D where double dip phenomenon disappears. Therefore, the wake prediction based on the wake model is possible after 3.66 D. In
Table 4, the section labeled ‘N/A’ denotes the point where wake prediction using the wake model is impossible.
As shown in
Table 4, the wind tunnel test shows that the wake model with the highest wake prediction accuracy is different for each measurement distance. The Jensen wake model considered in this study is a classical wake model with a top hat profile and predicts the wind speed profile in the wake region most unrealistically. For this reason, the Jensen wake model has the largest RMSE in all sections irrespective of the thrust coefficient change and distance from the wind turbine.
The Larsen wake model predicts the wind speed profile in the region close to the rotor based on the LES analysis. For this reason, it can be seen that the Larsen wake model has a relatively small error in the region close to the rotor, regardless of the thrust coefficient change of the wind turbine. However, as the distance increases, the Larsen wake model has higher errors than the other wake models.
As mentioned above, the Ainslie eddy viscosity wake model uses the wind speed initial profile derived from the wind tunnel test for the wake calculation. For this reason, it is confirmed that the wake prediction accuracy is higher than other wake models at the point where the initial wind speed profile is applied. However, the wake prediction accuracy is reduced compared to other wake models in the section where wake calculation is performed based on simplified parabolic RANS. This error can be thought to be caused by the inability to simulate the turbulence momentum change in the wake region in the process of simplifying the parabolized RANS.
The BP wake model and the Ishihara wake model are wind tunnel test-based wake models developed to predict the wake in the far wake region. However, although the two wake models are based on the scaled wind turbine model, unlike the Ishihara wake model, the BP wake model was developed based on the scaled wind turbine model which has a much lower tip speed ratio than the actual wind turbines. As a result, the wake prediction accuracy was relatively low compared to the Ishihara and the modified Ainslie eddy viscosity wake models in the below and above rated wind speed regions.
The Modified Ainslie eddy viscosity wake model showed similar wake prediction accuracy to the Ishihara wake model when the input wind speed was lower than the rated wind speed. However, it showed the best wake prediction accuracy in terms of RMSE among the wake models when the wind speed was higher than the rated wind speed in this study.
5. Conclusions
A modified Ainslie eddy viscosity wake model was proposed in this study. The initial wind speed profile and filter function of the original wake model were modified by wind tunnel test results. The modification of the wake model was performed using a scaled wind turbine model designed to have a tip speed ratio similar to that of a modern megawatt wind turbine. The scaled wind turbine model used in the wind tunnel test is controlled by a PLC with a C code to achieve the power control which is similar to the modern megawatt wind turbines. The wake measurement was performed using the hot-wire anemometry and a traversing system to obtain the two-dimensional wake field at the hub height of the scaled wind turbine model.
The wake prediction accuracy of the wake model was verified by comparing the mean wind speed deficit and the wake width. By comparing the measured data with the predictions from the wake models, it was found that the Larsen wake model has the highest wake prediction accuracy in the relatively close range (within 3.14 D). Such a wake prediction accuracy is considered due to the second term applied to the Larsen wake model which was obtained from CFD simulation. However, in the Larsen wake model, there is no term to consider the influence of the momentum change of the turbulence in the wake region on the wake. Therefore, it was found that the wake prediction accuracy decreases as the distance in the downstream direction increases.
For the original Ainslie eddy viscosity wake model, it was found that the wake prediction near the location where the empirical equation is applied is more accurate than the predictions from the other wake models. However, the original Ainslie eddy viscosity wake model is based on simplified parabolic RANS, and the wake prediction error increased as the downstream distance increased, as in the Larsen wake model.
Most wake models showed similar wake prediction accuracy in the region after 3.66 D. Among the wake models considered in this study, the proposed modified Ainslie eddy viscosity wake model and the Ishihara wake model showed the highest prediction accuracy. The BP wake model developed based on the LES results showed a better prediction accuracy compared with the other kinematic wake models, but was found to be less accurate compared with the modified Ainslie eddy viscosity or the insufficient result compared with the Ishihara wake model. The modified Ainslie eddy viscosity wake model showed the best wake prediction accuracy when the input wind speed was higher than the rated wind speed, that is, when the thrust coefficient of the wind turbine was reduced.
This wake prediction error, especially in the above rated wind speed region, can be considered to be the difference between the scaled wind turbine model used in this study and by others although further studies are required to clearly understand the cause. In addition, this study has limitations because ambient turbulent intensity and the thrust coefficients considered for the wake model development and validation were very limited. Therefore, to generalize the result, more development and validation must be carried in various conditions including the ambient turbulence intensity, thrust coefficients, and the yaw orientation of the upstream wind turbine.