LES Investigation of Terrain-Induced Turbulence in Complex Terrain and Economic Effects of Wind Turbine Control

In the present study, numerical wind simulation was conducted by reproducing the realistic topography near wind turbine sites with high spatial resolutions and using the Large-Eddy Simulation (LES) technique. The topography near wind turbine sites serves as the origin of the terrain-induced turbulence. The obtained numerical simulation results showed that the terrain-induced turbulence is generated at a small terrain feature located upstream of the wind turbine. The generated terrain-induced turbulence affects the wind turbine directly. The wind speed and wind direction at the wind turbine site are significantly changed with time. In the present study, a combination of the series of wind simulation results and on-site operation experience led to a decision to adopt an “automatic shutdown program”. Here, an “automatic shutdown program” means the automatic suspension of wind turbine operation based on the wind speed and wind direction meeting the conditions associated with significant effects of terrain-induced turbulence at a wind turbine site. The adoption of the “automatic shutdown program” has successfully led to a large reduction in the number of occurrences of wind turbine damage, thus, creating major positive economic effects.

The numerical wind prediction technique, RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain), which has been developed by the authors' research group, has a potential to resolve the above-mentioned issue all at once [7][8][9][10][11][12][13][14][15]. The core technology of the RIAM-COMPACT is under continuous development at the Research Institute for Applied Mechanics (RIAM), Kyushu University.
Recently, it has been reported that the availability factors of wind turbines on wind farms situated on complex terrain fall short of those originally projected; that is, reports of damage and breakage of the exteriors and interiors of wind turbines (e.g., breakage of yaw motors and yaw gears, and cracks on wind turbine blades), as well as wind turbines with, notably, low power output have surfaced. The main cause of these problems is that small variations in the topographical relief in the vicinity of wind turbines serve as the origin of turbulence (terrain-induced turbulence). The present paper will discuss a specific example of a wind-risk (terrain-induced turbulence) diagnosis using the RIAM-COMPACT natural terrain version software (Large-Eddy Simulation (LES) modeling).

Summary of the RIAM-COMPACT Natural Terrain Version Software
In the present study, to predict airflows over real complex terrain with high accuracy, the RIAM-COMPACT natural terrain version software (LES modeling) is used [7][8][9][10][11][12][13][14][15]. The RIAM-COMPACT natural terrain version software uses collocated grids in a general curvilinear coordinate system to avoid numerical instability. In these collocated grids, the physical velocity components and pressure are defined at the cell centers and variables, which result from the contravariant velocity components multiplied by the Jacobian, are defined at the cell faces. As for the simulation technique, the Finite-Difference Method (FDM) is adopted and a Large-Eddy Simulation (LES) model is used for the turbulence model. In LES, a spatial filter is applied to the flow field to separate eddies of various scales into Grid Scale (GS) components, which are larger than the computational grid cells, and Sub-Grid Scale (SGS) components, which are smaller than the computational grid cells. Large-scale eddies, i.e., the GS components of turbulence eddies, are directly numerically simulated without relying on the use of a physically simplified model. On the other hand, the main effect of small-scale eddies, i.e., the SGS components, is to dissipate energy and this dissipation is modeled based on the physical considerations of the SGS stress.
For the governing equations of the flow, a spatially-filtered continuity equation for incompressible fluid (Equation (1)) and a spatially-filtered Navier-Stokes equation (Equation (2)) are used: Supporting equations are given in Equations (3)- (8): ν SGS = (C s f s ∆) 2 S (4) f s = 1 − exp −z + /25 (7) The computational algorithm and the time marching method are based on the Fractional-Step (FS) method [16] and the Euler explicit method, respectively. The Poisson's equation for pressure is solved by the Successive Over-Relaxation (SOR) method. For discretization of all spatial terms in Equation (2), except for the convective term, a second-order central difference scheme is applied. For the convective term, a third-order upwind difference scheme is applied. An interpolation technique based on four-point differencing and four-point interpolation by Kajishima [17] is used for the fourth-order central differencing that appears in the discretized form of the convective term. For the weighting of the numerical diffusion term in the convective term discretized by third-order upwind differencing, α = 3.0 is commonly applied in the Kawamura-Kuwahara scheme [18]. However, α = 0.5 is used in the present study to minimize the influence of numerical diffusion. For LES subgrid-scale modeling, the commonly used Smagorinsky model [19] is adopted. A wall-damping function is used with a model coefficient of 0.1.

