Next Article in Journal
Performance Analysis of Heat Pump Dryer with Unit-Room in Cold Climate Regions
Next Article in Special Issue
Robust LFC Strategy for Wind Integrated Time-Delay Power System Using EID Compensation
Previous Article in Journal
Marginal Agricultural Land Low-Input Systems for Biomass Production
Previous Article in Special Issue
Fault-Ride Trough Validation of IEC 61400-27-1 Type 3 and Type 4 Models of Different Wind Turbine Manufacturers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigations on the Unsteady Aerodynamic Characteristics of a Horizontal-Axis Wind Turbine during Dynamic Yaw Processes

1
Key Laboratory of Condition Monitoring and Control for Power Plant Equipment, Ministry of Education, North China Electric Power University, Beijing 102206, China
2
China Huaneng Clean Energy Research Institute, Beijing 102209, China
3
Department of Aerospace Engineering, Iowa State University, Ames, IA 50011, USA
*
Author to whom correspondence should be addressed.
Energies 2019, 12(16), 3124; https://doi.org/10.3390/en12163124
Submission received: 11 June 2019 / Revised: 29 July 2019 / Accepted: 7 August 2019 / Published: 14 August 2019
(This article belongs to the Special Issue Modeling of Wind Turbines and Wind Farms)

Abstract

:
Wind turbines inevitably experience yawed flows, resulting in fluctuations of the angle of attack (AOA) of airfoils, which can considerably impact the aerodynamic characteristics of the turbine blades. In this paper, a horizontal-axis wind turbine (HAWT) was modeled using a structured grid with multiple blocks. Then, the aerodynamic characteristics of the wind turbine were investigated under static and dynamic yawed conditions using the Unsteady Reynolds Averaged Navier-Stokes (URANS) method. In addition, start-stop yawing rotations at two different velocities were studied. The results suggest that AOA fluctuation under yawing conditions is caused by two separate effects: blade advancing & retreating and upwind & downwind yawing. At a positive yaw angle, the blade advancing & retreating effect causes a maximum AOA at an azimuth angle of 0°. Moreover, the effect is more dominant in inboard airfoils compared to outboard airfoils. The upwind & downwind yawing effect occurs when the wind turbine experiences dynamic yawing motion. The effect increases the AOA when the blade is yawing upwind and vice versa. The phenomena become more dominant with the increase of yawing rate. The torque of the blade in the forward yawing condition is much higher than in backward yawing, owing to the reversal of the yaw velocity.

1. Introduction

