• International Journal of Automotive Technology
• /
• v.8 no.6
• /
• pp.687-696
• /
• 2007

To increase the reliability of auto-ignition in CAI engines, the thermodynamic properties of intake flow is often controlled using recycled exhaust gases, called internal EGR. Because of the internal EGR influence on the overall thermodynamic properties and mixing quality of the gases that affect the subsequent combustion behavior, optimizing the intake and exhaust valve timing for the EGR is important to achieve the reliable auto-ignition and high thermal efficiency. In the present study, fully 3D numerical simulations were carried out to predict the mixing characteristics and flow field inside the cylinder as a function of valve timing. The 3D unsteady Eulerian-Lagrangian two-phase model was used to account for the interaction between the intake air and remaining internal EGR during the under-lap operation while varying three major parameters: the intake valve(IV) and exhaust valve(EV) timings and intake valve lift(IVL). Computational results showed that the largest EVC retardation, as in A6, yielded the optimal mixing of both EGR and fuel. The IV timing had little effect on the mixing quality. However, the IV timing variation caused backflow from the cylinder to the intake port. With respect to reduction of heat loss due to backflow, the case in B6 was considered to present the optimal operating condition. With the variation of the intake valve lift, the A1 case yielded the minimum amount of backflow. The best mixing was delivered when the lift height was at a minimum of 2 mm.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1327-1339
• /
• 1999

Characteristics of two-phase flow and heat transfer were numerically investigated in a submerged gas Injection system. Effects of both the gas flow rate and bubble size were investigated. In addition, heat transfer characteristic and effects of heat transfer were investigated when temperature of the injected gas was different from that of the liquid. The Eulerian approach was used for the formulation of both the continuous and the dispersed phases. The turbulence in the liquid phase was modeled by the use of the standard $k-{\varepsilon}$ turbulence model. The interphase friction and heat transfer coefficient were calculated by means of correlations available in the literature. The turbulent dispersion of the phases was modeled by introducing a "dispersion Prandtl number". The plume region and the axial velocities are increased with increases in the gas flow rate and with decreases in the bubble diameter. The turbulent flow field grows stronger with the increases in the gas flow rate and with the decreases in the bubble diameter. In case that the heat transfer between the liquid and the gas is considered, the axial and the radial velocities are decreased in comparison with the case that there is no temperature difference between the liquid and the gas when the temperature of the injected gas is higher than the mean liquid temperature. The results in the present research are of interest in the design and the operation of a wide variety of material and chemical processes.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.6
• /
• pp.347-355
• /
• 2016

This study presents a Lattice Boltzmann Method (LBM) coupled with a momentum-exchange approach/fictitious domain (MEA/FD) method for the simulation of particle suspensions. The method combines the advantages of the LB and the FD methods by using two unrelated meshes, namely, a Eulerian mesh for the flow domain and a Lagrangian mesh for the solid domain. The rigid body conditions are enforced by the momentum-exchange scheme in which the desired value of velocity is imposed directly in the particle inner domain by introducing a pseudo body force to satisfy the constraint of rigid body motion, which is the key idea of a fictitious domain (FD) method. The LB-MEA/FD method has been validated by simulating two different cases, and the results have been compared with those through other methods. The numerical evidence illustrated the capability and robustness of the present method for simulating particle suspensions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.290-297
• /
• 2016

A two-phase Large Eddy Simulation(LES) has been conducted to investigate breakup and atomization of a liquid jet in a cross turbulent flow in a rectangular duct. Gas-droplet two-phase flow was solved by a coupled Eulerian-Lagrangian method which tracks every individual particles. Effects of liquid breakup models, sub-grid scale models, and a order of spatial discretization was investigated. The penetration depth in cross flow was comparable with experimental data by varying breakup model and LES scheme. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

• Journal of Radiation Protection and Research
• /
• v.48 no.1
• /
• pp.28-43
• /
• 2023