Wind Synopsis Analysis with the RIAM-COMPACT Natural Terrain Version Software
As described earlier, the availability factors of wind turbines in wind farms on complex terrain often fall below the originally projected values. In other words, problems, such as wind turbines with notably low power output and damage and breakage of the exteriors and interiors of wind turbines, have surfaced. The effects of turbulence (terrain-induced turbulence) have also been made note of for the Atsumi Wind Farm in the Aichi Prefecture, Japan, the site investigated in the present study. In the present paper, an LES turbulence simulation is conducted with high spatial resolution using the RIAM-COMPACT natural terrain version software. The main purpose of this simulation is to reproduce the relative comparisons of wind conditions among wind turbine locations.  Table 1). subgrid-scale modeling, the commonly used Smagorinsky model [19] is adopted. A wall-damping function is used with a model coefficient of 0.1.

Wind Synopsis Analysis with the RIAM-COMPACT Natural Terrain Version Software
As described earlier, the availability factors of wind turbines in wind farms on complex terrain often fall below the originally projected values. In other words, problems, such as wind turbines with notably low power output and damage and breakage of the exteriors and interiors of wind turbines, have surfaced. The effects of turbulence (terrain-induced turbulence) have also been made note of for the Atsumi Wind Farm in the Aichi Prefecture, Japan, the site investigated in the present study. In the present paper, an LES turbulence simulation is conducted with high spatial resolution using the RIAM-COMPACT natural terrain version software. The main purpose of this simulation is to reproduce the relative comparisons of wind conditions among wind turbine locations.  Table 1).    subgrid-scale modeling, the commonly used Smagorinsky model [19] is adopted. A wall-damping function is used with a model coefficient of 0.1.

Wind Synopsis Analysis with the RIAM-COMPACT Natural Terrain Version Software
As described earlier, the availability factors of wind turbines in wind farms on complex terrain often fall below the originally projected values. In other words, problems, such as wind turbines with notably low power output and damage and breakage of the exteriors and interiors of wind turbines, have surfaced. The effects of turbulence (terrain-induced turbulence) have also been made note of for the Atsumi Wind Farm in the Aichi Prefecture, Japan, the site investigated in the present study. In the present paper, an LES turbulence simulation is conducted with high spatial resolution using the RIAM-COMPACT natural terrain version software. The main purpose of this simulation is to reproduce the relative comparisons of wind conditions among wind turbine locations.  Table 1).

Simulation Set-Up
As shown in Figure 4, the dimensions of the computational domain are 2.0 (x) × 1.5 (y) × 1.0 (z) km, where x is the streamwise, y is the spanwise, and z is the vertical direction, respectively. A buffer zone is established in the upstream end of the computational domain in which the terrain irregularities are reduced by 97% to form flat terrain. Similarly, a buffer zone is added to the downstream end of the domain. The terrain elevation data were created with a spatial resolution of 5 m using contour lines from the basis map information of the Geospatial Information Authority of Japan (map scale 1:25,000). The total number of computational grid points, including those in the buffer zones in the upstream and downstream ends of the computational domain, are 391 (x) × 301 (y) × 41 (z), that is, approximately 4.8 million ( Figure 5). The grid spacing is uniformly 5 m in both the x-and y-directions. The vertical grid spacing decreases smoothly down to 0.35 m at the ground surface.
The wind direction set in the present simulation is south-easterly, which is a wind direction associated with confirmed abnormalities of wind turbines at this site. The wind velocity profile applied at the inflow boundary is based on a commonly used empirical power law ( Figure 6). A power law index is set to 5. Furthermore, at the side and upper boundaries, slip boundary conditions are applied. The convective outflow condition is applied at the outflow boundary. Non-slip boundary conditions are applied at the ground surface. The present LES is assumed to reproduce the wind tunnel testing. Therefore, the non-dimensional parameter, Re, in Equation (2) is the Reynolds number (=Uin h/ν) and is set to 10 4 . Figure 7 illustrates the characteristic scales relevant in the present simulation: h is the difference between the minimum and maximum ground surface elevations within the computational domain, Uin is the wind speed at the inflow boundary at the maximum ground surface elevation within the computational domain, and ν is the coefficient of the dynamic viscosity of air. The time step is set to Ät = 2 × 10 −3 h/Uin. In the present simulation, the effect of vertical thermal stratification (density stratification), which is generally present in the atmosphere, is

