1. Introduction
With the increase in the severity of the energy crisis and environmental pollution, wind energy, as one of the most commercialized forms of energy production, has developed rapidly. The height of turbines and the size of their blades have been increased to improve the utilization rate of wind energy. However, the wake effect between turbines reduces the wind speed and increases the turbulence at the downstream turbine, lowering its power generation efficiency and increasing its fatigue load [
1,
2]. Research has demonstrated that the wake distribution of wind turbines is mainly affected by the atmospheric boundary layer, wind turbine type, wind turbine yaw angle, and thrust coefficient [
3]. Some studies have shown that the yaw operation of the upstream turbine can change the distribution area of the wake, reducing the impact of the wake on the downstream turbine and significantly improving its output power. The wake effects of turbines have become an important research topic [
4].
Due to their low cost and high flexibility, wind tunnel tests are often used to study the wake characteristics of wind turbines. The wakes of a horizontal-axis wind turbine (HAWT) were extensively studied and comprehensively reviewed in [
5]. The researchers experimentally investigated the distribution characteristics of turbine wakes under operating conditions through wind tunnel tests using small-scale models. The experimental results facilitated the understanding of the wake evolutionary characteristics, as well as the validation and development of numerical models (such as large-eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS)) and analytical wake models [
6,
7]. Currently, most studies focus on the geometrical and kinematic similarities of scaled wind turbine models in wind tunnel tests, and there is a lack of comprehensive studies on the dynamic similarity. Thus, using wind tunnel tests to design scaled models that demonstrate wake characteristics with a reasonable degree of similarity to those of the corresponding prototypes is important. Generally speaking, to ensure that the wake distribution of a small-scale model wind turbine, as measured by a wind tunnel experiment, reflects that of a full-scale wind turbine, the model must be accurately scaled and designed to comply with the similarity principle. For wind turbine blades, the real wake characteristics can only be simulated when the tip speed ratio and Reynolds number of the model are the same as those of the full-scale wind turbine. In wind tunnel tests, it is easy to meet the requirements of the blade tip speed ratio. The authors of [
8] noted that if the blade tip speed ratio of the full-scale wind turbine is matched, the blade tip vortex of the wind turbine model presents the same spiral structure.
However, when the wind turbine rotates, the aerodynamic phenomenon of the blade depends on the Reynolds number of the flow. As it is usually limited by the scale ratio in a wind tunnel experiment, the Reynolds number of most scaled wind turbine models is generally in the range of Re ≈ 10
4. Researchers have demonstrated that the lift, drag, and stall behaviors of thick airfoils in wind turbines are highly sensitive when the Reynolds number is in the range of Re ≈ 10
4, while thick airfoils in geometrically similar models perform poorly in this low-Reynolds-number flow state, resulting in lower aerodynamic performance [
9]. The Reynolds number of a full-scale wind turbine is 2~3 orders of magnitude higher. It is impossible to maintain the Reynolds number of a full-scale wind turbine in a wind tunnel experiment, because the high speed required to achieve Reynolds similarity would compromise the similarity of the blade tip speed ratio between the full-scale wind turbine and the scaled model.
The lift, drag, thrust, and torque are the main aerodynamic parameters of an airfoil and the key factors affecting the operation of a wind turbine. The authors of [
10] showed that the wake expansion of a wind turbine model is related to the thrust coefficient. Scholars have also carried out research to determine the specific relationship between the thrust of a turbine model blade and the change in the Reynolds number. The authors of [
11] developed a 1:50 scale model wind turbine to simulate the Froude number needed to analyze the overall motion and response of offshore floating wind turbines. The model used a flat-plate airfoil (naca-ag24 airfoil) with a low Reynolds number that was able to produce a thrust coefficient equivalent to that of the full-scale wind turbine, thus achieving the experimental purpose. Unfortunately, this study did not investigate the wake characteristics. In developing a wake model for a single wind turbine, the authors of [
12] matched the thrust coefficient of a full-scale 2 MW wind turbine by adjusting the pitch angle of the blade with a low Reynolds number. It remains unclear whether a thrust coefficient as high as that of a large-scale wind turbine is suitable for a wind tunnel test. In a previous study, it was shown that, although the limitations of the wind tunnel experiment meant that the Reynolds number failed to match that of the full-scale turbine, the scaled model could reproduce the blade tip vortex, root vortex, and wake rotations observed in the wake development of the full-scale turbine [
13]. The above research shows that, in order to reproduce the turbine wake in a wind tunnel, the Reynolds number similarity can be relaxed to some extent, as long as the thrust of the model wind turbine is consistent with that of the full-scale wind turbine.
In [
14], it was suggested that, in contrast to the near wake region, the flow characteristics in the far wake region are greatly weakened by the aerodynamic characteristics of the wind turbine blades and are transformed into environmental turbulence, which plays a major role in the wake flow in this region. Generally, it is believed that there are three main turbulence sources in the far wake: atmospheric turbulence, caused by surface roughness and surface thermal convection; mechanical turbulence, caused by the wind turbine and tower; and wake turbulence, caused by tip vortex breaking. For instance, by comparing the particle image velocimetry (PIV) data of a scaled model with those of a full-scale turbine, one study came to the interesting conclusion that the velocity recovery rate in the far wake region mainly depends on atmospheric turbulence [
15]. This conclusion was also confirmed by two experimental investigations [
16,
17]. It is commonly recognized that the higher the atmospheric turbulence intensity, the faster the wake velocity deficit recovery. The reason for this is that the mechanical turbulence caused by the nacelle and tower of the wind turbine is greatly weakened. This study will mainly focus on the effects of the oncoming flow conditions and the aerodynamic performance of the blade airfoil based on dynamic similarity.
For a large-scale wind turbine, the aerodynamic characteristics of an airfoil geometrically scaled from a high-Reynolds-number full-scale wind turbine airfoil are reduced in the wind tunnel test, resulting in a low Reynolds number, which causes the obvious differences in wake distribution between the scaled wind turbine model and the full-scale turbine. The purpose of this study was twofold: firstly, to clarify the influence of dynamic similarity on wake modeling using a scaled wind turbine in a wind tunnel test; and secondly, to confirm the necessity of achieving thrust similarity in a wind tunnel test of turbine wake characteristics. In this study, a wind turbine model with thrust-optimized blades was designed to mimic the thrust coefficient of a full-scale wind turbine by varying the blade airfoil, chord length, and twist angle. We elucidated the problem of the insufficient thrust of geometrically similar blades through experimental measurements and the analysis of aerodynamic principles. Additionally, the specific aerodynamic performance parameters that most substantially affected the wake characteristics were quantitatively identified. We carried out a systematic study to compare the differences in wake characteristics between the wind turbine model designed according to geometric similarity and motion similarity (i.e., all characteristic dimensions of the wind turbine model were consistent with the full-scale wind turbine) and the wind turbine model designed according to thrust similarity.
2. Experimental Setup
The experiment was conducted in the open-type wind tunnel laboratory of Chongqing University (CQU). The test area was 15.12 m (length) × 2.4 m (width) × 1.8 m (height), and the maximum wind speed of the wind tunnel was 35 m·s
−1. Passive devices including baffles, spires, and roughness cells were implemented to generate the atmospheric boundary layer (ABL) flow, as shown in
Figure 1a. The inflow conditions are shown in
Table 1. The power law adopted for the upstream flow conditions was as follows:
where
is the average wind speed at the hub height
;
is the average wind speed at the height above ground z; and α is the exponential value of terrain exposure type A, stipulated in the wind load code of China. The experimental tests in this study were conducted to investigate the effects of the tip speed ratio and yaw angle on the aeromechanical performance and wake evolution of the wind turbine. The scale ratio of the wind turbine model was 1:250, it had a diameter of 608 mm (D) and a hub height of 550 mm (
), and the rotor had three blades, as shown in
Figure 1b. A comparison of the parameters between the scaled wind turbine model and the full-scale turbine is presented in
Table 2. The model was installed in the center of the turntable, and its position was about 8 m downstream from the inlet of the test area.
Figure 2 shows the profiles of the mean wind speed and turbulence intensity and the wind velocity spectrum simulated in the wind tunnel tests. Additionally, we measured the flow characteristics at the wind turbine model installation position (2
D) and behind the wind turbine model (4
D) to check the horizontal homogeneity of the upstream flows in the test region. The maximum deviation of these velocity profiles from the target ABL flow was no higher than 3%, demonstrating that the inflow had fully developed before reaching the wind turbine. Meanwhile, the blockage ratio of the wind turbine model in the wind tunnel test was about 4.47%. It has been found that when the blockage ratio of an experimental model is less than 10%, the wake expansion is similar to that in free flow [
9]. Therefore, the influence of the side wall and ceiling of the wind tunnel on the wind turbine wake could be ignored when ensuring the full development of the wake in the test area.
The geometrically similar wind turbine model was designed by directly scaling a working 5 MW wind turbine and maintaining strict geometric similarity, while the airfoil of the thrust-optimized wind turbine model was redesigned using the optimization methodology proposed in our previous research. The redesigned blade airfoils could appropriately mimic the thrust performance under low-Reynolds-number conditions. The geometric dimensions of the wind turbine tower and nacelle were similar to those of the prototype. Both types of blade were made of aluminum alloy by three-dimensional printing, and the surface was polished and painted to ensure smoothness. The rotational frequency of the turbine model was controlled by the servo motor placed inside the turbine nacelle to ensure kinematic similarity, which meant that the blade tip speed ratio (
) of the turbine model was consistent with that of the full-scale wind turbine. Under the incoming flow conditions of 6 m·s
−1, as shown in
Table 1, the rotary speed of the wind turbine was 1942 rpm, while under the incoming flow conditions of 9 m·s
−1, the rotor speed of the wind turbine was 2835 rpm. The full-scale wind turbine usually operates at a blade tip speed ratio ranging from 6 to 8. Therefore, our model was designed to operate at the optimum tip speed ratio. The Reynolds number based on the blade tip velocity and chord length was defined by Equation (2). The total length of a 5 WM full-scale wind turbine blade is 73.5 m, and an airfoil with a thickness-to-chord ratio of 23.95 and radial lengths ranging from 23 m to 53 m is used. The airfoil Reynolds number (Rec) of the geometrically similar blades was 3.1 × 10
4, while that of the thrust-optimized blades was 4 × 10
4, revealing that the Reynolds number of the blade airfoil was two orders of magnitude lower than that of the full-scale wind turbine operating under atmospheric flow.
where
is the tip velocity,
is the chord length, and
is the kinematic viscosity of the air.
The tower of the wind turbine model was made of a stainless-steel tube to ensure its stiffness and stability under the high wind and rotary speeds. The model was almost rigid, so the influence of wind turbine vibration on the wake could be ignored. In the present study, a high-accuracy multi-axis force/torque sensor (ATI-DELTA) was mounted at the tower base of the wind turbine model to measure the dynamic wind loads acting on the wind turbine model. This sensor can provide instantaneous measurements of the aerodynamic forces and moment (torque) along all three axes. The measurement range of the ATI-DELTA multi-axis force/torque sensor is 165 N for
x-axis forces and 15 N·m for torque, as shown in
Table 3. During the experiments, the sampling rate of the wind load data was 1 kHz. For each test case, the sampling time was set to 60 s to obtain sufficient samples of thrust forces. The force/torque sensor was deliberately calibrated before the wind tunnel tests, and the accuracy of the load balance was verified within ±1%. The drag force of the tower was excluded from the measurements by separately testing the parked wind turbine model; thus, the net thrust forces acting on the swept area of the turbine blades could be deducted.
A TFI Cobra Probe anemometry system was installed in the three-dimensional traversing mechanism in the wind tunnel to measure the turbulence in the far wake regions 4
D, 5
D, and 6
D behind the wind turbine model. The experimental setup is shown in
Figure 3, in which the streamwise direction is the
x-axis, the spanwise direction is the
y-axis, and the vertical direction is the
z-axis. By using this Cobra Probe, all three components of the velocity vector could be measured instantaneously. During the experiments, the maximum flow velocity measured in the wind tunnel was approximately 10 m·s
−1. Both a pitot-static tube and a temperature probe were used to calibrate the inflow velocity at the hub height. The accuracies of all the other velocity values measured in the present study were within ±0.15 m·s
−1. The sampling rate of the instantaneous velocity measurement was 500 Hz, with a sampling period of 60 s at each point of interest.
5. Conclusions
In this paper, the wake characteristics of geometrically similar and thrust-optimized models, designed by considering only geometrical similarity and motion similarity, respectively, were investigated in depth. Subsequently, we carried out further study of the wake characteristics of the thrust-optimized model. The major observations, results, and conclusions are listed below:
The thrust coefficient of the geometrically similar model was only 35%–40% of that of the full-scale wind turbine. The wake characteristics of the geometrically similar model differed significantly from the thrust-optimized model in terms of the normalized mean velocity deficit and the recovery rate. This indicates that obvious deviations in reproducing the wake characteristics of a full-scale wind turbine occur if only geometric and kinematic similarity are considered.
The wake distribution of the thrust-optimized model showed good agreement with the Gaussian wake model in the far wake region. Furthermore, it was found that the height of the center of the wake’s minimum velocity in the wind tunnel test decreased gradually with increasing downstream distance; a similar lowering of the wake center was observed for changes in the turbulence intensity. Furthermore, it is worth mentioning that the turbulence intensity distribution of the wake became asymmetric when the wind turbine was under yawed conditions. This led to an uneven recovery rate on both sides of the wake and an asymmetric wake velocity deficit profile in the spreading direction.
The normalized mean velocity profile and turbulence intensity profiles under the two inflow velocity conditions were identical, and the results for self-similarity, recovery rate, and wake distribution range were consistent. This indicates that the variation in the velocity of the inflow was not a key factor affecting the wake development.
Compared with the wake development under non-yawed conditions, the wakes under 10°, 20°, and 30° yaw angles were significantly deflected, as expected. The larger the yaw angle, the more obvious the deflection effect. The maximum lateral deflection distance downstream was close to 0.7D. The variation in the lateral normalized mean velocity profile according to the downstream distance for the thrust-optimized model under yawed conditions was verified.
The results of this study show that the problem of the mismatched performance between subscaled model blades and prototype blades was resolved. This experimental study combined theoretical analysis and experimental data in an attempt to solve the difficulties of aerodynamic evaluation in wind tunnel experiments involving geometrically similar models designed based on geometric and kinematic similarity with full-scale wind turbines. We revealed the mechanisms of the aerodynamic performance differences between geometrically similar models and full-scale wind turbines. Addtionally, we proposed a method for improving the aerodynamic performance of blades by using optimized airfoil shapes applicable to low-Reynolds-number conditions. The outcomes of this study could provide a useful reference and reasonable guidance for future wind tunnel tests of turbine wakes and supply the benchmark data for validating and verifying numerical simulations of wind turbine wakes.