• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.204-212
• /
• 2003

As the first step to investigate the nonlinear interactions between turbulence and marine structures inside a viscous NWT, a LES technique was applied to the turbulent channel flow for $Re_{T}=150$, in this paper. The employed models were 4 types, such as the Smagorinsky model, the Dynamic SGS model, the Structure Function model and the Generalized Normal Stress model. The simulated data in time-series for the LESs were averaged in both time and space and performed statistical analysis. And results of the LESs were compared with those of a DNS developed in the present study and two spectral methods by Yoon et al.(2003) & Kim et al.(1987). It seems to be quite difficult to evaluate their performances to the present problem, but is seen that the accuracy of LESs are still related to the number of grids(or fine grid size).

• Journal of the Korean Society of Propulsion Engineers
• /
• v.18 no.2
• /
• pp.1-9
• /
• 2014

To investigate the flame dynamics in an annular combustor with single swirl injector, a 3D large-eddy simulation (LES) and a level-set flamelet turbulent combustion model have been implemented. The LM6000 developed by GEAE has been used as the combustor of concern and boundary conditions are based on experimental data. The strong central toroidal recirculation zone induced by the volume expansion of the combustion gas and the vortex breakdown continuously occurred through the procession of the vortex with decreasing strength, are observed.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.1
• /
• pp.60-67
• /
• 2010

Recently, there has been an increase in the construction of various types of coastal structures for efficient wave dissipation, seawater exchange, and so on. Among these, a perforated slit caisson has been widely used to reduce the reflected wave energy and the wave pressure on the structure. Therefore, many studies on the wave force on a caisson, as well as the wave reflection from it, have been carried out with laboratory experiments and numerical analyses, considering it as a 2-D problem. However, because a structure like a perforated slit caisson has a variable 3-D shape, waves forces should be considered as a 3-D problem. Therefore, in this paper, a fully-nonlinear 3-D numerical model (LES-WASS-3D) is proposed to examine the reflection characteristics of a perforated slit caisson with two chambers. The numerical model, LES-WASS-3D, was verified in a 3-D wave field by a comparison with existing experimental data for wave reflection coefficients. Then, using the numerical results, the reflection from a perforated slit caisson with two chambers was examined in relation to wave steepness, chamber width, and the shape/porosity of perforated slit.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.3
• /
• pp.181-190
• /
• 2010

In case of constructing submerged breakwaters, the circulation current is occurred around the open inlet because of mean water level difference between front and rear sides of them. The aim of this study is to investigate the flow control structure of new-type submerged breakwater which is able to reduce mean water level at rear side of it. At first, the numerical model (LES-WASS-3D) is validated by comparing with existing experimental data. And then, numerical simulation is carried out to examine wave height, mean water level and mean flow around the newtype submerged breakwater. From the numerical results, it can be pointed out that the new-type submerged breakwater with drainage system reduces the rip current around the open inlet.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.5B
• /
• pp.483-489
• /
• 2011

In this study, the coupled-numerical model has been newly developed to investigate numerically scouring and deposition around a coastal structure like a submerged breakwater using a numerical wave model and a lagrangian particle model for sand transport. As a numerical wave model, LES-WASS-2D (Hur and Choi, 2008) is adopted. The model is able to consider the flow through a porous midium with inertial, laminar and turbulent resistance term and determine the eddy viscosity with LES turbulence model. Distinct element method (Cundall and Strack, 1979), which is able to apply to many dynamical analysis of particulate media, as a lagrangian particle model for sand transport is newly coupled to the numerical wave model. The numerical simulation has been carried out to examine the scour problem in front of an impermeable submerged breakwater using the newly coupled-numerical model. The numerical results has been compared qualitatively with an existing experimental data and then its applicability has been discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.3
• /
• pp.39-47
• /
• 1985

In this paper, a FEM analysis model was developed to solve the consolidation phenomena of embankment on the soft foundation. The developed FEM model was based on the Biot's consolidation equation which was coupled with one of three stress-strain constitutive relationships. In order to check the validity of the newly developed FEM model, the program input data were used by a test embankment which had been already constructed at Cubzac-les-ponts in France by Magnan et al. The FEM results compared to the experimental and analytical results which were obtained by the Magnan's group at Cubzac-les-ponts. The results compared showed that the consolidation phenomena were well explained by the author's FEM model which results were more accurate than the others. As for the pore water pressure, Christian-Boehmer's method used in this paper was considered preferable to Sandhu-Wilson's used by Magnan.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2011.02a
• /
• pp.333-333
• /
• 2011

We studied the H/D exchange kinetics of pure and acid dopped water-ice film by using the techniques of reactive ions scattering (RIS) and low energy sputtering (LES) with low kinetic energy cesium ion beam (<35 eV). From RIS, neutral water isotopomers were detected in the form of cesium-molecule ion clusters, $CsX^+$ (X= $H_2O$, HDO, $D_2O$). Ionic species, like $H_3O^+$, $DH_2O^+$, $D_2HO^+$, $D_3O^+$, adsorbed on the surface were ejected via LES process. Those techniques allowed us to trace the isotopomeric populations of water-ice film. To show the catalytic effect of excess proton in the H/D exchange reaction, our study was conducted with two types of water-ice films. In film 1, about 0.5 BL of $H_2O$ was adsorbed on HCl (0.1 ML) dopped $D_2O$ (8 BL) film. In film 2, similar amount of $H_2O$ used in film 1 was adsorbed on pure $D_2O$ film. Kinetic data were obtained from each film type for 90-110 K (film 1) and 110-130 K (film 2) and fitted with numerically integrated lines. Through the Arrhenius plot of kinetic coefficient deduced from fitting of the H/D exchange reaction, the activation energy of film 1 and 2 were estimated to be $10{\pm}3kJmol^{-1}$ and $17{\pm}4kJmol^{-1}$. This activation barrier difference could be understood from detailed pictures of H/D exchange. In film 2, both the formation of ion pair, $H_3O^+$ and OH. and proton transfer were needed for the H/D exchange. However, in film 1, only proton transfer was necessary but ion pair formation was not, so this might reduce the activation energy.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.449-452
• /
• 2008

Flow difference between reacting and non-reacting case in a swirl stabilized annular combustor is investigated using 3D Large Eddy Simulation with flamelet turbulent combustion model. The combustor of concern is the LM6000, lean premixed dry low-NOx annular combustor, developed by GEAE. Boundary conditions are based on experimental data. Heat release as a result of combustion put the dilatation of density in primary combustion zone highly increased so that the main swirl stream behind of a swirl cup stretched further downstream than that of non-reacting case. The oval shape of core flow in cross-section to flow direction, which clearly observed in non-reacting case, tends to be circle, and small vorticities in wide range in non-reacting case disappears, but the size of iso-vorticity increase in reacting case.

• 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.

• Proceedings of the KSR Conference
• /
• 2011.10a
• /
• pp.2885-2890
• /
• 2011

Simulation study were performed for fire location effect on the smoke spread in the deeply-underground subway station(DUSS). In this research, Shingumho station (The line # 5, Depth: 46m) has been selected as case-study for the analysis of smoke-spread effect with the different fire location. Field test data measured for actual fan in DUSS was applied as a condition of a simulation. The whole station was covered in this analysis and 4 million grids were generated for this simulation. The fire driven flow was analyzed case by case to compare the smoke-spread effect according to the fire location. In order to enhance the efficiency of calculation, parallel processing by MPI was employed and LES(large eddy simulation) method in FDS code was adopted.