Simulation Set-Up
As shown in Figure 4, the dimensions of the computational domain are 2.0 (x) × 1.5 (y) × 1.0 (z) km, where x is the streamwise, y is the spanwise, and z is the vertical direction, respectively. A buffer zone is established in the upstream end of the computational domain in which the terrain irregularities are reduced by 97% to form flat terrain. Similarly, a buffer zone is added to the downstream end of the domain. The terrain elevation data were created with a spatial resolution of 5 m using contour lines from the basis map information of the Geospatial Information Authority of Japan (map scale 1:25,000). The total number of computational grid points, including those in the buffer zones in the upstream and downstream ends of the computational domain, are 391 (x) × 301 (y) × 41 (z), that is, approximately 4.8 million ( Figure 5). The grid spacing is uniformly 5 m in both the xand y-directions. The vertical grid spacing decreases smoothly down to 0.35 m at the ground surface.
The wind direction set in the present simulation is south-easterly, which is a wind direction associated with confirmed abnormalities of wind turbines at this site. The wind velocity profile applied at the inflow boundary is based on a commonly used empirical power law ( Figure 6). A power law index is set to 5. Furthermore, at the side and upper boundaries, slip boundary conditions are applied. The convective outflow conditions are applied at the outflow boundary. Non-slip boundary conditions are applied at the ground surface. The present LES is assumed to reproduce the wind tunnel testing. Therefore, the non-dimensional parameter, Re, in Equation (2) is the Reynolds number (=U in h/ν) and is set to 10 4 . Figure 7 illustrates the characteristic scales relevant in the present simulation: h is the difference between the minimum and maximum ground surface elevations within the computational domain, U in is the wind speed at the inflow boundary at the maximum ground surface elevation within the computational domain, and ν is the coefficient of the dynamic viscosity of air. The time step is set to ∆t = 2 × 10 −3 h/U in . In the present simulation, the effect of vertical thermal stratification (density stratification), which is generally present in the atmosphere, is neglected. Furthermore, the effect of surface roughness is included by reconstructing the irregularities of the terrain surface in high resolution.

Simulation Results and Discussions
In the present study, discussions will focus on comparisons between the airflows at WT1 (Wind Turbine 1), which had the best operational performance of all investigated wind turbines, and the airflows at WT2 (Wind Turbine 2), which is presumed to be affected significantly by terrain-induced turbulence in the case of south-easterly wind. (Wind Turbine 1-Wind Turbine 2 will be abbreviated as WT1-WT2, hereafter). Figure 8 shows the wind velocity vectors at the hub height of the wind turbines (78 m above the ground surface). Specifically, the wind velocity vectors are those from the time-averaged wind field. The time interval used for averaging is 100 to 200 in non-dimensional time. Figure 9 shows wind velocity vectors at both of the investigated wind turbine locations from the same time period used for Figures 8 and 10. The results in Figures 8 and 9 are equivalent to those that would be found by a Reynolds-averaged modeling (RANS). An examination of both results reveals that there exist neither significant wind velocity fluctuations, indicating the presence of terrain-induced turbulence, nor severe deficits in the wind velocity profiles at the locations of the wind turbines. That is, these results would lead to the conclusion that the wind fields at and around the wind turbines are suitable for wind power generation.
However, breakage occurred on the yaw motor and yaw gear of WT2. To accurately simulate the behavior of terrain-induced turbulence, such as that investigated in the present study, the use of the RIAM-COMPACT natural terrain version software, which incorporates a LES model for the unsteady turbulence model, is highly effective. The RIAM-COMPACT software has the ability to output the temporal behaviors of various physical quantities and not just their time-averaged values. Figure 10 shows the temporal variation of (1) the fluctuating component (deviation from the mean value) of the streamwise (x) wind velocity at the hub height of the wind turbines, and (2)