Background: High-fidelity meteorological data is a prerequisite for the realistic simulation of atmospheric dispersion of radioactive materials near nuclear power plants (NPPs). However, many meteorological models frequently overestimate near-surface wind speeds, failing to represent local meteorological conditions near NPPs. This study presents a new high-resolution (approximately 1 km) meteorological downscaling method for modeling short-range (< 100 km) atmospheric dispersion of accidental NPP plumes. Materials and Methods: Six considerations from literature reviews have been suggested for a new dynamic downscaling method. The dynamic downscaling method is developed based on the Weather Research and Forecasting (WRF) model version 3.6.1, applying high-resolution land-use and topography data. In addition, a new subgrid-scale topographic drag parameterization has been implemented for a realistic representation of the atmospheric surface-layer momentum transfer. Finally, a year-long simulation for the Kori and Wolsong NPPs, located in southeastern coastal areas, has been made for 2016 and evaluated against operational surface meteorological measurements and the NPPs' on-site weather stations. Results and Discussion: The new dynamic downscaling method can represent multiscale atmospheric motions from the synoptic to the boundary-layer scales and produce three-dimensional local meteorological fields near the NPPs with a 1.2 km grid resolution. Comparing the year-long simulation against the measurements showed a salient improvement in simulating near-surface wind fields by reducing the root mean square error of approximately 1 m/s. Furthermore, the improved wind field simulation led to a better agreement in the Eulerian estimate of the local atmospheric dispersion. The new subgrid-scale topographic drag parameterization was essential for improved performance, suggesting the importance of the subgrid-scale momentum interactions in the atmospheric surface layer. Conclusion: A new dynamic downscaling method has been developed to produce high-resolution local meteorological fields around the Kori and Wolsong NPPs, which can be used in short-range atmospheric dispersion modeling near the NPPs.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.149-149
• /
• 2020

하천에서 물 흐름이 보와 댐과 같은 수공구조물을 지날 때 일반적으로 흐름상태에 다양하고 급진적인 변화가 발생한다. 특히 흐름이 구조물을 지나면서 사류(supercritical flow)로 변하고 다시 상류(subcritical flow)로 복원되면서 일어나는 도수(hydraulic jump) 현상은 수위의 급변화, 흐름 에너지 소산, 변동성이 강한 압력 분포 등이 특징이다. 이러한 흐름 특성들은 보나 여수로와 같은 수공구조물 자체의 성능뿐만 아니라 이들 수공구조물의 하류에서 발생하는 국부세굴로 인해 구조물의 안정성에 부정적인 영향을 줄 수 있다. 따라서 수공구조물을 설계할 때는 이들 구조물을 통과하는 흐름의 비정상 난류 흐름 특성을 정확하게 해석하여 반영하여야 한다. 이 연구에서는 k-omega SST 난류 모형과 자유수면의 급격한 변동을 해석하기 위한 하이브리드-VOF(hybrid volume of fluid)기법을 이용하여 도수현상을 수치적으로 재현하고자 한다. 기존 CFD(computational fluid Dynamics) 모델링에서는 자유수면 변동의 영향을 고려하기 위해 VOF 기법을 많이 사용하였다. 하지면 전통적인 VOF 기법은 다상흐름(multiphase flow)을 오직 부피분율(volume fraction)의 함수로만 고려하며 모의함으로써 강한 수면변동뿐만 아니라 공기연행(air entrainment)를 동반하는 난류 흐름을 모의하는데는 한계가 있다. 이 연구에서 이용하는 Eulerian 기법인 하이브리드 VOF 기법은 물과 공기의 각 상에 대하여 흐름 특성들을 개별적으로 계산하기 때문에 공기연행을 포함한 급변류 흐름에서 전통적인 VOF 기법보다 적용성이 우수하다. 이와 같은 난류모형과 자유수면 포착기법을 이용하여 3차원 비정상 난류 흐름 수치모형을 구축하여 수공구조물 주변에서 발생하는 강한 공기연행과 난류 특성를 보이는 급변류를 수치적으로 재현한다. 이 연구는 계산된 수치해석 결과를 기존 수리실험 결과와 비교하여 수치모형의 적용성을 평가하고 도수 현상에서 발생하는 독특한 흐름 특성을 제시한다.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.48.4-48.4
• /
• 2019

