• Journal of Advanced Marine Engineering and Technology
• /
• v.37 no.7
• /
• pp.729-736
• /
• 2013

The turbulent properties in a mixed statistically stationary flow were investigated experimentally by a pseudo stereoscopic PIV. In order to validate the experimental results, the profiles of the turbulent kinetic energy were evaluated with the flow features. A mechanical agitator having 6 blades was installed at the bottom of the mixing tank (D=60cm, H=60cm). The agitator was rotated with 80rpm clockwise and counter-clockwise. For the measurements, three cameras were used and all were synchronized. The images captured by one of the three cameras was used for the measurement of rotational speed, and the images captured by the other two cameras were used to measure three dimensional components of velocity vectors. All vectors captured at the same rotational angle were phase averaged to construct three-dimensional vector fields to reconstruct the spatial distribution of the flow properties. It was seen that the jet scrolling along the tank was the main source of mixing.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.2141-2145
• /
• 2004

As part of study on thermal hydraulic behavior in the reactor cavity under external vessel cooling in the APR (Advanced Power Reactor) 1400, one dimensional two phase flow of steady state in the reactor cavity have been analyzed to investigate a coolant circulation mass flow rate in the annulus region between the reactor vessel and the insulation material using the RELAP5/MOD3 computer code. The RELAP5/MOD3 results have shown that a two phase natural circulation flow of 300 - 600 kg/s is generated in the annulus region between the reactor vessel and the insulation material when the external vessel cooling has been applied in the APR 1400. An increase in the heat flux of the inner vessel leads to an increase of the coolant mass flow rate. An increase in the coolant outlet area leads to an increase in the coolant circulation mass flow rate, but the coolant inlet area does not effective on the coolant circulation mass flow rate. The change of the lower coolant outlet to a lower position affects the coolant circulation mass flow rate, but the variation trend is not consistent.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.126-134
• /
• 2008

We developed a three dimensional Navier-Stokes code based on the level set method to simulate two phase flows with high density ratio. The Navier-Stokes equations with consideration of the surface tension effects are solved by using SIMPLE algorithm on a non-staggered grid. The present code is validated by simulating two test problems. First one is to simulate a rising bubble inside a cube. The thickness of the interface of the bubble is shown to affect the pressure distribution around the interface. As the thickness decreases, the pressure field around the interface becomes more oscillatory. As the bubble rises, a ring vortex is shown to form around the interface and the bubble eventually develops into an ellipsoidal shape. Merge of two bubbles inside a container is secondly tested to show the robustness of the present code for two phase flow simulation. Numerical results show stable and reliable behavior during the process of merging of two bubbles. The velocity and pressure fields around the interface of bubbles are shown oscillation free during the merging of two bubbles.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Journal of KSNVE
• /
• v.6 no.2
• /
• pp.163-177
• /
• 1996

An investigation on thermohydraulic stability of flow oscillations in the CANada Deuterium Uranium-600(CANDU-6) heat transport system has been conducted. Flow oscillations in reactor coolant loops, comprising two heat sources and two heat sinks in series, are possibly caused by the response of the pressure to extraction of fluid in two-phase region. This response consists of two contributions, one arising from mass and another from enthalpy change in the two-phase region. The system computer code used in the investigation os SOPHT, which is capable of simulating steady states as well as transients with varying boundary conditions. The model was derived by linearizing and solving one-dimensional, homogeneous single- and two-phase flow conservation equations. The mass, energy and momentum equations with boundary conditions are set up throughout the system in matrix form based on a node-link structure. Loop stability was studied under full power conditions with interconnecting the two compressible two phase regions in the figure-of-eight circuit. The dominant function of the interconnecting pipe is the transfer of mass between the two-phase regions. Parametric survey of loop stability characteristics, i. e., damping ratio and period, has been made as a function of geometrical parameters of the interconnection line such as diameter, length, height and orifice flow coefficient. The stability characteristics with interconnection line has been clarified to provide a simple criterion to be used as a guide in scaling of the pipe.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.6
• /
• pp.38-46
• /
• 2011

In order to simulate the ignition-gas injection in the interior ballistics, a two-dimensional analytic code for two-phase flows has been developed. The Eulerian-Lagrangian approach and the low-dissipation simple high-resolution upwind scheme(LSHUS) have been adopted in the numerical code for the propellant combustion of the gun propelling charges. The ghost-cell extrapolation method has been used for the moving boundary in the chamber with the projectile movement. The calculation results of the developed code have been compared and verified through those of the dimensionless IBHVG2 code and the previous one-dimensional code. In comparison with the two-phase flows according to the drag models, the numerical analysis of the muzzle velocity has been affected by the drag model.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2626-2631
• /
• 2007

Air-water two phase natural circulation flow in the T-HERMES (Thermo-Hydraulic Evaluation of Reactor cooling Mechanism by External Self-induced flow)-1D experiment has been evaluated to verify and evaluate the experimental results by using the RELAP5/MOD3 computer code. The RELAP5 results have shown that an increase in the coolant inlet area leads to an increase in the water circulation mass flow rate. However, the water outlet area does not effective on the water circulation mass flow rate. As the coolant outlet moves to a lower position, the water circulation mass flow rate decreases. The water level is not effective on the water circulation mass flow rate. As the height increases in the air injection part, the void fraction increases. However, the void fraction in the upper part of the air injector maintains a constant value. An increase in the air injection mass flow rate leads to an increase in the local void fraction, but it is not effective on the local pressure.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2696-2701
• /
• 2007

For the analysis of a two-phase flow, the interaction between two phases such as the interfacial momentum or heat transfer is proportional to the interfacial area. So the interfacial area concentration (IAC) is one of the most important parameters governing the behavior of each phase. This study focuses on the development of a computational fluid dynamics (CFD) code for investigating a boiling flow with a one-group IAC transport equation. It was based on the two-fluid model and governing equations were calculated by SMAC algorithm. For checking the robustness of the developed code, the experiment of a subcooled boiling in a vertical annulus channel was analyzed to validate the capability of the IAC transport equation. As the results, the developed code was confirmed to have the capability in predicting multi-dimensional phenomena of vapor generation and propagation in a subcooled boiling.