Simulation Results and Discussions
In the present study, discussions will focus on comparisons between the airflows at WT1 (Wind Turbine 1), which had the best operational performance of all investigated wind turbines, and the airflows at WT2 (Wind Turbine 2), which is presumed to be affected significantly by terrain-induced turbulence in the case of south-easterly wind. (Wind Turbine 1-Wind Turbine 2 will be abbreviated as WT1-WT2, hereafter). Figure 8 shows the wind velocity vectors at the hub height of the wind turbines (78 m above the ground surface). Specifically, the wind velocity vectors are those from the time-averaged wind field. The time interval used for averaging is 100 to 200 in non-dimensional time. Figure 9 shows wind velocity vectors at both of the investigated wind turbine locations from the same time period used for Figure 8. The results in Figures 8 and 9 are equivalent to those that would be found by a Reynolds-averaged modeling (RANS). An examination of both results reveals that there exist neither significant wind velocity fluctuations, indicating the presence of terrain-induced turbulence, nor severe deficits in the wind velocity profiles at the locations of the wind turbines. That is, these results would lead to the conclusion that the wind fields at and around the wind turbines are suitable for wind power generation.
However, breakage occurred on the yaw motor and yaw gear of WT2. To accurately simulate the behavior of terrain-induced turbulence, such as that investigated in the present study, the use of the RIAM-COMPACT natural terrain version software, which incorporates a LES model for the unsteady turbulence model, is highly effective. The RIAM-COMPACT software has the ability to output the temporal behaviors of various physical quantities and not just their time-averaged values. Figure 10 shows the temporal variation of (1) the fluctuating component (deviation from the mean value) of the streamwise (x) wind velocity at the hub height of the wind turbines, and (2) the angle of the wind on horizontal and vertical cross-sections. In the figure, the horizontal axes indicate Subsequently, the fluctuating component (deviation from the mean value) of the streamwise (x) wind velocity will be examined (top panel in Figure 10). In contrast to the very small temporal fluctuations of the streamwise wind velocity at WT1, the presence of spike-like temporal fluctuations can be confirmed at WT2. Although it is beyond the scope of the present paper, further examination of temporal wind velocity changes at the lower, upper, left, and right ends of the swept area allows an examination of the moments exerted on the rotor.
An investigation of the temporal variation in the angle of the wind on the horizontal cross-section (middle panel in Figure 10) reveals the following: As in the case of the wind velocity, the amplitude of the fluctuations of the angle of the wind on the horizontal cross-section is quite small at WT1. On the other hand, changes in the wind direction exceeding 25.0 • are frequently seen at WT2 (see also Figure 11). Furthermore, the fact that many spikes occur in the direction of positive values suggests that easterly winds are generated on very short time scales due to the topographic effects.  Figure 10 also shows the definitions of the angles of the wind on the horizontal and vertical cross-sections. An inspection of the undulating patterns of the time series in Figure 10 allows a clear understanding of the change in the unsteady wind conditions occurring near the wind turbines. Subsequently, the fluctuating component (deviation from the mean value) of the streamwise (x) wind velocity will be examined (top panel in Figure 10). In contrast to the very small temporal fluctuations of the streamwise wind velocity at WT1, the presence of spike-like temporal fluctuations can be confirmed at WT2. Although it is beyond the scope of the present paper, further examination of temporal wind velocity changes at the lower, upper, left, and right ends of the swept area allows an examination of the moments exerted on the rotor.
An investigation of the temporal variation in the angle of the wind on the horizontal cross-section (middle panel in Figure 10) reveals the following: As in the case of the wind velocity, the amplitude of the fluctuations of the angle of the wind on the horizontal cross-section is quite small at WT1. On the other hand, changes in the wind direction exceeding 25.0° are frequently seen at WT2 (see also Figure 11). Furthermore, the fact that many spikes occur in the direction of positive values suggests that easterly winds are generated on very short time scales due to the topographic effects.   Figure 10).
Finally, the temporal variation in the angle of the wind on the vertical cross-section was examined (bottom panel in Figure 10). This figure also suggests that wind blowing upward and downward, exceeding 25.0° from the horizontal, is generated. It can be speculated that the minor landforms located upwind of the wind turbine contribute to these local wind direction changes. Figures 11 and 12 illustrate the wind field (wind velocity vectors) from times characterized by the large wind direction changes, indicated by "A" and "B" in Figure 10. An examination of Figure  11 reveals that the direction of the wind aloft at the location of WT2 (easterly wind) deviates significantly from the direction of the prevailing wind. Furthermore, an inspection of Figure 12, which shows the wind field from a different time from that in Figure 11, shows that large upward flow has been generated near the center of the hub of WT2. In addition, the vertical profiles of the streamwise wind velocity in Figure 12 do not follow the so-called wind profile power law; a large velocity deficit can be seen between the hub center and the lower end of the swept area.
In general, the power curve of a wind turbine (the catalog value) is specified using the velocity of the wind flowing into the hub center of the wind turbine on flat terrain, assuming the absence of the turbine itself. Furthermore, the vertical profile of the mean horizontal wind speed, which is used to evaluate the wind shear, is assumed to follow a power law, with a coefficient of approximately 5 to 7. Therefore, a large reduction in the generated electric power is expected in the presence of wind shear, which deviates significantly from the wind shear predicted using the profile based on the power law. Furthermore, such anomalous wind shear will likely be an increasingly important topic of research in the future in connection with the issues of wind turbine tower vibration and the fatigue strength of yaw gears. Figure 13 shows the vertical profiles of selected physical quantities at each wind turbine location. Figure 13a shows the vertical profile of the mean streamwise wind velocity. Figure 13b shows the vertical profile of the standard deviation of the streamwise velocity (the horizontal axis indicates the standard derivation normalized by the inflow wind velocity aloft, Uin). Figure 13c shows the turbulence intensity (%) evaluated by dividing the standard deviation of the streamwise velocity by the mean streamwise wind velocity, <u>, at each location. In all these figures, the vertical axes indicate the height above the ground (m). In these figures, the vertical range of the wind turbine swept area (the rotor diameter) is also illustrated. In the vertical profile of the mean streamwise wind velocity in Figure 13a, no severe wind velocity deficit is identified at any of the wind turbines, except for WT2. An inspection of the profiles in the vertical range of the wind turbines in Figures 13b,c reveals that the effects of terrain-induced turbulence are almost similar in magnitude at WT3 and WT4, and that, at WT2, it is extremely large due to the presence of the minor landforms upstream of these turbines. The effect of the terrain-induced turbulence is notably smaller at WT1 than at the other three wind turbines.  Figure 10).
Finally, the temporal variation in the angle of the wind on the vertical cross-section was examined (bottom panel in Figure 10). This figure also suggests that wind blowing upward and downward, exceeding 25.0 • from the horizontal, is generated. It can be speculated that the minor landforms located upwind of the wind turbine contribute to these local wind direction changes. Figures 11 and 12 illustrate the wind field (wind velocity vectors) from times characterized by the large wind direction changes, indicated by "A" and "B" in Figure 10. An examination of Figure 11 reveals that the direction of the wind aloft at the location of WT2 (easterly wind) deviates significantly from the direction of the prevailing wind. Furthermore, an inspection of Figure 12, which shows the wind field from a different time from that in Figure 11, shows that large upward flow has been generated near the center of the hub of WT2. In addition, the vertical profiles of the streamwise wind velocity in Figure 12 do not follow the so-called wind profile power law; a large velocity deficit can be seen between the hub center and the lower end of the swept area.
In general, the power curve of a wind turbine (the catalog value) is specified using the velocity of the wind flowing into the hub center of the wind turbine on flat terrain, assuming the absence of the turbine itself. Furthermore, the vertical profile of the mean horizontal wind speed, which is used to evaluate the wind shear, is assumed to follow a power law, with a coefficient of approximately 5 to 7. Therefore, a large reduction in the generated electric power is expected in the presence of wind shear, which deviates significantly from the wind shear predicted using the profile based on the power law. Furthermore, such anomalous wind shear will likely be an increasingly important topic of research in the future in connection with the issues of wind turbine tower vibration and the fatigue strength of yaw gears. Figure 13 shows the vertical profiles of selected physical quantities at each wind turbine location. Figure 13a shows the vertical profile of the mean streamwise wind velocity. Figure 13b-d shows the vertical profiles of the standard deviations of the three wind velocity components. In all these figures, the vertical axes indicate the height above the ground (m), and the horizontal axis is normalized by the inflow wind velocity aloft, U in . In addition, the vertical range of the wind turbine swept area (the rotor diameter) is also illustrated. In the vertical profile of the mean streamwise wind velocity in Figure 13a, no severe wind velocity deficit is identified at any of the wind turbines, except for WT2. An inspection of the profiles in the vertical range of the wind turbines in Figure 13b-d reveals that the effects of terrain-induced turbulence are almost similar in magnitude at WT3 and WT4, and that, at WT2, it is extremely large due to the presence of the minor landforms upstream of these turbines. The effect of the terrain-induced turbulence is notably smaller at WT1 than at the other three wind turbines.     To assess three-dimensional wind conditions, it is effective to release virtual particles and observe their trajectories (pathlines). Figures 14 and 15 show the results obtained using this method. In Figure 14, virtual particles are released between the right and left edges of the swept area at the height of the center of the wind turbine hub. On the other hand, virtual particles are released between the upper and lower edges of the swept area at the y-coordinate of the center of the wind turbine hub, as shown in Figure 15. As discussed thus far, in the vicinity of WT2, the airflow meanders significantly and three dimensionally in the horizontal and vertical directions due to the effects of the minor landforms located upstream of this wind turbine (Figures 14 and 15).
Finally, it should be recalled that the present study examines only the turbulence generated by the irregularities of the terrain surface and does not take into account the effects of the turbulence caused by gusts. To assess three-dimensional wind conditions, it is effective to release virtual particles and observe their trajectories (pathlines). Figures 14 and 15 show the results obtained using this method. In Figure 14, virtual particles are released between the right and left edges of the swept area at the height of the center of the wind turbine hub. On the other hand, virtual particles are released between the upper and lower edges of the swept area at the y-coordinate of the center of the wind turbine hub, as shown in Figure 15. As discussed thus far, in the vicinity of WT2, the airflow meanders significantly and three dimensionally in the horizontal and vertical directions due to the effects of the minor landforms located upstream of this wind turbine (Figures 14 and 15).
Finally, it should be recalled that the present study examines only the turbulence generated by the irregularities of the terrain surface and does not take into account the effects of the turbulence caused by gusts.

