• Wind and Structures
• /
• v.38 no.4
• /
• pp.261-275
• /
• 2024

Before performing an experimental study on the downburst-generated wave, it is necessary to examine the scale effects and corresponding corrections or compensations. Analysis of similarity is conducted to conclude the non-dimensional force ratios that account for the dynamic similarity in the interaction of downburst with wave between the prototype and the scale model, along with the corresponding scale factors. The fractional volume of fluid (VOF) method in association with the impinging jet model is employed to explore the characteristics of the downburst-generated wave numerically, and the validity of the proposed scaling method is verified. The study shows that the location of the maximum radial wind velocity in a downburst-wave field is a little higher than that identified in a downburst over the land, which might be attributed to the presence of the wave which changes the roughness of the underlying surface of the downburst. The impinging airflow would generate a concavity in the free surface of the water around the stagnation point of the downburst, with a diameter of about two times the jet diameter (D_jet). The maximum wave height appears at the location of 1.5D_jet from the stagnation point. Reynolds number has an insignificant influence on the scale effects, in accordance with the numerical investigation of the 30 scale models with the Reynolds number varying from 3.85 × 10⁴ to 7.30 × 10⁹. The ratio of the inertial force of air to the gravitational force of water, which is denoted by G, is found to be the most significant factor that would affect the interaction of downburst with wave. For the correction or compensation of the scale effects, fitting curves for the measures of the downburst-wave field (e.g., wind profile, significant wave height), along with the corresponding equations, are presented as a function of the parameter G.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.275-279
• /
• 2005

The vortex convection and induced flow field behind the KARI 3m x 4m LSWT gust generator was computed by using Computational Fluid Dynamics. For the accurate simulation of vortex convection, inviscid, laminar, Spalart-Allmars k-e and k-w turbulence models were tested with the NAL gust generator configuration and Spalart-Allmaras turbulence model was selected for the prediction of induced flow field behind the KARI LSWT gust generator. The wind tunnel test was also carried out at KARI LSWT and the results were compared with CFD prediction.

• Wind and Structures
• /
• v.23 no.6
• /
• pp.559-575
• /
• 2016

In order to investigate the influence of different blade positions on aerodynamic load and wind loads and load-effects of large scale wind turbine tower under the halt state, we take a certain 3 MW large scale horizontal axis three-blade wind turbine as the example for analysis. First of all, numerical simulation was conducted for wind turbine flow field and aerodynamic characteristics under different halt states (8 calculating conditions in total) based on LES (large eddy simulation) method. The influence of different halt states on the average and fluctuating wind pressure coefficients of turbine tower surface, total lift force and resistance coefficient, circular flow and wake flow characteristics was compared and analysed. Then on this basis, the time-domain analysis of wind loads and load-effects was performed for the wind turbine tower structure under different halt states by making use of the finite element method. The main conclusions of this paper are as follows: The halt positions of wind blade could have a big impact on tower circular flow and aerodynamic distribution, in which Condition 5 is the most unfavourable while Condition 1 is the most beneficial condition. The wind loads and load-effects of disturbed region of tower is obviously affected by different halt positions of wind blades, especially the large fluctuating displacement mean square deviation at both windward and leeward sides, among which the maximum response occurs in $350^{\circ}$ to the tower top under Condition 8; the maximum bending moment of tower bottom occurs in $330^{\circ}$ under Condition 2. The extreme displacement of blade top all exceeds 2.5 m under Condition 5, and the maximum value of windward displacement response for the tip of Blade 3 under Condition 8 could reach 3.35 m. All these results indicate that the influence of halt positions of different blades should be taken into consideration carefully when making wind-resistance design for large scale wind turbine tower.

• Wind and Structures
• /
• v.21 no.5
• /
• pp.461-487
• /
• 2015

In recent years, the Graphics Processing Unit (GPU) has become a competitive computing technology in comparison with the standard Central Processing Unit (CPU) technology due to reduced unit cost, energy and computing time. This paper describes the derivation and implementation of GPU-based algorithms for the analysis of wind loading uncertainty on high-rise systems, in line with the research field of probability-based wind engineering. The study begins by presenting an application of the GPU technology to basic linear algebra problems to demonstrate advantages and limitations. Subsequently, Monte-Carlo integration and synthetic generation of wind turbulence are examined. Finally, the GPU architecture is used for the dynamic analysis of three high-rise structural systems under uncertain wind loads. In the first example the fragility analysis of a single degree-of-freedom structure is illustrated. Since fragility analysis employs sampling-based Monte Carlo simulation, it is feasible to distribute the evaluation of different random parameters among different GPU threads and to compute the results in parallel. In the second case the fragility analysis is carried out on a continuum structure, i.e., a tall building, in which double integration is required to evaluate the generalized turbulent wind load and the dynamic response in the frequency domain. The third example examines the computation of the generalized coupled wind load and response on a tall building in both along-wind and cross-wind directions. It is concluded that the GPU can perform computational tasks on average 10 times faster than the CPU.

• Wind and Structures
• /
• v.10 no.2
• /
• pp.99-119
• /
• 2007

Estimations of wind flow over terrain are often needed for applications such as pollutant dispersion, transport safety or wind farm location. Whilst field studies offer very detailed information regarding the wind potential over a small region, the cost of instrumenting a natural fetch alone is prohibitive. Wind tunnels offer one alternative although wind tunnel simulations can suffer from scale effects and high costs as well. Computational Fluid Dynamics (CFD) offers a second alternative which is increasingly seen as a viable one by wind engineers. There are two issues associated with CFD however, that of accuracy of the predictions and set-up and simulation times. This paper aims to address the two issues by demonstrating, by way of an investigation of wind potential for the Askervein Hill, that a good level of accuracy can be obtained with CFD (10% for the speed up ratio) and that it is possible to automate the simulations in order to compute a full wind rose efficiently. The paper shows how a combination of script and session files can be written to drive and automate CFD simulations based on commercial software. It proposes a general methodology for the automation of CFD applied to the computation of wind flow over a region of interest.

• Journal of Environmental Science International
• /
• v.26 no.12
• /
• pp.1307-1319
• /
• 2017

In this study, we investigated the impact of different initial data on atmospheric modeling results using the Weather Research and Forecast (WRF) model. Four WRF simulations were conducted with different initialization in March 2015, which showed the highest monthly mean $PM_{10}$ concentration in the recent ten years (2006-2015). The results of WRF simulations using NCEP-FNL and ERA-Interim were compared with observed surface temperature and wind speed data, and the difference of grid nudging effect on WRF simulation between the two data were also analyzed. The FNL simulation showed better accuracy in the simulated temperature and wind speed than the Interim simulation, and the difference was clear in the coastal area. The grid nudging effect on the Interim simulation was larger than that of the FNL simulation. Despite of the higher spatial resolution of ERA-Interim data compared to NCEP-FNL data, the Interim simulation showed slightly worse accuracy than those of the FNL simulation. It was due to uncertainties associated with the Sea Surface Temperature (SST) field in the ERA-Interim data. The results from the Interim simulation with different SST data showed significantly improved accuracy than the standard Interim simulation. It means that the SST field in the ERA-Interim data need to be optimized for the better WRF simulation. In conclusion, although the WRF simulation with ERA-Interim data does not show reasonable accuracy compared to those with NCEP-FNL data, it would be able to be Improved by optimizing the SST variable.

• Bulletin of the Korean Space Science Society
• /
• 2008.10a
• /
• pp.25.2-25.2
• /
• 2008

Polar rain is a spatially uniform precipitation of electrons with energies around 100eV that penetrate into the polar cap region where geomagnetic field lines are connected to the Interplanetary Magnetic Fields (IMF). Since their occurrences depend on the IMF sector polarity, they are believed to originate from the field aligned component of the solar wind. However, statistically direct correlation between polar rain and solar wind has not been shown. In this presentation, we examined specifically the IMF strength influence on the polar rain flux variation by classifying of IMF sector polarities. For this study, we employed the polar rain flux data measured by STSAT-1 and compared them with the solar wind parameters obtained from the WIND and ACE satellites. We found the direct mutuality between polar rain flux and IMF strength with correlation coefficient above 0.5. This proportional tendency appears stronger when the northern hemisphere is in the away sector of the IMF, which could be associated with a favorable geometry for magnetic reconnection. Simple particle trajectory simulation clearly shows why polar rain intensity depends on the IMF sector polarity. These results are consistent with the direct entry model of Fairfield et al.(1985), while low correlation coefficient with solar wind density, the similarity between slops of both energy spectra shows that transport process occur without acceleration.

• KIEAE Journal
• /
• v.11 no.2
• /
• pp.75-82
• /
• 2011

In Korea, the ratio of population in urban areas used to be only 50.1% in 1970, but with the value risen to 90.8% in 2009, urbanization is going on rapidly. Urbanization, which occurs by the rampantly planted buildings, has become major source of raising building density, changing wind direction and reducing wind amount, and such reductions are affecting even inside the building. In each year, among the total energy consumption in Korea, residential portion takes up significant ratio, and specifically the ratio of apartment house is shown to be highest. In order to solve such problem, many studies are being conducted for the improvement of natural ventilation performance. The natural ventilation performance of apartment house are significantly determined by the characteristics of external and internal structure, but in macroscopic perspective, the performance is established fundamentally by the layout characteristics of the main building of the apartment house in preparation for wind conditions. So far researches on raising the thermal comfort through elevation of ventilation performance have been conducted actively, but many of them propose only theoretical concepts deduced through wind path analysis, and do not include any indicator to measure ventilation performance simply only with area data from layout planning stage. Therefore, in this study, gap ratio a wind field measuring indicator was developed, and after the ventilation characteristics by layout types and main building uniformity were identified, the scope of gap ratio efficient for ventilation and that of uniformity were clarified, followed by verification through simulation.

• Structural Engineering and Mechanics
• /
• v.70 no.4
• /
• pp.479-497
• /
• 2019

In the current code design, the use of a uniform internal pressure coefficient of cooling towers as internal suction cannot reflect the 3D characteristics of flow field inside the tower body with different ventilation rate of shutters. Moreover, extreme weather such as heavy rain also has a direct impact on aerodynamic force on the internal surface and changes the turbulence effect of pulsating wind. In this study, the world's tallest cooling tower under construction, which stands 210m, is taken as the research object. The algorithm for two-way coupling between wind and rain is adopted. Simulation of wind field and raindrops is performed iteratively using continuous phase and discrete phase models, respectively, under the general principles of computational fluid dynamics (CFD). Firstly, the rule of influence of 9 combinations of wind speed and rainfall intensity on the volume of wind-driven rain, additional action force of raindrops and equivalent internal pressure coefficient of the tower body is analyzed. The combination of wind velocity and rainfall intensity that is most unfavorable to the cooling tower in terms of distribution of internal pressure coefficient is identified. On this basis, the wind/rain loads, distribution of aerodynamic force and working mechanism of internal pressures of the cooling tower under the most unfavorable working condition are compared between the four ventilation rates of shutters (0%, 15%, 30% and 100%). The results show that the amount of raindrops captured by the internal surface of the tower decreases as the wind velocity increases, and increases along with the rainfall intensity and ventilation rate of the shutters. The maximum value of rain-induced pressure coefficient is 0.013. The research findings lay the basis for determining the precise values of internal surface loads of cooling tower under extreme weather conditions.

• Wind and Structures
• /
• v.28 no.2
• /
• pp.71-87
• /
• 2019

The yaw and interference effects of blades affect aerodynamic performance of large wind turbine system significantly, thus influencing wind-induced response and stability performance of the tower-blade system. In this study, the 5MW wind turbine which was developed by Nanjing University of Aeronautics and Astronautics (NUAA) was chosen as the research object. Large eddy simulation on flow field and aerodynamics of its wind turbine system with different yaw angles($0^{\circ}$, $5^{\circ}$, $10^{\circ}$, $20^{\circ}$, $30^{\circ}$ and $45^{\circ}$) under the most unfavorable blade position was carried out. Results were compared with codes and measurement results at home and abroad, which verified validity of large eddy simulation. On this basis, effects of yaw angle on average wind pressure, fluctuating wind pressure, lift coefficient, resistance coefficient,streaming and wake characteristics on different interference zone of tower of wind turbine were analyzed. Next, the blade-cabin-tower-foundation integrated coupling model of the large wind turbine was constructed based on finite element method. Dynamic characteristics, wind-induced response and stability performance of the wind turbine structural system under different yaw angle were analyzed systematically. Research results demonstrate that with the increase of yaw angle, the maximum negative pressure and extreme negative pressure of the significant interference zone of the tower present a V-shaped variation trend, whereas the layer resistance coefficient increases gradually. By contrast, the maximum negative pressure, extreme negative pressure and layer resistance coefficient of the non-interference zone remain basically same. Effects of streaming and wake weaken gradually. When the yaw angle increases to $45^{\circ}$, aerodynamic force of the tower is close with that when there's no blade yaw and interference. As the height of significant interference zone increases, layer resistance coefficient decreases firstly and then increases under different yaw angles. Maximum means and mean square error (MSE) of radial displacement under different yaw angles all occur at circumferential $0^{\circ}$ and $180^{\circ}$ of the tower. The maximum bending moment at tower bottom is at circumferential $20^{\circ}$. When the yaw angle is $0^{\circ}$, the maximum downwind displacement responses of different blades are higher than 2.7 m. With the increase of yaw angle, MSEs of radial displacement at tower top, downwind displacement of blades, internal force at blade roots all decrease gradually, while the critical wind speed decreases firstly and then increases and finally decreases. The comprehensive analysis shows that the worst aerodynamic performance and wind-induced response of the wind turbine system are achieved when the yaw angle is $0^{\circ}$, whereas the worst stability performance and ultimate bearing capacity are achieved when the yaw angle is $45^{\circ}$.