The use of wind power has grown rapidly over the past few decades. The aerodynamics of the wind turbine are a core subject of wind power generation. Owing to the intermittency of wind, the orientation of a wind turbine must be frequently adjusted, typically using the yaw mechanism, to ensure the rotor disk is perpendicular to the direction of the wind. Yaw control is one of the most important strategies for improving the efficiency of wind energy conversion. Under the yawed condition, velocity component of the upstream wind is parallel to the plane of rotation. Thus, the angle of attack of airfoil along blade spanwise exhibits periodic variation in one revolution [1]. The aerodynamic performance of a wind turbine under yaw, as well as the surrounding flow, show highly unsteady characteristics, resulting in high fatigue loads that can dramatically reduce the operational life of a wind turbine. Therefore, understanding the unsteady aerodynamic characteristics of wind turbines under yaw is important for improving their design [2]. In particular, accurate predictions of unsteady aerodynamic loads under yawed conditions are important for blade design and operational control.
Knowledge of complex three-dimensional (3D) flow phenomena over wind turbine blades and near wake under the yawing condition are still lacking. Thus, improved numerical models that consider the combined effects of stall delay, unsteady flow, and other significant effects such as the retreating & advancing effect and upwind & downwind effect are needed. The aim of present paper is to further elucidate the aerodynamic and 3D flow characteristics of yawed and yawing turbines.
The aerodynamics of wind turbines under yaw have been extensively investigated using both experiments and numerical simulations. Detailed aerodynamic performance measurements along the blade in combination with 3D flow field measurements have helped further our understanding of the flow mechanism. Wind tunnel experiments performed by Hand et al. [3] on a two-blade wind turbine provide a benchmark for numerical simulations. In their study, Hand and colleagues obtained the pressure and aerodynamic force coefficients for an azimuthal distribution at different radial locations under five yawed cases: 0°, 10°, 30°, 45°, and 45°. The angle of attack was found to vary approximately with the blade azimuth angle according to a sinusoidal function. Further to this, Schepers et al. [4] performed wind tunnel experiments on a three-bladed wind turbine, the MEXICO rotor, under three different yaw angles: 15°, 30°, and 45°. They obtained sectional aerodynamic forces and torques for different azimuth angles; however, the sinusoidal variation rule did not correlate well along the inner span of the blade owing to the velocities induced by root vortices. Sant et al. [5] investigated the TUDelft reference rotor in an open jet wind tunnel under both axial and yawed conditions and assessed the blade element moment (BEM) method, which was found to be limited to the skewed wake under yawed conditions.
Based on the work of Schepers et al. [6], Micallef et al., performed stereo particle image velocimetry measurements to study the blade aerodynamic performance and near wake development. Under yaw, flow in the windward region exhibited inboard motion due to rapid motion of the tip vortices away from the blade. In field tests, Ven et al. [7] applied in-field measurements of a turbine to investigate the three dimensional stall phenomena under yawed inflow conditions. Imbalances in the crossflow fraction along the radial direction of the blade were observed between azimuth angles of 90° and 270°. The phenomenon is caused by yawed inflow leading to in-wash and outwash on the blade surface. Furthermore, Dai et al. [8] provided a detailed investigation of the yaw effect based on Supervisory Control and Data Acquisition (SCADA), and presented characteristics of the power coefficient and rotor torque under yaw. Interestingly, many yaw error control solutions based on the smart prediction of the wind direction had been used to maximum power extraction of wind power [9,10]. Further to this, Saenz-Aguirre et al. [11] claim that the drawbacks of the yaw control are the very large time constant and the strict yaw angle rate owing to the high mechanical loads, resulting that yaw control is less studied in the literature than that pitch and speed control. Furthermore, Munters et al. [12] performed wind-farm control research with a new dynamic yaw control based on the large-eddy simulation, and found that with the consideration of the unsteady-control interacting with flow dynamics, dynamic yaw control for the given wind farm setup is better than induction control in the process of power extraction. In fact, the yawing operation of wind turbine is very common, and the aerodynamic characteristics of wind turbine under yaw had been extensively studied in the literature with numerical simulation.
Numerical methods to evaluate the aerodynamic performance of wind turbines mainly include BEM, vortex theory, and computational fluid dynamics (CFD). The most commonly used method is BEM, which is fast with low computational costs. Ryu et al. [13] developed a successive under-relaxation technique based on the classical BEM to compute specific radial locations of the blade and found that yawed inflow leads to periodic changes in the angle of attack on different spanwise sections. Jimenez et al. [14] firstly established an analytical wake model with a simple formula to predicit the wake skew effect based on the momentum conservation and top-hat model offered by Jensen [15] for the velocity deficit effect. However, the experimental validation was not sufficient as mentioned by the author and the model overestimate the wake deflection of wind turbine, owing to the inaccurate top-hat model for velocity deficit. With the goal of optimize the yaw angles of the wind turbines, Gebraad et al. [16] presented a new flow redirection and induction in steady-state model based on internal parametric model for wake effects, which demonstrate the optimization control for increasing the energy production of the wind plant. Recently, Bastankhah and Porte-Agel [17] developed a new analytical model on the assumption of the Gaussian distribution to predict the wake deflection and far wake velocity deficit, which showed a good agreement with the experimental data. However, some parameters of that model have not specified which is needed to be researched. Based on the works of them, Qian Guo-Wei et al. [18] presented a new analytical wake model with Gaussian distribution function of velocity deficit and momentum conservation in the lateral direction, considering the ambient turbulence intensity, thrust coefficient and yawed effect, which enables the model to have good application under various conditions. Yaw not only influences the wake flow behind the wind turbine, but also flows on the surface of the blade and 3D aerodynamic characteristics of the rotor [19]. Castellani [20] analyzed a small three-blade wind turbine and investigated the effects of large yaw angles, using numerical and experimental methods, and found that the 3P blade frequency is clearly visible. In addition, non-uniform flow around the blades was modeled using a simplified wake model modified by Prandtl’s tip loss factor.
The time constant used in dynamic inflow models results in an unrealistic drop at the tip region, and therefore, cannot accurately model for the dynamic conditions of wind turbines. For instance, both the tip loss model and dynamic inflow model have several shortcomings when used for non-design load cases but can possibly be improved by free vortex methods (FVM), which combine the blade aerodynamics with wake flow computational models. Qiu et al. [21] investigated the dynamic variation of loads on a wind turbine blade during the yaw process using an improved lift line method and proposed a wake model comprised of the vortex sheet model and tip vortex model. Further, the yaw rate and dynamic wake were shown to significantly influence the shaft thrust and torque of the wind turbine. Micallef et al. [22] performed both experimental and numerical investigations based on an FVM to investigate the tip vortex generation of a horizontal-axis wind turbine (HAWT) under yaw. Vorticity on the suction side under non-yawed conditions was shown to be more concentrated than for the yawed condition, in which vorticity spreads over a small region at the tip of the turbine blade.
Several aerodynamic correction models have been proposed and widely used, including Pitt and Peter’s yaw correction [23,24], Suzuki’s dynamic inflow model [25], Prandtl’s tip loss function [26], Du and Selig’s 3D stall delayed model [27], and Buhl’s wake correction [28]. Nonetheless, classical BEM and FWM are still not sufficiently reliable for predicting the aerodynamic load distributions on wind turbine blades, especially for the yawed and stalled rotor conditions. Investigations were performed on the complicated yaw aerodynamic problem of wind turbines using more accurate CFD simulations using detailed visualization of the 3D flow [29]. Tongchitpakdee et al. [30] investigated the aerodynamic characteristics of the National Renewable Energy Laboratory (NREL) Phase VI wind turbine modeled by an unstructured grid under different wind speeds and yaw angles using the BL model. At low wind speeds (<7 m/s) and under low yaw angles, the flow remains attached over most regions of the rotor. Further to this, Yu et al. [31] performed time-accurate aerodynamic simulations of the NREL Phase VI wind turbine under yawed flow conditions, which was also modeled with an unstructured mesh, along with the k-ω shear stress transport (SST) turbulent model. They found that the retreating side shown aerodynamic loads (Cn, Ct) higher magnitudes of loads than on the advancing side. However, the three-dimensional unsteady stall effect at different yaw angles was not studied. Schulz et al. [32] modeled and analyzed a generic 2.4-MW wind turbine with over-set structured grids and yaw angles ranging from −50° to 50° using the FLOWer CFD solver. The simulation was performed using a detached eddy simulation method and revealed azimuthal non-uniformity of the load variation along spanwise sections of the blade due to the retreating & advancing effects of the blade. Furthermore, the deformation characteristics of the blade were considered in the yawed case of the wind turbine simulation. Jeong et al. [33] investigated the effects of yaw errors on the aerodynamic and aeroelastic behaviors of NREL 5-MW HAWT blades, and showed how yaw misalignment adversely affects the dynamic aeroelastic stability of the blade. Dai et al. [34] performed an aeroelastic analysis on the Tjæreborg wind turbine under yaw, and found that fluid-solid coupling results in higher averaged power and thrust, as well as violent oscillation amplitudes.
Yaw is a continuous rotation process with a constant or variable velocity of rotation. Qiu et al. [21] investigated the shaft torque of the blade and instantaneous wake flow under the dynamic yawing process with a yaw rate between 5 and 20°/s for the NREL Phase VI wind turbine. For larger wind turbines, such as the NREL 5-MW et al., these yaw rate seem to be impossible owing to the larger moment of inertia of the rotor-nacelle subsystem due to enlargement of the blade and tower length. Leble et al. [35] carried out a series of CFD simulations to investigate the performance of the DTU 10-MW wind turbine and modeled the turbine as a structured grid with moving boundaries. Simulations were performed with a yaw rate of 0.3°/s and 3° yaw angle and the results showed that overall power under the dynamic yawing condition was much larger than for the yawed case. Then, Wen et al. [36] investigated the NREL 5-MW yawing dynamic sinuous motion with an averaged yaw rate of 1.2°/s and 2.4°/s for the case of f = 0.1 Hz and f = 0.2 Hz, respectively. The dynamic yaw motion induced the upwind & downwind yawing effect, which considerably influenced the AOA of blade sections. As the blade is yawing upwind, the AOA increases and vice versa. In summary, the unsteady aerodynamic characteristics of wind turbine blades under yaw considering full 3D rotational effects are still unclear. In the present paper, unsteady dynamic CFD simulations were performed to investigate the aerodynamic characteristics of a wind turbine blade during the rotational revolutions using a full 3D wind turbine model. With the assumption of the rigid body, the main contribution of this paper is the proposal of a new grid methodology for analyzing the overall performance using URANS simulation, aerodynamic loads, and flow field of wind turbines in the yawed and yawing processes. The method can be used to analyze interactions between transient aerodynamic phenomena associated with the wind turbine control system. The other main objective of this study was to investigate the effects of yaw on the dynamic output power, rotor thrust, and the blade sectional aerodynamic characteristics caused by continuous changes in the yaw angle. The main novelty of the current work is the inclusion of the dynamic yawing simulation of wind turbine with the multiple structured domain sliding mesh. The rest of this paper is organized as follows: Section 2 describes the geometry model and computational method. In Section 3, simulation results under several yaw angles are compared to experimental data to validate the model. The discussion mainly focuses on the sectional azimuthal variations of aerodynamic force coefficients, as well as the pressure distribution on the blade. Finally, the conclusions are summarized in Section 4.

2. Numerical Model

The 5-MW reference wind turbine (RWT) designed by the NREL was used in this study [37]. The blade uses the Delft University (DU) airfoil family with a relative thickness ranging from 21 to 40%. The blade chord, relative thickness, and twist angle all have non-linear distributions. The rotor diameter is 126 m and the wind turbine operates at a wind speed of 11.4 m/s with a rotational velocity of 12.1 rpm, resulting in the tip speed ratio of λ = 7. The definitions of yaw angle and azimuth angle are shown in Figure 1.

2.1. Computational Domain and Boundary Conditions