The Economic Effects
At the Atsumi Wind Farm in the Aichi Prefecture, Japan (start date of operation: 17 January 2007), the terrain-induced turbulence affected the wind turbines in the form of worn brake pads for yaw motors and fractured planetary gears both due to frequent yaw controls and, also, in the form of hydraulic system failures due to frequent pitch controls.
In the present study, the numerical wind simulation was conducted by reproducing the topography near the wind turbine sites in high resolution and using the Large-Eddy Simulation (LES) technique. The topography near the wind turbine sites serves as the origin of the terrain-induced turbulence. The obtained numerical results showed that the terrain-induced turbulence is generated at a small terrain feature located upstream of WT2 in the case of south-easterly wind. The generated turbulence affects the wind turbine directly. The wind speed and wind direction at the wind turbine site are significantly changed with time.
At the Atsumi Wind Farm, a combination of the series of wind simulation results and on-site operation experience led to a decision to adopt an "automatic shutdown program" for WT1 and WT2 (introduced on 23 May 2011; modified on 26 July 2011). Here, "automatic shutdown program" means the automatic suspension of wind turbine operation based on the wind speed and wind direction meeting the conditions associated with significant effects of terrain-induced turbulence at a wind turbine site. Specifically, after modifications, the program was configured such that WT1

The Economic Effects
At the Atsumi Wind Farm in the Aichi Prefecture, Japan (start date of operation: 17 January 2007), the terrain-induced turbulence affected the wind turbines in the form of worn brake pads for yaw motors and fractured planetary gears both due to frequent yaw controls and, also, in the form of hydraulic system failures due to frequent pitch controls.
In the present study, the numerical wind simulation was conducted by reproducing the topography near the wind turbine sites in high resolution and using the Large-Eddy Simulation (LES) technique. The topography near the wind turbine sites serves as the origin of the terrain-induced turbulence. The obtained numerical results showed that the terrain-induced turbulence is generated at a small terrain feature located upstream of WT2 in the case of south-easterly wind. The generated turbulence affects the wind turbine directly. The wind speed and wind direction at the wind turbine site are significantly changed with time.
At the Atsumi Wind Farm, a combination of the series of wind simulation results and on-site operation experience led to a decision to adopt an "automatic shutdown program" for WT1 and WT2 (introduced on 23 May 2011; modified on 26 July 2011). Here, "automatic shutdown program" means the automatic suspension of wind turbine operation based on the wind speed and wind direction meeting the conditions associated with significant effects of terrain-induced turbulence at a wind turbine site. Specifically, after modifications, the program was configured such that WT1