We studied the metal-distribution of isolated Milky-way mass galaxy using various hydrodynamic solvers and investigated the difference of the result between AMR and SPH codes. In particle-based codes, physical quantities like mass or metallicity defined in each particle are conserved unless being injected explicitly by the effect of the supernova, whereas in the Eulerian codes the diffusion is simply accomplished by hydro-equation. Therefore, without including explicit physics of diffusion on the SPH- codes, the metal mixing in the galaxy or CGM only can be accomplished by the direct motion of the particles, however, the standard-SPH codes depress the instability of the turbulent fluid mixing. In this work, we simulated under common initial conditions, common gas-physics like cooling-heating models, and star-formation feedback using ENZO(AMR) GIZMO and GADGET-2 codes. We additionally included a metal-diffusion algorithm on the SPH-codes, which follows the subgrid-turbulent mixing model investigated by Shen et al. (2010) and compared the effect of the metal-outflow on the halo region of the galaxy in different hydro-solvers. We also found that for the implementation of the diffusion scheme in the SPH-codes, the existence of a sufficient number of the gas-particles, which is the carrier of the metals, is necessary. So we tested a new initial condition for proper implementation of the diffusion scheme on the SPH simulations. By comparing the metal-contamination of the circumgalactic medium with different hydrodynamics models, we quantify the diffusion strength of AMR codes using diffusion parameterization of the SPH codes and also suggest the calibration solutions in the different behavior of codes in metal-outflow.

• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.216-227
• /
• 2004

A new cavitating model by using bubble size distribution based on bubbles-mass has been proposed. Both liquid and vapor phases are treated with Eulerian framework as a mixture containing minute cavitating bubbles. In addition vapor phase consists of various sizes of vapor bubbles, which are distributed to classes based on their mass. The bubble number-density for each class was solved by considering the change of the bubble-mass due to phase change as well as generation of new bubbles due to heterogeneous nucleation. In this method, the bubble-mass is treated as an independent variable, and the other dependent variables are solved in spatial coordinates and bubble-mass coordinate. Firstly, we employed this method to calculate bubble nucleation and growth in stationary super-heated liquid nitrogen, and bubble collapse in stationary sub-cooled one. In the case of bubble growth in super-heated liquid, bubble number-density of the smallest class based on its mass is increased due to the nucleation. These new bubbles grow with time, and the bubbles shift to larger class. Therefore void fraction of each class is increased due to the growth in the whole class. On the other hand, in the case of bubble collapse in sub-cooled liquid, the existing bubbles are contracted, and then they shift to smaller class. It finally becomes extinct at the smallest one. Secondly, the present method is applied to a cavitating flow around NACA00l5 foil. Liquid nitrogen and liquid oxygen are employed as working fluids. Cavitation number, $\sigma$, is fixed at 0.15, inlet velocities are changed at 5, 10, 20 and 50m/s. Inlet temperatures are 90K in case of liquid nitrogen, and 90K and 1l0K in case of liquid oxygen. 110K of oxygen is corresponding to the 90K of nitrogen because of the same relative temperature to the critical one, $T_{r}$=$T/T_c^{+}$. Cavitating flow around the NACA0015 foils was properly analyzed by using bubble size distribution. Finally, the method is applied to a cavitating flow in an inducer of the LE-7A hydrogen turbo-pump. This inducer has 3 spiral foils. However, for simplicity, 2D calculation was carried out in an unrolled channel at 0.9R cross-section. The channel moves against the fluid at a peripheral velocity corresponding to the inducer revolutions. Total inlet pressure, $Pt_{in}$, is set at l00KPa, because cavitation is not generated at a design point, $Pt_{in}$=260KPa. The bubbles occur upstream of the foils and collapse between them. Cavitating flow in the inducer was successfully predicted by using the bubble size distribution.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.1-10
• /
• 2020

The application of multiphase flows is increasingly being applied to analyze phenomena such as single phase flows where the fluid boundary changes continuously over time or the problem of mixing a liquid phase and a gas phase. In particular, multiphase flow models that take into account incompressible Newtonian fluids for liquid and gas are often applied to solve the problems of the free water surface such as wave fields. In general, multi-phase flow models require time-based the surface tracking of each fluid's phase boundary, which determines the accuracy of the final calculation of the model. This study evaluates the advection performance of representative VOF-type boundary tracking techniques applied to various CFD numerical codes. The effectiveness of the FCT method to control the numerical flux to minimize the numerical diffusion in the conventional VOF-type boundary tracking method and advection calculation was mainly evaluated. In addition, the possibility of tracking performance of free surface using CIP method (Yabe and Aoki, 1991) was also investigated. Numerical results show that the FCT-VOF method introducing an anti-diffusive flux to precent excessive diffusion is superior to other methods under the confined conditions in this study. The results from this study are expected to be used as an important basic data in selecting free surface tracking techniques applied to various numerical codes.