• Nuclear Engineering and Technology
• /
• v.54 no.5
• /
• pp.1914-1928
• /
• 2022

To improve the heat transfer efficiency of the reactor fuel assembly, it is necessary to accurately calculate the two-phase flow boiling characteristics and the critical heat flux (CHF) in the fuel assembly. In this paper, a Eulerian two-fluid model combined with the extended wall boiling model was used to numerically simulate the 5 × 5 fuel rod bundle with spacer grids (four sets of mixing vane grids and four sets of simple support grids without mixing vanes). We calculated and analyzed 11 experimental conditions under different pressure, inlet temperature, and mass flux. After comparing the CHF and the location of departure from the nucleate boiling obtained by the numerical simulation with the experimental results, we confirmed the reliability of computational fluid dynamic analysis for the prediction of the CHF of the rod bundle and the boiling characteristics of the two-phase flow. Subsequently, we analyzed the influence of the spacer grid and mixing vanes on the void fraction, liquid temperature, and secondary flow distribution. The research in this article provides theoretical support for the design of fuel assemblies.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.121-128
• /
• 2009

During a hypothetical high-pressure accident in a nuclear power plant (NPP), molten corium can be ejected through a breach of a reactor pressure vessel (RPV) and dispersed by a following jet of a high-pressure steam in the RPV. The dispersed corium is fragmented into smaller droplets in a reactor cavity of the NPP by the steam jet and released into other compartments of the NPP by a overpressure in the cavity. The fragments of the corium transfer thermal energy to the ambient air in the containment or interact chemically with steam and generate hydrogen which may be burnt in the containment. The thermal loads from the ejected molten corium on the containment which is called direct containment heating (DCH) can threaten the integrity of the containment. DCH in a NPP containment is related to many physical phenomena such as multi-phase hydrodynamics, thermodynamics and chemical process. In the evaluation of the DCH load, the melt dispersion rates depending on the RPV pressure are the most important parameter. Mostly, DCH was evaluated by using lumped-analysis codes with some correlations obtained from experiments for the dispersion rates. In this study, MC3D code was used to evaluate the dispersion rates in the APR1400 NPP during the high-pressure accidents. MC3D is a two-phase analysis code based on Eulerian four-fields for melt jet, melt droplets, gas and water. The dispersion rates of the corium melt depending on the RPV pressure were obtained from the MC3D analyses and the values specific to the APR1400 cavity geometry were compared to a currently available correlation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.05a
• /
• pp.99-102
• /
• 2009

LES(Large eddy simulation) of breakup and atomization of a liquid jet into cross turbulent flow was performed. Two phase flow between a gas phase and a liquid phase was modeled by a mixed numerical scheme of both Eulerian and Lagrangian methods for gas and liquid phases respectively. The first and second breakup of liquid column was observed. The penetration depth in cross flow was comparable with experimental data for several variant of a liquid-gas momentum flux ratio by varying liquid injection velocities. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• Journal of Mechanical Science and Technology
• /
• v.14 no.5
• /
• pp.553-561
• /
• 2000

This paper presents a numerical model of multi phase flow of the mixtures of molten material-liquid-vapor, particularly in thermal nonequilibrium. It is a two-dimensional, transient, three-fluid model in Eulerian coordinates. The equations are solved numerically using the finite difference method that implicitly couples the rates of phase changes, momentum, and energy exchange to determine the pressure, density, and velocity fields. To examine the model's ability to predict an experimental data, calculations have been performed for tests of pouring hot particles and molten material into a water pool. The predictions show good agreement with the experimental data. It appears, however, that the interfacial heat transfer and breakup of molten material need improved models that can be applied to such high temperature, high pressure, multi phase flow conditions.

• Clean Technology
• /
• v.25 no.1
• /
• pp.19-32
• /
• 2019

