| 内田 孝紀（うちだ たかのり） | データ更新日：2024.04.18 |
教授 ／
応用力学研究所
附属再生可能流体エネルギー研究センター
| 1. | Koichiro Shibuya, Takanori Uchida, Wake asymmetry of yaw state wind turbines induced by interference with wind towers, Energy, https://doi.org/10.1016/j.energy.2023.128091, 280, 2023.10, Wind turbine wakes are known to cause significant reductions in power generation and increased load on downwind wind turbines. Therefore, controlling wind turbine wakes is crucial, and yaw steering is considered one of the most effective methods of doing so. In this study, wind tunnel experiments using a wind turbine model and Large-eddy simulations (LES) were conducted to gain a better understanding of the wake of a yawed wind turbine. The results showed a clear difference in the horizontal distribution of wake velocity at hub height for positive and negative yaw angles, as well as both lateral and vertical wake deflections. Further investigation revealed that wind towers have a significant effect on the wakes of yawing wind turbines. For positive yaw angles, the velocity deficit is larger above the rotor, resulting in a vertical wake shape caused by the interference between the blade wakes and detached flow from the tower. Conversely, for negative yaw angles, the velocity deficit is larger below the rotor, resulting in a horizontal wake shape. These findings will facilitate the development of more accurate wake models and control methods..  |
| 2. | ＠内田 孝紀, 最新版リアムコンパクトソフトウェアの紹介と将来展望, 日本風力エネルギー学会誌, 47, 3, 465-470, 2023.10. |
| 3. | ＠Koichi Watanabe, Megumi Matsumoto, Thandar Nwe, Yuji Ohya, Takashi Karasudani, Takanori Uchida , Power Output Enhancement of Straight-Bladed Vertical-Axis Wind Turbines with Surrounding Structures, Energies, https://doi.org/10.3390/en16186719, 16, 18, 6719, 2023.09, Wind tunnel experiments were conducted by installing wind-acceleration structures on both sides of a straight-bladed vertical-axis wind turbine (VAWT) to improve the output performance of the turbine. In the case of Venturi-shape structures, a curved shape with a large outlet opening produced a higher power output than straight or brimmed Venturi shapes. More importantly, two simple flat plates installed upstream of the wind turbine achieved the highest power enhancement of 2.4 times the power of the bare wind turbine. From the analysis of the flow visualization results, the power enhancement was attributed to the increase in lift force on the blades in the upstream region due to the acceleration of the gap flow between the flat plates, and the decrease in drag force on the blades toward the upstream region due to stagnation of the flow behind the plates.. |
| 4. | ＠村上 礼雄,＠山崎 将英, ＠肥髙 邦彦, ＠飯野 哲平, ＠馬場 好孝, ＠内田 孝紀, ＠高桑 晋, ＠渡邊 慶一郎, ＠長谷川 聖矢, ＠藤川 凛太郎, 洋上風力発電所の利用率向上に向けた風車ウエイクの観測と予測モデルの開発, 日本エネルギー学会大会講演要旨集/ 第32回日本エネルギー学会大会, https://doi.org/10.20550/jietaikaiyoushi.32.0_68, 68-69, 2023.08, 抄録 In offshore wind firm, wind turbines are densely placed in a narrow area, and a wind speed attenuation and strong turbulence generated behind the turbines, known as wakes, cause decrease power with other wind turbines. In this study, to clarify the wind speed attenuation as a function of the distance from the wind turbine, we measured the wakes of the 3.3 MW wind turbine in commercial operation and the small model wind turbine installed in a wind tunnel. Furthermore, we applied this wake data to machine learning to develop new wake prediction models, which showed wind speed profiles within a relative error of 5%.. |
| 5. | ＠Tran The Hung, ＃Masato Hijikuro, ＠Masayuki Anyoji, ＠Takanori Uchida, ＠Takuji Nakashima, ＠Keigo Shimizu, Surface flow and aerodynamic drag of Ahmed body with deflectors, Experimental Thermal and Fluid Science, 145, 110887(18pages), 2023.07, This study investigated the surface flow and the aerodynamic characteristics of a 7/10 scale Ahmed model equipped with deflectors at several angles. The experiments were conducted at Reynolds number of 2.45 × 105 based on the model height. A global luminescent oil film skin-friction measurement technique was applied for the flow fields on the slant and base surfaces. The aerodynamic forces of the model were measured to understand the effect of the deflection angle on the aerodynamic characteristics. The pressure was measured and pressure drags acting on the slant and base surfaces were also calculated. The results showed that the drag and lift of the model suddenly dropped at the deflection angle above −5°. The pressure distribution on the slant and base surfaces became flat with a maximum drag reduction of 8 for the model with a deflection angle of 0°. The skin-friction measurement indicated that the drag reduction was connected to the breakdown of the separation bubble and longitudinal vortices on the slant. At the fully separated flow state, the skin-friction results illustrated two small vortices and a counter-clockwise vortex on the slant. Additionally, the wake vortices shifted from horizontal to vertical structures. The relation among the flow fields on the slant surface, flow on the symmetric plan, pressure distribution, and drag was illustrated for the model with different deflection angles.. |
| 6. | ＠渡邉 康一, ＃柳川 泰我, ＠内田 孝紀, 風力による排気塔ファン消費電力の省エネルギー化, 2023年度年次研究発表会梗概集・日本風工学会, https://doi.org/10.14887/jaweam.2023.0_151, 151-152, 2023.07. |
| 7. | ＠Takanori Uchida, Koichiro Shibuya, Gustavo Richmond-Navarro, Williams R. Calderón-Muñoz, A wind tunnel investigation of yawed wind turbine wake impacts on downwind wind turbine performances and wind loads, Wind Engineering, https://doi.org/10.1177/0309524X221150219, 47, 3, 0309524X221150219, 2023.06, In the current work we experimentally explored yawed wind turbine wake impacts on downwind wind turbine performances and wind loads. The lab-scale wind turbine model with a rotor diameter (D) of 0.442 m and a height of 1 m (=2.26D) was installed in a closed-circuit boundary layer wind tunnel (test section: 15 m long × 3.6 m wide × 2.0 m high) of the Research Institute for Applied Mechanics (RIAM) of Kyushu University. Power performance tests were initially conducted with a single rotor in isolation in order to characterize a rotor’s power output in stand-alone conditions. A detailed comparison of the tests revealed that the power output decreased rapidly as the yaw angle (γ) increased. It is presumed that the power output decrease in yawed cases is mainly due to the decrease in the effective rotor area and the change in the angle of the incoming wind flow with respect to the wind turbine blade. Next, using two wind turbine models aligned with the dominant inflow direction, the merging wakes behaviors caused by three different lateral separation distances were tested: (a) Case 1 (y = 0), (b) Case 2 (y = 0.5D), and (c) Case 3 (y = 1D). Here, the separation distance between the two wind turbine models was fixed at 6D in all cases. Extremely large power output deficits of 46%–76% were seen in the Case 1 configuration. This is mainly due to the significant wake velocity deficits induced by the upwind wind turbine model. In the Case 2 configuration with γ values of 20° and 30°, a significant increase in the power output of the downwind wind turbines was observed. Similar to Case 1 configuration, these results are considered to be mainly due to the upwind turbine-induced wake velocity deficits and wake deflection. Finally, in the Case 3 configuration, no significant difference was found in all of the results, and the tendency was almost the same. We show that the wake velocity deficits induced by the upwind wind turbine model had almost no effect on the power output of the downwind wind turbine model. We evaluated the total power output of the two turbines. As a result, in the Case 2 configuration with 20° yaw angle, the total power output of the two wind turbine models was the highest due to the increase in the power output of the downwind wind turbine model. In order to investigate the main cause of the significant increase in the power output of the downwind wind turbine model at 20° and 30° yaw angles in the Case 2 configuration, we measured the lateral wind speed distribution at the 6D position on the downwind side of the upwind wind turbine model by using the ultrasonic anemometer. As a results, it was clarified that the peak of the wake velocity deficits induced by the upwind wind turbine model is clearly shifted away from the downwind turbine such that it experiences a smaller deficit due to wake steering. Also, with wake steering the upwind turbine-induced wake velocity deficits may be smaller due to the reduction in rotor area. Finally, it is extremely important to understand the wind load acting on the downwind wind turbine model operating within the wake region induced by the upwind yawed wind turbine model when the maximum power output is generated. It can be seen that as the yaw angle of the upwind wind turbine model increased, the power output generated by the downwind wind turbine model and the streamwise wind load acting on it also increased. However, it was also clarified that the streamwise wind load acting on the downwind wind turbine model in this situation did not exceed the stand-alone value.. |
| 8. | Takanori UCHIDA, Teppei TANAKA, Ryuta SHIZUI, Hiroto ICHIKAWA, Ryo TAKAYAMA, Kazuomi YAHAGI, Ryoya OKUBO, GPU simulation of wake effects at the Horns Rev 1 offshore wind farm using the CFD porous disk wake model, Wind Engineering, https://doi.org/10.1177/0309524X221132003, 47, 2, 408-421, 0309524X221132003, 2023.04, To verify the effectiveness of the GPU simulation of wake effects at a large-scale offshore wind farm, we ran an in-house large-eddy simulation (LES) solver with a CFD porous disk wake model for the Horns Rev 1 wind farm. For this numerical research, we prepared the latest workstation equipped with a Xeon W-2265 CPU and an NVIDIA RTX A6000 GPU. We clarified that the calculation speed of the single GPU of the NVIDIA RTX A6000 is approximately 10 times faster than the calculation speed of the Xeon W-2265. Careful data analysis and visualization of the unsteady turbulent flow fields obtained in the current LES study suggest that the mutual interference of the wakes developed by wind turbines may frequently form a local speed-up region around wind turbines, located on the downstream side of large offshore wind farms.. |
| 9. | ＠静居 竜大, ＠内田 孝紀, ＠田中 鉄平, ＠市川 弘人, ＠高山 糧, ＠矢作 和臣, 日本域における大気安定度分布の地域特性, 風力エネルギー学会 論文集, https://doi.org/10.11333/jwearonbun.46.4_38, 46, 4, 38-49, 2023.03, To improve the forecast accuracy of the annual energy production from wind power, the characteristics of the atmospheric stability distribution are summarized using data from wind observation towers at ten sites in Japan. The bulk Richardson numbers indicate that stable conditions prevailed at seven out of ten locations. An analysis of each met mast indicated that the prevailing stability varied by season and differed not only according to the area but also according to the wind farm scale. These differences are thought to result from the distance to the sea and the wind shear. In addition, we compared the results according to the observational data and the bulk Richardson number calculated using the Weather Research and Forecasting model（WRF） during each season. The prevailing atmospheric stability at the observation points is affected by the atmospheric stability of the land and the surrounding ocean. The results of previous research and this study indicate that considering the atmospheric stability enables more correct estimations of the annual wind power energy production. If possible, seasonal and local differences in the atmospheric stability should be considered to better estimate the wind power generation.. |
| 10. | ＠Tran The Hung, ＃Masato Hijikuro, ＠Masayuki Anyoji, ＠Takanori Uchida, ＠Takuji Nakashima, ＠Keigo Shimizu, Effect of a short, bio-mimetic control device on aerodynamic drag of Ahmed body, Journal of Fluid Enginerring, https://doi.org/10.1115/1.4056341, 145, 3, 031206 (11 pages), 2023.03, The effects of short automatic moving deflectors (AMD) on the aerodynamic characteristics of Ahmed body are considered in this study. AMDs, known as biomimetic control devices, were added to the leading edge of the slanted surface. Its position was automatically adjusted in a separation flow. The aerodynamic drag, the pressure, and the skin-friction distribution on the slanted surface were measured for the model with three deflectors with lengths of 9%, 18%, and 30% of the slant. Particle image velocimetry was also utilized to assess the flow on the vertical symmetric plane. The Reynolds number based on the height of the model is between 1.44 × 105 and 2.80 × 105. The results showed that at a low Reynolds number, a short deflector increases the drag of the model. The effectiveness of the deflector in reducing the drag arises at a high velocity, where a maximum drag reduction of 11% was observed. The deflectors also reduced the lift coefficient by as much as 89%. Global luminescent oil-film skin-friction measurements showed that in the low drag state, the structure of the longitudinal vortexes and the separation bubble disappear on the surface. A complex flow structure is classified for the baseline model and the model with deflectors. The relationship between the surface flow, pressure distribution, and flow on the symmetric vertical plane is discussed in detail.. |
| 11. | Takanori UCHIDA, Susumu TAKAKUWA, Keiichiro WATANABE, Seiya HASEGAWA, Yoshitaka BABA Reo MURAKAMI, Masahide YAMASAKI, Kunihiko HIDAKA, Numerical visualization of wind turbine wakes using passive scalar advection-diffusion equation and its application for wake management, Wind Engineering, 10.1177/0309524X221113011, 46, 6, 1870-1887, 0309524X221113011, 2022.12. |
| 12. | ＠Tran The Hung, ＃Masato Hijikuro, ＠Masayuki Anyoji, ＠Takanori Uchida, ＠Takuji Nakashima, ＠Keigo Shimizu, Deflector effect on flow behavior and drag of an Ahmed body under crosswind conditions, Journal of Wind Engineering and Industrial Aerodynamics, https://doi.org/10.1016/j.jweia.2022.105238, 231, 105238, 2022.12, This study investigates the effect of a deflector on Ahmed model drag and the flow on the slant during steady crosswind conditions by experimental methods at a based-height Reynolds number of 2.45 × 105. The length deflector is 0.09 times the length of the slant. The deflector was fixed at the leading edge of the slant with an upward angle of 5° to the horizontal axis. Force, pressure, and global skin-friction measurements were conducted to understand the relation between drag, pressure distribution, and flow fields on the slant surface. The results showed that the deflector reduces aerodynamic drag at low yaw angles. This result is due to the breakdown structure of longitudinal vortexes and the separation bubble, which leads to fully separated flow on the leading edge of the slant. However, at yaw angles above 8°, a large reversed flow region forms on the slant and the drag of the model increases. A similar flow structure on the slant is observed for yaw angles of 12° and 15°. Conversely, the model lift is reduced with the deflector at all yaw angles tested. Detailed pressure distribution and complex flow structure on the slant provide insight in this study.. |
| 13. | ＠鈴木 英之, ＠内田 孝紀, ＠谷山 賀浩, 洋上風力発電の展望と風況解析技術, 電気学会, https://doi.org/10.1541/ieejjournal.142.519, 142, 8, 519-522, 2022.08. |
| 14. | Susumu TAKAKUWA, Takanori UCHIDA, Improvement of Airflow Simulation by Refining the Inflow Wind Direction and Applying Atmospheric Stability for Onshore and Offshore Wind Farms Affected by Topography, Energies, 10.3390/en15145050, 15, 14, 5050, 2022.07. |
| 15. | 内田 孝紀, 高桑 晋, 渡邊 慶一郎, 長谷川 聖矢, 馬場 好孝, 村上 礼雄, 山崎 将英, 肥高 邦彦, UAVを用いた大型商用風車のウェイク計測の試み, 一般社団法人 日本風工学会 2022年度年次研究発表会梗概集, https://doi.org/10.14887/jaweam.2022.0_87, 2022.07, 本研究では，特別な観測機器をドローンに搭載することなく，ドローン自身を風速センサとして活用する手法を試みる．最初に，大型風洞設備を用い，ドローンを風速センサとして活用するための校正試験を行った．次に，風洞スケールおよび実機スケールでの風車ウェイクの気流計測を試みた．. |
| 16. | Gustavo Richmond-Navarro, Takanori Uchida, Williams R. Calderón-Muñoz, Shrouded wind turbine performance in yawed turbulent flow conditions, Wind Engineering, 10.1177/0309524X211036041, 46, 2, 518-528, Aritical Number 0309524X211036041, 2022.04, Wind turbines represent a growing energy source worldwide. In many cases, operating in turbulent and changing wind direction spots. In this work, we use a wind tunnel to analyze the effect of the turbulence in a wind turbine provided with a Wind Lens flow concentrator, under yaw conditions, for turbulence intensity values of 10% and 15%. Measurements are made of the power coefficient as a function of the Tip Speed Ratio using two types of Wind Lens, CiiB5 and CiiB10, at yaw angles from 0° to 30°. In general, for the turbine with Wind Lens, an increase of the yaw angle causes a reduction of the power coefficient. If the turbine operates with the CiiB10, the stronger the turbulence, the greater performance is obtained. In conclusion, for the case of turbulent flow and yaw = 20° or less, the Wind Lens turbine offers better performance than without the flow concentrator.. |
| 17. | ＠内田 孝紀、＠小野 謙二、＠飯田 明由、＠吉村 忍、＠加藤 千幸、＠山出 吉伸、＠今村 博、＠植田 祐子, スパコン版リアムコンパクトによる風車ウエイクの相互干渉に関する大規模数値シミュレーション―ウエイク領域内の平均風速欠損量の予測精度検証―, 風力エネルギー学会 論文集, https://doi.org/10.11333/jwearonbun.45.4_71, 45, 4, 71-82, 2022.03, In the current project, the first author and the second author play a central role in conducting high-speed tuning and wind turbine wake analysis of the supercomputer version of RIAM-COMPACT. In this report, we applied the supercomputer version RIAM-COMPACT to the wind turbine wake simulation from the wind tunnel scale to the utility-scale wind turbine wakes. As a result, we clarified unsteady wake aerodynamics of wind turbines including multiple wake interactions with high precision and high reality.. |
| 18. | ＠Takanori UCHIDA, ＠Yves GAGNON, Effects of continuously changing inlet wind direction on near-to-far wake characteristics behind wind turbines over flat terrain, Journal of Wind Engineering and Industrial Aerodynamics, 10.1016/j.jweia.2021.104869, 220, 104869, 1-17, 2022.01, The wake characteristics of a utility-scale wind turbine under realistic atmospheric boundary layer conditions are affected by the continuously changing wind direction arriving at the wind turbine. In the present study, the effects of continuous changes in the incoming wind direction were studied for a wind turbine on flat terrain, with the objective of understanding the wake characteristics of the wind turbine. Thus, understanding the effects of continuously changing incoming wind direction on the wake characteristics of wind turbines over flat terrain is important in the design of wind farm layouts, including in the design of offshore wind power plants. For this purpose, a computational fluid dynamics (CFD) approach using large-eddy simulations (LES) was adopted in the present study. An in-house LES-solver based on the actuator line (AL) aerodynamics technique was constructed in order to successfully capture the wake structure behind the wind turbine. First, experimental investigations on both a blade-only wind turbine scale model and a full 3D wind turbine scale model (isolated wind turbine) were conducted for a fixed inlet wind condition, the latter including the nacelle and the tower. Through a detailed comparison of the wind tunnel experimental and numerical results, the prediction accuracy of the in-house LES-solver was verified and validated for fixed inlet wind conditions. On the basis of the validation results obtained, and using the full 3D wind turbine scale model, the effects of the continuously changing inlet wind conditions on the wake characteristics in the near- and far-wake regions were numerically investigated. In addition, the effects of the wind turbine nacelle and tower on the wake characteristics were also investigated. The numerical results show that the most significant impact of the effects of the continuously changing wind direction was the rapid recovery of the mean velocity deficits in the wind turbine wake region. Further, at the x = 10D position (D is the rotor diameter) downstream of the wind turbine, the non-dimensional streamwise mean velocity was 0.93, which nearly matches the approaching flow speed, under an optimal tip speed ratio of 4.0, compared to the fixed wind direction scenario.. |
| 19. | ＠内田 孝紀、＠高桑 晋, 複雑地形における風力資源の数値予測に対する大気安定度の影響, 風力エネルギー学会 論文集, https://doi.org/10.11333/jwearonbun.45.1_1, 45 （2021）, 1, 1-9, 2021.12. |
| 20. | ＠Tadasuke YOSHIDA, ＠Takanori UCHIDA, ＠Masaki INUI, ＠Yoshihiro TANIYAMA, ＠Yuki FUKATANI, Wind turbine wake evaluation using a vertical profiling Lidar and new CFD porous disk model, 第15回動力エネルギー際会議(ICOPE-2021) The Japan Society of Mechanical Engineers, https://doi.org/10.1299/jsmeicope.2021.15.2021-0175, セッションID: 2021-0175, 2021.11, Wind turbine wakes reduce the power generation and life of the downstream wind turbines. Accurately predicting the impact of wind turbine wakes is very important for evaluating the feasibility of large wind farms. We have recently proposed a computational fluid dynamics (CFD) porous disk (PD) wake model as an intermediate method between engineering wake models and CFD wake models in order to predict accurately the time-averaged wind speed deficits in the wind turbine wakes. In this study, to further evaluate the validity of the CFD PD wake model, we additionally measured the wind turbine wakes of 2 MW‐class downwind turbines installed on the coastal area using a vertical profiling lidar (ZephIR ZX300), and considered them including the data in the previous report to clarify its airflow characteristics in detail. Based on the measurement results by the lidar, we reproduced the wind turbine wake using the CFD PD wake model. The simulated vertical distribution of wind speed by the CFD PD wake model corresponded with the measurement result of the lidar.. |
| 21. | Yuki FUKATANI, Takanori UCHIDA, Yoshihiro TANIYAMA, Zhiren BAI, Tadasuke YOSHIDA, Masaki INU, Wind turbine wake evaluation considering atmospheric stability using vertical profiling LiDARs, 第15回動力エネルギー際会議(ICOPE-2021) The Japan Society of Mechanical Engineers, https://doi.org/10.1299/jsmeicope.2021.15.2021-0175, セッションID: 2021-0175, 2021.11, In wind farms, an upstream turbine causes the turbine wake which leads to a decrease in the power output and an unsteady loading to the downstream turbine due to the velocity deficit and the increased level of turbulence. Therefore, accurate prediction of the wake flow is important for evaluating a feasibility of large scale wind farms. In this research, we report results of analysis of seaward and landward wake data measured using three vertical profiling LiDARs. In addition, the effect of atmospheric stability on the wake is investigated using atmospheric reanalysis data. The results show that considering atmospheric stability for seaward wake is important especially when a wide range of atmospheric stability can be observed.. |
| 22. | ＠浜口 慶一郎, ＠髙田 稜一, ＠内田 孝紀, ＠杉谷 賢一郎, ＠谷川 博哉, ＠野口 尚史, ＠平田 勝哉, 極低レイノルズ数における基本的翼型の空力特性と後流特性, 一般社団法人 日本機械学会 2021年度 年次大会, https://doi.org/10.1299/jsmemecj.2021.S051-01, 2021.09. |
| 23. | ＃吉岡 達也, ＠澁谷 光一郎, ＠内田 孝紀, ヨーステアリング手法を適用した大規模洋上ウィンドファームの最適化に関する研究, 一般社団法人 日本機械学会 2021年度 年次大会, https://doi.org/10.1299/jsmemecj.2021.S055-05, 2021.09. |
| 24. | ＃吉岡 達也, ＠村本 幸央, ＠内田 孝紀, UAVの飛行データを活用した新たな風況観測手法の開発および風車後流計測への適用, 一般社団法人 日本機械学会 2021年度 年次大会, https://doi.org/10.1299/jsmemecj.2021.S055-05, 2021.09. |
| 25. | Takanori Uchida, Ryo Araya, Applications of the Atmospheric Transport and Diffusion of LES Modeling to the Spread and Dissipation of COVID-19 Aerosol Particles inside and outside the Japan National Stadium (Tokyo Olympic Stadium), Modelling and Simulation in Engineering, 10.1155/2021/8822548, 2021, Article ID 8822548, 1-19, 2021.06, In this paper, we use an analysis function for gas diffusion known as the Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain (RIAM-COMPACT), which was developed for complex terrain, in Airflow Analyst software, and apply it to the spread and dissipation of a fluid layer (assuming the fluid layer contains COVID-19 particles). First, to verify the prediction accuracy of the gas diffusion using RIAM-COMPACT, comparisons with past wind tunnel test results conducted on simple and complex terrains are presented under neutral atmospheric stability. The results of the numerical simulations carried out in this study show good agreement with the wind tunnel experiments for both simple and complex terrains. Next, a model of the Japan National Stadium (Tokyo Olympic Stadium) was constructed using 3D detailed topographic Advanced World 3D Map (AW3D) data generated by combining high-resolution satellite images. We tried to reproduce the hypothetical spread and dissipation of the fluid layer (assuming the fluid layer contains COVID-19 particles) inside and outside of the Japan National Stadium using Airflow Analyst implemented with the RIAM-COMPACT analysis function for gas diffusion. We paid special attention to the effect of wind ventilation driven by natural wind. The numerical results under various scenarios show that ventilation driven by natural wind is very effective for the Japan National Stadium.. |
| 26. | ＠内田 孝紀、＠小野 謙二、＠飯田 明由、＠吉村 忍、＠加藤 千幸、＠山出 吉伸、＠今村 博、＠植田 祐子, スパコン版リアムコンパクトによる風車ウエイクの相互干渉シミュレーション, 日本風工学会 2021年度年次研究発表会梗概集, https://doi.org/10.14887/jaweam.2021.0_101, 101-102, 2021.06. |
| 27. | 内田 孝紀, 高桑 晋, 複雑地形における風力資源の数値予測に対する大気安定度の影響, 日本風力エネルギー学会論文集, 45, 1(通巻137), 1-9, 2021.05. |
| 28. | ＠Takanori Uchida, Tadasuke Yoshida, Masaki Inui, Yoshihiro Taniyama, Doppler Lidar Investigations of Wind Turbine Near-Wakes and LES Modeling with New Porous Disc Approach, Energies, 10.3390/en14082101, 14, 8, 1-33, Article Number 21012101, 2021.04. |
| 29. | 内田孝紀, 最新の数値風況シミュレーション技術リアムコンパクトが実現するバーチャルウィンドファームー開発の歴史と将来展望ー, 日本風力エネルギー学会誌, 44, 4, 666-671, 2021.02. |
| 30. | 内田孝紀, 地形性乱流が風車ブレードに与える危険度を判定するための新しい評価スケールの提案, 風力エネルギー学会誌, 44, 3, 381-384, 2020.11, 本研究では，鹿児島県いちき串木野市羽島地区に位置する串木野れいめい風力発電所(平成24年11月より運転開始)の風車10号機を対象とし，風況データおよび風車ブレード歪データの実計測を行い，地形性乱流が風車ブレードに与える影響について詳細に考察した．. |
| 31. | 内田孝紀, 大規模洋上ウィンドファームのバンカビリティ（融資適格性）評価に資する汎用的使用を目的とした新しい風車ウエイクモデルの開発とその応用, 月刊『 クリーンエネルギー 』, 11月号, 42-49, 2020.11. |
| 32. | Amr M. Halawa, Takanori Uchida, Koichi Watanabe, Yuji Ohya, Validation Study of Multi-Rotor Systems Using Two Shrouded Wind Turbines, Journal of Physics: Conference Series, 10.1088/1742-6596/1618/3/032017, 1618, 1-15, The Science of Making Torque from Wind (TORQUE 2020), 2020.09, Multi rotor systems (MRS) have shown a great potential as a future application of wind energy. In this study, the aim is simulating the aerodynamic performance of a an MRS using fully-resolved shrouded wind turbine blades then validating with experimental data. MRS for wind turbine configurations have been studied using both numerical and experimental approaches. Different case studies have been studied and the power output comparison have been reported. The wind lens turbine (WLT) obviously shows a large increase in power output compared to the bare rotor case. Besides, the twin side-by-side (SBS) WLTs shows even larger power increase compared to the case of single WLT by around 20% for computational fluid dynamics (CFD) calculations at the optimum tip speed ratio s/D of around 0.2. The increase in power coefficient in close proximity can be explained by flow interference and gap flow behaviors. Previously, we used simplified models for blade modeling like actuator line method (ALM) and actuator disk method (ADM). However, currently we are improving the accuracy using CFD with full-scale blades with higher grid resolutions. As the number of units for an MRS is increased, the increase in power output becomes larger and larger. This is because that the gap flows between brimmed diffuser-augmented wind turbines (B-DAWT) in a MRS are accelerated and cause lowered pressure regions due to vortex interaction behind the brimmed diffusers. Thus, an MRS with more B-DAWTs can draw more wind into turbines showing higher power output.. |
| 33. | Takanori Uchida, Effects of Inflow Shear on Wake Characteristics of Wind-Turbines over Flat Terrain, Energies, 10.3390/en13143745, 13, 14, 3745, 2020.07, [URL], The scope of the present study was to understand the wake characteristics of wind-turbines under various inflow shears. First, in order to verify the prediction accuracy of the in-house large-eddy simulation (LES) solver, called RIAM-COMPACT, based on a Cartesian staggered grid, we conducted a wind-tunnel experiment using a wind-turbine scale model and compared the numerical and experimental results. The total number of grid points in the computational domain was about 235 million. Parallel computation based on a hybrid LES/actuator line (AL) model approach was performed with a new SX-Aurora TSUBASA vector supercomputer. The comparison between wind-tunnel experiment and high-resolution LES results showed that the AL model implemented in the in-house LES solver in this study could accurately reproduce both performances of the wind-turbine scale model and flow characteristics in the wake region. Next, with the LES solver developed in-house, flow past the entire wind-turbine, including the nacelle and the tower, was simulated for a tip-speed ratio (TSR) of 4, the optimal TSR. Three types of inflow shear, N = 4, N = 10, and uniform flow, were set at the inflow boundary. In these calculations, the calculation domain in the streamwise direction was very long, 30.0 D (D being the wind-turbine rotor diameter) from the center of the wind-turbine hub. Long-term integration of t = 0 to 400 R/Uin was performed. Various turbulence statistics were calculated at t = 200 to 400 R/Uin. Here, R is the wind-turbine rotor radius, and Uin is the wind speed at the hub-center height. On the basis of the obtained results, we numerically investigated the effects of inflow shear on the wake characteristics of wind-turbines over a flat terrain. Focusing on the center of the wind-turbine hub, all results showed almost the same behavior regardless of the difference in the three types of inflow shear.. |
| 34. | Takanori UCHIDA, Susumu TAKAKUWA, Numerical Investigation of Stable Stratification Effects on Wind Resource Assessment in Complex Terrain, Energies, 10.3390/en13153896, 6638, 2020.07. |
| 35. | Takanori Uchida, Yoshihiro Taniyama, Yuki Fukatani, Michiko Nakano, Zhiren Bai, Tadasuke Yoshida, Masaki Inui, A new wind turbine CFD modeling method based on a porous disk approach for practical wind farm design, Energies, 10.3390/en13123197, 13, 12, 3197, 2020.06, [URL], In this study, the new computational fluid dynamics (CFD) porous disk (PD) wake model was proposed in order to accurately predict the time-averaged wind speed deficits in the wind turbine wake region formed on the downstream side by the 2-MW wind turbine operating at a wind speed of 10 m/s. We use the concept of forest canopy model as a new CFD PD wake model, which has many research results in the meteorological field. In the forest canopy model, an aerodynamic resistance is added as an external force term to all governing equations (Navier-Stokes equations) in the streamwise, spanwise, and vertical directions. Therefore, like the forest model, the aerodynamic resistance is added to the governing equations in the three directions as an external force term in the CFD PD wake model. In addition, we have positioned the newly proposed the LES using the CFD PD wake model approach as an intermediate method between the engineering wake model (empirical/analytical wake model) and the LES combined with actuator disk (AD) or actuator line (AL) models. The newly proposed model is intended for use in large-scale offshore wind farms (WFs) consisting of multiple wind turbines. In order to verify the validity of the new method, the optimal model parameter CRC was estimated by comparison with the time-averaged wind speed database in the wind turbine wake region with fully resolved geometries, combined with unsteady Reynolds-averaged Navier-Stokes (RANS) equations, implemented using the ANSYS(R) CFX(R) software. Here, product names (mentioned herein) may be trademarks of their respective companies. As a result, in the range from x = 5D of the near wake region to x = 10D of the far wake region, by selecting model parameter CRC, it was clarified that it is possible to accurately evaluate the time-averaged wind speed deficits at those separation distances. We also examined the effect of the spatial grid resolution using the CFD PD wake model that is proposed in the present study, clarifying that the spatial grid resolution has little effect on the simulation results shown here.. |
| 36. | Koichi WATANABE, Sho FUKUTOMI, Yuji OHYA, Takanori UCHIDA, An Ignored Wind Generates More Electricity A Solar Updraft Tower to a Wind Solar Tower, International Journal of Photoenergy, 10.1155/2020/4065359, 2020, 2020.01, [URL], A solar updraft tower is one of the wind power generation plants which utilizes solar energy. The purpose of this study was to ascertain whether the tower was also able to utilize crosswind energy. Wind tunnel experiments and numerical simulations were conducted simulating the crosswind. The results showed that suctioned updraft speed in the tower was proportional to the crosswind speed, and its conversion rate depended on the tower configuration. A diffuser-shaped tower with a vortex generator achieved to produce the updraft whose speed exceeded the crosswind speed. It was due to the low pressure created by the vortex atop the tower and to the diffuser effect. The crosswind utilization enables the simple power generation device to generate electricity during the night, and the hybrid utilization of renewable energies contributes to the increasing wind energy market.. |
| 37. | Takanori Uchida, Kenichiro Sugitani, Numerical and Experimental Study of Topographic Speed-Up Effects in Complex Terrain, Energies, 13, 15, 3896, 2020.01, Our research group is developing computational fluid dynamics (CFD)-based software for wind resource and energy production assessments in complex terrain called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University (RIAM)-Computational Prediction of Airflow over Complex Terrain), based on large eddy simulation (LES). In order to verify the prediction accuracy of RIAM-COMPACT, we conduct a wind tunnel experiment that uses a two-dimensional steep ridge model with a smooth surface. In the wind tunnel experiments, airflow measurements are performed using an I-type hot-wire probe and a split film probe that can detect forward and reverse flows. The results of the numerical simulation by LES are in better agreement with the wind tunnel experiment using the split film probe than the results of the wind tunnel experiment using the I-type hot wire probe. Furthermore, we calculate that the two-dimensional ridge model by changing the length in the spanwise direction, and discussed the instantaneous flow field and the time-averaged flow field for the three-dimensional structure of the flow behind the model. It was shown that the eddies in the downwind flow-separated region formed behind the two-dimensional ridge model were almost the same size in all cases, regardless of the difference in the length in the spanwise direction. In this study, we also perform a calculation with a varying inflow shear at the inflow boundary. It was clear that the size in the vortex region behind the model was almost the same in all the calculation results, regardless of the difference in the inflow shear. Next, we conduct wind tunnel …. |
| 38. | Takanori UCHIDA, Keiji ARAKI, Reproduction of Local Strong Wind Area Induced in the Downstream of Small-scale Terrain by Computational Fluid Dynamic (CFD) Approach, Modelling and Simulation in Engineering, 10.1155/2019/7378315, Volume 2019, Article ID 7378315, 12 pages, 2019.12, In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation In this research, the computational fluid dynamic (CFD) approach that has been used in wind power generation. |
| 39. | Takanori UCHIDA, Numerical Investigation of Terrain-Induced Turbulence in Complex Terrain Using High-Resolution Elevation Data and Surface Roughness Data Constructed with a Drone, Energies, 10.3390/en12193766, 12(19), 3766, 2019.10, Using the method based on unmanned aerial vehicle (UAV) imagery, two kinds of data can be obtained: the digital elevation model (DEM) for the digital expression of terrain, and the digital surface model (DSM) for the digital expression of the surface of the ground, including trees. In this research, a 3D topography model with a horizontal spatial resolution of 1 m was reproduced using DEM. In addition, using the differences between the DEM and DSM data, we were able to obtain further detailed information, such as the heights of trees covering the surface of the ground and their spatial distribution. Therefore, the surface roughness model and the UAV imagery data were directly linked. Based on the above data as input data, a high-resolution 3D numerical flow simulation was conducted. By using the numerical results obtained, we discussed the effect of the existence of surface roughness on the wind speed at the height of the hub of the wind turbine. We also discussed the effect of the differences in the spatial resolution in the horizontal direction of the computational grid on the reproductive precision of terrain-induced turbulence. As a result, the existence and the vortex structure of terrain-induced turbulence occurring near the target wind turbine was clearly revealed. It was shown that a horizontal grid resolution of about 5 m was required to reproduce terrain-induced turbulence formed from topography with an altitude of about 127 m. By the simulation using the surface roughness model, turbulence intensity higher than class A in the International Electrotechnical Commission (IEC) turbulence category was confirmed at the present study site, as well as the measured data.. |
| 40. | Takanori UCHIDA, Ryo ARAYA, Practical Applications of the Large-Eddy Simulation Technique for Wind Environment Assessment around New National Stadium, Japan (Tokyo Olympic Stadium), Open Journal of Fluid Dynamics, 10.4236/ojfd.2019.94018, 09, 04, Article ID:96061, 2019.10. |
| 41. | 内田 孝紀, 川島 泰史, 地形性乱流が風車ブレードに与える危険度を判定するための新しい定量化指標の提案, 九州大学応用力学研究所所報, 156, 1-22, 2019.09, 本研究では，鹿児島県いちき串木野市羽島地区に位置する串木野れいめい風力発電所(平成24年11月より運転開始)の風車10号機を対象とし，風況データおよび風車ブレード歪データの実計測を行い，地形性乱流が風車ブレードに与える影響について詳細に考察した．. |
| 42. | Koichi WATANABE, Yuji OHYA, Takanori UCHIDA, Power Output Enhancement of a Ducted Wind Turbine by Stabilizing Vortices around the Duct, Energies, 10.3390/en12163171, 12(16), 3171, 2019.08, [URL], brimmed-diffuser augmented turbine (called a wind lens turbine: WLT) actively uses vortices around the brim to enhance its power output. However, the vortices are usually unstable and asymmetric. This study attempted to stabilize the vortices to enhance the power output of a WLT. Then, we investigated new approaches using vortex stabilization plates and polygonal brims in wind tunnel experiments and numerical simulations. Both approaches achieved a 1.5–3.8% increase in power output compared with a standard WLT. Our numerical simulations revealed a periodicity existing in a fluctuating vortex structure on the circular brim. Importantly, vortex stabilization plates and polygonal brims must be the same periodic scale to suppress the vortex fluctuation and stabilize the vortices effectively. In addition, a larger brim tended to enhance the stabilizing effects. We believe that this discovery rovides an easy way to increase the power output of existing wind turbines. It is particularly important in light of advances in wind energy technology and the increasing wind energy market.. |
| 43. | Takanori UCHIDA, Numerical Prediction of Strong Wind Induced by Topographic Effect, Open Journal of Fluid Dynamics, 10.4236/ojfd.2019.93015, 09, 03, ISSN Online: 2165-3860, ISSN Print: 2165-3852, 2019.08. |
| 44. | Takanori Uchida, Yasushi Kawashima, New assessment scales for evaluating the degree of risk of wind turbine blade damage caused by terrain-induced turbulence, Energies, 10.3390/en12132624, 12, 13, 2019.07, [URL], The present study scrutinized the impacts of terrain-induced turbulence on wind turbine blades, examining measurement data regarding wind conditions and the strains of wind turbine blades. Furthermore, we performed a high-resolution large-eddy simulation (LES) and identified the three-dimensional airflow structures of terrain-induced turbulence. Based on the LES results, we defined the Uchida-Kawashima Scale_1 (the U-K scale_1), which is a turbulence evaluation index, and clarified the existence of the terrain-induced turbulence quantitatively. The threshold value of the U-K scale_1 was determined as 0.2, and this index was confirmed to not be dependent on the inflow profile, the influence of the horizontal grid resolution, and the influence of the computed azimuth. In addition, we defined the Uchida-Kawashima Scale_2 (the U-K scale_2), which is a fatigue damage evaluation index based on the measurement data and the design value obtained by DNV GL’s Bladed. DNV GL (Det Norske Veritas Germanischer Lloyed) is a third party certification body in Norway, and Bladed has been the industry standard aero-elastic wind turbine modeling software. Using the U-K scale_2, the following results were revealed: the U-K scale_2 was 0.86 1.0 (exceeding the designed value) in the case of easterly wind. As a result, it was revealed that the blades of the target wind turbine were directly and strongly affected by terrain-induced turbulence when easterly winds occurred.. |
| 45. | Takanori Uchida, Susumu Takakuwa, A Large-Eddy Simulation-Based Assessment of the Risk of Wind Turbine Failures Due to Terrain -Induced Turbulence over a Wind Farim inComplex Terrain, Energies, 10.3390/en12101925, 19pages, 2019.05, The first part of the present study investigated the relationship among the number of yaw gear and motor failures and turbulence intensity (TI) at all the wind turbines under investigation with the use of in situ data. The investigation revealed that wind turbine #7 (T7), which experienced a large number of failures, was affected by terrain-induced turbulence with TI that exceeded the TI presumed for the wind turbine design class to which T7 belongs. Subsequently, a computational fluid dynamics (CFD) simulation was performed to examine if the abovementioned observed wind flow characteristics could be successfully simulated. The CFD software package that was used in the present study was RIAM-COMPACT, which was developed by the first author of the present paper. RIAM-COMPACT is a nonlinear, unsteady wind prediction model that uses large-eddy simulation (LES) for the turbulence model. RIAM-COMPACT is capable of simulating flow collision, separation, and reattachment and also various unsteady turbulence–eddy phenomena that are caused by flow collision, separation, and reattachment. A close examination of computer animations of the streamwise (x) wind velocity revealed the following findings: As we predicted, wind flow that was separated from a micro-topographical feature (micro-scale terrain undulations) upstream of T7 generated large vortices. These vortices were shed downstream in a nearly periodic manner, which in turn generated terrain-induced turbulence, affecting T7 directly. Finally, the temporal change of the streamwise (x) wind velocity (a non-dimensional quantity) at the hub-height of T7 in the period from 600 to 800 in non-dimensional time was re-scaled in such a way that the average value of the streamwise (x) wind velocity for this period was 8.0 m/s, and the results of the analysis of the re-scaled data were discussed. With the re-scaled full-scale streamwise wind velocity (m/s) data (total number of data points: approximately 50,000; time interval: 0.3 s), the time-averaged streamwise (x) wind velocity and TI were evaluated using a common statistical processing procedure adopted for in situ data. Specifically, 10-min moving averaging (number of sample data points: 1932) was performed on the re-scaled data. Comparisons of the evaluated TI values to the TI values from the normal turbulence model in IEC61400-1 Ed.3 (2005) revealed the following: Although the evaluated TI values were not as large as those observed in situ, some of the evaluated TI values exceeded the values for turbulence class A, suggesting that the influence of terrain-induced turbulence on the wind turbine was well simulated.. |
| 46. | Kenji Ono, Takanori Uchida, High-Performance Parallel Simulation of Airflow for Complex Terrain Surface, Modelling and Simulation in Engineering, 10.1155/2019/5231839, 2019, Article ID 5231839, 10 pages, 2019.02, [URL], It is important to develop a reliable and high-throughput simulation method for predicting airflows in the installation planning phase of windmill power plants. This study proposes a two-stage mesh generation approach to reduce the meshing cost and introduces a hybrid parallelization scheme for atmospheric fluid simulations. The meshing approach splits mesh generation into two stages: in the first stage, the meshing parameters that uniquely determine the mesh distribution are extracted, and in the second stage, a mesh system is generated in parallel via an in situ approach using the parameters obtained in the initialization phase of the simulation. The proposed two-stage approach is flexible since an arbitrary number of processes can be selected at run time. An efficient OpenMP-MPI hybrid parallelization scheme using a middleware that provides a framework of parallel codes based on the domain decomposition method is also developed. The preliminary results of the meshing and computing performance show excellent scalability in the strong scaling test.. |
| 47. | Takanori Uchida, Large-Eddy Simulation of Airflow over a Steep,Three-Dimensional Isolated Hill with Multi-GPUs Computing, Open Journal of Fluid Dynamics, 10.4236/ojfd.2018.84027, 8, 4, 416-434, 2018.11, The present research attempted a Large-Eddy Simulation (LES) of airflow over a steep, three-dimensional isolated hill by using the latest multi-cores multi-CPUs systems. As a result, it was found that 1) turbulence simulations using approximately 50 million grid points are feasible and 2) the use of this system resulted in the achievement of a high computation speed, which exceeded the speed of parallel computation attained by a single CPU on one of the latest supercomputers. Furthermore, LES was conducted by using the multi-GPUs systems. The results of these simulations revealed the following findings: 1) the multi-GPUs environment which used the NVDIA Tesla M2090 or the M2075 could simulate turbulence in a model with as many as approximately 50 million grid points. 2) The computation speed achieved by the multi-GPUs environments exceeded that by parallel computation which used four to six CPUs of one of the latest supercomputers.. |
| 48. | Takanori Uchida, Numerical Investigation of Terrain-induced Turbulence in Complex Terrain by Large-eddy Simulation (LES) Technique, Energies, 10.3390/en11061530, 11(10), 2638, 15pages, 2018.10, In the present study, field observation wind data from the time of the wind turbine blade damage accident on Shiratakiyama Wind Farm were analyzed in detail. In parallel, high-resolution large-eddy simulation (LES) turbulence simulations were performed in order to examine the model’s ability to numerically reproduce terrain-induced turbulence (turbulence intensity) under strong wind conditions (8.0–9.0 m/s at wind turbine hub height). Since the wind velocity and time acquired from the numerical simulation are dimensionless, they are converted to full scale. As a consequence, both the standard deviation of the horizontal wind speed (m/s) and turbulence intensity evaluated from the field observation and simulated wind data are successfully in close agreement. To investigate the cause of the wind turbine blade damage accident on Shiratakiyama Wind Farm, a power spectral analysis was performed on the fluctuating components of the observed time series data of wind speed (1 s average values) for a 10 min period (total of 600 data) by using a fast Fourier transform (FFT). It was suggested that the terrain-induced turbulence which caused the wind turbine blade damage accident on Shiratakiyama Wind Farm was attributable to rapid wind speed and direction fluctuations which were caused by vortex shedding from Tenjogadake (elevation: 691.1 m) located upstream of the wind farm.. |
| 49. | Takanori UCHIDA, Graham Li, Comparison of RANS and LES in the Prediction of Airflow Field over Steep Complex Terrain, Open Journal of Fluid Dynamics, 10.4236/ojfd.2018.83018, 08, 286-307, Article ID:87086,22 pages, 2018.09, The present study compared the prediction accuracy of the three CFD software packages for simulating airflow around a three-dimensional, isolated hill with a steep slope: 1) WindSim (turbulence model: RNG k-ε RANS), 2) Meteodyn WT (turbulence model: k-L RANS), which are the leading commercially available CFD software packages in the wind power industry and 3) RIAM-COMPACT (turbulence model: standard Smagorinsky LES), which has been developed by the lead author of the present paper. Distinct differences in the airflow patterns were identified in the vicinity of the isolated hill (especially downstream of the hill) between the RANS results and the LES results. No reverse flow region (vortex region) characterized by negative wind velocities was identified downstream of the isolated hill in the result from the simulation with WindSim (RNG k-ε RANS) and Meteodyn WT (k-L RANS). In the case of the simulation with RIAM-COMPACT natural terrain version (standard Smagorinsky LES), a reverse flow region (vortex region) characterized by negative wind velocities clearly forms. Next, an example of wind risk (terrain-induced turbulence) diagnostics was presented for a large-scale wind farm in China. The vertical profiles of the streamwise (x) wind velocity do not follow the so-called power law wind profile; a large velocity deficit can be seen between the hub center and the lower end of the swept area in the case of the LES calculation (RIAM-COMPACT).. |
| 50. | 髙橋 江梨香, 大西 一希, 谷川 博哉, 内田 孝紀, 杉谷 賢一郎, 平田 勝哉, 極低レイノルズ数領域での様々な二次元翼型の空力特性と流れ, 一般社団法人 日本機械学会 2018年度 年次大会, 10.1299/jsmemecj.2018.S0510205, 2018.09, We have been requiring more precise knowledge about the aerodynamic characteristics of the airfoil especially in low Reynolds-number ranges less than 106, because of the recent miniaturisation of robots such as unmanned aerial vehicles known as UAVs or drones and micro air vehicles known as MAVs, in addition to the importance of insect/bird flight dynamics, small-scale machines like micro fluid machineries and micro combustion engines and so on. In the present study, we investigate the relationship between various aerodynamic characteristics and attack angle α for flat plate NACA0015 and iNACA0015 (the NACA0015 placed back to front) by numerical analysis and water-tank experiment at Re = 1.0×102 － 8.0×102. Then, we investigate the effects of attack angle α upon various aerodynamic characteristics such as the lift coefficient CL, drag coefficient CD and the lift-to-drag ratio CL/CD . In order to discuss the revealed α effects, we visualise the flow around the airfoils using the Q value, the helicity, streamlines’ distribution around an airfoils at α = 0 – 24 deg., together with PIV analysis based on water-tank experiment.. |
| 51. | Takanori Uchida, Computational Investigation of the Causes of Wind Turbine Blade Damage at Japan's Wind Farm in Complex Terrain, Journal of Flow Control, Measurement & Visualization, 10.4236/jfcmv.2018.63013, 6, 3, 152-167, 2018.07, During the passage of Typhoon 0918 (Melor) over southern Honshu in Japan on 7 and 8 October 2009, strong winds with extremely high turbulence fluctuations were observed over Shirataki Mountain and the surrounding mountains in Shimonoseki, Yamaguchi Prefecture, Japan. These strong winds caused damage to wind turbine blades at the Shiratakiyama Wind Farm owned by Kinden Corporation. In order to investigate the causes of the blade damage, the airflow characteristics from the time of the incidences are first simulated in detail with the combined use of the WRF-ARW mesoscale meteorological model and the RIAM-COMPACT LES turbulence model (CFD model). Subsequently, in order to evaluate the wind pressure acting on the wind turbine blades, an airflow analysis is separately performed for the vicinity of the blades with the RANS turbulence model. Finally, the stress on the blades is investigated using the FEM with the RANS analysis results as the boundary conditions.. |
| 52. | Takanori Uchida, Designed Wind Speed Evaluation Technique in Wind Turbine Installation Point by Using the Meteorological Model and CFD Model, Journal of Flow Control, Measurement & Visualization, 10.4236/jfcmv.2018.63014, 06, 03(2018), Article ID:85916,17 pages, 2018.07, It is highly important in Japan to choose a good site for wind turbines, becausethe spatial distribution of wind speed is quite complicated over steepcomplex terrain. We have been developing the unsteady numerical model called the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain).The RIAM-COMPACT is based on the LES (Large-Eddy Simulation).The object domain of the RIAM-COMPACT is from several m to several km,and can predict the airflow and gas diffusion over complex terrain with highprecision. In the present paper, the design wind speed evaluation technique inwind turbine installation point by using the mesoscale meteorological modeland RIAM-COMPACT CFD model was proposed. The design wind speed tobe used for designing WTGs can be calculated by multiplying the ratio of the mean wind speed at the hub-height to the mean upper-air wind speed at the inflow boundary, i.e. , the fractional increase of the mean hub-height wind speed, by the reduction ratio, R. The fractional increase of the mean hub-height wind speed was evaluated using the CFD simulation results. This method was proposed as Approach 1 in the present paper. A value of 61.9 m/s was obtained for the final design wind speed, Uh, in Approach 1. In the evaluation procedure of the design wind speed in Approach 2, neither the above-mentioned reduction rate, R, nor an upper-air wind speed of 1.7 Vo, where Vo is the reference wind speed, was used. Instead, the value of the maximum wind speed which was obtained from the typhoon simulation for each of the investigated wind directions was adopted. When the design wind speed was evaluated using the 50-year recurrence value, the design wind speed was 48.3 m/s. When a somewhat conservative safety factor was applied, that is, when the 100 year recurrence value was used instead, the design wind speed was 52.9 m/s.. |
| 53. | Takanori Uchida, LES Investigation of Terrain-Induced Turbulence in Complex Terrain and Economic Effects of Wind Turbine Control, Energies, 10.3390/en11061530, 11(6), 1530, 15pages, 2018.06, 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.. |
| 54. | Takanori Uchida, Computational Fluid Dynamics Approach to Predict the Actual Wind Speed over Complex Terrain, Energies, 10.3390/en11071694, 11(7), 1694, 13pages, 2018.06, This paper proposes a procedure for predicting the actual wind speed for flow over complex terrain with CFD. It converts a time-series of wind speed data acquired from field observations into a time-series data of actual scalar wind speed by using non-dimensional wind speed parameters, which are determined beforehand with the use of CFD output. The accuracy and reproducibility of the prediction procedure were investigated by simulating the flow with CFD with the use of high spatial resolution (5 m) surface elevation data for the Noma Wind Park in Kagoshima Prefecture, Japan. The errors of the predicted average monthly wind speeds relative to the observed values were less than approximately 20%.. |
| 55. | Takanori Uchida, Computational Fluid Dynamics (CFD) Investigation of Wind Turbine Nacelle Separation Accident over Complex Terrain in Japan, Energies, 10.3390/en11061485, 11(6), 1485, 13pages, 2018.06, We have developed an unsteady and non-linear wind synopsis simulator called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain) to simulate the airflow on a micro scale, i.e., a few tens of km or less. In RIAM-COMPACT, the large-eddy simulation (LES) has been adopted for turbulence modeling. LES is a technique in which the structures of relatively large eddies are directly simulated and smaller eddies are modeled using a sub-grid scale model. In the present study, we conducted numerical wind diagnoses for the Taikoyama Wind Farm nacelle separation accident in Japan. The simulation results suggest that all six wind turbines at Taikoyama Wind Farm are subject to significant influence from separated flow (terrain-induced turbulence) which is generated due to the topographic irregularities in the vicinity of the wind turbines. A proposal was also made on reconstruction of the wind farm.. |
| 56. | Takanori UCHIDA, A New Proposal for Vertical Extrapolation of Offshore Wind Speed and an Assessment of Offshore Wind Energy Potential for the Hibikinada Area, Kitakyushu, Japan, Energy and Power Engineering, 10.4236/epe.2018.104011, 10, 04, 11pages, 2018.04, The author’s research group has been conducting research on applications of various meteorological Grid Point Value (GPV) data offered by the Japan Meteorological Agency (JMA) to the field of wind power generation. In particular, the group’s research has been focusing on the following areas: 1) the use of GPV data from the JMA Meso-Scale Model (MSM-S; horizontal resolution: 5 km) and the JMA Local Forecast Model (LFM-S; horizontal resolution: 2 km), and 2) examinations of the prediction accuracy of local wind assessment with the use of these data. In both the MSM-S and the LFM-S, grid points are fixed at 10 m above the sea (ground) surface. The purpose of the present study is to establish a method in which the values of the MSM-S and LFM-S wind speed data from the 10 m height are used as the reference wind speed and are, using a power law, vertically extrapolated to 80 to 90 m heights, typical hub-heights of offshore wind turbines. For this purpose, the present study examined time-averaged vertical profiles of wind speed over the ocean based on the MSM-S data and in-situ data in the Hibikinada area, Kitakyushu City, Fukuoka Prefecture, Japan. As a result, it was revealed that a strong wind shear existed close to the sea surface, between the 10 and 30 m heights. In order to address the above-mentioned wind shear, a two-step vertical extrapolation method was proposed in the present study. In this method, two values of N, specifically for low and high altitudes (below and above approximately 30 m, respectively), were calculated and used. The data were created for the five years between 2012 and 2016. Similarly to previous analyses, the analysis of the created data set indicated that the potential of offshore wind power generation in the Hibikinada area, Kitakyushu City is quite high.. |
| 57. | Takanori UCHIDA, Large-Eddy Simulation and Wind Tunnel Experiment of Airflow over Bolund Hill, Open Journal of Fluid Dynamics, 10.4236/ojfd.2018.81003, 8, 30-43, 2018.03. |
| 58. | Yasushi KAWASHIMA, Takanori UCHIDA, Effects of Terrain-Induced Turbulence on Wind Turbine Blade Fatigue Loads, Energy and Power Engineering, 10.4236/epe.2017.913053, 9, 843-857, 2017.12. |
| 59. | Takanori UCHIDA, Three-Dimensional Numerical Simulation of Stably Stratified Flows over a Two-Dimensional Hill, Open Journal of Fluid Dynamics, 10.4236/ojfd.2017.74039, 7, 579-595, 2017.12. |
| 60. | Takanori UCHIDA, High-Resolution LES of Terrain-Induced Turbulence around Wind Turbine Generators, by Using Turbulent Inflow Boundary Conditions, Open Journal of Fluid Dynamics, 10.4236/ojfd.2017.74035, 7, 511-524, 2017.12. |
| 61. | Takanori UCHIDA, CFD Prediction of the Airflow at a Large-Scale Wind Farm above a Steep, Three-Dimensional Escarpment, Energy and Power Engineering, 10.4236/epe.2017.913052, 9, 829-842, 2017.12. |
| 62. | Takanori UCHIDA, High-Resolution Micro-Siting Technique for Large Scale Wind Farm Outside of Japan Using LES Turbulence Model, Energy and Power Engineering, 10.4236/epe.2017.912050, 9, 802-813, 2017.12. |
| 63. | 川島泰史, 内田孝紀, 清木荘一郎, 近藤勝俊, 非定常乱流モデルLESによる地形性乱流診断と風車安全管理への応用, 九州大学応用力学研究所所報, 153, 96-103, 2017.10. |
| 64. | 谷山賀浩, 内田孝紀, 松島啓二, 中野三知子, 猪俣麻子, 田中翔一, 大迫俊樹, 実風車模型を用いた風車後流風速の定量的評価手法の開発, 日本風力エネルギー学会論文集, 153, 104-108, 2017.10. |
| 65. | 川島 泰史, 内田 孝紀, 清木 荘一郎, 近藤 勝俊, 地形性乱流が風車構造強度に与える影響に関する研究（非定常乱流モデルLESによる地形性乱流診断）, 日本風力エネルギー学会論文集, 41, 2, 17-24, 2017.08, 本研究は，山岳地形に建設された風車を研究対象に取り上げ，風車運転データとブレード歪みデータの同時計測を行い，風車ブレードに与える地形性乱流の相関性を定量的に明らかにするとともに高解像度数値風況シミュレーションによってこれらの地形性乱流が形成されるメカニズムを明らかにした．さらに，風車ブレードの疲労蓄積割合を評価する手法を明らかにし，多くの重要な新しい知見を得た．. |
| 66. | 内田孝紀, 風車ウエイクのラージ・エディ・シミュレーション（LES）, 計算工学, 22, 3, 3613-3617, 2017.05. |
| 67. | Koichi Watanabe, Yuji Ohya, Takanori Uchida, Tomoyuki Nagai, Numerical Prediction and Field Verification Test of Wind-Power Generation Potential in Nearshore Area Using a Moored Floating Platform, Journal of Flow Control, Measurement & Visualization, 10.4236/jfcmv.2017.52002, 5, 2, 2017.04. |
| 68. | Yuji Ohya, Masaki Wataka, Koichi Watanabe, Takanori Uchida, Laboratory experiment and numerical analysis of a new type of solar tower efficiently generating a thermal updraft, energies, 10.3390/en9121077, 9, 12, 2016.12, [URL], A new type of solar tower was developed through laboratory experiments and numerical analyses. The solar tower mainly consists of three components. The transparent collector area is an aboveground glass roof, with increasing height toward the center. Attached to the center of the inside of the collector is a vertical tower within which a wind turbine is mounted at the lower entry to the tower. When solar radiation heats the ground through the glass roof, ascending warm air is guided to the center and into the tower. A solar tower that can generate electricity using a simple structure that enables easy and less costly maintenance has considerable advantages. However, conversion efficiency from sunshine energy to mechanical turbine energy is very low. Aiming to improve this efficiency, the research project developed a diffuser-type tower instead of a cylindrical tower, and investigated a suitable diffuser shape for practical use. After changing the tower height and diffuser open angle, with a temperature difference between the ambient air aloft and within the collector, various diffuser tower shapes were tested by laboratory experiments and numerical analyses. As a result, it was found that a diffuser tower with a semi-open angle of 4° is an optimal shape, producing the fastest updraft at each temperature difference in both the laboratory experiments and numerical analyses. The relationships between thermal updraft speed and temperature difference and/or tower height were confirmed. It was found that the thermal updraft velocity is proportional to the square root of the tower height and/or temperature difference.. |
| 69. | Yuji Ohya, Masaki Wataka, Koichi Watanabe, Takanori Uchida, Laboratory Experiment and Numerical Analysis of a New Type of Solar Tower Efficiently Generating a Thermal Updraft, ENERGIES, 10.3390/en9121077, 9, 12, 2016.12, A new type of solar tower was developed through laboratory experiments and numerical analyses. The solar tower mainly consists of three components. The transparent collector area is an aboveground glass roof, with increasing height toward the center. Attached to the center of the inside of the collector is a vertical tower within which a wind turbine is mounted at the lower entry to the tower. When solar radiation heats the ground through the glass roof, ascending warm air is guided to the center and into the tower. A solar tower that can generate electricity using a simple structure that enables easy and less costly maintenance has considerable advantages. However, conversion efficiency from sunshine energy to mechanical turbine energy is very low. Aiming to improve this efficiency, the research project developed a diffuser-type tower instead of a cylindrical tower, and investigated a suitable diffuser shape for practical use. After changing the tower height and diffuser open angle, with a temperature difference between the ambient air aloft and within the collector, various diffuser tower shapes were tested by laboratory experiments and numerical analyses. As a result, it was found that a diffuser tower with a semi-open angle of 4 degrees is an optimal shape, producing the fastest updraft at each temperature difference in both the laboratory experiments and numerical analyses. The relationships between thermal updraft speed and temperature difference and/or tower height were confirmed. It was found that the thermal updraft velocity is proportional to the square root of the tower height and/or temperature difference.. |
| 70. | 川島 泰史, 内田 孝紀, 風車構造強度条件を考慮した風車最適配置手法の確立を目的としたリアムコンパクトによる数値風況診断ー串木野れいめい風力発電所を例としてー（その１：計測データによる風況特性分析）, 九州大学応用力学研究所所報, 151, 24-34, 2016.09. |
| 71. | 内田 孝紀, 洋上風況調査における風速の高度補正に関する新しい提案 その１, 九州大学応用力学研究所所報, 151, 4-11, 2016.09. |
| 72. | 内田 孝紀, 洋上風況調査における風速の高度補正に関する新しい提案 その2, 九州大学応用力学研究所所報, 151, 12-23, 2016.09. |
| 73. | 大城善郎, 伊藤芳樹, 内田孝紀, 高桑 晋, 勝呂幸男, 松澤幸一, 相原雅彦, 太田健一郎, 複雑地形における超音波風向風速計を用いた高精度風況解析(青森県下北地方岩屋ウィンドファームにおける2高度観測), 日本風力エネルギー学会論文集, 40, 2, 7-12, 2016.08. |
| 74. | 内田孝紀, LESによる陸上および洋上ウィンドファームの数値風況予測, ターボ機械, 44, 7, 15-22, 2016.07. |
| 75. | Masahiro Konishi, Shigeyuki Sako, Takanori Uchida, Ryou Araya, Koui Kim, Yuzuru Yoshii, Mamoru Doi, Kotaro Kohno, Takashi Miyata, Kentaro Motohara, Masuo Tanaka, Takeo Minezaki, Tomoki Morokuma, Yoichi Tamura, Toshihiko Tanabé, Natsuko Kato, Takafumi Kamizuka, Hidenori Takahashi, Tsutomu Aoki, Takao Soyano, Ken'ichi Tarusawa, The University of Tokyo Atacama Observatory 6.5m Telescope: enclosure design and wind analysis, Ground-based and Airborne Telescopes VI, https://doi.org/10.1117/12.2231385, 9906, 2016.07, We present results on the computational fluid dynamics (CFD) numerical simulations as well as the wind tunnel experiments for the observation facilities of the University of Tokyo Atacama Observatory 6.5m Telescope being constructed at the summit of Co. Chajnantor in northern Chile. Main purpose of this study starting with the baseline design reported in 2014 is to analyze topographic effect on the wind behavior, and to evaluate the wind pressure, the air turbulence, and the air change (ventilation) efficiency in the enclosure. The wind velocity is found to be accelerated by a factor of ~ 1.2 to reach the summit (78 m sec-1 expected at a maximum), and the resulting wind pressure (3,750 N m-2) is used for the framework design of the facilities. The CFD data reveals that the open space below the floor of the facilities works efficiently to drift away the air turbulence near the ground level which could significantly affect the dome seeing. From comparisons of the wind velocity field obtained from the CFD simulation for three configurations of the ventilation windows, we find that the windows at a level of the telescope secondary mirror have less efficiency of the air change than those at lower levels. Considering the construction and maintenance costs, and operation procedures, we finally decide to allocate 13 windows at a level of the observing floor, 12 at a level of the primary mirror, and 2 at the level of the secondary mirror. The opening area by those windows accounts for about 14% of the total interior surface of the enclosure. Typical air change rate of 20-30 per hour is expected at the wind velocity of 1 m sec-1.. |
| 76. | Takanori Uchida, Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software ―Report 1: For the Case of a Three-Dimensional Isolated-Hill With Steep Slopes―, Reports of Research Institute for Applied Mechanics, Kyushu University, 150, 47-59, 2016.03. |
| 77. | Takanori Uchida, Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software ―Report 3: For the Case of a Three-dimensional Cube―, Reports of Research Institute for Applied Mechanics, Kyushu University, 150, 71-83, 2016.03. |
| 78. | 内田 孝紀, 風車ウエイクの「ゆらぎ」に関するアクチュエータラインモデルを用いた高解像度LES, 九州大学応用力学研究所所報, 150, 25-33, 2016.03. |
| 79. | 内田孝紀, 青柳達郎, 渡邊文人, 見上伸, 市販CFDソフトウエアによる実地形を対象とした年間平均風速および年間発電電量の予測, 九州大学応用力学研究所所報, 150, 34-39, 2016.03. |
| 80. | 内田 孝紀, 「みなと100年公園」において2014年12月1日に発生した小型風車の破損事故に関する風況調査, 九州大学応用力学研究所所報, 150, 40-46, 2016.03. |
| 81. | 内田 孝紀, 非定常乱流モデルLESによる地形性乱流の数値的再現性, 日本風力エネルギー学会論文集, 39, 4, 53-60, 2016.02. |
| 82. | 内田 孝紀, LESによる数値風況診断に基づいた風車制御とその経済効果, 日本風力エネルギー学会論文集, 39, 4, 61-68, 2016.02. |
| 83. | Shinsuke Okada, Takanori Uchida, Takashi Karasudani, Yuji Ohya, Improvement in Solar Chimney Power Generation by Using a Diffuser Tower, Journal of Solar Energy Engineering, Transactions of the ASME, 138, 1, 2016.02. |
| 84. | 川島 泰史, 内田 孝紀, 複雑地形における気象庁局地数値予報モデルデータ(LFM)を用いた簡易風況推定法の試み－串木野れいめい風力発電所を例として－, 九州大学応用力学研究所所報, 149, 51-63, 2015.09. |
| 85. | 内田 孝紀, 風車および小規模地形の周辺流れに対する温度成層の効果－その1：流れ場の可視化－, 九州大学応用力学研究所所報, 149, 85-90, 2015.09. |
| 86. | 内田 孝紀, 福岡市博多湾を対象にした気象GPVデータによる洋上風況解析, 九州大学応用力学研究所所報, 149, 64-71, 2015.09. |
| 87. | 内田 孝紀, 渡邊 文人, 見上 伸, 市販CFDソフトウエアによる急峻な3次元孤立峰を対象とした気流場解析(第2報), 九州大学応用力学研究所所報, 149, 72-77, 2015.09. |
| 88. | 内田 孝紀, 鵜沢 憲, 大規模な崖状地形に建設された大型ウインドファームを対象とした気流場解析, 九州大学応用力学研究所所報, 149, 78-84, 2015.09. |
| 89. | Takanori Uchida, Fumihito Watanabe, Shin Mikami, Analysis of the Airflow Field around a Steep, Three-dimensional Isolated Hill with Commercially Available CFD Software, Reports of Research Institute for Applied Mechanics, Kyushu University, 149, 91-98, 2015.09. |
| 90. | Shinsuke Okada, Takanori Uchida, Takashi Karasudani, Yuji Ohya, Improvement in Solar Chimney Power Generation by Using a Diffuser Tower, JOURNAL OF SOLAR ENERGY ENGINEERING-TRANSACTIONS OF THE ASME, 10.1115/1.4029377, 137, 3, 2015.06, The solar chimney prototype, operated in Spain from 1982 to 1989, verified the concept of the solar chimney. The power generation mechanism in this system is to turn the wind turbine placed inside a high rise cylindrical hollow tower by an induced thermal updraft. As long as the thermal updraft is induced inside the tower by the solar radiation, this system can produce electricity. The disadvantage of this system is the low power generation efficiency compared to other solar energy power generation systems. To overcome this disadvantage, we improved the mechanism in order to augment the velocity of the air which flows into the wind turbine. By applying a diffuser tower instead of a cylindrical one, the efficiency of the systems power generation is increased. The mechanism that we investigated was the effect of the diffuser on the solar chimney structure. The inner diameter of the tower expands as the height increases so that the static pressure recovery effect of the diffuser causes a low static pressure region to form at the bottom of the tower. This effect induces greater airflow within the tower. The laboratory experiment, as does the computational fluid dynamics (CFD) analysis of the laboratory sized model, shows that the proposed diffuser type tower induces a velocity approximately 1.38-1.44 times greater than the conventional cylindrical type. The wind power generation output is proportional to the cube of the incoming wind velocity into the wind turbine; therefore, approximately 2.6-3.0 times greater power output can be expected from using the diffuser type tower.. |
| 91. | 内田 孝紀, 太鼓山風力発電所のナセル落下事故に対する数値流体力学的アプローチによる一考察, 日本風力エネルギー学会論文集, 39, 1, 6-13, 2015.05, [URL], Because a significant portion of the topography in Japan is characterized by steep, complex terrain, which results in a complex spatial distribution of wind speed, great care is necessary for selecting a site for the construction of Wind Turbine Generators (WTGs). We have developed a Computational Fluid Dynamics (CFD) model for unsteady flow called Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain (RIAM-COMPACT®). The RIAM-COMPACT® CFD model is based on Large-Eddy Simulation (LES) technique. In this paper, the numerical wind simulation over the Taikoyama wind farm was executed using the high resolution elevation data. As a result, the numerical results obtained suggest that all the six WTGs in the Taikoyama wind farm are subject to significant influence from separated flow (terrain-induced turbulence) which is due to the topographical irregularity in front of WTGs. The new reproduction proposal was also done.. |
| 92. | Kenji Kitagawa, Shintaro Iwama, Sho Fukui, Yuuki Sunaoka, Hayato Yazawa, Atsushi Usami, Masaaki Naramoto, Takanori Uchida, Satoshi Saito, Hiromi Mizunaga, Effects of Components of the Leaf Area Distribution on Drag Relations for Cryptomeria Japonica and Chamaecyparis Obtusa, European Journal of Forest Research, 10.1007/s10342-014-0859-6, 134, 3, 403-414, 2015.05, The objectives of this study were to clarify the effects of components of the leaf area distribution on the drag coefficient of crowns and streamlining (e.g., leaf area index; LAI, outline of the crown shape, and clumpiness) and to contribute to the accumulation of data on drag relations by quantifying data for Chamaecyparis obtusa and Cryptomeria japonica. We conducted drag experiments while simultaneously capturing dynamic crown images for 28 Ch. obtusa crowns and 13 Cr. japonica crowns to analyze the relationships between the leaf area distribution components and drag coefficient or streamlining. The static drag coefficient increased with the LAI for Ch. obtusa and with decreasing clumpiness for Cr. japonica. The reduction rate of the static drag coefficient decreased with increasing clumpiness for Ch. obtusa and with a combination of increasing LAI and decreasing clumpiness for Cr. japonica. The reduction rate of the static drag coefficient had a clear relationship with the decreasing rate of the dynamic crown projected area of obstacles (foliage elements, branches, and stems) calculated from captured video images under windy conditions for Cr. japonica, while Ch. obtusa did not show clear relationship between them. The drag coefficients assuming non-porous crown; Cmax estimated by simple model combining LAI and clumpiness were approximately 1.0 in Ch. obtusa and 0.5 in Cr. japonica and were equivalent to the dynamic drag coefficients from video image under windy condition. The combination of LAI and clumpiness provided simple estimation for drag relations and enable to link crown structure to wind damage easier.. |
| 93. | 川島泰史, 内田孝紀, 風車の運転方法の検討を目的としたリアムコンパクトによる数値風況診断－串木野れいめい風力発電所を例として－, 九州大学応用力学研究所所報, 148.0, 59.0-65.0, 2015.03, [URL], Recently to become, in wind farm, which was built on the complex terrain, the operation rate is less than the initial expectations, ie, power generation output and significantly bad windmill, failure of windmill inside and outside (for example, failure of the yaw motor and Yogia, windmill blades cracks, etc.) of the problem of is actualized. The main cause is a change in the windmill most recent slight terrain relief becomes origin, I considered to be a disorder of the wind (terrain turbulence) generated from there. In response to this situation, in our group, precise numerical wind diagnosis by real terrain version RIAM-COMPACT® software (the window risk assessment) was performed, has conducted studies safe windmill operating method 1) . In this paper, under the cooperation of Kyudenko new energy Co., Ltd., Numerical wind conditions diagnosis for the purpose of examination of windmill 10 Unit method of operation of the kushikinoreimei wind farm (start the operation from November 2012) since the was carried out, we report on the results.. |
| 94. | 内田孝紀, 渡邊文人, 見上 伸, 市販CFDソフトウエアによる急峻な3次元孤立峰を対象とした気流場解析, 九州大学応用力学研究所所報, 148.0, 35.0-41.0, 2015.03, [URL], The present study compared the prediction accuracy of two CFD software packages for simulating a flow around a three-dimensional, isolated hill with a steep slope: 1) Meteodyn WT (turbulence model: k-L RANS), which is one of the leading commercially available CFD software packages in the wind power industry and 2) RIAM-COMPACT® (turbulence model: the standard Smagorinsky LES), which has been developed by the lead author of the present paper. Although the Reynolds number settings differed between the simulation with Meteodyn WT and that with RIAM-COMPACT®, distinct differences in the flow patterns were identified in the vicinity of the isolated hill (especially downstream of the isolated hill) between the flows simulated by the two CFD software packages.. |
| 95. | 内田孝紀, 福岡市内における陸上と洋上の風況特性, 九州大学応用力学研究所所報, 148.0, 51.0-58.0, 2015.03, [URL], In this report, we explain the outline of the weather GPV data such as MSM. Next, the onshore and offshore wind characteristics in Fukuoka-shi are reported through comparison between the actual measurement data (observed value) and the weather GPV data.. |
| 96. | 内田孝紀, 福岡市内における冬季の風況特性, 九州大学応用力学研究所所報, 148.0, 43.0-49.0, 2015.03, [URL], In this report, we explain the outline of the weather GPV data such as MSM. Next, the flow characteristics of winter in Fukuoka-shi are reported through comparison between the actual measurement data (observed value) and the weather GPV data.. |
| 97. | Takanori Uchida, An Examination of the Taikoyama Wind Farm Nacelle Separation Accident Using a CFD Approach, Reports of Research Institute for Applied Mechanics, Kyushu University, 148.0, 15.0-24.0, 2015.03, [URL], Because a significant portion of the topography of Japan is characterized by steep, complex terrain, which results in a complex spatial distribution of wind speed, great care is necessary for selecting a site for the construction of wind turbines. We have developed a computational fluid dynamics (CFD) model for unsteady flow called Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain (RIAM-COMPACT®). The RIAM-COMPACT® CFD model is based on the large-eddy simulation (LES) technique. In this paper, a numerical wind simulation for the Taikoyama Wind Farm is performed using high-resolution terrain elevation data. The results suggest that all six wind turbines at the Taikoyama Wind Farm are subject to significant influence from separated flow (terrain-induced turbulence) which is generated due to the topographical irregularities in the vicinity of the wind turbines. A proposal has been also made on reconstruction of the wind farm.. |
| 98. | Shinsuke Okada, Takanori Uchida, Takashi Karasudani, Yuji Ohya, Improvement in solar chimney power generation by using a diffuser tower, Journal of Solar Energy Engineering, Transactions of the ASME, 10.1115/1.4029377, 137, 3, 2015.01, [URL], The solar chimney prototype, operated in Spain from 1982 to 1989, verified the concept of the solar chimney. The power generation mechanism in this system is to turn the wind turbine placed inside a high rise cylindrical hollow tower by an induced thermal updraft. As long as the thermal updraft is induced inside the tower by the solar radiation, this system can produce electricity. The disadvantage of this system is the low power generation efficiency compared to other solar energy power generation systems. To overcome this disadvantage, we improved the mechanism in order to augment the velocity of the air which flows into the wind turbine. By applying a diffuser tower instead of a cylindrical one, the efficiency of the systems power generation is increased. The mechanism that we investigated was the effect of the diffuser on the solar chimney structure. The inner diameter of the tower expands as the height increases so that the static pressure recovery effect of the diffuser causes a low static pressure region to form at the bottom of the tower. This effect induces greater airflow within the tower. The laboratory experiment, as does the computational fluid dynamics (CFD) analysis of the laboratory sized model, shows that the proposed diffuser type tower induces a velocity approximately 1.38-1.44 times greater than the conventional cylindrical type. The wind power generation output is proportional to the cube of the incoming wind velocity into the wind turbine; therefore, approximately 2.6-3.0 times greater power output can be expected from using the diffuser type tower.. |
| 99. | Fumihito Watanabe, Takanori Uchida, Micro-siting of wind turbine in complex terrain Simplified fatigue life prediction of main bearing in direct drive wind turbines, Wind Engineering, 10.1260/0309-524X.39.4.349, 39, 4, 349-368, 2015.01, [URL], In Japan, 1,516 wind turbine accidents have been recorded between year 2004 and 2012, and 84% of them were for turbines in complex terrains. The longest downtime was associated with damage to main shafts or bearings with an average downtime of 5.7 months. Careful micro-siting in complex terrains can prevent these accidents from happening. The objective of the paper is to provide an intermediate step that allows consultants, developers, and wind farm owners to further evaluate micro-siting of wind turbines in complex terrains, prior to load simulations by manufacturers. The author developed a simplified method to predict fatigue life of a main rear bearing in direct drive wind turbines. The method uses hub-height 10 minutes wind data as an input. The validation with an actual accident showed practically good agreement of 12.0 years of the predicted life against 12.7 years of the actual life. The method was also applied to quantify the effect of a curtailment. The proposed curtailment increased the predicted life to 35.5 years with a relatively small range of wind speed at a direction of only I % frequency distribution. With the proposed method, it is possible to layout turbines where fatigue life of a main rear bearing is longer than its design life. The method can also be applied to existing turbines in order to spot turbines that require careful maintenance, and to consider an installation of condition monitoring system on a bearing.. |
| 100. | Takanori Uchida, Validation Testing of the Prediction Accuracy of the Numerical Wind Synopsis Prediction Technique RIAM-COMPACT for the Case of the Bolund Experiment-Comparison against a Wind-Tunnel Experiment-, Reports of Research Institute for Applied Mechanics, Kyushu University, 147, 7-14, 2014.09. |
| 101. | 内田 孝紀, 川島 泰史, 沿岸部における気象GPVデータを用いた簡易風況推定法の試み－鹿児島県内の風力発電所を例として－, 九州大学応用力学研究所所報, 第147号, pp.15-29, 2014.09. |
| 102. | 内田 孝紀, 川島 泰史, 山間部における気象GPVデータを用いた簡易風況推定法の試み－阿蘇車帰風力発電所を例として－, 九州大学応用力学研究所所報, 第147号, pp.31-43, 2014.09. |
| 103. | 内田 孝紀, 京都府太鼓山風力発電所で起きたナセル落下事故に対する数値風況面からの一考察, 九州大学応用力学研究所所報, 第147号, pp.45-50, 2014.09. |
| 104. | 内田孝紀; 丸山敬; 西嶋一欽, 技術および最新情報 複雑地形上の風車サイトを対象にした風況の品質管理と風車の運転制御, 風力エネルギー, 38.0, 3.0, 344.0-349.0, 2014.04, [URL]. |
| 105. | 内田 孝紀, 市販CFDソフトウェアによる複雑乱流場の再現性－第2報 急峻な傾斜角度を有する2次元峰地形を対象として－, 九州大学応用力学研究所所報, 第146号, pp.147-155, 2014.03. |
| 106. | 内田 孝紀, 市販CFDソフトウェアによる複雑乱流場の再現性－第3報 3次元立方体を対象として－, 九州大学応用力学研究所所報, 第146号, pp.157-167, 2014.03. |
| 107. | Yuji Ohya, Takanori Uchida, Tomoyuki Nagai, Near Wake of a Horizontal Circular Cylinder in Stably Stratified Flows, Open Journal of Fluid Dynamics, Vol.3, 311-320, 2013.12. |
| 108. | Kana Kamimura, Satoshi Saito, Hiroko Kinoshita, Kenji Kitagawa, Takanori Uchida, Hiromi Mizunaga, Analysis of wind damage caused by multiple tropical storm events in Japanese Cryptomeria japonica forests, Forestry, 10.1093/forestry/cpt011, 86, 4, 411-420, 2013.10, [URL], This study analyzed wind damage caused by tropical storms from 1991 to 2007 to Japanese forests mainly consisting of Cryptomeria japonica. Statistical analyses based on logistic regression and Cox regression models were conducted in relation to conditions at the forest and stand level. Known damage to forests managed by Kyusyu Rinsan Corporation (KR forests), located on the Kyushu Island, was analyzed at the forest level, using tropical storm characteristics such as air pressure, precipitation and periods when the forests were within the storm zone as predictors. Wind damage was also examined at the stand level (150 analysis points) using Cox regression models, according to stand age, site index, terrain conditions, management practices and wind velocity indicators (horizontal and vertical velocity vectors). The results indicated that at the forest level, higher maximum hourly wind speed and longer periods of >15 m s-1 of wind speed were significantly correlated to damage occurrence. At the stand level, indicators of upward vertical velocity, thinning treatments and site index were positively associated with the probability of wind damage. For instance, stands receiving higher upward vertical velocities and thinning treatment within 2 years were more likely to have reduced stability against tropical storms. Stands with higher and lower site index than average also showed lower stability in our analysis.. |
| 109. | 内田 孝紀, Graham Li, 市販CFDソフトウェアによる複雑乱流場の再現性－急峻な傾斜角度を有する3次元孤立地形を対象として－, 九州大学応用力学研究所所報, 第145号, pp.107-119, 2013.09. |
| 110. | 齊藤 俊哉, 内田 孝紀, 荒屋 亮, ロケット発射時における地上風の数値風況予測の試行, 九州大学応用力学研究所所報, 第145号, pp.121-126, 2013.09. |
| 111. | Kana Kamimura, Satoshi Saito, Hiroko Kinoshita, Kenji Kitagawa, 内田 孝紀, Hiromi Mizunaga, Analysis of wind damage caused by multiple tropical storm events in Japanese Cryptomeria japonica forests, Forestry, 10.1093/forestry/cpt011, 0, pp.1-10, 2013.04. |
| 112. | 内田 孝紀, 大屋 裕二, 複雑地形上の風車ウエイクのラージ・エディ・シミュレーション, JWE : 日本風工学研究会誌 : Journal of Wind Engineering, 10.14887/jaweam.2013.0.135.0, 135, 135-136, 2013.04, 乱流モデルLESとアクチュエータラインモデルを用いて，複雑地形上の風車ウエイクの再現を試みた．その結果，風車背後にはTip Vortexの存在が確認されるものの，地形からの剥離流の影響により，その構造は急速に崩れていることが示された．. |
| 113. | 内田 孝紀, 川島 泰史, 荒屋 亮, 気象GPVデータの風力発電分野への活用に関する検討, 九州大学応用力学研究所所報, 第144号, pp.33-40, 2013.03. |
| 114. | 内田 孝紀, 辰己 賢一, 川島 泰史, 荒屋 亮, メソ気象モデルWRF-ARWを用いた複雑地形上の数値風況予測, 九州大学応用力学研究所所報, 第144号, pp.41-47, 2013.03. |
| 115. | 内田 孝紀, 高橋 幸平, 大屋 裕二, レンズ風車の相互干渉に関する数値的研究, 風工学シンポジウム論文集, pp.401-406, 2012.12. |
| 116. | 高橋 周平, 上野 祥彦, 大屋 裕二, 烏谷 隆, 内田 孝紀, 鍔付きディフューザ風車の性能に対する支持構造の後流の影響, 風工学シンポジウム論文集, pp.407-412, 2012.12. |
| 117. | Yuji Ohya, Takanori Uchida, Takashi Karasudani, Masaru Hasegawa, Hiroshi Kume, Numerical studies of flow around a wind turbine equipped with a flanged-diffuser shroud using an actuator-disk model, Wind Engineering, 10.1260/0309-524X.36.4.455, 36, 4, 455-472, 2012.08, [URL], Unsteady 3-D direct numerical simulations based on the finite-difference method (FDM) were performed for flow fields around a wind turbine equipped with a flanged-diffuser shroud. Generally, it is difficult to numerically simulate the flow around rotating bodies such as the blades of wind turbines because of the unsteady flow generated by moving bodies with complex geometry. Therefore, we have devised an actuator-disk model for a wind turbine for simulating the drag and rotational forces exerted on the fluid by the wind turbine. By incorporating the body forces derived from the actuator-disk model into the external terms in the Navier-Stokes equations, the unsteady flow around a wind turbine can be simulated. The results of numerical simulations were compared with the results from wind tunnel experiments and showed good agreement for the velocity and pressure fields.. |
| 118. | 内田 孝紀, 荒屋 亮, 福岡市街地を対象とした大規模数値風況予測, 九州大学応用力学研究所所報, 第142号, pp.55-62, 2012.03. |
| 119. | 内田 孝紀, 荒屋 亮, 実在市街地における浮力効果を考慮したガス拡散の大規模数値予測, 九州大学応用力学研究所所報, 第142号, pp.63-70, 2012.03. |
| 120. | 岡林 一木, 永山 圭憲, 原 智宏, 堀 英一, 大屋 裕二, 内田 孝紀, 複雑地形へ適用可能な非定常乱流拡散数値シミュレーション技術の開発 －風洞実験代替に向けて－, 三菱重工技報, Vol.49 No.1, pp.75-81, 2012.01. |
| 121. | Shuhei Takahashi, Yuya Hata, Yuji Ohya, Takashi Karasudani, Takanori Uchida, Behavior of the blade tip vortices of a wind turbine equipped with a brimmed-diffuser shroud, Energies, 10.3390/en5125229, 5, 12, 5229-5242, 2012.01, [URL], To clarify the behavior of the blade tip vortices of a wind turbine equipped with a brimmed-diffuser shroud, called a "Wind-Lens turbine", we conducted a three-dimensional numerical simulation using a large eddy simulation (LES). Since this unique wind turbine consists of not only rotating blades but also a diffuser shroud with a broad-ring brim at the exit periphery, the flow field around the turbine is highly complex and unsteady. Previously, our research group conducted numerical simulations using an actuator-disc approximation, in which the rotating blades were simply modeled as an externalforce on the fluid. Therefore, the detailed flow patterns around the rotating blades and the shroud, including the blade tip vortices, could not be simulated. Instead of an actuator-disc approximation, we used a moving boundary technique in the present CFD simulation to simulate the flow around a rotating blade in order to focus especially on blade tip vortices. The simulation results showed a pair ofvortices consisting of a blade tip vortex and a counter-rotating vortex which was generated between the blade tip and the inner surface of the diffuser. Since these vortices interacted with each other, the blade tip vortex was weakened by the counter-rotating vortex. The results showed good agreement with past wind tunnel experiments.. |
| 122. | 内田 孝紀，丸山 敬，大屋 裕二, 流体工学CFDモデルを用いた連続的な風向変化の再現性について, 日本風力エネルギー学会論文集, Vol.35，通巻99, pp.7-13, 2011.11. |
| 123. | Takanori UCHIDA, Takashi MARUYAMA, Tetsuya TAKEMI, Yuichiro OKU, Yuji OHYA and Graham Li, Proposal of Designed Wind Speed Evaluation Technique in WTG Installation Point by Using the Meteorological Model and CFD Model, 九州大学応用力学研究所所報, 第141号, pp.1-12, 2011.10. |
| 124. | Takanori UCHIDA, Takashi MARUYAMA, Hirohiko ISHIKAWA, Masaru ZAKO and Akira DEGUCHI, Investigation of the Causes of Wind Turbine Blade Damage at Shiratakiyama Wind Farm in Japan-A Computer Simulation Based Approach-, 九州大学応用力学研究所所報, 第141号, pp.13-25, 2011.10. |
| 125. | 内田孝紀, 大屋裕二, LES技術を用いたウインドファーム風況診断―熊本県阿蘇車帰風力発電所を例として―, 土木学会論文集A2（応用力学）Vol.67 特集号, 2011.09. |
| 126. | Takanori Uchida, Takashi Maruyama, Yuji Ohya, New evaluation technique for WTG design wind speed using a CFD-model-based unsteady flow simulation with wind direction changes, Modelling and Simulation in Engineering, 10.1155/2011/941870, 2011, 2011.09, [URL], Because a significant portion of the topography in Japan is characterized by steep, complex terrain, which results in a complex spatial distribution of wind speed, great care is necessary for selecting a site for the construction of wind turbine generators (WTG). We have developed a CFD model for unsteady flow called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, computational prediction of airflow over complex terrain). The RIAM-COMPACT CFD model is based on large eddy simulation (LES). The computational domain of RIAM-COMPACT can extend from several meters to several kilometers, and RIAM-COMPACT can predict airflow and gas diffusion over complex terrain with high accuracy. The present paper proposes a technique for evaluating the deployment location of a WTG. The proposed technique employs the RIAM-COMPACT CFD model and simulates a continuous wind direction change over 360 degrees.. |
| 127. | Takanori Uchida, Yuji Ohya, Kenichiro Sugitani, Comparisons between the wake of a wind turbine generator operated at optimal tip speed ratio and the wake of a stationary disk, Modelling and Simulation in Engineering, 10.1155/2011/749421, 2011, 1-7, 2011.04, [URL], The wake of a wind turbine generator (WTG) operated at the optimal tip speed ratio is compared to the wake of a WTG with its rotor replaced by a stationary disk. Numerical simulations are conducted with a large eddy simulation (LES) model using a nonuniform staggered Cartesian grid. The results from the numerical simulations are compared to those from wind-tunnel experiments. The characteristics of the wake of the stationary disk are significantly different from those of the WTG. The velocity deficit at a downstream distance of 10 D (D: rotor diameter) behind the WTG is approximately 30 to 40 of the inflow velocity. In contrast, flow separation is observed immediately behind the stationary disk (≤2D), and the velocity deficit in the far wake (10 D) of the stationary disk is smaller than that of the WTG. Copyright © 2011 Takanori Uchida et al.. |
| 128. | Takanori Uchida, Takashi Maruyama and Yuji Ohya, New Evaluation Technique for WTG Design Wind Speed using a CFD-model-based Unsteady Flow Simulation with Wind Direction Changes, Modelling and Simulation in Engineering, Volume 2011 (2011), 2011.03. |
| 129. | Takanori Uchida, Yuji Ohya and Kenichiro Sugitani, Comparisons Between The Wake Of A Wind Turbine Generator Operated At Optimal Tip Speed Ratio And The Wake Of A Stationary Disk, Modelling and Simulation in Engineering, Volume 2011 (2011), 2011.03. |
| 130. | 内田 孝紀, 丸山 敬, 大屋 裕二, 連続的な風向変化を考慮した非定常数値風況予測による 風車設置地点における設計風速評価手法の提案, 風力エネルギー協会誌, Vol.34，通巻96, pp.129-134, 2011.02. |
| 131. | Takanori Uchida, Yuji Ohya, Latest developments in numerical wind synopsis prediction using the RIAM-COMPACT® CFD model-design wind speed evaluation and wind risk (terrain-induced turbulence) diagnostics in Japan, Energies, 10.3390/en4030458, 4, 3, 458-474, 2011.01, [URL], Because a significant portion of the topography in Japan is characterized by steep, complex terrain, which results in a complex spatial distribution of wind speed, great care is necessary for selecting a site for the construction of Wind Turbine Generators (WTGs). We have developed a CFD model for unsteady flow called Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain (RIAM-COMPACT®). The RIAM-COMPACT® CFD model is based on Large-Eddy Simulation (LES) technique. The computational domain of RIAM-COMPACT® can extend from several meters to several kilometers, and RIAM-COMPACT® can predict airflow and gas diffusion over complex terrains with high accuracy. First, the present paper proposes a technique for evaluating the deployment location of WTGs. Next, wind simulation of an actual wind farm was executed using the high resolution elevation data. As a result, an appropriate point and an inappropriate point for locating WTGs were shown based on the numerical results obtained. This cause was found to be a topographical irregularity in front of WTGs.. |
| 132. | 内田孝紀, 水永博己, 齊藤哲, 上村佳奈, 木下裕子, 丸山敬, 森林・林業分野における風害シミュレーション-数値風況予測技術(CFD)によるアプローチ-, 森林立地学会誌，森林立地, 52(2)，pp.67-77, 2010.12. |
| 133. | 秦祐也, 上野祥彦, 大屋裕二, 烏谷隆, 内田孝紀, 鍔つきディフューザ風車まわりの流れ場の渦構造に関する数値的研究, 第21回風工学シンポジウム論文集, pp.227-232, 2010.12. |
| 134. | 内田孝紀, 大屋裕二, 最適周速比における風車後流のラージ・エディ・シミュレーション -第2報 格子乱流により生成した流入変動風の影響-, 第21回風工学シンポジウム論文集, pp.233-238, 2010.12. |
| 135. | 秦祐也, 永井修平, 大屋裕二, 辻美奈子, 内田孝紀, 烏谷隆, ポーラス状の外側部をもつ円柱まわりの流れの数値シミュレーション, 第21回風工学シンポジウム論文集, pp.257-262, 2010.12. |
| 136. | 丸山敬, 石川裕彦, 内田孝紀, 出口啓, メソスケールモデルとLESを用いたウインドファーム周辺の気流解析, 第21回風工学シンポジウム論文集, pp.209-214, 2010.12. |
| 137. | 川島 泰史, 内田 孝紀, 荒屋 亮, 藤本 弘明, 海外サイトにおける標高データ生成と 流体工学モデルRIAM-COMPACT®を用いた風力発電のマイクロサイティング, 風力エネルギー協会誌, Vol.34，通巻.95, pp.97-101, 2010.12. |
| 138. | 内田孝紀, 大屋裕二, Three-Dimensional Numerical Simulation of Stably Stratified Flows Over a Two-Dimensional Hill -Effect of Stratification on the Non-Periodic Separation and Reattachment of the Flow-, 九州大学応用力学研究所所報, 第139号, pp.57-68, 2010.09. |
| 139. | 内田孝紀, 烏谷隆, 大屋裕二, 九州大学伊都キャンパスにおける流体工学モデルRIAM-COMPACT®を用いた100kW級風レンズ風車のマイクロサイティング -第2報 土地造成を考慮した場合-, 九州大学応用力学研究所所報, 第139号, pp.131-139, 2010.09. |
| 140. | Takanori UCHIDA, Yuji OHYA, Challenge to Huge Computation of Airflow around Urban Area by using RIAM-COMPACT® CFD Model, Proceedings of EAEP2010/The 4th International Symposium on the Asian Environmental Problems, pp.191-194, 2010.09. |
| 141. | Hirotaka HANO, Takanori UCHIDA, Yuji OHYA, Wake Structure Behind Wind Turbine Generator in Turbulent Boundary Layer, Proceedings of EAEP2010/The 4th International Symposium on the Asian Environmental Problems, pp.195-200, 2010.09. |
| 142. | 内田 孝紀，丸山 敬，竹見 哲也，奥 勇一郎，大屋 裕二，李 貫行, 気象モデルと流体工学モデルを用いた 風車設置地点における設計風速評価手法の提案, 風力エネルギー協会誌, 第34巻, 第2号, pp.118-124, Vol.34，通巻94, pp.118-124, 2010.08, 本研究では，最新の研究成果を取り入れ，より現実に近い強風性状を取り込み，かつ，安全性を確保したものとして風車設置地点における耐風設計用の基準風速を求める一手法を提案したものである．. |
| 143. | Takanori Uchida and Yuji Ohya, Large-Eddy Simulation of Topography-Induced Turbulence around WTG by using the RIAM-COMPACT® CFD Model, Proceedings of RENEWABLE ENERGY 2010 (RE2010), 2010.06. |
| 144. | 内田 孝紀, 大屋 裕二, 李 貫行, 風車立地点近傍に発生する地形乱流の高解像度LES, 風力エネルギー協会誌, Vol.34，通巻.93, pp.121-126, 2010.05. |
| 145. | Yuji Ohya, Takanori Uchida, Near Wake of a Horizontal Circular Cylinder in Stably Stratified Flows, 九州大学応用力学研究所所報, 第138号, pp.13-22, 2010.03. |
| 146. | 内田孝紀, 烏谷隆, 大屋裕二, 九州大学伊都キャンパスにおける流体工学モデルRIAM-COMPACT®を用いた100kW級風レンズ風車のマイクロサイティング, 九州大学応用力学研究所所報, 第138号, pp.41-54, 2010.03. |
| 147. | Takanori Uchida and Yuji Ohya, HIGH RESOLUTION LES OF TURBULENT AIRFLOW OVER COMPLEX TERRAIN, Proceedings of Seventh Asia-Pacific Conference on Wind Engineering (APCWE-VII), pp.405-408, 2009.11. |
| 148. | Tomohiro Hara, Yuji Ohya, Takanori Uchida, Ryohji Ohba, Wind-tunnel and numerical simulations of the coastal thermal internal boundary layer, Boundary-Layer Meteorology, 10.1007/s10546-008-9343-5, 130, 3, 365-381, 2009.02, [URL], Wind-tunnel experiments in a thermally stratified wind tunnel and direct numerical simulations were performed to simulate the thermal internal boundary layer (TIBL) that developed over a coastal area in a sea-breeze flow. The results of the simulations were analyzed to investigate turbulence structure in the TIBL. To study the effects of the atmospheric stability over the sea on the TIBL, two vertical profiles of temperature were created in the upstream portion of the wind-tunnel experiment and the direct numerical simulation. Turbulence statistics of the TIBL changed significantly according to the temperature profile over the sea, indicating that the stability of the flow over the sea has a significant effect on the structure and turbulence characteristics of the TIBL. Furthermore, the TIBL heights were estimated from the vertical profiles of the local Richardson number. The estimated TIBL heights agreed with those predicted by a pre-existing relation, suggesting that both the wind-tunnel experiment and the direct numerical simulation accurately reproduced the growth of the TIBL.. |
| 149. | 内田 孝紀, 丸山 敬, 竹見 哲也, 大屋 裕二, 道下 和明, 複雑地形上の風車ハブ高さ風速に与える流入気流性状と標高データの影響, 第20回風工学シンポジウム論文集, pp.139-144, 2008.12. |
| 150. | 内田孝紀, 大屋裕二, 非定常・非線形風況シミュレータRIAM-COMPACTによる野間ウィンドパークを対象とした年間発電電力量の推定精度検証, 風力エネルギー協会誌, Vol.32，No.3, pp.122-131, 2008.11. |
| 151. | Takanori Uchida, Yuji Ohya, The wind risk management in the wind farm by using the RIAM-COMPACT CFD code, Proceedings of China Wind Power 2008 & Global Wind Power 2008, 2008.10. |
| 152. | Yuji Ohya, Takanori Uchida, Laboratory and numerical studies of the atmospheric stable boundary layers, Journal of Industrial Aerodynamics, 10.1016/j.jweia.2008.02.037, 96, 10-11, 2150-2160, 2008.10, [URL], We have investigated the turbulence phenomena of SBLs for a wide range of stability, particularly focusing on the effects of strong stratification on turbulent boundary layers, using a thermally stratified wind tunnel. In parallel with wind tunnel experiments, to understand the turbulence features and fluid dynamics in detail, we have also performed numerical simulations of SBLs under the boundary conditions similar to those in the wind tunnel experiments. The numerical studies based on a finite-difference method (FDM) are direct Navier-Stokes simulations without any turbulence model (DNS). Under the Boussinesq approximation, the governing equations consist of the Navier-Stokes, continuity and energy equations for 3D incompressible stratified flows. Stable stratification rapidly suppresses the fluctuations of streamwise velocity and temperature as well as the vertical velocity fluctuation. Momentum and heat fluxes are also significantly decreased with increasing stability and become nearly zero over the whole boundary depth of the boundary layer with very strong stability. From the flow visualization in both wind tunnel experiment and DNS, wave-like motions driven by buoyancy and waves due to the Kelvin-Helmholtz instability can be observed locally and intermittently in a SBL flow with strong stability.. |
| 153. | Takanori Uchida, Yuji Ohya, Micro-siting technique for wind turbine generators by using large-eddy simulation, Journal of Industrial Aerodynamics, 10.1016/j.jweia.2008.02.047, 96, 10-11, 2121-2138, 2008.10, [URL], It is highly important in Japan to choose a good site for wind turbines, because the spatial distribution of wind speed is quite complicated over steep complex terrain. We are developing the unsteady numerical model called the Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain (RIAM-COMPACT). The RIAM-COMPACT is based on the large-eddy simulation (LES). The object domain of the RIAM-COMPACT is from several m to several km. First, to test the accuracy of the RIAM-COMPACT we have performed an experimental and a numerical simulation of a non-stratified airflow past a two-dimensional ridge and a three-dimensional hill in a uniform flow. Attention is focused on airflow characteristics in a wake region. For this purpose, the velocity components were measured with a split-film probe (SFP) in the wind tunnel experiment. Through comparison of the experimental and numerical results, they showed a good agreement. The accuracy of both of the wind tunnel experiment by the SFP and also numerical simulation by the RIAM-COMPACT were confirmed. Next, we have applied the RIAM-COMPACT to the airflow over real complex terrain. The numerical results obtained by the RIAM-COMPACT demonstrated that the changes induced on the wind field by the topographic effect, such as the local wind acceleration and the flow separation, were successfully simulated. Furthermore, the forecast accuracy of an actual wind speed is examined.. |
| 154. | Yuji Ohya, Takanori Uchida, Laboratory and numerical studies of the atmospheric stable boundary layers, JOURNAL OF WIND ENGINEERING AND INDUSTRIAL AERODYNAMICS, 10.1016/j.jweia.2008.02.037, 96, 10-11, 2150-2160, 2008.10, We have investigated the turbulence phenomena of SBLs for a wide range of stability, particularly focusing on the effects of strong stratification on turbulent boundary layers, using a thermally stratified wind tunnel. In parallel with wind tunnel experiments, to understand the turbulence features and fluid dynamics in detail, we have also performed numerical simulations of SBLs under the boundary conditions similar to those in the wind tunnel experiments. The numerical studies based on a finite-difference method (FDM) are direct Navier-Stokes simulations without any turbulence model (DNS). Under the Boussinesq approximation, the governing equations consist of the Navier-Stokes, continuity and energy equations for 3D incompressible stratified flows. Stable stratification rapidly suppresses the fluctuations of streamwise velocity and temperature as well as the vertical velocity fluctuation. Momentum and heat fluxes are also significantly decreased with increasing stability and become nearly zero over the whole boundary depth of the boundary layer with very strong stability. From the flow visualization in both wind tunnel experiment and DNS, wave-like motions driven by buoyancy and waves due to the Kelvin-Helmholtz instability can be observed locally and intermittently in a SBL flow with strong stability. (C) 2008 Elsevier Ltd. All rights reserved.. |
| 155. | 内田孝紀, 杉谷賢一郎, 大屋裕二, 縮尺模型を用いた風洞実験における幾つかの知見, 応用力学論文集, Vol.11，pp.789-798, 2008.08. |
| 156. | Takanori Uchida, Yuji Ohya, Numerical Simulation of Airflow around Urban Area by using the RIAM-COMPACT CFD Model, Proceedings of Sino-Japan International Symposium on The East Asian Environmental Problems (EAEP2008), pp.51-53, 2008.08. |
| 157. | Yuji Ohya, Takanori Uchida, Laboratory and numerical studies of the atmospheric stable boundary layers, Journal of Wind Engineering & Industrial Aerodynamics, Vol.96, pp.2150-2160, 2008.07. |
| 158. | Takanori Uchida and Yuji Ohya, Micro-siting Technique for Wind Turbine Generators by Using Large-Eddy Simulation, Journal of Wind Engineering & Industrial Aerodynamics, Vol.96, pp.2121-2138, 2008.07, 風力業界で未解決課題であった風車に対する風の乱れ(ウィンドリスク)に対して，NEDO技術開発機構の産業技術研究助成事業(若手研究グラント)に採択され，研究代表者として3年間の研究開発を実施した．本研究を通じ，先端的数値風況予測モデル「RIAM-COMPACT®(リアムコンパクト)」を駆使し，世界で初めてウィンドリスクの存在を視覚的に特定することに成功し，その力学的機構を解明した．一連の成果に対し，2010年科学技術分野の文部科学大臣表彰･若手科学者賞を受賞した．. |
| 159. | T.Uchida and Y.Ohya, Verification of the Prediction Accuracy of Annual Energy Output at Noma Wind Park by the Non-Stationary and Non-Linear Wind Synopsis Simulator, RIAM-COMPACT, Journal of Fluid Science and Technology, Vol.3, No.3, pp.344-358, 2008.06. |
| 160. | Yuji Ohya, Reina Nakamura, Takanori Uchida, Intermittent bursting of turbulence in a stable boundary layer with low-level jet, Boundary-Layer Meteorology, 10.1007/s10546-007-9245-y, 126, 3, 349-363, 2008.03, [URL], The atmospheric stable boundary layer (SBL) with a low-level jet is simulated experimentally using a thermally stratified wind tunnel. The turbulence structure and flow characteristics are investigated by simultaneous measurements of velocity and temperature fluctuations and by flow visualization. Attention is focused on the effect of strong wind shear due to a low-level jet on stratified boundary layers with strong stability. Occasional bursting of turbulence in the lower portion of the boundary layer can be found in the SBL with strong stability. This bursting originates aloft away from the surface and transports fluid with relatively low velocity and temperature upward and fluid with relatively high velocity and temperature downward. Furthermore, the relationship between the occurrence of turbulence bursting and the local gradient Richardson number (Ri) is investigated. The Ri becomes larger than the critical Ri, Ricr = 0.25, in quiescent periods. On the other hand, the Ri number becomes smaller than Ricr during bursting events.. |
| 161. | Yuji Ohya, Reina Nakamura, Takanori Uchida, Intermittent Bursting of Turbulence in a Stable Boundary Layer with Low-Level Jet,Boundary-Layer Meteorology, Boundary-Layer Meteorology, vol.26, No.3, pp.349-363, 2008.01. |
| 162. | Takanori Uchida, High Resolution LES of Airflow over Complex Terrain, Proceedings of APCOM'07-EPMESC XI, 2007.12. |
| 163. | 内田孝紀, 大屋裕二, CFDを用いた複雑地形上の実風速推定法の提案, 応用力学論文集, Vol.10，pp.733-740, 2007.08. |
| 164. | 内田孝紀, 杉谷賢一郎, 大屋裕二, 九州大学伊都キャンパスにおける風観測 —第二報 乱流強度分布と発電電力量の試算について—, 九州大学応用力学研究所所報, 第131号，pp.13-24, 2006.12. |
| 165. | Takanori Uchida, Yuji Ohya, Application of LES technique to diagnosis of wind farm by using high resolution elevation data, JSME International Journal, Series B: Fluids and Thermal Engineering, 10.1299/jsmeb.49.567, 49, 3, 567-575, 2006.12, [URL], We are developing the numerical model called the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain). The object domain of this numerical model is from several m to several km, and can predict the airflow and the gas diffusion over complex terrain with high precision. The RIAM-COMPACT has already been marketed by certain tie-up companies. The estimation of the annual electrical power output is also possible now based on the field observation data. In the present study, wind simulation of an actual wind farm was executed using the high resolution elevation data. As a result, an appropriate point and an inappropriate point for locating a wind turbine generator were shown based on the numerical results obtained. This cause was found to be a topographical irregularity in front of the wind turbine generator.. |
| 166. | 友清衣利子, 内田孝紀, 前田潤滋, 台風0423号通過時における佐賀県小城市周辺風況のCFD解析—風況に及ぼす大気安定度の影響—, 第19回風工学シンポジウム論文集, pp.205-210, 2006.11. |
| 167. | 岡田臣右, 大屋裕二, 内田孝紀, 烏谷隆, 杉谷賢一郎, 垂直型集風構造体の集風性能評価, 第19回風工学シンポジウム論文集, pp.145-150, 2006.11. |
| 168. | 内田孝紀, 大屋裕二, 杉谷賢一郎, 最適周速比における風車後流と静止円盤後流の比較, 第19回風工学シンポジウム論文集, pp.187-192, 2006.11. |
| 169. | Takanori Uchida, Yuji Ohya, Diagnosis of Airflow Characteristics in Wind Farm by Using the Unsteady Numerical Model RIAM-COMPACT, Proceedings of Renewable Energy 2006, 2006.10. |
| 170. | 内田孝紀，大屋裕二，諏訪部哲也，李貫行, 非定常・非線形風況シミュレータRIAM-COMPACTによるウィンドファーム風況診断の提案, 日本風力エネルギー協会誌, Vol.30，No.2，通巻78，pp.101-108, 2006.09. |
| 171. | T.Uchida and Y.Ohya, Application of LES Technique to Diagnosis of Wind Farm by Using High Resolution Elevation Data, JSME International Journal, 「Environmental Flows」, Series B, Vol.49, No.3, pp.567-575, 2006.09. |
| 172. | 内田孝紀, 大屋裕二, 友清衣利子, 前田潤滋, 地形性強風の数値予測と格子解像度の影響, 応用力学論文集, Vol.9，pp.795-802, 2006.08. |
| 173. | Takanori Uchida, Yuji Ohya, Application of LES technique to diagnosis of wind farm by using high resolution elevation data, JSME INTERNATIONAL JOURNAL SERIES B-FLUIDS AND THERMAL ENGINEERING, 49, 3, 567-575, 2006.08, We are developing the numerical model called the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain). The object domain of this numerical model is from several m to several km, and can predict the airflow and the gas diffusion over complex terrain with high precision. The RIAM-COMPACT has already been marketed by certain tie-up companies. The estimation of the annual electrical power output is also possible now based on the field observation data. In the present study, wind simulation of an actual wind farm was executed using the high resolution elevation data. As a result, an appropriate point and an inappropriate point for locating a wind turbine generator were shown based on the numerical results obtained. This cause was found to be a topographical irregularity in front of the wind turbine generator.. |
| 174. | 丸山 敬，石川裕彦，内田孝紀，河井宏允，大屋裕二, 台風0418号通過時の宮島周辺の強風場に関する数値シミュレーション, 日本風工学会論文集, Vol.31，No.3，pp.95-104, 2006.07. |
| 175. | 友清衣利子, 喜多村美保, 内田孝紀, 前田潤滋, 台風0423号通過時における佐賀県小城市周辺風況のCFD解析, 九州大学大学院人間環境学研究院紀要, 第10号，pp.49-56, 2006.07. |
| 176. | 内田孝紀, 大屋裕二, 安定成層流体中の地形効果とその可視化, 可視化情報学会誌, 2006.07. |
| 177. | Takashi Maruyama, Hirohiko Ishikawa, Takanori Uchida, Hiromasa Kawai, Yuji Ohya, Numerical simulation of strong wind fields around Miyajima during Typhoon Songda in 2004, Journal of Wind Engineering, 10.5359/jwe.31.95, 31, 3, 95-104, 2006.07, [URL], The strong wind fields around Miyajima in Hiroshima Bay were simulated numerically during Typhoon Songda in 2004. The PSU/MM5 model was used for the calculation of meso-scale regions. The strong wind that passed over the sea and blew into the Hiroshima Bay was simulated. The predicted wind fields were examined comparing with the observed records and the field investigation of damage to the houses and buildings. The calculated maximum wind speed map was well correlated to the distribution of damage rate by strong wind in Hiroshima Prefecture. Unsteady wind fields were also computed by Large Eddy Simulation in the fine region around Miyajima. The local wind characteristics caused by the topography around Itsukushima Shrine were investigated by using the calculated results. The down flow along the valley to the south of the shrine was simulated in the neutral atmospheric condition. The local wind flows around the shrine were discussed.. |
| 178. | 丸山敬, 石川裕彦, 内田孝紀, 河井宏允, 大屋裕二, 台風0418 号通過時の宮島周辺の強風場に関する数値シミュレーション, 日本風工学会誌, 10.5359/jwe.31.95, 108, 95-104, 2006.07. |
| 179. | 内田孝紀, PC1台で動く，風が見える，風力発電適地選定ソフトRIAM-COMPACT(リアムコンパクト)の紹介, 日本機械学会論文集, Vol.109，No.1050，pp.411(トピックス記事), 2006.05. |
| 180. | 丸山敬，石川裕彦，内田孝紀，河井宏允，大屋裕二, 台風0418号通過時の厳島神社周辺における地形性の強風を再現する試み, 京都大学防災研究所年報, 第49号B, 2006.04. |
| 181. | 内田孝紀, 杉谷賢一郎, 大屋裕二, 一様流中に置かれた急峻な単純地形まわりの気流性状の評価—3次元孤立峰モデルの場合—, 日本風工学会論文集, Vol.31，No.2，pp.63-74, 2006.04. |
| 182. | Takanori Uchida, Kenichirou Sugitani, Yuji Ohya, Evaluation on wind characteristics around a steep simple terrain in a uniform flow-case of a three-dimensional isolated-hill model-, Journal of Wind Engineering, 10.5359/jwe.31.63, 31, 2, 63-74, 2006.04, [URL], The purpose of this research is to construct a database of a non-stratified airflow past a steep simple terrain under an imposition of a uniform flow, and, in addition, is to do the accuracy inspection of the numerical model under development at present. This numerical model is referred to as the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain), and is for the purpose of the prediction of airflow over complex terrain with several m to several km. This paper describes the experimental and numerical study of a non-stratified airflow past a three-dimensional isolated-hill in a uniform flow as the second phase. The Reynolds number, based on the uniform flow and the height of the hill, is about 104. Airflows around the hill include the unsteady vortex shedding. Attention is focused on airflow characteristics in a wake region. For this purpose, the velocity components in the streamwise direction were measured with a SFP (Split-Film Probe) in the wind tunnel experiment. In addition, the flow visualization was performed by using the smoke-wire technique. Through comparison of the experimental and numerical results, they showed a good agreement. The accuracy of both of the wind tunnel experiment by the SFP and also numerical simulation by the RIAM-COMPACT were confirmed as the result.. |
| 183. | 内田孝紀，大屋裕二, ナセル搭載の風向・風速計で測定された複雑地形上の風況特性—岬の場合—, 九州大学応用力学研究所所報, 第130号，pp.35-52, 2006.03. |
| 184. | 内田孝紀, 大屋裕二, 荒屋亮, 田辺正孝, 川島泰史, 非定常・非線形風況シミュレータRIAM-COMPACTを用いたウィンドファーム風況診断, 九州大学応用力学研究所所報, 第130号，pp.53-60, 2006.03. |
| 185. | Toshiyuki SANADA, Masatoshi FUJINO, Daisuke MATSUSHITA, Takanori UCHIDA, Hikaru MATSUMIYA, Masao WATANABE, Yoshinori HARA, Minori SHIROTA, Numerical Site Calibration on a Complex Terrain and its Application for Wind Turbine Performance Measurements, 2006 European Wind Energy Conference Proceedings, 2006.02. |
| 186. | 内田孝紀, 大屋裕二, 種々の安定成層場における三宅島火山ガスの挙動の数値シミュレーション, 第19回数値流体力学シンポジウム講演論文集(CD-ROM), 2005.12. |
| 187. | 弓本桂也, 吉田保衡, 鵜野伊津志, 内田孝紀, 大屋裕二, メソスケール気象モデルを用いた高解像度風況シミュレーション —野間岬風力発電サイトを対象として—, 日本気象学会 機関誌「天気」, VOL.52，NO.11，pp.11-17, 2005.11. |
| 188. | 内田孝紀, 大屋裕二, 風車単体後流の渦構造解明に向けた基礎的研究, 九州大学応用力学研究所所報, 第129号，pp.123-128, 2005.09. |
| 189. | 内田孝紀, 大屋裕二, 複雑地形上の非定常風況シミュレーションにおける流出境界断面の取扱いについて, 九州大学応用力学研究所所報, 第129号，pp.129-133, 2005.09. |
| 190. | 内田孝紀，大屋裕二，荒屋亮，田辺正孝，川島泰史, 風況シミュレーションのための紙地図からの高解像度地形データの構築, 九州大学応用力学研究所所報, 第129号，pp.135-141, 2005.09. |
| 191. | T. Uchida and Y. Ohya, Micro-siting Technique for Wind Turbine Generator by Using One PC -Introduction of the RIAM-COMPACT based on LES-, EXPO WCWRF 2005(CD-ROM), 2005.09. |
| 192. | 内田孝紀, 大屋裕二, 風車単体後流の渦構造解明に向けた基礎的研究, 日本流体力学会年会2005講演論文集(CD-ROM), 2005.09. |
| 193. | 内田孝紀, 大屋裕二, 一様流中の2次元崖状地形まわりの気流性状に関する数値的研究, 応用力学論文集, Vol.8，pp.831-838, 2005.08. |
| 194. | 内田孝紀, 杉谷賢一郎, 大屋裕二, 九州大学元岡キャンパスにおける風観測 —第一報 平均的な風況特性について—, 九州大学応用力学研究所所報, 第128号，pp.37-51, 2005.03. |
| 195. | 弓本桂也, 吉田保衡, 鵜野伊津志, 内田孝紀, 大屋裕二, 野間岬風力発電サイトを対象とした風況シミュレーション, 九州大学大学院総合理工学報告, Vol.26，No.1，pp.9-14, 2004.12. |
| 196. | 内田孝紀, 大屋裕二, 鵜野伊津志, 急峻な岬まわりの気流場シミュレーション —気象力学モデルCSU-RAMSと流体力学モデルRIAM-COMPACTの適用性について—, 第18回風工学シンポジウム論文集, pp.35-40, 2004.12. |
| 197. | 河内昭紀, 岡林一木, 吉門洋, 小林恵三, 北林興二, 大屋裕二, 内田孝紀, 井手靖雄, 複雑地形へ適用可能な拡散モデル開発の概要, 環境管理, Vol.40，No.10，pp.50-61, 2004.10. |
| 198. | 内田孝紀, 大屋裕二, 人工的に生成された流入変動風を用いた建物周辺流れのラージ・エディ・シミュレーション, 九州大学応用力学研究所所報, 第127号，pp.95-107, 2004.09. |
| 199. | 内田孝紀, 大屋裕二, 非定常風況・拡散シミュレータRIAM-COMPACTを用いた Intel Itanium2プロセッサとAMD Opteronプロセッサの性能評価 —自動並列計算の有効性について—, 九州大学応用力学研究所所報, 第127号，pp.87-94, 2004.09. |
| 200. | 内田孝紀, 大屋裕二, 田辺正孝, 葛生和人, 橋本晃慈, 実地形を対象に非定常風況・拡散シミュレータRIAM-COMPACTの紹介 —最新のインテルPentium4プロセッサ搭載Windows PCの製作と性能評価—, 九州大学応用力学研究所所報, 第127号，pp.75-86, 2004.09. |
| 201. | 内田孝紀, 大屋裕二, 田辺正孝, 葛生和人, 橋本晃慈, 風力タービン適地選定のためのコンピュータによるマイクロサイティング技術 —実地形を対象にした非定常風況・拡散シミュレータRIAM-COMPACTの紹介—, 日本風力エネルギー協会誌, Vol.28，No.2，pp.46-51, 2004.09. |
| 202. | 内田孝紀, 大屋裕二, 田辺正孝, 葛生和人, 橋本晃慈, 風力タービン適地選定のためのコンピュータによるマイクロサイティング技術 —実地形を対象にした非定常風況・拡散シミュレータRIAM-COMPACT—, 土木施工, Vol.45，No.8，pp.49-55, 2004.08. |
| 203. | Yuji Ohya, Takanori Uchida, Laboratory and numerical studies of the convective boundary layer capped by a strong inversion, Boundary-Layer Meteorology, 10.1023/B:BOUN.0000027913.22130.73, 112, 2, 223-240, 2004.08, [URL], The convective boundary layer (CBL) with a wide range of stability is simulated experimentally using a thermally stratified wind tunnel, and numerically by direct numerical simulation (DNS). The turbulence structures and flow characteristics of various CBL flows, capped by a strong temperature inversion and affected by surface shear, are investigated. The various vertical profiles of turbulence statistics similar to those from the observed CBL in the field are successfully simulated in both the wind-tunnel experiment and in DNS. The comparison of the wind-tunnel data and DNS results with those of atmospheric observations and water-tank studies shows the crucial dependence of the turbulence statistics in the upper part of the layer on the strength of the inversion layer, as well as the modification of the CBL turbulence regime by the surface shear.. |
| 204. | Yuji Ohya, Takanori Uchida, Laboratory and numerical studies of the convective boundary layer capped by a strong inversion, Boundary-Layer Meteorology, 10.1023/B:BOUN.0000027913.22130.73, 112, 2, 223-240, 2004.08, The convective boundary layer (CBL) with a wide range of stability is simulated experimentally using a thermally stratified wind tunnel, and numerically by direct numerical simulation (DNS). The turbulence structures and flow characteristics of various CBL flows, capped by a strong temperature inversion and affected by surface shear, are investigated. The various vertical profiles of turbulence statistics similar to those from the observed CBL in the field are successfully simulated in both the wind-tunnel experiment and in DNS. The comparison of the wind-tunnel data and DNS results with those of atmospheric observations and water-tank studies shows the crucial dependence of the turbulence statistics in the upper part of the layer on the strength of the inversion layer, as well as the modification of the CBL turbulence regime by the surface shear. © 2004 Kluwer Academic Publishers.. |
| 205. | Takanori Uchida, Kenichirou Sugitani, Yuji Ohya, Evaluation on wind characteristics around a steep simple terrain in a uniform flow-case of a two-dimensional ridge model, Journal of Wind Engineering, 10.5359/jwe.29.35, 29, 3, 35-43, 2004.07, [URL], The purpose of this research is to construct a database of a non-stratified airflow past a steep simple terrain under an imposition of a uniform flow, and, in addition, is to do the accuracy inspection of the numerical model under development at present. This numerical model is referred to as the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain), and is for the purpose of the prediction of airflow over complex terrain with several m to several km. This paper describes the experimental and numerical study of a non-stratified airflow past a two-dimensional ridge in a uniform flow as the first phase. The Reynolds number, based on the uniform flow and the height of the ridge, is about 10 4. Airflows around the ridge include the unsteady vortex shedding. Attention is focused on airflow characteristics in a wake region. For this purpose, the velocity components in the streamwise direction were measured with a SFP (Split-Film Probe) in the wind tunnel experiment. In addition, the flow visualization was performed by using the smoke-wire technique. Through comparison of the experimental and numerical results, they showed a good agreement. The accuracy of both of the wind tunnel experiment by the SFP and also numerical simulation by the RIAM-COMPACT were confirmed as the result.. |
| 206. | 内田孝紀, 杉谷賢一郎, 大屋裕二, 一様流中に置かれた急峻な単純地形まわりの気流性状の評価—2次元尾根モデルの場合—, 日本風工学会論文集, No.100，pp.35-43, 2004.06. |
| 207. | 内田孝紀, 大屋裕二, 風況シミュレータRIAM-COMPACTによる拡散場シミュレーション —その2．安定時の場合—, 九州大学応用力学研究所所報, 第126号，pp.17-22, 2004.03. |
| 208. | 内田孝紀, 大屋裕二, 風況シミュレータRIAM-COMPACTによる拡散場シミュレーション —その1．不安定時の場合—, 九州大学応用力学研究所所報, 第126号，pp.9-15, 2004.03. |
| 209. | 内田孝紀, 大屋裕二, PCクラスタを用いた風況予測シミュレータRIAM-COMPACTの開発 —その2．スカラー並列計算機における経過時間の比較—, 九州大学応用力学研究所所報, 第126号，pp.1-8, 2004.03. |
| 210. | 内田孝紀, 大屋裕二, 風況シミュレータRIAM-COMPACTによる インテルプロセッサ搭載ハイパフォーマンスコンピュータの性能評価, 九州大学情報基盤センター広報, Vol.4，No.1，pp.1-14, 2004.03. |
| 211. | 内田孝紀, 大屋裕二, 小規模なスカラー並列計算機を用いた風況予測シミュレータRIAM-COMPACTの実用化へ向けた開発, 九州大学情報基盤センター年報, 第4号，pp.19-30, 2004.03. |
| 212. | 内田孝紀, 大屋裕二, PCクラスタを用いた風況予測シミュレータRIAM-COMPACTの実用化へ向けた検討, 日本風力エネルギー協会誌, Vol.27，No.4，pp.8-11, 2003.10. |
| 213. | 内田孝紀, 大屋裕二, 風況予測シミュレータRIAM-COMPACTの開発 —風況精査とリアルタイムシミュレーション—, 日本流体力学会誌「ながれ」, Vol.22，pp.417-428, 2003.10. |
| 214. | 内田孝紀, 大屋裕二, PCクラスタを用いた風況予測シミュレータRIAM-COMPACTの開発 —その1．種々の計算機におけるCPU時間の比較—, 九州大学応用力学研究所所報, 第125号，pp.5-20, 2003.09. |
| 215. | 内田孝紀, 大屋裕二, 乱流境界層に埋没した急峻な孤立峰まわりの風況シミュレーション —大きな乱れ強さを伴う流入風の場合—, 応用力学論文集, Vol.6，pp.819-826, 2003.08. |
| 216. | 内田孝紀，大屋裕二, パッシブ粒子法による三宅島火山ガス挙動の可視化, 可視化情報学会論文集, Vol.23, No.7, pp.58-65, 2003.07. |
| 217. | 内田孝紀, ブラフボディまわりの数値流体シミュレーションとその可視化 —マイクロスケールからメソスケールまで—, 九州大学情報基盤センター広報, Vol.3, No.2，pp.135-148, 2003.07. |
| 218. | Yuji Ohya, Takanori Uchida, Turbulence structure of stable boundary layers with a near-linear temperature profile, Boundary-Layer Meteorology, 10.1023/A:1023069316164, 108, 1, 19-38, 2003.07, [URL], By using a thermally stratified wind tunnel, we have successfully simulated stably stratified boundary layers (SBL), in which the mean temperature increases upward almost linearly. We have investigated the flow structure and the effects of near-linear stable stratification on the transfer of momentum and heat. The vertical profiles of turbulence quantities exhibit different behaviour in two distinct stability regimes of the SBL flows with weak and strong stability. For weak stability cases, the turbulent transfer of momentum and heat is basically similar to that for neutral turbulent boundary layers, although it is weakened with increasing stability. For strong stability cases, on the other hand, the time-mean transfer is almost zero over the whole boundary-layer depth. However, the instantaneous turbulent transfer frequently occurs in both gradient and counter-gradient directions in the lower part of the boundary layer. This is due to the Kelvin-Helmholtz (K-H) shear instability and the rolling up and breaking of K-H waves. Moreover, the internal gravity waves are observed in the middle and upper parts of all stable boundary layers.. |
| 219. | 内田孝紀, 杉谷賢一郎, 大屋裕二, 一様流中の2次元崖状地形まわりの気流性状に関する実験的研究, 日本風工学会論文集, No.95，pp.233-244, 2003.04. |
| 220. | Takanori Uchida, Kenichiro Sugitani, Yuji Ohya, Experimental study on wind characteristics around a two-dimension escapement in a uniform flow, Journal of Wind Engineering, 28, 2, 233-244, 2003.04, This paper describes the experimental study of a non-stratified airflow past a two-dimensional escarpment in a uniform flow. The Reynolds number, based on the uniform flow and the height of the escarpment, is about 104. The slope gradient of the escarpment is 25, 35 and 45 degree. Airflows around the escarpment include the unsteady separated and reattaching flow (hereafter called a "separation bubble"), where the separation occurs from a sharp comer. Attention is focused on an influence of a surface roughness on airflow characteristics in a wake region. For this purpose, the velocity components were measured with the X-wire probe and the split-film probe. In addition, the flow visualization was performed by using the smoke-wire technique. Through comparison between the experimental results with the surface roughness and those without it, the significant difference in the airflow characteristics is confirmed in the separation bubble. This is mainly due to the size in the separation bubble. The size in the separation bubble with the surface roughness is much larger than the one without it. In the case under an imposition of the surface roughness, the velocity is strongly defected near the slope surface. As a result, the production of the vorticity in the separated shear layer is also inhibited, leading to the elongation in the separation bubble.. |
| 221. | 内田孝紀, 杉谷賢一郎，大屋裕二, 3次元数値モデルによる九大新キャンパスの風況予測シミュレーション —第2報建物群まわりの風環境予測—, 九州大学情報基盤センター年報, 第3号，pp.57-66, 2003.03. |
| 222. | T. Uchida, Y. Ohya, Stable stratification effect on the separated and reattaching flow behind two-dimensional topography, 九州大学応用力学研究所所報, 第124号，pp.17-24, 2003.03. |
| 223. | 内田孝紀, 大屋裕二, 安定成層流体中の地形効果に関する数値的研究 —急峻な孤立峰の場合—, 九州大学応用力学研究所所報, 第124号，pp.9-16, 2003.03. |
| 224. | 内田孝紀, 大屋裕二, 安定成層場における山越え気流の三次元数値シミュレーション —非定常な剥離—再付着流れに対する安定成層の効果—, 日本流体力学会誌「ながれ」, Vol.22，pp.65-78, 2003.02. |
| 225. | Takanori Uchida, Yuji Ohya, Large-eddy simulation of turbulent airflow over complex terrain, Journal of Industrial Aerodynamics, 10.1016/S0167-6105(02)00347-1, 91, 1-2, 219-229, 2003.01, [URL], In order to develop an overall efficient and accurate method of predicting an unsteady three-dimensional airflow over a complex terrain with characteristic length scales on the order of kilometers, we recently developed the CFD codes referred to as the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain). In this paper, we carried out the calculation of turbulent airflow over a real complex terrain in a horizontal region of 9.5km × 5km with a relatively fine spatial resolution of 50 m. This area covers the new campus of Kyushu University. In order to generate instantaneous velocity fluctuations in an approaching flow, an unsteady flow field in the driver unit is directly calculated. The numerical results obtained by RIAM-COMPACT demonstrated that the changes induced on the wind field by the topographic effect, such as the local wind acceleration and the flow separation, were successfully simulated. We emphasized the following features as regards the wind field over the new campus area under the influence of the west wind. The wind field in the central region of the new campus area is strongly influenced by the wake region generated behind Mt. Hiyama (244m). The airflows moving around the side of Mt. Hiyama exhibit relatively small fluctuations. When these winds reach Mt. Ishigateke (99m), they are locally accelerated at the peak. Consequently, a local speed-up effect is confirmed in the southern part of the new campus area.. |
| 226. | Takanori Uchida, Yuji Ohya, Large-eddy simulation of turbulent airflow over complex terrain, Journal of Wind Engineering and Industrial Aerodynamics, 10.1016/S0167-6105(02)00347-1, 91, 1-2, 219-229, 2003.01, In order to develop an overall efficient and accurate method of predicting an unsteady three-dimensional airflow over a complex terrain with characteristic length scales on the order of kilometers, we recently developed the CFD codes referred to as the RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, Computational Prediction of Airflow over Complex Terrain). In this paper, we carried out the calculation of turbulent airflow over a real complex terrain in a horizontal region of 9.5km × 5km with a relatively fine spatial resolution of 50 m. This area covers the new campus of Kyushu University. In order to generate instantaneous velocity fluctuations in an approaching flow, an unsteady flow field in the driver unit is directly calculated. The numerical results obtained by RIAM-COMPACT demonstrated that the changes induced on the wind field by the topographic effect, such as the local wind acceleration and the flow separation, were successfully simulated. We emphasized the following features as regards the wind field over the new campus area under the influence of the west wind. The wind field in the central region of the new campus area is strongly influenced by the wake region generated behind Mt. Hiyama (244m). The airflows moving around the side of Mt. Hiyama exhibit relatively small fluctuations. When these winds reach Mt. Ishigateke (99m), they are locally accelerated at the peak. Consequently, a local speed-up effect is confirmed in the southern part of the new campus area. © 2002 Elsevier Science Ltd. All rights reserved.. |
| 227. | 内田孝紀, 大屋裕二, 急峻な孤立地形まわりの複雑乱流場に対するLES —流入気流の乱れの影響について—, 第17回風工学シンポジウム論文集, pp.137-142, 2002.12. |
| 228. | 内田孝紀, 大屋裕二, 急峻な2次元山からの非定常な渦放出に対する数値的考察, 応用力学論文集, Vol.5，pp.735-742, 2002.08. |
| 229. | R. Ohba, T. Hara, S. Nakamura, Y. Ohya, Takanori Uchida, Gas diffusion over an isolated hill under neutral, stable and unstable conditions, Atmospheric Environment, 10.1016/S1352-2310(02)00642-8, 36, 36-37, 5697-5707, 2002.12, [URL], Wind tunnel experiments of gas diffusion were performed over flat terrain and over an isolated three-dimensional hill under neutral, stable and unstable (sea breeze) conditions. Conditions of airflow in the wind tunnel were determined so as to satisfy the similarity rule for the bulk Richardson number, by controlling temperature profiles and wind velocity of the thermally stratified wind tunnel. Typical characteristics were observed under each condition of atmospheric stability; reversed vortex behind the hill in neutral condition, downward slope wind in stable one and convective motion in unstable one. We compared these experiments with the results of a Direct Numerical Simulation (DNS) model for the wind velocity over the hill under neutral conditions. The numerical results showed good agreement with the experimental results.. |
| 230. | 内田孝紀, 大屋裕二, ネストグリッドを用いた複雑地形上の風況予測シミュレーション, 日本風工学会論文集, No.92，pp.135-144, 2002.07. |
| 231. | Takanori Uchida, Yuji Ohya, Numerical simulation on wind fields over complex terrain using nested grids, Journal of Wind Engineering, 27, 3, 135-144, 2002.07, In order to develop an overall efficient and accurate model of simulating an unsteady three-dimensional airflow over complex terrain with characteristic length scales on the order of kilometers, we have been examining the large-eddy simulation (LES) technique using a finite-difference method (FDM). These LES codes are referred to as the RIAM-COMPACT (Research institute for Applied Mechanics, Kyushu University, Computational JPrediction of Airflow over Complex Terrain), and are based on two grid systems and corresponding variable arrangement: one is an orthogonal staggered grid; the other is a generalized curvilinear collocated grid. In this paper, using the RIAM-COMPACT based on a generalized curvilinear collocated grid, we have performed the calculations of a non-stratified airflow over real complex terrain. This area covers the new campus of Kyushu University. To investigate more clearly the influence of topography on the airflow over the new campus area, we employed nested grids: the outer grid domain is 9.75km × 4.65km × 1.46km with horizontal grid spacing of 50m and the inner grid domain is 5km × 2.7km × 1.46km with horizontal grid spacing of 25m. The numerical results were compared with the wind tunnel experiment, and airflow characteristics due to the topographic effect, such as the wind speed-up and the separated flow, were successively simulated.. |
| 232. | 内田孝紀, 大屋裕二, 3次元数値モデルによる九大新キャンパスの風況シミュレーション, 九州大学情報基盤センター年報, 第2号，pp.99-106, 2002.03. |
| 233. | Takanori Uchida, Yuji Ohya, Numerical study of stably stratified flows over a two-dimensional hill in a channel of finite depth, Fluid Dynamics Research, 10.1016/S0169-5983(01)00025-9 , 29, 4, 227-250, 2001.10, [URL], Stably stratified flows over a two-dimensional hill in a channel of finite depth are investigated numerically at a Reynolds number of 2000, which is based on the uniform upstream velocity U and the hill height h. As a first step, we assume a free-slip condition on the ground, both upstream and downstream of the hill, and impose a no-slip condition only on the hill surface. Such a configuration corresponds to that of the previous towing tank experiments and numerical studies. For strong stratification (1 d on the hill, and emphasize the following features, where N is the buoyancy frequency and H is the channel depth. For 1.1 ≤ K ≤ 1.7, columnar disturbances with mode n = 1 are dominant so that the flow around the hill shows a persistent periodic unsteadiness. This flow unsteadiness is mainly due to the periodic shedding of upstream advancing columnar disturbances with mode n = 1 with a clockwise circulation. For 1.8 ≤ K ≤ 2, as columnar disturbances with mode n = 2 become dominant, the flow around the hill rapidly reaches an almost steady state. In addition, through the calculations with the blockage ratio H/h = 6, 10 and corresponding Re=20, 100 and 2000, it is found that the normalized periods of Cd oscillations have a strong dependence on both the H/h and Re. As a next step, to investigate the flow around the hill under real atmospheric situations, we have performed calculations under imposition of a no-slip condition on the ground, particularly focusing on the effect of stable stratifications on the unsteady separated-reattached flow behind the hill. The flow around the hill exhibits different behavior over the whole range 0 ≤ K ≤ 2, corresponding to the difference in the boundary condition on the ground. For 0 ≤ K ≤ 0.9, the vortex shedding from the separation bubble behind the hill occurs. For K = 1 and 1.1, the vortex shedding is strongly suppressed so that the flow around the hill rapidly reaches an almost steady state. Under strong stratification (1 |
| 234. | 内田孝紀, 大屋裕二, 九州大学新キャンパス移転地上の風況場に関する数値シミュレーションと風洞実験, 九州大学応用力学研究所所報, 第121号，pp.59-72, 2001.09. |
| 235. | 内田孝紀, 大屋裕二, ラージ・エディ・シミュレーションによる局地的風況予測モデルRIAM-COMPACTの評価 —その2．壁面せん断乱流場による計算精度の検証と実地形上の風況場解析への適用—, 九州大学応用力学研究所所報, 第121号，pp.73-85, 2001.09. |
| 236. | 内田孝紀, 大屋裕二, 安定成層場で風下波の砕波により発生する局地強風の数値的考察, 応用力学論文集, Vol.4，pp.655-664, 2001.08. |
| 237. | 内田孝紀, 大屋裕二, ラージ・エディ・シミュレーションによる局地的風況予測モデルの評価, 九州大学応用力学研究所所報, 第120号，pp.29-34, 2001.02. |
| 238. | 内田孝紀, 大屋裕二, 局地的風況予測に対するLESの適用, 第16回風工学シンポジウム論文集, pp.59-64, 2000.11. |
| 239. | 内田孝紀, 大屋裕二, 風工学における数値シミュレーションのための流入変動風の生成法に関する研究 その1. 一定圧力勾配により駆動される十分に発達した境界層流のLES, 九州大学応用力学研究所所報, 第119号，pp.75-81, 2000.09. |
| 240. | 内田孝紀, 大屋裕二, 丘陵地形を過ぎる安定成層流の2次元数値シミュレーション —格子分解能と数値計算法の影響について—, 応用力学論文集, Vol.3，pp.729-736, 2000.08. |
| 241. | 内田孝紀, 大屋裕二, 有限深さ流路内の2次元丘陵地形を過ぎる安定成層流の非定常性 —第2報 非定常な剥離・再付着流れに対する安定成層の効果—, 日本流体力学会誌「ながれ」, Vol.18，pp.308-320, 1999.11. |
| 242. | 内田孝紀, 大屋裕二, 定成層場で地形効果により発生する局地強風の数値的検討, 応用力学論文集, Vol.2，pp.583-592, 1999.08. |
| 243. | Takanori Uchida, Yuji Ohya, Numerical simulation of atmospheric flow over complex terrain, Journal of Industrial Aerodynamics, 10.1016/S0167-6105(99)00024-0, 81, 283-293, 1999.05, In order to develop an overall efficient and accurate method of simulating an unsteady three-dimensional atmospheric flow over topography, we examined two grid systems and corresponding variable arrangements: one is a body-fitted coordinate (BFC) grid system based on a collocated variable arrangement; the other is an orthogonal grid system based on a staggered variable arrangement. Using these codes, we calculated the wind system over topography such as an isolated hill and real complex terrain. Both codes remarkably removed the numerical difficulties such as the convergence of the SOR method in solving the pressure Poisson equation, resulting in numerical results with much higher accuracy. Despite the differences in the grid system and in variable arrangement, no significant differences in the flow pattern between the both numerical results were found. Compared with the previous studies, the numerical results obtained are very satisfactory in the sense that overall characteristic flows are successfully simulated irrespective of the simulation codes.. |
| 244. | 内田孝紀, 大屋裕二, 安定成層した非粘性流体中に地形効果により励起された内部重力波の数値解析, 九州大学応用力学研究所所報, 第85号，pp.13-22, 1999.02. |
| 245. | 内田孝紀, 大屋裕二, 山越え成層流の差分解析—地形形状と格子分解能の影響—, 九州大学応用力学研究所所報, 第85号，pp.1-12, 1999.02. |
| 246. | 内田孝紀, 大屋裕二, 一般曲線座標系・コロケート格子を用いた複雑地形上の安定成層流に関する風況予測, 第15回風工学シンポジウム論文集, pp.131-136, 1998.12. |
| 247. | 内田孝紀, 大屋裕二, 差分法における幾つかの知見, 九州大学応用力学研究所所報, 第84号，pp.29-43, 1998.10. |
| 248. | 内田孝紀, 大屋裕二, 有限流路内の2次元丘陵地形を過ぎる安定成層流のLES, 応用力学論文集, Vol.1，pp.615-624, 1998.08. |
| 249. | 内田孝紀，大屋裕二, 有限流路内の2次元丘陵地形を過ぎる安定成層流の非定常性, 日本流体力学会誌「ながれ」, Vol.17，pp.45-56, 1998.02. |
| 250. | 内田孝紀, 大屋裕二, DNSによる複雑地形を過ぎる安定成層流の風況予測 —矩形格子系と一般曲線格子系との比較—, 九州大学応用力学研究所所報, 第82号，pp.113-126, 1997.11. |
| 251. | Takanori Uchida, Y. Ohya, A numerical study of stably stratified flows over a two-dimensional hill - Part I. Free-slip condition on the ground, Journal of Industrial Aerodynamics, 10.1016/S0167-6105(97)00096-2, 67-68, 493-506, 1997.01, [URL], Stably stratified flows over a two-dimensional hill in a channel of finite depth are analyzed numerically by using a newly-developed multi-directional finite-difference method at a Reynolds number Re = 2000. To simplify the phenomena occurring in the flow around the hill, the free-slip condition for the velocity is assumed on the ground, and the nonslip condition is imposed only on the hill surface. Attention is focused on the unsteadiness in the flow around the hill for the cases of K( = NH/πU) > 1 where N and U are the buoyancy frequency and free-stream velocity and H is the domain depth. The flow unsteadiness is discussed, being associated with shedding of the upstream advancing columnar disturbance.. |
| 252. | 内田孝紀，大屋裕二，小園茂平, 3次元複雑地形を過ぎる安定成層流の数値シミュレーション, 第14回風工学シンポジウム論文集, pp.91-96, 1996.12. |
本データベースの内容を無断転載することを禁止します。