The author et al. have been developing a nonstationary and nonlinear wind-synopsis simulator (large-eddy simulation (LES) turbulence model) applicable to analyses of wind synopsis and the diffusion fields of any complex and steep terrain in the world which is called the Research Institute for Applied Mechanics, Kyushu University, COM putational Prediction of Airflow over Complex Terrain (RIAM-COMPACT)  -  . The analysis space is focused on a narrow range from several hundred meters to several kilometers, and it is possible to numerically predict the topographic effects of winds with high accuracy, such as flow separation and the formation of reverse-flow regions caused by that, local increase of flow speed, and the reattachment of a separated shear layer.
Typhoon Tokage (0423) approached Kyushu, Japan, on 20 October 2004. Even though it did not make landfall in Kyushu, the range of its strong winds covered a wide area, and Ogimachi (Ogi City at present, hereafter referred to as Ogimachi) in Saga Prefecture suffered extensive damage (Figure 1). The direction of the wind of Typhoon 0423 at the time when it caused the damage was from northwest to north-northwest, and if the geographical features surrounding the area are considered, the influence of the Tenzan Mountains located on the north side of Ogimachi is supposed to have been great. In other words, it is supposed that the strong winds of the typhoon were further strengthened due to local geographical features and, as a result, the damage from the strong winds was greatly amplified  -  . In order to clarify the mechanism of the occurrence of strong winds accompanied by the effects of terrain, the use of the numerical simulation is very effective, in which conditions can be set without the influence of ground roughness and temperature distribution.
In this research, two tasks were investigated and discussed using the RIAM-COMPACT method described above. A task was to conduct a numerical simulation in a broad range targeting the north-northwest wind direction in order to clarify the mechanism of the strong-wind phenomenon that occurred around Ogimachi. Another task was to investigate the reproducibility of regional-scale terrains, in other words, the effect of horizontal grid resolution on calculation results.
The structure of this paper is as follows. Chapter 2 explains the numerical calculation method. Chapter 3 explains wind-velocity data record when Typhoon 0423 was passing through. Chapter 4 explains calculation results and discussion. Chapter 5 explains the effect of horizontal grid resolution on calculation results. Finally, Chapter 6 is the conclusion.
2. LES by Finite-Difference Method Based on Collocated Grid of Generalized Curvilinear Coordinate System
The numerical calculation method, calculation areas, and boundary conditions are described here. In order to numerically predict wind flows over a complex terrain with a high degree of accuracy, avoiding unstable oscillations, a collocated grid of generalized curvilinear coordinates was used. The collocated grid here is a grid system in which physical velocity components and pressure were defined in the cell centers of the calculation grid, and the variables of contravariant velocity components multiplied by Jacobian are defined in the cell interface. The numerical calculation method was conducted based on the finite-difference method (FDM), and the large-eddy simulation (LES) was conducted. For the governing equations of the flow, a filtered continuity equation for incompressible fluid (Equation (1)) and a filtered Navier-Stokes equation (Equation (2)) were used.
Figure 1. Disasters caused by strong winds targeted by this research.
The calculation algorithm was based on the fractional step (F-S) method  , and the time-marching method was based on the Euler explicit method. Poisson’s equation was solved by the successive over-relaxation (SOR) method. All discretizations of the terms of space were based on a second-order central-difference scheme except for the convection term of Equation (2). The term of convection was based on third-order upwind difference scheme by Kajishima  . For the numerical diffusion coefficient due to third-order upwind difference scheme, 0.5 is used as opposed to 3.0 from the Kawamura-Kuwahara scheme  in order to minimize the influence of numerical diffusion. The standard Smagorinsky model was used for the SGS model  . The model coefficient was assumed to be 0.1 by using a wall-damping function.
The calculation area and boundary conditions are shown in Figure 2. The horizontal calculation area is 16 km2 from the center of Ogimachi and the vertical calculation area is 5 km. The topographic geometry data with 100 m horizontal spatial resolution was constructed based on the Geospatial Information Authority of Japan (GSI) with 50 m horizontal spatial resolution. The number of grid points to each direction is 161 × 161 × 51 points in the main-flow-direction (x), spanwise direction (y), and vertical direction (z), respectively. The grid widths of directions x and y were approximately uniform intervals, with a horizontal resolution of approximately 100 m. The grid width was nonuniform in the z direction so that the density of grid points increased smoothly toward the ground surface (Δzmin = 3 m). The targeted wind direction was north–northwest. The inflow boundary surface was given with a velocity profile following the 1/7th power law. The upper boundary surface and the side boundary surface were given with slip conditions, and the outflow boundary surface was given with convective outflow conditions. The ground surface was given with no-slip conditions. The nondimensional parameter Re of Equation (2) is a Reynolds number (= Uin h/ν). 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. In addition, as in Uchida   , the effects of the surface roughness were taken into consideration by reconstructing surface irregularities in high resolution. The comparison between RANS results and the present LES results are summarized in the latest article  , and the prediction accuracy of the present LES approach by comparison with wind tunnel experiments is discussed in the article  . The nondimensional parameter Re was assumed to be 104. Here, h is the altitude difference in the calculation domain (h = 1029 m), Uin is the wind velocity at the inflow boundary surface at the height of the maximum terrain elevation within the calculation domain and ν is a kinematic coefficient of viscosity. The nondimensional time increment was assumed to be t = 2 × 10−3 h/Uin.
The Ogimachi and Network for Wind Measurement in Kyushu (NeWMeK) observation sites (Nos. 22 and 27) and Saga University (hereafter referred as to Saga) are shown in Figure 2. NeWMeK is a wide-area high-density wind-observation system of Kyushu Electric Co., Ltd. This system can capture the distribution characteristics of strong winds, such as the process of generation and dissipation of local winds in a wide area in real time. This system has 123 observation stations using electrical-transmission steel towers in all areas of Kyushu and can continuously and simultaneously record data in one-second
Figure 2. Calculation area, boundary conditions and computational grid.
intervals. The No. 22 and 27 NeWNeK observation sites observed strong winds when a typhoon passed through. This is described later. Saga was used as the site for comparison with Ogimachi.
3. Wind-Velocity Data Record When Typhoon 0423 Was Passing Through
Table 1 shows wind-velocity data recorded at NeWMeK observation sites 22 and 27 at the time when Typhoon 0423 was passing through for 10 minutes, from 16:30 to 16:40, 20 October 2004. If it is taken into consideration that the anemometer of No. 27 was installed 18 m above the ground, it can be expected that strong winds extensively occurred near the ground during this time period.
Table 1. Wind-velocity data recorded at the Network for Wind Measurement in Kyushu (NeWMeK) observation sites.
4. Calculation Results and Discussion
4.1. Visualization of Flow and Mechanism for Strong-Wind Generation
In this section, flow visualization in the neighborhood of Ogimachi is shown, and the mechanism of strong-wind generation in the neighborhood of Ogimachi is discussed. Velocity-vector diagrams (instantaneous flow field) at 20 m above the ground are shown in Figure 3. Figure 3(a) and Figure 3(b) shows numerical results at different times. Figure 4 shows the contour diagram (instantaneous flow field) of the main-flow-direction (x) velocity components corresponding to the time of Figure 3. The range of contour here is 0.6 ≤ (u/Uin) ≤ 1.3, and this range is shown by being divided into 30. Therefore, the areas in which contours are shown mean that they locally have stronger winds than the surrounding areas. Pay attention to Ogimachi and Saga in Figure 3. Almost no strong-wind areas were generated in the areas surrounding Saga. Contrarily, it can be found that Ogimachi had strong winds, the same as those of the areas located high above sea level on the upstream side. Furthermore, this can be clearly understood if you scrutinize it together with Figure 4, which shows that speed-increasing areas always generate in the neighborhood of Ogimachi. This was confirmed in the flow animation. What is interesting is that the local increasing-speed area extended in a stripe shape from the side of the upper stream, as shown in Figure 4 (the area surrounded with the dotted line). This suggests that speed-increasing flow is generated in the upper stream of Ogimachi, and it flows down to Ogimachi without decreasing speed.
It was clarified that the generation sources of speed-increasing areas in the surroundings of Ogimachi are in the areas indicated with the solid lines in Figure 3 and Figure 4. Therefore, this part is enlarged and shown in Figure 5. The result of superposing a streamline on the contour diagram of the main-flow-direction (x) velocity components (instantaneous flow field), shown in Figure 4, is shown in Figure 5. Pay attention to the streamline in which virtual particles are released from a height of 20 m above the ground and the particle trajectories are drawn. In the mountain terrain of the upper stream of Ogimachi, there is a small basin in the neighborhood of the area surrounded with solid lines. As a result, the flows converge can be captured. In other words, the strong winds in the neighborhood of Ogimachi are presumably caused because the flows of wind converge at the small basin in the mountain terrain in the
Figure 3. Velocity vector diagrams (instantaneous flow field) at 20 m above the ground. (a) Nondimensional time = T1; (b) Nondimensional time = T1 + Δt.
Figure 4. Contour diagram of main-stream-direction (x) velocity components 20 m above ground. 0.6 ≤ (u/Uin) ≤ 1.3 is divided into 30. (a) Nondimensional time = T1; (b) Nondimensional time = T1 + Δt.
upper stream, which leads to gap flows, so to speak. Furthermore, as shown in Figure 6, it is considered that the reason why the speed-increasing flow did not slow down so much is that the terrain (wind path) is easy for wind to go through.
4.2. Time-History Waveform and Vertical Profile of Various Quantities of Turbulent Airflow
In this section, more quantitative discussions are outlined. Figure 7 shows the
Figure 5. Streamline and contour diagram of main-stream-direction (x) velocity components drawn from 20 m above ground. 0.6 ≤ (u/Uin) ≤ 1.3 is divided into 30.
Figure 6. Mechanism of local strong winds generating in the surrounding areas of Ogimachi.
time variations (20 and 30 m above ground) of real-scale velocity components of main-flow-direction (x). For the purpose of comparison, the time-history waveform of Saga (20 and 30 m above ground) is also shown with it. The way to convert the nondimensional wind-velocity value that was output from the numerical simulation to a realistic scale is explained. In this calculation, wind velocities and times are normalized, as shown in Equations (3) and (4), based on the reference scales.
Normalization of Wind Velocities and Times
Figure 7. Time-history waveforms in real scale of main-stream-direction (x) velocity components in the points of Ogimachi and Saga.
Superscript (−) to which the spatial filter was applied is omitted as a matter of convenience. In order to convert wind velocity and time to real-scale values, you only have to substitute the values of specific reference scale h (m) and Uin (m/s). As for wind velocity, Uin (m/s) was set so that the nondimensional average wind-velocity value corresponding to the No. 22 (51 m above ground) NeWMeK observation site coincides with the average wind velocity (18.1 m/s) for ten minutes from 16:30 to 16:40, 20 October 2004, shown in Table 1. This was caused because it was supposed that the effects of surrounding terrains were small because the No. 22 NeWMeK observation site is located at the side of the upper stream of the calculation area and relatively high above ground at 51 m.
As shown in Figure 7, no high wind velocity emerged at the Saga test points from 20 to 30 m above ground during the weak-wind condition. On the contrary, average wind velocity from 20 to 30 m high above ground at the point of Ogimachi was greater than Saga. It was also observed that large wind-velocity values often emerged in the shape of a spike. Average wind velocity obtained from the time-history waveform of Figure 7 is shown in Table 2.
The average speed profile of main-flow-direction (x) and the vertical profile of standard deviation are shown in Figure 8 and Figure 9, respectively. The vertical axis is z* (m) above ground in real scale, and the horizontal axis is a nondimensional value normalized by inflow wind velocity Uin. In Figure 9, it can be seen that the standard deviation for each mesh point is not the turbulence intensity divided by average wind velocity of that point. In both Figure 8 and Figure 9, the position of 20 m above ground is shown. If you pay attention to the average velocity profile shown in Figure 8, you can notice that winds about three times stronger occur around Ogimachi than in Saga at 20 m above the ground. In the
Table 2. Average values in real scale of main-flow-direction (x) velocity components in the points of Ogimachi and Saga.
Figure 8. Average velocity profile in the points of Ogimachi and Saga, main-flow (x) velocity components.
Figure 9. Vertical profile of standard deviation in the points of Ogimachi and Saga, main-flow (x) velocity components.
vertical profile of the standard deviation shown in Figure 9, values are clearly larger in Ogimachi than those in Saga. Therefore, as mentioned in Figure 7, it became clear that greatly fluctuating and very strong winds occurred in Ogimachi.
5. Grid-Resolution Effects
It was shown in the previous section that subtle changes of topographic reliefs become a source of topographic wind generation (local strong wind, gap flow in this research). Therefore, in order to reproduce such flow phenomena using a numerical simulation, it is very important to clarify how accurately small terrains should be reproduced, in other words, to what extent is horizontal grid resolution required. In this research, three cases of numerical simulations with different grid dissolutions were conducted, and the above problems are discussed.
The calculation area is the area surrounded by a solid line, as shown in Figure 1. This area has a horizontal space of 10 km to the main-stream-direction (x), a horizontal space of 6 km to the orthogonal direction to the main-stream-direction (y), and a space of 4 km was set up in the vertical direction (z). Three cases of dissolution, shown in Table 3, are discussed. The number of meshes in the vertical direction (z) of each of the three cases is 51, and their grid dissolutions are also the same. The first cell from the ground, the minimum mesh width, is approximately 2.4 - 2.5 m.
Figure 10 shows the comparison of the reproducibility of surrounding areas of topographic-strong-wind (local-strong-wind) generation sources. No significant difference was observed between the 50 m horizontal grid resolution of Figure 10(a) and the 100 m horizontal grid resolution of Figure 10(b) (resolution by the calculation of the wide area up to the previous section). Meanwhile, the result of the 200 m horizontal grid resolution of Figure 10(c) shows noticeable differences in topographic reliefs compared with those of Figure 10(a) and Figure 10(b). The obtained calculation result by each grid resolution is shown below.
Figure 11 shows the contour diagram (instantaneous flow field) of the main-stream-direction (x) velocity component 20 m above ground. The range of the contour is shown by dividing 0.6 ≤ (u/Uin) ≤ 1.3 into 30 in the same way as before. The speed-increasing flow, from a small basin of mountain terrain in the upper stream of Ogimachi, can be observed in the resolution shown in
Table 3. Grid resolutions of three discussed cases.
Figure 10. Comparison of terrain reproducibility. (a) Case 1, Horizontal grid resolution 50 m; (b) Case 2, Horizontal grid resolution 100 m; (c) Case 3, Horizontal grid resolution 200 m.
Figure 11. Contour diagram of main-stream-direction (x) velocity component 20 m above ground. 0.6 ≤ (u/Uin) ≤ 1.3 is divided into 30. (a) Case 1, Horizontal grid resolution 50 m; (b) Case 2, Horizontal grid resolution 100 m; (c) Case 3, Horizontal grid resolution 200 m.
Figure 11(a) (horizontal grid resolution; 50 m) and Figure 11(b) (horizontal grid resolution; 100 m), the same as Figure 4. Meanwhile, as the result of the horizontal grid resolution at 200 m shown in Figure 11(c), local winds from the valley were almost not reproduced.
Figure 12. Time-history waveform of main-flow-direction (x) velocity component, nondimensional scale. (a) Case 1, Horizontal grid resolution 50 m; (b) Case 2, Horizontal grid resolution 100 m; (c) Case 3, Horizontal grid resolution 200 m.
The effects of the grid resolution can be understood more clearly by showing the time-history waveform of the speed-variation component. Figure 12 shows the time changes of main-stream-direction (x) speed-variation components, u’ (= y – Uave) at 200 m above ground. Both the vertical and the horizontal axis are on a nondimensional scale. The time scale of the horizontal axis almost corresponds to the real scale of Figure 7. Large-fluctuation velocities (variable amplitudes), as shown in Figure 7, can be observed by the 50 m resolution shown in Figure 12(a), and the 100 m resolution shown in Figure 12(b). It is revealed that high-frequency components are reproduced more noticeably in the 50 m resolution shown in Figure 12(a) in particular, along with the increase of grid resolution. Contrarily, as a result of the 200 m resolution shown in Figure 12(c), almost all of the speed variations due to the passing of speed-increasing flow are not captured, but it looks like they are smoothed.
The mechanism of the generation of strong winds generated in Ogimachi, Saga Prefecture caused by Typhoon 0423 was investigated using a nonstationary and nonlinear wind-synopsis simulator, RIAMM-COMPACT. Furthermore, a discussion was conducted about how the difference of terrain reproducibility, grid resolution, which is required for the reproduction of these topographic strong winds, affects the calculation result. The findings obtained in this research are as follows:
1) The direction of the wind at the time when the damage was caused by Typhoon 0423 ranged from northeast to north-northwest, and it was shown that the Tenzan Mountains located in the north side of Ogimachi had an influence on that.
2) It was observed that there is a small basin in the mountain terrain in the upper stream of Ogimachi and all flows were converging there. In other words, it became clear that it is highly possible that the strong winds that were generated in the areas surrounding Ogimachi were caused by wind flows converging to a small basin of mountain terrain in the upper stream and formed a gap flow.
3) It was shown that a horizontal resolution of 50 - 100 m is required for the reproducibility of local-scale terrain.
This work is an update of the previous conference abstract  . The author expresses appreciation to the person concerned at that time.
This work was supported by JSPS KAKENHI Grant Number 17H02053.
 Uchida, T. (2018) Large-Eddy Simulation of Airflow over a Steep, Three-Dimensional Isolated Hill with Multi-GPUs Computing. Open Journal of Fluid Dynamics, 8, 416-434.
 Uchida, T. and Li, G. (2018) Comparison of RANS and LES in the Prediction of Airflow Field over Steep Complex Terrain. Open Journal of Fluid Dynamics, 8, 286-307.
 Uchida, T. (2018) Computational Fluid Dynamics (CFD) Investigation of Wind Turbine Nacelle Separation Accident over Complex Terrain in Japan. Energies, 11, 1485.
 Uchida, T. (2017) Three-Dimensional Numerical Simulation of Stably Stratified Flows over a Two-Dimensional Hill. Open Journal of Fluid Dynamics, 7, 579-595.
 Uchida, T. (2017) High-Resolution LES of Terrain-Induced Turbulence around Wind Turbine Generators by Using Turbulent Inflow Boundary Conditions. Open Journal of Fluid Dynamics, 7, 511-524.
 Uchida, T. (2017) CFD Prediction of the Airflow at a Large-Scale Wind Farm above a Steep, Three-Dimensional Escarpment. Energy and Power Engineering, 9, 829-842.
 Uchida, T. (2017) High-Resolution Micro-Siting Technique for Large Scale Wind Farm Outside of Japan Using LES Turbulence Model. Energy and Power Engineering, 9, 802-813.
 Uchida, T. and Ohya, Y. (2008) Micro-Siting Technique for Wind Turbine Generators by Using Large-Eddy Simulation. Journal of Wind Engineering & Industrial Aerodynamics, 96, 2121-2138.
 Uchida, T. and Ohya, Y. (2006) Application of LES Technique to Diagnosis of Wind Farm by Using High Resolution Elevation Data. JSME International Journal Series B, 49, 567-575.
 Uchida, T. and Ohya, Y. (2003) Large-Eddy Simulation of Turbulent Airflow over Complex Terrain. Journal of Wind Engineering & Industrial Aerodynamics, 91, 219-229.
 Uchida, T. and Ohya, Y. (1999) Numerical Simulation of Atmospheric Flow over Complex Terrain. Journal of Wind Engineering & Industrial Aerodynamics, 81, 283-293.
 Yang, Z., Calderer, A., He, S., et al. (2019) Measurement-Based Numerical Study of the Effects of Realistic Land Topography and Stratification on the Coastal Marine Atmospheric Surface Layer. Boundary-Layer Meteorology, 171, 289-314.
 Yang, Z., Calderer, A., He, S., Sotiropoulos, F., Doyle, J.D., Flagg, D.D., MacMahan, J., Wang, Q., Haus, B.K., Graber, H.C. and Shen, L. (2018) Numerical Study on the Effect of Air-Sea-Land Interaction on the Atmospheric Boundary Layer in Coastal Area. Atmosphere, 9, 51.
 Gargallo-Peiró, A., Avila, M., Owen, H., Prieto-Godino, L. and Folch, A. (2018) Mesh Generation, Sizing and Convergence for Onshore and Offshore Wind Farm Atmospheric Boundary Layer Flow Simulation with Actuator Discs. Journal of Computational Physics, 375, 209-227.
 Misu, Y. and Ishihara, T. (2018) Prediction of Frequency Distribution of Strong Crosswind in a Control Section for Train Operations by Using Onsite Measurement and Numerical Simulation. Journal of Wind Engineering and Industrial Aerodynamics, 174, 69-79.
 Fukuhara, T., Tanimoto, S. and Araki, K. (2016) Methods to Estimate Spatial Distribution of Local Meteorological Conditions along Railway Line. Quarterly Report of RTRI, 57, 268-274.
 Araki, K., Fukuhara, T., Shimamura, T. and Imai, T. (2011) Method for Detecting Railway Line Sections Exposed to Strong Winds Using Numerical Simulations. Quarterly Report of RTRI, 52, 27-33.
 Burlando, M., Freda, A., Ratto, C.F. and Solari, G. (2010) A Pilot Study of the Wind Speed along the Rome-Naples HS/HC Railway Line. Part 1—Numerical Modelling and Wind Simulations. Journal of Wind Engineering and Industrial Aerodynamics, 98, 392-403.
 Freda, A. and Solari, G. (2010) A Pilot Study of the Wind Speed along the Rome-Naples HS/HC Railway Line: Part 2—Probabilistic Analyses and Methodology Assessment. Journal of Wind Engineering and Industrial Aerodynamics, 98, 404-416.
 Takemi, T., Kusunoki, K., Araki, K., Imai, T., Bessho, K., Hoshino, S. and Hayashi, S. (2010) Representation and Localization of Gusty Winds Induced by Misocyclones with a High-Resolution Meteorological Modeling. Theoretical and Applied Mechanics Japan, 58, 121-130.
 Kim, J. and Moin, P. (1985) Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations. Journal of Computational Physics, 59, 308-323.
 Kajishima, T. (1994) Upstream-Shifted Interpolation Method for Numerical Simulation of Incompressible Flows. Transactions of the Japan Society of Mechanical Engineers Series B, 60, 3319-3326. (In Japanese)
 Kawamura, T., Takami, H. and Kuwahara, K. (1986) Computation of High Reynolds Number Flow around a Circular Cylinder with Surface Roughness. Fluid Dynamics Research, 1, 145-162.
 Smagorinsky, J. (1963) General Circulation Experiments with the Primitive Equations. Part 1, Basic Experiments. Monthly Weather Review, 91, 99-164.
 Uchida, T. and Ohya, Y. (2007) Numerical Simulation of Local Strong Wind Induced by Topographic Effect. The Proceedings of the Fluids Engineering Conference, Hiroshima, 17 November 2007. (In Japanese)