As regional air pollution gets worse by the sulfur oxides emitted from various types of vessels passing through the many countries, the International Maritime Organization establishes the emission control areas and regulates sulfur dioxide in those areas. In order to satisfy these regional regulations, the fuel selection method and the exhaust gas post-treatment device are applied to the ships. Due to the economic reasons, the post-treatment method of exhaust gas for reducing the amount of sulfur oxides discharged is mainly preferred. The scrubber which is dominantly used in the ships are the spray type system where the sprayed liquid drops used for capturing the soluble sulfur dioxides in the exhaust gas. The performance of the spray type system depends on the size distribution of the sprayed droplets. In order to evaluate this performance, we designed counterflow type scrubber and cyclone scrubber and evaluated the desulfurization efficiency and the amount of droplet evaporation according to the size of each droplet by using computational fluid dynamics. The Eulerian-Eulerian analysis method was used because the scrubber had a gas-liquid two-phase flow inside the scrubber. When the diameter of the droplet was $100{\mu}m$, $300{\mu}m$, $500{\mu}m$ and $700{\mu}m$. As a result, both of scrubbers showed high desulfurization efficiency and low evaporation amount at $500{\mu}m$ and $700{\mu}m$.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.171-175
• /
• 2008

In the present study, a numerical simulations are investigated in a stirred solid/liquid system by using Eulerian multi-phase model. The transient flow field of liquid phase and distribution of solid particles are predicted in stirred tanks consisting of 4-pitched paddles impeller and baffles. The effects of number and width of baffles on the mixing time and the quality of solid suspension in a stirred tank are presented numerically. The result shows that the mixing time decreases as the width and number of baffles increase.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.330-340
• /
• 2005

Compressible two-phase flow is analyzed based on the arbitrary Lagrangian-Eulerian (ALE) formulation. For water, Tamman type stiffened equation of state is used. Numerical fluxes are calculated using the ALE two-phase Godunov scheme which assumes only that the speed of sound and pressure can be provided whenever density and internal energy are given. Effects of the approximations of a material interface speed are Investigated h method Is suggested to assign a rigid body boundary condition effectively To validate the developed code, several well-known problems are calculated and the results are compared with analytic or other numerical solutions including a single material Sod shock tube problem and a gas/water shock tube problem The code is applied to analyze the refraction and transmission of shock waves which are impacting on a water-gas interface from gas or water medium.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.29-35
• /
• 2006

A code is developed to analyze a spherically symmetric underwater explosion. The arbitrary Lagrangian-Eulerian(ALE) Godunov scheme for two-phase flow is used to calculate numerical fluxes through moving control surfaces. For detonation gas of TNT and liquid water, the Jones-Wilkins-Lee(JWL) equation of states and the isentropic Tait relation are used respectively. It is suggested to use the Godunov variable to estimate the velocity of a material interface. The code is validated through comparisons with other results on the gas-water shock tube problem. It is shown that the code can handle generation of discontinuity and recovering of continuity in the normal velocity near the material interface during shock waves interact with the material interface. The developed code is applied to analyze a spherically symmetric underwater explosion. Repeated transmissions of shock waves are clearly captured. The calculated period and maximum radius of detonation gas bubble show good agreements with experimental and other numerical results.

• Nuclear Engineering and Technology
• /
• v.40 no.6
• /
• pp.477-488
• /
• 2008

The objective of the paper is to analyze the thermally induced density wave oscillations in water cooled boiling water reactors. A transient thermal hydraulic model is developed with a characteristics-based implicit finite-difference scheme to solve the nonlinear mass, momentum and energy conservation equations in a time-domain. A two-phase flow was simulated with a one-dimensional homogeneous equilibrium model. The model treats the boundary conditions naturally and takes into account the compressibility effect of the two-phase flow. The axial variation of the heat flux profile can also be handled with the model. Unlike the method of characteristics analysis, the present numerical model is computationally inexpensive in terms of time and works in a Eulerian coordinate system without the loss of accuracy. The model was validated against available benchmarks. The model was extended for the purpose of studying the flow-induced density wave oscillations in forced circulation and natural circulation boiling water reactors. Various parametric studies were undertaken to evaluate the model's performance under different operating conditions. Marginal stability boundaries were drawn for type-I and type-II instabilities in a dimensionless parameter space. The significance of adiabatic riser sections in different boiling reactors was analyzed in detail. The effect of the axial heat flux profile was also investigated for different boiling reactors.