The computational domain, as shown in Figure 2a, is a rectangular zone that is 330D in length and 240D in width, where D is the rotor diameter. The computational domain is made up of two main zones, the far-field zone and internal zone, as depicted in Figure 2b. The far-field is used to set the inflow conditions. The internal zone can be further subdivided into two parts: (1) cylinder yaw zone with a radius of 6D and a height of 7.5D; (2) rotor roll zone with 1.5R radius and 1.6-m height. The cylinder zone was used to configure the dynamic yaw motion. The distance between the inlet and rotational cylinder was 120D. The rotor is located at the center of the rotor roll zone.
Figure 3 shows the boundary conditions used in the simulations. A uniform wind speed of 11.4 m/s was set at the inlet of the domain with a turbulence intensity of 5%. Boundary conditions for the up and bottom planes in the domain are free-slip and no-slip walls. The blade surface was set as a no-slip wall. The pressure outlet condition was assigned to the outlet of the domain. To simplify the static yaw simulation, the rotational axis was fixed, while the inlet velocity direction was varied to generate different angles mimicking different yaw condition. In Fluent, it is possible to select a mesh that can process non-conformal interfaces, namely boundaries interfaces between cell zones, for which the mesh node locations are not identical. In the current simulation case, the matching mesh interface option was applied to three interfaces, for example, the interface between the far-field zone and yawing zone, which is more accurate than the mapped interface option.

2.2. Grid Setup

Building multiple blocks is recommended for studying the yawing dynamics of coupled effects within two axes (yawing axis and rotational axis) since the relative rotation between zones can be investigated. Figure 4a,b show grids in the three different zones of the turbine (far-field, yaw cylinder, and rotation) generated using the ICEM software. The NREL 5-MW rotor was modeled including the hub (without the tower) [31]. The grid number for the far-field was 3.35 million, which is enough to transition from the inflow condition to the internal zone. Dynamic yaw motion was set, along with the rotational velocity in the cylinder yaw zone, with 3.25 million grids. Figure 4a illustrates the grids of the rotor zone, generated using AutoGrid tool in the NUMECA software. The grid number for this zone was 3.76 million, the first layer wall normal distance is about 10−4 m. The total grid number was about 10.27 million. The work of Tran and Kim [38] employed a 6 million-cell grid for a 5-MW wind turbine. The blade surface was resolved with 93 cells along the span. In order to resolve the boundary layer, the grids around the blades were refined to keep y+ less than 5 on the blade surface, which satisfies the needs of the T-SST turbulent model.

2.3. Numerical Methods