The Economic Effects
At the Atsumi Wind Farm in the Aichi Prefecture, Japan (start date of operation: 17 January 2007), the terrain-induced turbulence affected the wind turbines in the form of worn brake pads for yaw motors and fractured planetary gears both due to frequent yaw controls and, also, in the form of hydraulic system failures due to frequent pitch controls.
In the present study, the numerical wind simulation was conducted by reproducing the topography near the wind turbine sites in high resolution and using the Large-Eddy Simulation (LES) technique. The topography near the wind turbine sites serves as the origin of the terrain-induced turbulence. The obtained numerical results showed that the terrain-induced turbulence is generated at a small terrain feature located upstream of WT2 in the case of south-easterly wind. The generated turbulence affects the wind turbine directly. The wind speed and wind direction at the wind turbine site are significantly changed with time.
At the Atsumi Wind Farm, a combination of the series of wind simulation results and on-site operation experience led to a decision to adopt an "automatic shutdown program" for WT1 and WT2 (introduced on 23 May 2011; modified on 26 July 2011). Here, "automatic shutdown program" means the automatic suspension of wind turbine operation based on the wind speed and wind direction meeting the conditions associated with significant effects of terrain-induced turbulence at a wind turbine site. Specifically, after modifications, the program was configured such that WT1 would stop operation for the wind direction range of 130 • to 190 • , with wind speeds of 9.0 m/s or higher, and that WT2 would stop operation for the wind direction range of 90 • to 130 • , with wind speeds of 10.0 m/s or higher (see Figure 16). The adoption of the "automatic shutdown program" has successfully led to a large reduction in the number of occurrences of wind turbine damage, thus, creating major positive economic effects. To evaluate the economic effects, the performance of the wind turbines from the two years preceding the adoption of the "automatic shutdown program" was used as a reference and changes in the capacity and availability factors were calculated from the reduction in downtime of the wind turbines from the two years following the adoption of the "automatic shutdown program." The following major economic effects have been realized based on converting the performance of the two wind turbines combined from the two years preceding and the two years following the adoption of the "automatic shutdown program" into the performance of a single wind turbine per year: (1) A reduction in the repair costs by 9.322 million yen per year per wind turbine, (2) An increase in the availability factor by 8.05% (87.3%→95.4%), and (3) An increase in the capacity factor by 1.7% (21.2%→22.9%).
The number of occurrences of damage to wind turbines can be reduced significantly by combining "quality control of local wind conditions" by RIAM-COMPACT, that is, numerical wind diagnoses (analyses of terrain-induced turbulence), with appropriate wind turbine control based on the results of the diagnoses. At the same time, this combination of approaches also results in significant positive economic effects. In the future development of the wind power industry in Japan, it is essential to adopt "quality control of local wind conditions" for each of the country's land-based wind turbines. would stop operation for the wind direction range of 130° to 190°, with wind speeds of 9.0 m/s or higher, and that WT2 would stop operation for the wind direction range of 90° to 130°, with wind speeds of 10.0 m/s or higher (see Figure 16). The adoption of the "automatic shutdown program" has successfully led to a large reduction in the number of occurrences of wind turbine damage, in addition to, creating positive economic effects. To evaluate the economic effects, the performance of the wind turbines from the two years preceding the adoption of the "automatic shutdown program" was used as a reference and changes in the capacity and availability factors were calculated from the reduction in downtime of the wind turbines from the two years following the adoption of the "automatic shutdown program." The following major economic effects have been realized based on converting the performance of the two wind turbines combined from the two years preceding and the two years following the adoption of the "automatic shutdown program" into the performance of a single wind turbine per year: (1) A reduction in the repair costs by 9.322 million yen per year per wind turbine, (2) An increase in the availability factor by 8.05% (87.3%→95.4%), and (3) An increase in the capacity factor by 1.7% (21.2%→22.9%).
The number of occurrences of damage to wind turbines can be reduced significantly by combining suitability inspections of local wind conditions" by RIAM-COMPACT, that is, numerical wind simulations (analyses of terrain-induced turbulence), with appropriate wind turbine control based on the results of the diagnoses. At the same time, this combination of approaches also results in significant positive economic effects. In the future development of the wind power industry in Japan, it is essential to adopt suitability inspections of local wind conditions" for each of the country's land-based wind turbines.

