In order to closely examine the wind flow characteristics that have caused failures of the yaw gears and motors on T7, a CFD simulation was conducted. The CFD software used in the present study is RIAM-COMPACT, which was developed by the first author of the present paper. Because the details of the numerical simulation methods used in RIAM-COMPACT have been discussed in previous papers [
13,
14,
15,
16,
17,
18,
19,
20,
21], they will be omitted here. In the present study, the LES is assumed to reproduce the wind tunnel testing. Therefore, the effects of atmospheric stability associated with vertical thermal stratification of the atmosphere and inflow turbulence were neglected. Thus, TI calculated from numerical results is smaller than measured data. The standard Smagorinsky model was used for the subgrid-scale model (SGS) model [
22]. The model coefficient was assumed to be 0.1 by using a wall-damping function. In addition, as in [
13,
16,
17,
18], the effects of the surface roughness were taken into consideration by reconstructing surface irregularities in high resolution. A comparison between the Reynolds-averaged modeling (RANS) results and the present LES results is summarized in a recent article [
14], and the prediction accuracy of the present LES approach by comparison with wind tunnel experiments is discussed in [
21].
Figure 10 illustrates the computational domain and grid used for the present study. For the study, the complex terrain of the wind farm and its surroundings were numerically constructed using the 10-m resolution land surface digital elevation model (DEM) from the Geospatial Information Authority of Japan (GSI). The dimensions of the computational domain were 10.0 km × 4.0 km × 4.0 km in full scale in the streamwise (x), spanwise (y), and vertical (z) directions, respectively. The computational domain was set in such a way that T7 was located in the center of the x-y plane of the computational domain. The grid spacing for all directions varied so that grid density was high in the vicinity of T7. The minimum grid spacing for the x- and y-directions was set to approximately 8 m, and the minimum grid spacing for the z-direction was set to approximately 1.7 m. The number of grid points was 501 × 201 × 101 in the x-, y-, and z-directions, respectively, which resulted in a total number of approximately 10 million grid points. The wind direction considered for the simulation was west-north-westerly. (West-north-westerly wind is wind that flows from an angle of 292.5° clockwise from north, which was set to 0° in the wind direction coordinate system adopted in the present study.) Because inflow boundary was located over the ocean (
Figure 10), the vertical wind profile of the inflow streamwise wind velocity was set according to a power law (α = 0.1, where α is the power law exponent, i.e., N = 10, where N is the inverse of α). Other boundary conditions are detailed in [
13,
14,
15,
16,
17,
18,
19,
20,
21].
3.1. Results of Non-Dimensional Simulation and Discussions
The governing equations of the flow adopted in the present study (i.e., a filtered continuity equation for incompressible viscous flows and filtered Navier–Stokes equations for incompressible viscous flows) were non-dimensionalized with the use of 1) the difference between the minimum and maximum terrain elevations within the computational domain and 2) the inflow streamwise wind velocity at the height of the maximum surface elevation in the computational domain. Therefore, all the physical variables that are output by the model are non-dimensional quantities. In the current sub-section, the non-dimensional data of the three wind velocity components that were obtained from the simulation will be analyzed, and the results will be shown and discussed.
Figure 11 shows the temporal change of the fluctuating parts of the three wind velocity components at the hub height of T7 (60 m above the ground surface) in the non-dimensional time period from 600 to 800. Specifically, the data values shown for each wind velocity component are those that were obtained by subtracting the period-averaged wind velocity component from the original time series of the wind velocity component.
Figure 11 shows that the values of the streamwise (x), spanwise (y), and vertical (z) wind velocity components all fluctuated significantly in time.
Figure 12 shows vertical profiles of statistical quantities of the turbulent flow at the site of T7 from the non-dimensional time period from 600 to 800.
Figure 12a shows the vertical profiles of the streamwise wind velocity. Specifically, the red line in
Figure 12 indicates the vertical profile of the streamwise inflow velocity, and the blue line indicates the vertical profile of the mean streamwise wind velocity at the site of T7 from the time period under investigation.
Figure 12a also includes the values of the speed-up ratio at the bottom of the swept area (29 m above the ground surface), at the hub center (60 m above the ground surface), and at the top of the swept area (91 m above the ground surface), where the speed-up ratio is defined as the ratio of the streamwise wind velocity at a height of interest above the ground surface at the site of T7 to the inflow streamwise wind velocity at the height of interest. These results show that, due to terrain effects, the streamwise wind velocity increased locally at the wind turbine site, and additionally, there was no significant streamwise wind velocity deficit at the site. As can be presumed from
Figure 10 and
Figure 13, the locally increased streamwise wind velocity likely occurred as the wind flowed uphill along the terrain and into the wind turbine.
Figure 12b shows the vertical profiles of the standard deviations of the streamwise (x), spanwise (y), and vertical (z) wind velocity components. The values of the standard deviations for all three components are relatively large, reflecting the temporal change of the fluctuating parts of the wind velocity components in
Figure 11. Examinations of the values of the standard deviations at the hub center (60 m above the ground surface) in
Figure 12b reveal that the value of the standard deviation of the vertical (z) wind velocity component is large and that the ratio of the values of the standard deviations of the three wind velocity components at the hub center was σ
1:σ
2:σ
3 = 1.0:0.7:0.65, which clearly indicates that there was an influence of terrain-induced turbulence at this site. (This finding will be discussed again in
Section 3.2).
To investigate details of the flow field in the vicinity of T7,
Figure 13 illustrates the temporal change of the streamwise (x) wind velocity as contour plots. The visualized streamwise wind velocity field in
Figure 13 shows that, as a separation vortex (indicated by the arrows in
Figure 13) that was shed upstream of the wind turbine passed through the wind turbine, the wind velocity field surrounding the wind turbine changed significantly.
Figure 14 shows the vertical profiles of the streamwise (x) wind velocity at the site of T7 from the same times for which the cross-sectional views of the streamwise (x) wind velocity in
Figure 13 were created. Immediately before the separation vortex passed through the wind turbine (
Figure 13 (a)), the vertical profile of the streamwise (x) wind velocity showed a local increase of the velocity due to terrain effects and thus showed no significant wind velocity deficit with respect to the power law profile of the streamwise (x) wind velocity (
Figure 14a). As the separation vortex that had been located upwind of the turbine approached the turbine, a velocity deficit occurred in the layer between the hub center (60 m above the ground surface) and the bottom of the rotor (
Figure 14b). At the time at which the separation vortex arrived at the wind turbine (
Figure 14c), negative wind shear was evident between the hub center height (60 m above the ground surface) and heights that were slightly higher than the hub center height. As illustrated in
Figure 14d,e, after the passage of the separation vortex, the wind velocity recovered to values predicted by the power law.
A close examination of computer animations of the simulations results in
Figure 13 and
Figure 14 led to the following finding: The sequence of wind flow patterns described above, that is, large vortex shedding that originated from the micro-topographical features upstream of the wind turbine and its accompanying generation and disappearance of a separation vortex, occurred in a nearly periodic manner. A complex wind flow field with a vertical profile of the streamwise wind velocity such as the one in
Figure 14c, which is generally rare, periodically formed in the vicinity of the wind turbine. More specifically, this vertical profile of streamwise wind velocity deviated significantly from the power law and also had negative wind shear. It can be surmised that when complex wind flow with such a profile passes through a wind turbine, it causes a large wind load on the turbine. In addition, because the structure of the abovementioned complex wind flow is three-dimensional, wind loads on the left and right side of the swept area of a wind turbine in such a wind flow are expected to differ. Thus, it can be speculated that such wind loads would exert force on the wind turbine in such a way that they would forcibly rotate the nacelle of the wind turbine, which in turn would cause impact loads on both the yaw gears and motors and result in the failures of the yaw gears and motors in the end. Accordingly, it may be possible to make prior assessments of wind turbine failure risks due to terrain-induced turbulence by studying, with the use of CFD, wind velocity fluctuations in the vicinity of a wind turbine and the vertical profiles of statistical quantities of the three velocity components (i.e., the three-dimensional flow structure) within the swept area.
3.2. Re-Scaled Dimensional Simulation Results and Discussions
The temporal change of the streamwise (x) wind velocity (a non-dimensional quantity) in the period from 600 to 800 in non-dimensional time at the hub height of T7 was re-scaled such that the average value of the streamwise (x) wind velocity at the hub height of this turbine for this period became 8.0 m/s. The abovementioned re-scaling procedure can be summarized as follows:
- (1)
The average value of the streamwise (x) wind velocity (a non-dimensional quantity) at the hub height of T7 was calculated for the period from 600 to 800 in non-dimensional time. The calculated average value was 1.087 in the present study.
- (2)
A correction coefficient was calculated so that the average value of the streamwise (x) wind velocity at the hub height of T7 in the period under investigation was 8.0 m/s in full scale. Then, the non-dimensional wind velocity data from the entire simulation time period were multiplied by the calculated correction coefficient. The calculated correction coefficient was 7.36 (= 8.0/1.087) in the present study. With this procedure, the streamwise (x) non-dimensional wind velocity was converted to full-scale wind velocity (m/s).
- (3)
The time in the period from 600 to 800 in non-dimensional time was converted to full scale using the equation T = t (h/Uin), where T is full-scale time (s), t is non-dimensional time, h is the difference between the minimum and maximum terrain elevations within the computational domain (m), and Uin is the streamwise wind velocity (m/s) at the height of the maximum surface elevation in the computational domain at the inflow boundary. In the present study, the 200 non-dimensional time period, i.e., the period from 600 to 800 in non-dimensional time, was converted to approximately 15,500 s (approximately 4 h) in full scale. The time step was 0.3 s in full scale.
Figure 15 shows the temporal change of the full-scale streamwise (x) wind velocity (m/s) that was obtained from the rescaling procedure with the method described above (total data points: 50,000; time interval: 0.3 s). The green line in the figure indicates 8.0 m/s, which is the average streamwise (x) wind velocity that the re-scaling procedure was designed to attain for the hub height of T7 in the time period under investigation.
Figure 16 shows a histogram of the streamwise (x) wind velocity data from
Figure 15 with bin widths of 1 m/s. The average streamwise (x) wind velocity that the re-scaling procedure was designed to attain for the time period investigated was 8.0 m/s in the present study. As a result, the re-scaled streamwise (x) wind velocity ranged between approximately 4.0 and 11.0 m/s, and the occurrence frequency of the wind velocity class of 9.0 to 10.0 m/s in particular was large.
In the present study, following a common statistical processing procedure adopted for in situ data, a 10-min moving average filter (1932-point averaging filter) was applied to the time series of the re-scaled streamwise (x) wind velocity (m/s) in
Figure 15 (total data points: 50,000; time interval: 0.3 s) to evaluate the values of the moving-averaged wind velocity and the corresponding TI (
Figure 17) (48,068 data points). In
Figure 17a, the green line indicates 8.0 m/s, which is the average streamwise wind velocity (x) that the re-scaling procedure in this study was designed to attain for the time period under investigation.
Figure 17b shows the evaluated TI values. These values were obtained with Equation (1) below using a moving-averaged filter with a window length of 10-min and the wind velocity data within the window (number of sample data points: 1932).
where
where u(t) is the instantaneous streamwise wind velocity,
is the average of the instantaneous streamwise wind velocity within the 10-min moving-averaged window, and u’ is the fluctuating component of the streamwise wind velocity due to turbulence.
Figure 17a,b shows that the values of the average streamwise velocity and TI fluctuate in a correlated manner.
That is, as the average wind velocity in
Figure 17a increases, the TI in
Figure 17b decreases. Conversely, as the average wind velocity in
Figure 17a decreases, the TI in
Figure 17b increases. A further examination of the temporal change of the TI in
Figure 17b reveals that the TI changes in large amplitude with the increasing and decreasing average wind velocity. The average value of TI was 0.19 (the green line in
Figure 17b), which is relatively large. This result also indicates that T7 was strongly affected by terrain-induced turbulence.
Finally, the TI values from
Figure 17 were examined by comparing them with those from the NTM in IEC 61400-1 Ed.3 (2005) (
Figure 18). The NTM defines wind turbine classes as in
Table 1. V
ref in
Table 1 represents the 50-year return period values of 10-min average wind speed. I
ref is the expected value of TI for a wind speed of 15 m/s.
For NTM, the values of and
and
are calculated using the streamwise (x) wind velocity as
where I
ref is the expected value for the turbulence intensity for V = 15 m/s; TI is the turbulence intensity; U is the 10-min average streamwise wind velocity (m/s); σ is wind velocity standard deviation (m/s); and subscript 90q is the 0.9 quantile value.
Wind turbine designers design wind turbines in such a way that they meet both the wind turbine class and turbulence class requirements. Wind power providers are able to reduce their business risks by confirming that the values of TI at their wind turbine sites lie under the curve defined by Equation (4).
Figure 18 shows that the values of the streamwise (x) TI simulated with the abovementioned method were not as large as those observed in situ (
Figure 6). However, some of the simulated values exceeded the TI values for turbulence class A, suggesting that the influence of terrain-induced turbulence on the wind turbine was well simulated.
Based on turbulence spectral relationships, the spanwise (y) wind-velocity standard deviation, σ
2, and the vertical (z) wind-velocity standard deviation, σ
3, are given with respect to the streamwise (x) wind-velocity standard deviation, σ
1, as in Equations (5) and (6). Both of these equations were derived from turbulence spectra from wind flow over flat terrain.
In the present study, the value of the standard deviation of the vertical (z) wind velocity, σ
3, at the hub height of T7 in
Figure 12b was fairly large as discussed earlier and shown below: