• Wind and Structures
• /
• v.24 no.3
• /
• pp.205-221
• /
• 2017

The unsteady flow field disturbance between the blades and tower is one of the primary factors affecting the aerodynamic performance of wind turbine. Based on the research object of a 3MW horizontal axis wind turbine which was developed independently by Nanjing University of Aeronautics and Astronautics, numerical simulation on the aerodynamic performance of wind turbine system in halt state with blades in different position was conducted using large eddy simulation (LES) method. Based on the 3D unsteady numerical simulation results in a total of eight conditions (determined by the relative position with the tower during the complete rotation process of the blade), the characteristics of wind pressure distributions of the wind turbine system and action mechanism of surrounding flow field were analysed. The effect of different position of blades on the aerodynamic performance of wind turbine in halt state as well as the disturbance effect was evaluated. Results of the study showed that the halt position of blades had significant effect on the wind pressure distribution of the wind turbine system as well as the characteristics of flow around. Relevant conclusions from this study provided reference for the wind-resistant design of large scale wind turbine system in different halt states.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.842-845
• /
• 2005

A method for the numerical simulation of two-dimensional free-surface flow resulting from the propagation of regular gravity waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid finite-differences and spectral-method scheme. The formulation includes a boundary-fitted transformation, and is suitable for extension to incorporate large-eddy simulation (LES) and large-wave simulation (LWS) terms for turbulence and breaking wave modeling, respectively. Results are presented for the simulation of the free-surface flow over two different bottom topographies, with constant slope values of 1:10 and 1:20, two different inflow wave lengths and two different inflow wave heights. An absorption outflow zone is utilized and the results indicate minimum wave reflection from the outflow boundary. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion, while wave heights remain more or less unchanged. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3874-3897
• /
• 2023

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High-Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next-generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds-averaged Navier-Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

• Journal of the Korean Society for Railway
• /
• v.17 no.5
• /
• pp.328-335
• /
• 2014

Numerical prediction methods were applied to investigate the turbulent air currents and air-conditioning efficiency in an underground subway station, and the results compared to experimental data. The Shin-gumho Station($8^{th}$ floor underground and 43.6m in depth) in Seoul was selected for the analysis. The entire station was covered for simulation and the ventilation mode was ordinary. The ventilation diffusers were modeled as 95 square shapes of $0.6m{\times}0.6m$ in the lobby and as 222 square shapes in the platform. Cooling air of $47,316m^3/h$ was supplied and the returned air of $33,980m^3/h$ is exhausted in the lobby and the cooling air of $33,968m^3/h$ is supplied and the returned air of $76,190m^3/h$ was exhausted in the platform which is the same as the experimental data. The cases of the screen-door-closed and open were respectively investigated. A total of 7.5million grids were generated and the whole domain divided into 22 blocks for MPI efficiency of calculation. Large eddy simulation (LES) was applied to solve the momentum and energy equation.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.11 no.4
• /
• pp.133-140
• /
• 2008

The purpose of this study is to provide response methods to minimize the damage from chemical terrorism in a naturally ventilated indoor system using several types of dispersion simulations. Three chemical warfare agents such as sarin(GB), phosgene and chlorine gas which have high potential to be used in terror or to be involved with accidents were selected in this simulation. Fire dynamic simulation based on Large Eddy Simulation which is effective because of less computational effort and detailed expression of the dispersion flow was adopted to describe the dispersion behavior of these agents. When the vent speed is 0.005m/s, the heights of 0.1 agent mass fraction are 0.9m for sarin, 1.0m for phosgene and 1.1m for chlorine gas, and the maximum mass fraction are 0.27 for all three agents. However, when the vent speed is increased to 0.05m/s, the heights of 0.1 agent mass fraction become 1.6m for all three agents and maximum mass fraction inside the room increase to 0.70 for sarin, 0.58 for phosgene and 0.53 for chlorine gas. It is shown that molecular weight of the agents has an important role for dispersion, and it is important to install ventilation system with height less than 1.6m to minimize the damage from chemical toxicity.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1435-1442
• /
• 2014

This study presented numerical simulations of smooth-bed flows in the rectangular open-channel using the source code by OpenFOAM. For the analysis of the turbulent flow, Large Eddy Simulations were carried out and the dynamic sub-grid scale model proposed by Germano et al. (1991) is used to model the residual stress term. In order to analyze the coherent structure, the uw quadrant method proposed by Lu and Willmarth (1973) is used and the contribution rate and the fraction time of the instantaneous Reynolds stress are obtained in the Reynolds stress. The results by the present study are analyzed and compared with data from previous laboratory studies and direct numerical simulations. It is found that the contribution rate of the ejection events is larger than that of sweep events over the buffer layer in the open-channel flow over the smooth bed, however, the frequency of the sweep event is higher than that of the ejection events.

• Wind and Structures
• /
• v.21 no.3
• /
• pp.289-309
• /
• 2015

Vertical profiles of mean and fluctuating wind velocities over water waves were studied, by performing Large-Eddy Simulations (LES) on a fully developed turbulent boundary layer over simplified water waves. The water waves were simplified to two-dimensional, periodic and non-evolving. Different wave steepness defined by $a/{\lambda}$ (a : wave amplitude; ${\lambda}$ : wavelength) and wave age defined by $c/U_b$ (c: phase velocity of the wave; $U_b$ : bulk velocity of the air) were considered, in order to elaborate the characteristics of mean and fluctuating wind profiles. Results shows that, compared to a static wave, a moving wave plays a lesser aerodynamic role as roughness as it moves downstream slower or a little faster than air, and plays more aerodynamic roles when it moves downstream much faster than air or moves in the opposite direction to air. The changes of gradient height, power law index, roughness length and friction velocity with wave age and wave amplitude are presented, which shed light on the wind characteristics over real sea surfaces for wind engineering applications.

• The Journal of the Acoustical Society of Korea
• /
• v.43 no.1
• /
• pp.95-102
• /
• 2024

Due to technological advances, the cruising speed of high-speed trains is increasing, and aerodynamic noise generated from the flow outside the train has been an important consideration in the design stage. To accurately predict the flow-induced noise, high-resolution generation of sound sources in the near field and low-dissipation of sound propagation in the far field are required. This should be accompanied by a numerical grid and time resolution that can properly consider both temporal and spatial scales for each component of the real high-speed train. To overcome these challenges, this research simultaneously calculates the external flow and acoustic fields of five high-speed train cars of real-scale and at operational running speeds using a threedimensional unsteady Large Eddy Simulation technique. To verify the numerical analysis, the measurements of the wall pressure fluctuation and numerical results are compared. The Ffowcs Williams and Hawking equation is used to predict the acoustic power radiated from the high-speed train. This research is expected to contribute to noise reduction based on the analysis of the aerodynamic noise generation mechanism of high-speed trains.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Journal of Korea Water Resources Association
• /
• v.41 no.2
• /
• pp.197-212
• /
• 2008

Recently, stream flow analysis has been accomplished by one or two dimensional equations and was applied by simple momentum equations and fixed energy conservations which contain many condition limits. In this study, FLOW-3D using CFD (Computational Fluid Dynamics) was applied to stream flow analysis which can solve three dimensional RANS (Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behaviors and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as ${\kappa}-{\varepsilon}$, RNG (Renormalized Group) ${\kappa}-{\varepsilon}$ and LES (Large Eddy Simulation). Numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow, water level pressure and eddy flows around the piers and transverse overflow type structures. These results will be able to used by basis data that catch hold of effects on long-term bed elevation changes, sediment accumulations, scours and water aggravations by removal of obsolete transverse over flow type structures in urban stream.