Conclusions
The wind conditions at the Atsumi Wind Farm in the Aichi Prefecture, Japan (in operation since March 2007) were numerically predicted using the RIAM-COMPACT natural terrain version software. The simulation results showed that the minor landforms located upstream of WT2 serve as the origin of turbulence (terrain-induced turbulence), in the case of south-easterly wind, and that WT2 is strongly affected by this turbulence.

Conclusions
The wind conditions at the Atsumi Wind Farm in the Aichi Prefecture, Japan (in operation since March 2007) were numerically predicted using the RIAM-COMPACT natural terrain version software. The simulation results showed that the minor landforms located upstream of WT2 serve as the origin of turbulence (terrain-induced turbulence), in the case of south-easterly wind, and that WT2 is strongly affected by this turbulence.
The use of time-averaged RANS turbulence models creates severe difficulties in simulating the effects of terrain-induced turbulence generated locally due to minor landforms and small topographical irregularities distributed sparsely around wind turbines. As a result, the use of RANS models tends to miss the presence of such effects. Therefore, it is effective to adopt an approach based on an unsteady turbulence model (e.g., LES), such as the RIAM-COMPACT.
In addition, relative comparisons of the effects of terrain-induced turbulence among wind turbines, such as that presented in the present study, can be made appropriately using only detailed terrain data accurately representing the land formations at a site of interest, without the use of any meteorological data. If the results from such comparisons were to be effectively applied to wind turbine operation control, the number of internal and external breakdowns of wind turbines caused by terrain-induced turbulence would decrease dramatically, which would likely allow the availability factors of wind turbines to increase significantly. For example, meteorological conditions in which large wind direction and speed changes occur so rapidly that they cannot be adequately counteracted by the pitch angle control of the wind turbine blades (i.e., specific wind directions in which such meteorological conditions occur) can be predicted in advance for the position of each nacelle. When the actual meteorological conditions meet the previously identified meteorological conditions, measures, such as halting wind power generation, can be taken.
In the future, regardless of whether the terrain is flat or complex, it is necessary; (1) to conduct a detailed wind synopsis diagnosis, such as the one presented in this paper, for all wind directions around a wind turbine; (2) to accurately understand the three-dimensional local wind conditions at each wind turbine location and the temporal change of the vertical profiles of the wind velocities at wind turbine locations; and (3) to carefully consider the wave-like structures in the time series of the wind velocities and the turbulence intensities at the lower, upper, left, and right ends of the swept area.