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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.2
• /
• pp.162-169
• /
• 2007

This study is focused on the comparison of a 3 dimensional numerical and hydraulic model experiment for the flow phenomenon when a lock gate is opened. The lock gate is designed to discharge the flood flow rate at $218m^3/s$ of Solicheon at the Kun Jang national industry complex. The three dimensional ${\kappa}-{\epsilon}$ turbulent model of ANSYS CFX-10 of the computational fluid dynamics(CFD) program was used. The characteristics of CFX-10 are able to be simulated effectively for turbulent flow, especially the flow separation of the boundary layer of the two phase interface of air and water. The velocity and the flow pattern of the numerical model was showed to be similar to the results of the hydraulic model experiment.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.313-316
• /
• 2011

This paper investigates CFD investigation on single phase supersonic nozzle flow and 2-phase subson ic flow prior to rocket nozzle supersonic 2-phase flow with water injection within the flame deflector. Numerical results of supersonic nozzle single phase flow showed no notable unrealistic behavior as it captures the usual shock cell structures. Three-dimensional 2-phase flow analysis has also been performed to verify whether the approach can grab the droplet behavior during cooling by water injection. It is expected these basic studies will enhance the cooling problem analysis of supersonic 2-phase rocket plume in the future.

• Wind and Structures
• /
• v.30 no.5
• /
• pp.451-463
• /
• 2020

Two-dimensional numerical simulations are conducted at a low Reynolds number Re = 100 to investigate the near wake of three identical circular cylinders that are arranged in an equilateral triangular configuration. The incident angle of the three-cylinder configuration with respect to incoming flow is varied from θ = 0° to 60°, while the spacing between adjacent cylinders (L) covers a wide range of L/D = 1.25-7.0, where D is diameter of the cylinder. Typical flow structures in the near wake of the three-cylinder configuration are identified, including a single Karman vortex street, bistable flip-flopping near wake, anti-phase and/or in-phase vortex shedding, shear layer reattachment, and vortex impingement, depending on the configuration (L/D, θ). The behavior of Strouhal number (St) is discussed in detail, echoing the distinct structures of near wake. Furthermore, fluid forces on the individual cylinders are examined, which, though highly depending on (L/D, θ), exhibit a close correlation to the near wake behavior.

• Nuclear Engineering and Technology
• /
• v.32 no.5
• /
• pp.433-445
• /
• 2000

The interfacial area concentration (IAC) is one of the most important parameters in the two-fluid model for two-phase flow analysis. The IAC can be measured by a local conductivity probe method that uses the difference of conductivity between water and air/steam. The number of sensors in the conductivity probe may be differently chosen by considering the flow regime of two-phase flow. The four sensor conductivity probe method predicts the IAC without any assumptions of the bubble shape. The local IAC can be obtained by measuring the three dimensional velocity vector elements at the measuring point, and the directional cosines of the sensors. The five sensor conductivity probe method proposed in this study is based on the four sensor probe method. With the five sensor probe, the local IAC for a given referred measuring area of the probe can be predicted more exactly than the four sensor probe. In this paper, the mathematical approach of the five sensor probe method for measuring the IAC is described, and a numerical simulation is carried out for ideal cap bubbles of which the sizes and locations are determined by a random number generator.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.72-83
• /
• 2022

This paper presents a numerical investigation of two-phase natural circulation flows established when external reactor vessel cooling is applied to a severe accident of the APR1400 reactor for the in-vessel retention of the core melt. The coolability limit due to external reactor vessel cooling is associated with the natural circulation flow rate around the lower head of the reactor vessel. For an elaborate prediction of the natural circulation flow rate using a thermal-hydraulic system code, MARS-KS1.5, a three-dimensional computational fluid dynamics (CFD) simulation is conducted to estimate the flow rate and pressure distribution of a liquid-state coolant at the brink of significant void generation. The CFD calculation results are used to determine the loss coefficient at major flow junctions, where substantial pressure losses are expected, in the nodalization scheme of the MARS-KS code such that the single-phase flow rate is the same as that predicted via CFD simulations. Subsequently, the MARS-KS analysis is performed for the two-phase natural circulation regime, and the transient behavior of the main thermal-hydraulic variables is investigated.

• Journal of Energy Engineering
• /
• v.22 no.1
• /
• pp.8-16
• /
• 2013

A core catcher is considered as a promising engineered system to stabilize the molten corium in the containment during a postulated severe accident in a nuclear power plant. Conceptually, the core catcher consists of a carbon steel body, sacrificial material, protection material, and engineered cooling channel. The cooling capacity of the engineered cooling channel should be guaranteed to remove the decay heat of the molten corium. The flow in ex-vessel core catcher is a combined problem of a two-phase flow in the engineered cooling channel and a single-phase natural circulation in the whole core catcher system. In this study, the analysis of the test facility for the core catcher using the CUPID code, which is a three-dimensional thermal-hydraulic code for the simulation of two-phase flows, was carried out to evaluate its cooling capacity.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.2
• /
• pp.362-373
• /
• 2002

A stereoscopic PIV (SPIV) measurement system based on the translation configuration was developed and applied to the flow behind a forward-swept axial-fan. Measurement of three orthogonal velocity components is essential for flow analysis of three-dimensional flows such as flow around a fan or propeller. In this study, the translation configuration was adopted to calculate the out-of-plane velocity component from 2-D PIV data obtained from two CCD cameras. The error caused by the out-of-plane motion was estimated by direct comparison of the 2-D PIV and 3-D SPIV results that measured from the particle images captured simultaneously. The comparison shows that the error ratio is relatively high in the region of higher out-of-plane motion near the axial fan blade. The turbulence intensity measured by the 2-D PIV method is bigger by about 5.8% in maximum compared with that of the 3-D SPIV method. The phase-averaged velocity field results show that the wake behind an axial fan has a periodic flow structure with respect to the blade phase and the characteristic flow structure is shifted downstream in the next phase.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.264-267
• /
• 2010

This research focuses on the development of numerical code to deal with compressible two phase flow around three dimensional objects combined with cavitation model suggested by Weishyy et al. with k-e turbulent model. The cryogenic cavitation is carried out by considering the thermodynamic effect on physical properties of cryogenic fluids in physical point of view and implementing the temperature sensitivity in the energy equation of the government equations in numerical point of view, respectively. The formulation has been extensively validated for both liquid nitrogen and liquid hydrogen by simulating the experiments of Hord on hydrofoils. Then, simulations of cavitating turbopump inducers at their design flow rate are presented. Results over a broad range of Nss numbers extending from single-phase flow conditions through the critical head break down point are discussed. In particular, thermal depression effects arising from cavitation in cryogenic fluids are identified and their impact on the suction performance of the inducer quantified.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.233-234
• /
• 2003

Fly ash enters axial compressor when a turbomachinery is operated in an adverse environment. We have numerically investigated erosion of the blade and shroud in the turbulent compressor passage flow under the influence of gas-particle two-phase interaction. There have appeared quasi-three dimensional calculations on this subject but not the complete three-dimensional gas-particle interaction as done in the present work. Lagrangian particle tracing technique is used on the base of parallel processing for efficient calculation. Accuracy of the present code is tested using the benchmark lPL nozzle. In the DFVLR compressor blades, we have shown that a large number of particles passing through the tip clearance make impact on the blade tip and on the shroud. Higher degree of erosion is resulted by the heavier particles due to the centrifugal force.