The unsteady Reynolds-averaged Navier-Stokes equations (URANS) were solved in ANSYS Fluent. The k-ω transitional SST turbulent model was used for turbulence modeling considering both laminar and transitional effects at the blade surface. The sliding mesh technique was used for data exchange between the rotational zone and the stationary zone. To ensure time-accurate calculations, a dual-time implicit time integration algorithm based on linearized second-order Euler backward differencing was used. The time-step size was equivalent to a rotor azimuthal increment of 5° coupled with 20 pseudo-time sub-iterations. The residual of the continuity equation was reduced by at least three orders throughout the calculations. For comparison, BEM computations using the Fast software were also performed. Work of Tran and Kim [38] employed 2° increments of the azimuth angle. The yawed case was achieved by changing the inlet velocity component to implement yaw inflow.
The yaw velocity can be constant, shown as the black dashed line in Figure 5. In fact, yawing is a process that changes gradually, similar to the definition of 2-s (or 4-s) start or stop duration. The selected variation of angular velocity for studying the influence of start-stop duration on the dynamic aerodynamic characteristics under yawing is shown a solid black line in Figure 6, which can be separated into three part: 1. Start duration with sinusoidal variation law of yaw velocity, 2. the stage with constant yaw velocity of 0.3°/s, 3. Stop duration with sinusoidal law of yaw velocity.
In the present investigation, the yawing process of the NREL 5-MW turbine was calculated under two different start-stop durations, 2-s and 4-s. Variation of the yaw velocity is shown in Figure 6, and changes more gradually with the 4-s duration of the start-stop process compared to the 2-s duration. During the start-stop yawing process, the sinusoidal variation of yaw velocity was applied to achieve transient yawing of the wind turbine, as follows:
y a w V e l ( t ) = r p m × sin ( 2 π ( 1 / ( d u r a t i o n × 4 ) × t )
where rpm = 0.3°/s, according to the description of the NREL5 WM design manual [37], and duration is either 2-s or 4-s. Up to 0.3°/s, the wind turbine changes from transient to a yaw angle of approximately 20° with the maximum rpm.

3. Analysis of Results

The NREL 5-MW reference wind turbine (RWT) was studied [28]. The yawed and yawing wind turbines were considered. The parameters used in the simulation are listed in Table 1 for each case. Two of the cases assumed uniform inflow and the blades were assumed to be rigid in the model.
Physically, the flow-field around a rotating HAWT is significantly influenced by the existence of wind shear, turbulence, gust, and yaw motion of the nacelle. For a yawing wind turbine, flow characteristics are more complex than those of a static yawed wind turbine, and the additional wind contribution effects transmitted to the rotor due to the nacelle motion must be considered. Therefore, accurate prediction of the unsteady aerodynamics calculated by many conventional numerical approaches is still questionable for a yawing wind turbine. In this study, unsteady CFD simulations based on the sliding structure mesh technique were performed to analyze the yawing motion of the wind turbine due to the nacelle motion. Thus, to investigate the effects of vortex-wake-blade interactions on the aerodynamic performance of the wind turbine, the yawing motion of the rotating turbine blades due to the nacelle motion was considered.
The CFD method was first verification with grid independence, turbulent model studies and time step studies. Then the computational results are validated using both fixed yaw and yawing of the NREL 5-MW wind turbine. The unsteady aerodynamic loads of the yawing wind turbine were more sensitive to change as the yaw angle was varied. Nine simulation conditions were performed based on a fixed yaw angle of 0°, 5°, 10°, 15°, 20°, 25°, and 30°.
Next, the total aerodynamic performance, which was used to validate the results of the simulation, will be discussed. Thereafter, aerodynamic characteristics and three-dimensional stall characteristics at different span sections will be investigated. Finally, the results of two separate simulations of the yawing start-stop stage will be used to examine the mechanism of the yawing effect.

3.1. Verification of the Computational Results

3.1.1. Grid Independence

Table 2 shows the different torques calculated for five different numbers of cells. Every mesh was refined along the airfoil circumferential direction. The grid numbers were computed in the T-SST turbulent mode and simulations were performed with the following grid numbers: 1.6 million, 3.09 million, 3.67 million, 5.31 million, and 6.05 million. In addition, two grid numbers (1.6 million and 3.67 million) were investigated under different yaw angles (10° and 20°), as shown in Table 2. For non-yaw condition, the design torque at 11.4 m/s is about 4.08 × 106 N·m, which is taken as a reference value. The relative errors for the results using other four mesh are 1.9~0.49%. The relative error for mesh of 3.67 million cells is smaller than 1% and the computational cost using this mesh is moderate. Thus this mesh was selected for further investigation.

3.1.2. Turbulent Model Studies

In this section, we present results of the unsteady independence study based on four turbulent models: SST, k-kl- ω , Reynold stress, and T-SST. The rotor torque under axial flow obtained for each turbulent model is presented in Table 3. The computed torque was similar for each of the four turbulent models. The simulation cases using the SST model and Reynold stress model were used to interpret the fully developed flow mechanism without transient processes, resulting less torque compared to using Transient SST turbulent model. When the wind turbine operates under low wind speeds, transient phenomena occur on the suction side of the blade; therefore, the T-SST turbulent model was selected for static and yawing process studies.

3.1.3. Time Step Studies

In the unsteady numerical simulation, the time step could influence the overall performance of the model, therefore, it is necessary to investigate time-step independence. Table 4 lists the averaged rotor torques during one revolution for non-yawed case obtained using the T-SST turbulent model with three different time steps (72, 144, and 360). The comparisons show slight differences among the results (less than 3%). The result using 72 time steps in one revolution gave a better agreement comparing to the designed value. Considering the computational cost, a time step of 72 during one revolution was selected for subsequent simulations.

3.2. The Validation of Numerical Simulation Results for Yawed Wind Turbine

After the verification of the simulation, the subsections below will show some aerodynamic analysis of wind turbine under yawed and yawing simulation.

3.2.1. Overall Performance Analysis for Different Yaw Angles

Figure 7 shows the overall performance of the rotor under yaw. As shown in Figure 7a, only a small amount of variation in rotor power can be observed under small yaw angles (≤5°). When the yaw angle exceeds 5°, the rotor torque significantly decreases as the yaw angle increases. The torque value computed using the BEM method, with or without the Beddoes stall model, is higher than the value computed by CFD since the BEM computation does not consider flow separation and flow transient phenomena. Zhu [39] extensively investigated the combined effects of rotational augmentation and dynamic stall, and found that the hysteresis loop of aerodynamic load is much larger compared to 2D simulations. In fact, for the NREL 5-MW wind turbine and a wind speed of 11.4 m/s, flow separation and 3D radial flow mainly occur on the inner board region, resulting in a lower torque and lower thrust. The deviation in the power and thrust between CFD and BEM was 4.1% and 5.91%, respectively. Compared with the axial free inflow, the three functions on cos(γ), cos2(γ), cos3(γ) of the power and thrust are shown with dashed, dotted and dash-dotted line, respectively. The variation of averaged power in yaw conditions will decrease by cos2(γ); the averaged thrust agrees well with cos(γ).

3.2.2. Aerodynamic Load Analysis along the Span of the Blade

Theoretical analysis is necessary to determine the mechanism behind the yawed effect. Figure 8 shows the velocity diagram under the yawed condition. A velocity component exists in the blade revolution plane and can be projected into radial and chordwise components denoted Vr and Vc, as shown in Figure 8a. The retreating and advancing effects occur because of Vc, which causes the sectional load and AOA to fluctuate periodically. The maximum AOA and maximum load occur when the blade is in the 12 o’clock direction in the current simulation setup. Additionally, loads on the retreating side are much larger than on the advancing side. Figure 9b,c quantify the influence of Vc, which is calculated by:
V r e l = ( V o ( cos γ a ) + η V d y n ) 2 + ( ω r ( 1 + b ) β V 0 sin γ cos φ ) 2
(1) Time-averaged load analysis on spanwise section under yawed condition
Figure 9, Figure 10 and Figure 11 show the analysis of the variation of averaged AOA, Cn, and Ct under one rotor revolution. The value of Cn and Ct is calculated by:
C n = x c = 0 x c = 1 C p - u p d ( x c ) x c = 0 x c = 1 C p d o w n d ( x c )
C t = y 0 y 1 C p - u p d ( y c ) y 0 y 1 C p d o w n d ( y c )
where Cp-up and Cp-down denote the pressure coefficients in the upper and down side of the airfoil, respectively; y0 and y1 mean the leading edge position and trailing edge position of the airfoil.
The AOA distribution presented in Figure 9 takes into account the advancing and retreating effects, as previously shown by Castellani [20] in the yawed simulation of a HAWT using the BEM model. Comparing the simulation results for CFD and FAST, the same trends can be observed along the blade span, whereas the AOA value is much larger than the FAST result. Figure 10 and Figure 11 show the variation of aerodynamic load spanwise along the blade. For the axial flow, some large abnormal fluctuations in the AOA and aerodynamic loads can be observe d. Differences in the results of the CFD and FAST methods mainly occur along the inner board, owing to flow separation, which lead to smaller aerodynamic loads compared to those computed by the BEM model that does not consider flow separation.
(2) The spanwise section aerodynamic load analysis under yawed conditions
Figure 12 and Figure 13 show the distribution of the AOA with respect to the azimuth angle under a yaw angle of 15° and 30°. Fluctuations become larger as the yaw angle increases due to the advancing and retreating effect, which is partly caused by the inflow velocity component under yaw. After post-processing, the AOA using the combination of the CFD sectional airfoil aerodynamic loads and the BEM method, the maximum and minimum AOA during one revolution occur at an azimuth angle of 0° and 180°, respectively. Wen [36] found that under the impact of non-uniform effects (due to variations of the induced factor caused by the radial position and azimuth angle in the rotational plane), the maximum AOA tends to occur at 90° for inboard airfoils. The maximum and minimum aerodynamic loads occur at an azimuth angle of 90° and 270°, respectively. Thus, future computations of the AOA should consider non-uniform effects.
Figure 14 shows the variation of the averaged AOA with aerodynamic load at five typical spanwise sections with respect to yaw angle. Fluctuations of the aerodynamic performance are clearly observed and become larger as the yaw angle increases.
The mean aerodynamic performance remains relatively constant with less decrease due to the yaw effect. Variation of the AOA and aerodynamic loads exhibit much higher fluctuations in the inner board compared to the middle and outer board. With the increasing of the yaw angle, the circumferential loads and AOA also cause higher maximum loads.

3.3. Numerical Simulation of Dynamic Yawing Wind Turbine

This section presents results of the numerical analysis of the NREL 5-MW RWT under the dynamic yawing process.

3.3.1. Torque Characteristics of Wind Rotor

Figure 15 shows the distribution of the wind rotor torque under different start-stop yaw velocities. During the yawing process, fluctuation of the wind rotor torque is small. At the beginning of dynamic yaw, larger torques occur due to changes in the rotor position. As the yaw angle increases, the torque gradually decreases and is 5 × 105 N·m larger in the case of the 2-s duration compared to the 4-s duration. The 2-s duration may induce a larger rotor torque since the yaw angle changes more quickly than with 4-s duration. The reasons are given below. At the start and stop stage, the yaw velocity is changed with the sinusoidal variation law related to the frequency, the higher frequency of yaw velocity caused much larger of the power and thrust, which is similar to the variation of power under platform pitching [38,40]. When the wind turbine is yawing with the constant yaw velocity, the additional velocity inducing by yawing is the same, and the power is only decreased with the square cosine of yaw angle. When wind turbine yawed to the stage of stop, the power under the 2-s case decrease much faster than under 4-s cases, which is similar to the yaw start period.
The results of 58 rotor revolutions were analyzed. During each revolution, 12 torque measurements were collected. Then, the fast Fourier transform was used to obtain the rotor aerodynamic frequency. Figure 16 shows the torque power spectra of the rotor under 2-s and 4-s duration yawing process. The blade passing frequency of the NREL 5-MW turbine is about 0.2017 Hz (1P fluctuation), and is clearly the main frequency of the two yawing processes. Both are 0.6 Hz (which is approximately a 3P rotor fluctuation), which is similar to results presented by Castellani [20]. Figure 16a shows a secondary frequency of 0.2 Hz (1P fluctuation). Due to the less sampling data in the sinusoidal stage, the yawing start-stop frequency (2-s and 4-s duration with a corresponding main frequency of 1/8 Hz and 1/16 Hz, respectively) cannot be captured in the computation of the torque power structure.

3.3.2. Torque Characteristics of Blade

Figure 17 shows the variation of torque under different yaw rates in two periodic processes of start-stop yawing. Similarly, the torque of the blade under a 2-s duration is larger than the torque of the 4-s duration. Interestingly, the torque of the blade in the forward yaw stage (shown in Figure 18 within 0~66.8 min, yaw angle from 0° to 20°) is larger than that of the backward yaw stage (shown in Figure 18 within 66.8~133.6 min, yaw angle from 20° back to 0°) since the dynamic yawing effect, which generates the dynamic velocity, reduces the relative velocity.

3.3.3. Wake Flow Characteristics

Wake effects are important in analyzing wind turbine aerodynamics. For convenience, the dynamic yaw was classified into six cases, as illustrated in Figure 18. Cases 1,3,4,6 include the forward-start yaw, forward-stop yaw, backward-start yaw, and backward-stop yaw, respectively. For Cases 2 and 5, the yaw angle is 10° and they represent forward yaw and backward yaw, respectively. The grey line in the figure illustration the variation of yaw angular velocity during the simulation, and shows that except the Cases 2 and 5, other cases are all in the simulation about the variation of yaw angular velocity.
(1)
Forward yaw-start stage (Case 1). The velocity contours at t = 1/8T (T = 8 s or 16 s) at the beginning of dynamic yaw are shown in Figure 19. The near wake velocity flow structure shows the same status at the beginning of the yawing process for both the 2-s and 4-s start-stop duration. The different yaw velocities have very little influence on the velocity distribution in the 2D range downstream.
(2)
Yaw angle of 10° under forward yawing stage (Case 2). Sketches of the velocity streamlines of the two dynamic yawing processes (yaw angle of 10°) are shown in Figure 20a,b are similar to those obtained for the yawed case, as shown in Figure 20c. The dynamic yaw rotates with a fixed yaw velocity of 0.3°/s, and similar results are observed for the dynamic process. More energy intermediate effects can be observed between the wake zone and the main flow zone than in the yawed case. This may be due to the effects of dynamic stall.
(3)
Forward-yaw-stop stage (Case 3). Figure 21 shows the instantaneous velocity contours for the dynamic yawing and yawed cases. Both dynamic yawing processes result in a much larger wake zones than under the yawed case. Meanwhile, some deflection occurs in the velocity wake, which is similar to the wake deflection effect reported in the work of Qian [18], which took into account the velocity deficit and turbulent intensity using a Gaussian-based wake model.
(4)
Backward-yaw start stage (Case 4). Figure 22 illustrates the velocity wake in the backward-yaw start stage. The wake zone velocity field retains almost the same flow structure. To refine the simulation results, large eddy simulations can be used to improve the interpretation of the flow process.
(5)
Yaw angle of 10° under backward yawing stage (Case 5). When the wind turbine rotates about the Y-axis with a yaw angle of 10°, the velocity wake gradually become symmetrical and the 4-s duration simulation recovers to the symmetry state faster than 2-s duration simulation, as shown in Figure 23.
(6)
Backward-yaw stop stage (Case 6). As shown in Figure 24, the wind turbine returns to its initial state, the wind direction is normal to the rotor rotational plane, and the three-dimensional flow structure is symmetrical in streamwise.
In summary, dynamic yawing based on two start-stop durations shows the approximate flow structure, and the upwind and downwind effect induced by yawing process expands the wake zone earlier and much larger than the yawed condition.

3.3.4. Aerodynamic Characteristics along Blade Spanwise Section

To investigate the yawing effect, additional velocity induced by the yawing wind rotor should also be examined. Figure 25 illustrates the velocity triangle along the span of the blade under dynamic yawing. The dashed line of the yawing zone indicates positive yaw (negative yaw was not investigated in the present simulation). Four process variations can be extracted according to the dynamic yawing stage. Note that R = 0 means that initial position of the rotational axis of the wind turbine is the rotor hub position, which is different from the platform yawing process used for the offshore wind turbines. The direction of Vdyn is different on both sides of the yaw axis. In the process of yawing counterwise, Vdyn and the inflow wind speed create an acute angle in the right side of yaw axis, while Vdyn and the inflow wind speed create an obtuse angle in the left one of yaw axis, as shown in Figure 26a. The process of yawing clockwise is just the contrary to the case of yawing counterwise. The absolute formulation of Vdyn can be written as:
V dyn = | ω y a w r sin ( φ ) |
where, Vdyn become zero at the azimuth angle of 0° and 180° in the current setup of wind turbine.
The relative velocity of the section under dynamic yawing can be defined as:
V r e l = ( V o ( cos γ a ) + η V dyn ) 2 + ( ω r ( 1 + b ) β V 0 sin γ cos φ ) 2
where η = 1 if the direction of Vdyn and inflow wind velocity creates an acute angle and η = −1 when the direction of Vdyn and inflow wind velocity creates an blunt angle. In the present yawing simulation, if the wind turbine is rotating in the positive yaw direction, β = 1.
Figure 26 shows the variation of the AOA in three typical section (r/R = 0.21, 0.4, or 0.67), with respect to yaw angles 0–20–0–20–0°. The upper abscissa is the wind rotor yaw angle, while the bottom abscissa is the wind turbine rotational cycle. The vertical coordinate indicates the AOA distribution along the radial direction of the blade. The three sections are indicated by solid lines of different thickness of airfoil. Inner board airfoils are the thickest and outer board one are thinnest.
The AOA gradually decreases along the span of the blade. Due to the effect of dynamic yawing, the AOA oscillates during the rotational period of the rotor, with the combined effects of both retreating & advancing and upwind & downwind. Figure 26a shows the calculated results of 2-s start-stop yaw rate. In this case, the start-stop velocity is faster, and the overall fluctuation are larger than those in the 4-s scenario. The start-stop effect mainly affects the AOA near the outer board but has less influence on other areas. In addition, the start-stop process leads to changes in the AOA.
Figure 27 show the distribution of the normal force coefficient at three typical sections (r/R = 0.21, 0.4, 0.67) within the dynamic process, including two forward yawing and backward yawing states with a yaw angle of 20°. Variation of Cn under dynamic yawing with a 2-s duration causes larger differences at the radial position (r/R = 0.4) than the 4-s duration. Under the yawed case, the average normal force coefficient of the r/R = 0.21 section with a yaw angle of 20° is 1.1, and the maximum and minimum load coefficients are 1.7 and 0.5, respectively (see Figure 14). In the current yawing case, the maximum and minimum load factors are 1.9 and 1.0, which are due to the downwind and upwind effects of the yaw dynamics, similar to the horizontal wind shear effect, but more pronounced than typical horizontal wind shear effects. In the start and stop duration of yawing, the yawing velocity of 2-s case has much higher frequency than 4-s case, resulting to much higher additional velocity and aerodynamic loads. Thus the overall performance of 2-s yawing presents much larger fluctuation in the process of yawing.
Figure 28 shows how the tangential force coefficient varies with simulation time at three typical sections (r/R = 0.21, 0.4, and 0.67). The results suggest the aerodynamic load inside the inner board is more influenced under dynamic yawing with the 2-s during than the 4-s duration. Fast shifts of the yaw angle under the yawing start-stop stage influence the aerodynamic load along the inner and middle blade span.

4. Conclusions

This paper has presented the results of unsteady numerical simulations investigating the static and dynamic aerodynamic performances of a NREL 5-MW HAWT under yawed and yawing conditions. The blades were considered rigid for cases with a prescribed motion, and the tower was not included in the computational domain. The numerical simulation results based on the T-SST turbulent model were consistent with the BEM data. Fluctuating and alternative loads were observed for the yawed case, resulting in the retreating and advancing effect. As the yawing angle decreased, overall performance also decreased according to the cosine law. The results suggest larger variations in power for the dynamic yaw case, compared to yawed cases. The wind turbine suffers from the coupled effects of the rotational axes (yawing axis and rotating axis). For maximum loading of rotor, the combined effects of the yawed condition and yawing resulted in an azimuth around the first periodic, which became larger as the yaw angle increased.
In the simulations of two different yawing start-stop durations (2-s and 4-s), the faster shift of yaw angle (2-s duration) resulted in a larger torque than the 4-s duration, owing to the additional velocity induced by the yawing effect of different sinusoidal frequency and yaw velocity magnitude. Interestingly, the torque of the rotor and blade was much higher in forward yawing compared to backward yawing, caused by reversing the direction of the yaw velocity. The wake deflection occurs in the near wake flow structure and expands more than in the yawed case. The different start-stop durations under dynamic yawing have little influence on the near wake velocity contour. The next step will be to include deformation of the blade in the simulation of multiple axial angular motions.

Author Contributions

S.K. and H.H. contributed to the conception of the study. X.W. and Z.Y. contributed significantly to analysis and manuscript preparation; Z.Y. performed the data analyses and wrote the manuscript; X.W. helped perform the analysis with constructive discussions.

Funding

This research was supported by National Natural Science Foundation of China (51576065), and the Fundamental Research Funds for the Central Universities (2015RCY05, 2018ZD09).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BEMBlade Element Momentum method
FVMFree Vortex Method
VTVortex Theory
CFDComputational Fluid Dynamic
SPIVStereo Particle Image Velocimetry
SCADASupervisory Control And Data Acquisition
FASTFatigue Aerodynamic Structure Turbulence software
FFTFast Fourier Transform
URANSUnsteady Reynolds Averaged Navier-Stokes equations
T-SSTTransient Shear Stress turbulence model
HAWTHorizontal Axis Wind Turbine
AOAAngle of Attack

Nomenclature

γ Yaw Angle
φ Azimuth Angle
β Pitch Angle
α g e o m Geometric angle of attack
α e f f Effective angle of attack
V0Inflow velocity
VrelRelative velocity
ω Rotation velocity of wind turbine
TTime
f Shedding/Meandering frequency
S t Strouhal number
R e Reynold number
λ blade tip speed ratio
CnNormal forces coefficient on the local airfoil
CtTangential forces coefficient on the local airfoil
Cp-upPressure coefficient at the upper side of airfoil
Cp-downPressure coefficient at the down side of airfoil
ClLife forces coefficient
CdDrag forces coefficient
CXTangent forces coefficient on the rotational plane
CZAxial forces coefficient on the rotational plane
CpPressure coefficient
CpoPower coefficent
CtThrust coefficent
TWind rotor torque
FtWind rotor thrust

References

  1. Rahimi, H.; Schepers, J.G.; Shen, W.Z.; García, N.R.; Schneider, M.S.; Micallef, D.; Ferreira, C.J.S.; Jost, E.; Klein, L.; Herráez, I. Evaluation of different methods for determining the angle of attack on wind turbine blades with CFD results under axial inflow conditions. Renew. Energy 2018, 125, 866–876. [Google Scholar] [CrossRef] [Green Version]
  2. Chehouri, A.; Younes, R.; Ilinca, A.; Perron, J. Review of performance optimization techniques applied to wind turbines. Appl. Energy 2015, 142, 361–388. [Google Scholar] [CrossRef]
  3. Hand, M.M.; Simms, D.A. Unsteady Aerodynamics Experiment Phase VI: Wind Tunnel Test Configurations and Available Data Campaigns; NREL/TP-500-29955; National Renewable Energy Lab: Golden, CO, USA, 2001.
  4. Schepers, J.G.; Boorsma, K.; Cho, T. Analysis of Mexico Wind Tunnel Measurements: Final Report of IEA Task 29, Mexnext (Phase 1); Energy Research Centre of the Netherlands: Amsterdam, The Netherlands, 2012. [Google Scholar]
  5. Sant, T. Improving BEM-Based Aerodynamic Models in Wind Turbine Design Codes. Ph.D. Thesis, University of Malta, Msida, Malta, 2007; pp. 45–80. [Google Scholar]
  6. Micallef, D.; Bussel, G.V.; Ferreira, C.S.; Sant, T. An investigation of radial velocities for a horizontal axis wind turbine in axial and yawed flows. Wind Energy 2013, 16, 529–544. [Google Scholar] [CrossRef]
  7. Vey, S.; Lang, H.M.; Nayeri, C.N.; Paschereit, C.O.; Pechlivanoglou, G.; Weinzierl, G. Utility scale wind turbine yaw from a flow visualization view. In Proceedings of the ASME Turbo Expo 2015: Turbine Technical Conference and Exposition, Montreal, QC, Cananda, 15 June 2015. V009T46A021. [Google Scholar]
  8. Dai, J.; Yang, X.; Hu, W.; Wen, L.; Tan, Y. Effect investigation of yaw on wind turbine performance based on SCADA data. Energy 2018, 149, 684–696. [Google Scholar] [CrossRef]
  9. Song, D.; Yang, J.; Fan, X.; Liu, Y.; Liu, A.; Chen, G.; Joo, Y.H. Maximum power extraction for wind turbines through a novel yaw control solution using predicted wind directions. Energy Convers. Manag. 2018, 157, 587–599. [Google Scholar] [CrossRef]
  10. Song, D.; Fan, X.; Yang, J.; Liu, A.; Chen, S.; Joo, Y.H. Power extraction efficiency optimization of horizontal-axis wind turbines through optimizing control parameters of yaw control systems using an intelligent method. Appl. Energy 2018, 224, 267–279. [Google Scholar] [CrossRef]
  11. Saenz-Aguirre, A.; Zulueta, E.; Fernandez-Gamiz, U.; Lozano, J.; Lopez-Guede, J.M. Artificial Neural Network Based Reinforcement Learning for Wind Turbine Yaw Control. Energies 2019, 12, 436. [Google Scholar] [CrossRef]
  12. Munters, W.; Meyers, J. Dynamic Strategies for Yaw and Induction Control of Wind Farms Based on Large-Eddy Simulation and Optimization. Energies 2018, 11, 177. [Google Scholar] [CrossRef]
  13. Ryu, K.; Kang, S.H. Prediction of Aerodynamic Loads for NREL Phase VI Wind Turbine Blade in Yawed Condition. Int. J. Aeronaut. Space Sci. 2016, 17, 157–166. [Google Scholar] [CrossRef] [Green Version]
  14. Jiménez, Á.; Crespo, A.; Migoya, E. Application of a LES technique to characterize the wake deflection of a wind turbine in yaw. Wind Energy 2010, 13, 559–572. [Google Scholar] [CrossRef]
  15. Jensen, N.O. A Note on Wind Generator Interaction; Risø National Laboratory: Roskilde, Denmark, 1983; ISBN 87-550-0971-9. [Google Scholar]
  16. Gebraad, P.M.O.; Teeuwisse, F.W.; van Wingerden, J.W.; Fleming, P.A.; Ruben, S.D.; Marden, J.R.; Pao, L.Y. Wind plant power optimization through yaw control using a parametric model for wake effects-a CFD simulation study. Wind Energy 2014, 19, 95–114. [Google Scholar] [CrossRef]
  17. Bastankhah, M.; Porté-Agel, F. Experimental and theoretical study of wind turbine wakes in yawed conditions. J. Fluid Mech. 2016, 806, 506–541. [Google Scholar] [CrossRef]
  18. Qian, G.W.; Ishihara, T.A. New Analytical Wake Model for Yawed Wind Turbines. Energies 2018, 11, 665. [Google Scholar] [CrossRef]
  19. Micallef, D.; Sant, T. Chapter 2: A review of wind turbine yaw aerodynamics. In Wind Turbines-Design, Control and Applications; InTech: Rijeka, Croatia, 2016; pp. 27–53. [Google Scholar]
  20. Castellani, F.; Astolfi, D.; Natili, F.; Marti, F. The Yawing Behavior of Horizontal-Axis Wind Turbines: A Numerical and Experimental Analysis. Machines 2019, 7, 15. [Google Scholar] [CrossRef]
  21. Qiu, Y.X.; Wang, X.D.; Kang, S.; Zhao, M.; Liang, J.Y. Predictions of unsteady HAWT aerodynamics in yawing and pitching using the free vortex method. Renew. Energy 2014, 70, 93–106. [Google Scholar] [CrossRef]
  22. Micallef, D.; Ferreira, C.S.; Sant, T.; Bussel, G.V. Experimental and numerical investigation of tip vortex generation and evolution on horizontal axis wind turbines. Wind Energy 2016, 19, 1485–1501. [Google Scholar] [CrossRef]
  23. Peters, D.A.; Haquang, N. Dynamic inflow for practical applications. J. Am. Helicopter Soc. 1988, 33, 64–68. [Google Scholar] [CrossRef]
  24. Gaonkar, G.H.; Peters, D. Review of dynamic inflow modeling for rotorcraft flight dynamics. In 27th Structures, Structural Dynamics and Materials Conference; VERTICA: Massachusetts, MA, USA, 1988; Volume 12, pp. 213–242. [Google Scholar]
  25. Suzuki, A. Application of Dynamic Inflow Theory to Wind Turbine Rotors. Ph.D. Thesis, The University of Utah, Salt Lake City, UT, USA, 2000. [Google Scholar]
  26. Glauert, H. Airplane propellers. In Aerodynamic Theory; Springer: Berlin, Germany, 1935; pp. 169–360. [Google Scholar]
  27. Du, Z.; Selig, M. A 3D Stall-Delay Model for Horizontal Axis Wind Turbine Performance Predictions; AIAA-98-0021; ASME: Reno, NV, USA, 1998; pp. 9–19. [Google Scholar]
  28. Buhl, M.L. A new empirical relationship between thrust coefficient and induction factor for the turbulent windmill state. Natl. Renew. Energy Lab. 2005, 47, 777–780. [Google Scholar]
  29. Yan, S.; Shi, S.P.; Chen, X.M.; Wang, X.D.; Mao, L.Z.; Liu, X.J. Numerical simulations of flow interactions between steep hill terrain and large scale wind turbine. Energy 2018, 151, 740–747. [Google Scholar] [CrossRef]
  30. Tongchitpakdee, C.; Benjanirat, S.; Sankar, L. Numerical Simulation of the Aerodynamics of Horizontal Axis Wind Turbines Under Yawed Flow Conditions. J. Sol. Energy Eng. 2005, 127, 464–474. [Google Scholar] [CrossRef]
  31. Yu, D.O.; You, J.Y.; Kwon, O.J. Numerical investigation of unsteady aerodynamics of a Horizontal-axis wind turbine under yawed flow conditions. Wind Energy 2013, 16, 711–727. [Google Scholar] [CrossRef]
  32. Schulz, C.; Letzgus, P.; Lutz, T.; Kramer, E. CFD study on the impact of yawed inflow on loads, power and near wake of a generic wind turbine. Wind Energy 2017, 20, 253–268. [Google Scholar] [CrossRef]
  33. Jeong, M.S.; Kim, S.W.; Lee, I.; Yoo, S.J.; Park, K.C. The impact of yaw error on aeroelastic characteristics of a horizontal axis wind turbine blade. Renew. Energy 2013, 60, 256–268. [Google Scholar] [CrossRef]
  34. Dai, L.; Zhou, Q.; Zhang, Y.W.; Yao, S.G.; Kang, S.; Wang, X.D. Analysis of wind turbine blades aeroelastic performance under yaw conditions. J. Wind Eng. Ind. Aerodyn. 2017, 171, 273–287. [Google Scholar] [CrossRef]
  35. Leble, V.; Barakos, G. 10-MW wind turbine performance under pitching and yawing motion. J. Sol. Energy Eng. 2017, 139, 041003. [Google Scholar] [CrossRef]
  36. Wen, B.R.; Tian, X.L.; Dong, X.J.; Peng, Z.K.; Zhang, W.M.; Wei, K.X. A numerical study on the angle of attack to the blade of a horizontal-axis offshore floating wind turbine under static and dynamic yawed conditions. Energy 2019, 168, 1138–1156. [Google Scholar] [CrossRef]
  37. Jonkman, J.; Butterfield, S.; Musial, W.; Scott, G. Definition of a 5-MW Reference Wind Turbine for Offshore System Development; National Renewable Energy Lab. (NREL): Golden, CO, USA, 2009.
  38. Tran, T.T.; Kim, D.H. The platform pitching motion of floating offshore wind turbine: A preliminary unsteady aerodynamic analysis. J. Wind Eng. Ind. Aerodyn. 2015, 142, 65–81. [Google Scholar] [CrossRef] [Green Version]
  39. Zhu, C.; Wang, T.; Zhong, W. Combined Effect of Rotational Augmentation and Dynamic Stall on a Horizontal Axis Wind Turbine. Energies 2019, 12, 1434. [Google Scholar] [CrossRef]
  40. Wen, B.R.; Dong, X.; Tian, X.; Peng, Z.; Zhang, W.; Wei, K. The power performance of an offshore floating wind turbine in platform pitching motion. Energy 2018, 154, 508–521. [Google Scholar] [CrossRef]
Figure 1. Sketch of the asymmetric relative velocity under yaw.
Figure 1. Sketch of the asymmetric relative velocity under yaw.
Energies 12 03124 g001
Figure 2. Computational Zones: (a) far-field (static); (b) yaw zone; (c) rotational zone.
Figure 2. Computational Zones: (a) far-field (static); (b) yaw zone; (c) rotational zone.
Energies 12 03124 g002
Figure 3. Boundary conditions.
Figure 3. Boundary conditions.
Energies 12 03124 g003
Figure 4. Computational mesh for NREL 5-MW wind turbine; (a) rotor cylinder; (b) blade surface mesh.
Figure 4. Computational mesh for NREL 5-MW wind turbine; (a) rotor cylinder; (b) blade surface mesh.
Energies 12 03124 g004
Figure 5. The sketch for start-stop yawing process.
Figure 5. The sketch for start-stop yawing process.
Energies 12 03124 g005
Figure 6. Variation of yaw velocity under different transient start-stop processes.
Figure 6. Variation of yaw velocity under different transient start-stop processes.
Energies 12 03124 g006
Figure 7. Variation of the power and thrust of the wind turbine under yaw.
Figure 7. Variation of the power and thrust of the wind turbine under yaw.
Energies 12 03124 g007
Figure 8. Section airfoil induced velocity Vc under stable yaw: (a) magnitude of radial and chordwise velocity; (b) direction of radial and chordwise velocity; (c) velocity diagram of blade element.
Figure 8. Section airfoil induced velocity Vc under stable yaw: (a) magnitude of radial and chordwise velocity; (b) direction of radial and chordwise velocity; (c) velocity diagram of blade element.
Energies 12 03124 g008
Figure 9. The profile of AOA along the blade span under yawed cases.
Figure 9. The profile of AOA along the blade span under yawed cases.
Energies 12 03124 g009
Figure 10. The profile of Cn along the blade span under yawed cases.
Figure 10. The profile of Cn along the blade span under yawed cases.
Energies 12 03124 g010
Figure 11. The profile of Ct along the blade span under yawed cases.
Figure 11. The profile of Ct along the blade span under yawed cases.
Energies 12 03124 g011
Figure 12. Azimuthal aerodynamic loads at different spanwise sections during one revolution under a yaw angle of 15°.
Figure 12. Azimuthal aerodynamic loads at different spanwise sections during one revolution under a yaw angle of 15°.
Energies 12 03124 g012
Figure 13. Azimuthal aerodynamic loads at different spanwise sections during one revolution under a yaw angle of 30°.
Figure 13. Azimuthal aerodynamic loads at different spanwise sections during one revolution under a yaw angle of 30°.
Energies 12 03124 g013
Figure 14. Variation of time-averaged AOA and aerodynamic load under yaw during one revolution.
Figure 14. Variation of time-averaged AOA and aerodynamic load under yaw during one revolution.
Energies 12 03124 g014
Figure 15. Torque characteristics of rotor under two different yawing rotational velocities.
Figure 15. Torque characteristics of rotor under two different yawing rotational velocities.
Energies 12 03124 g015
Figure 16. Frequency of rotor torque for a yawing start-stop process with a duration of 2-s and 4-s.
Figure 16. Frequency of rotor torque for a yawing start-stop process with a duration of 2-s and 4-s.
Energies 12 03124 g016
Figure 17. Torque characteristics of blades under two different yaw velocities.
Figure 17. Torque characteristics of blades under two different yaw velocities.
Energies 12 03124 g017
Figure 18. Six cases during the start-stop stage under dynamic yaw.
Figure 18. Six cases during the start-stop stage under dynamic yaw.
Energies 12 03124 g018
Figure 19. Velocity contours for Case 1.
Figure 19. Velocity contours for Case 1.
Energies 12 03124 g019
Figure 20. Velocity contours for Case 2.
Figure 20. Velocity contours for Case 2.
Energies 12 03124 g020
Figure 21. Velocity contours for Case 3.
Figure 21. Velocity contours for Case 3.
Energies 12 03124 g021
Figure 22. Velocity contours for Case 4.
Figure 22. Velocity contours for Case 4.
Energies 12 03124 g022
Figure 23. Velocity contours for Case 5.
Figure 23. Velocity contours for Case 5.
Energies 12 03124 g023
Figure 24. Velocity contours for Case 6.
Figure 24. Velocity contours for Case 6.
Energies 12 03124 g024
Figure 25. Velocity diagram for a blade section under dynamic yaw at radius r: (a) sketch of the yawing process for the velocity dynamic analysis; (b) velocity diagram along the span of the blade.
Figure 25. Velocity diagram for a blade section under dynamic yaw at radius r: (a) sketch of the yawing process for the velocity dynamic analysis; (b) velocity diagram along the span of the blade.
Energies 12 03124 g025
Figure 26. Comparison of AOA under different start-stop durations for three different spanwise sections under dynamic yawing.
Figure 26. Comparison of AOA under different start-stop durations for three different spanwise sections under dynamic yawing.
Energies 12 03124 g026
Figure 27. Comparison of Cn under different start-stop durations for three different spanwise sections under dynamic yawing.
Figure 27. Comparison of Cn under different start-stop durations for three different spanwise sections under dynamic yawing.
Energies 12 03124 g027
Figure 28. Comparison of Ct under different start-stop durations for three different spanwise sections under dynamic yawing.
Figure 28. Comparison of Ct under different start-stop durations for three different spanwise sections under dynamic yawing.
Energies 12 03124 g028
Table 1. Simulation cases.
Table 1. Simulation cases.
ConfigurationAngle/Amplitude (°)Unsteady Computation Time (s)
Aligned089.26
Fixed yaw5, 10, 15, 20, 25, 3089.25
Dynamic yawing 0–20–0126
Table 2. Computed torque of different grid size averaged in one revolution.
Table 2. Computed torque of different grid size averaged in one revolution.
Torque (N·m)Number of Cells (Million)
Yaw Angle (°) 1.63.093.675.316.05Designed Value
4.0 × 1064.05 × 1064.06 × 1064.11 × 1064.10 × 1064.08 × 106
relative error1.9%0.73%0.49%0.73%0.49%
Table 3. Time averaged rotor torque calculated using different unsteady turbulence models.
Table 3. Time averaged rotor torque calculated using different unsteady turbulence models.
Turbulent Modelsk-ω SSTk-kl-ωReynold StressT-SSTDesigned Value
Torque/(N·M)3.89 × 1064.17 × 1063.76 × 1064.06 × 1064.08 × 106
relative error4.65%2.2%7.8%0.49%
Table 4. Rotor torque calculated using different time step sizes averaged over one revolution.
Table 4. Rotor torque calculated using different time step sizes averaged over one revolution.
Time Step72144360Designed Value
Torque/(N·M)4.06 × 1064.16 × 1064.19 × 1064.08 × 106
relative error0.49%1.9%2.69%

Share and Cite

MDPI and ACS Style

Wang, X.; Ye, Z.; Kang, S.; Hu, H. Investigations on the Unsteady Aerodynamic Characteristics of a Horizontal-Axis Wind Turbine during Dynamic Yaw Processes. Energies 2019, 12, 3124. https://doi.org/10.3390/en12163124

AMA Style

Wang X, Ye Z, Kang S, Hu H. Investigations on the Unsteady Aerodynamic Characteristics of a Horizontal-Axis Wind Turbine during Dynamic Yaw Processes. Energies. 2019; 12(16):3124. https://doi.org/10.3390/en12163124

Chicago/Turabian Style

Wang, Xiaodong, Zhaoliang Ye, Shun Kang, and Hui Hu. 2019. "Investigations on the Unsteady Aerodynamic Characteristics of a Horizontal-Axis Wind Turbine during Dynamic Yaw Processes" Energies 12, no. 16: 3124. https://doi.org/10.3390/en12163124

APA Style

Wang, X., Ye, Z., Kang, S., & Hu, H. (2019). Investigations on the Unsteady Aerodynamic Characteristics of a Horizontal-Axis Wind Turbine during Dynamic Yaw Processes. Energies, 12(16), 3124. https://doi.org/10.3390/en12163124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop