• 한국전산유체공학회:학술대회논문집
• /
• 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 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.

• Journal of ILASS-Korea
• /
• v.27 no.3
• /
• pp.144-154
• /
• 2022

Gas injection is a technique applied to improve throttling in liquid rocket engines and atomization in effervescent injectors. When a gas is injected into a liquid, it creates a two-phase flow inside the injector. The changes (bubbly flow, slug flow, annular flow, etc.) in the two-phase flow affect the injector's spray characteristics. In this study, cold-flow tests were performed by using three injectors with different orifice diameters and four aerators with different gas injection hole diameters. The experiments were done by changing the thrust ratio (liquid mass flow rate ratio) and gas-liquid mass flow rate ratio. Two-phase flow transition, breakup length, and discharge coefficient according to the injector/aerator design and flow conditions were investigated in detail.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.15 no.1
• /
• pp.63-69
• /
• 2009

Gas volumetric fractions and pressure loss are very important parameters in understanding and predicting gas-liquid two-phase flows. They are also essential to design large heat exchanging system in many industries, boiler and refrigerating systems mounted at ships. This paper therefore presents a theoretical method of predicting the pressure loss and gas volumetric fractions in gas-liquid two-phase flows for the whole range of pipe inclinations. The theoretical analysis is based on the two-fluid stratified flow model. It also provides the results of the comparisons between this theoretical analysis results and previous experimental results.

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

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.73-80
• /
• 2010

In this paper, we present the exact Riemann solver for the compressible liquid-gas two-phase shock tube problems. We hereby consider both isentropic and non-isentropic two-phase flows. The shock tube has a diaphragm in the mid-section which separates the liquid medium on the left and the gas medium on the right. By rupturing the diaphragm, various waves are observed on the phasic field variables such as pressure, density, temperature and void fraction in the form of rarefaction wave, shock wave and material interface (contact discontinuity). Both phases are treated as compressible fluids using the linearized equation of state or the stiffened-gas equation of state. We solve several shock tube problems made of a high/low pressure in the liquid and a low/high pressure in the gas. The wave propagations are well resolved by the exact Riemann solutions.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1903-1908
• /
• 2004

A preconditioned numerical method for gas-liquid two-phase flows is applied to solve cavitating flow. The present method employs a finite-difference dual time-stepping integration procedure and the MUSCLTVD scheme. A homogeneous equilibrium cavitation model is used. The present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation, large density changes and incompressible flow characteristics at low Mach number. Some internal flows such as convergent-divergent nozzles are computed using this method. Comparisons of predicted and experimental results are provided and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.169-177
• /
• 2005

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.11a
• /
• pp.301-306
• /
• 1996

To predict the average liquid level under the condition of the countercurrent stratified two-phase flow in a pipe, an analytical model has been suggested. This is made by introducing the interfacial level gradient into the liquid-phase and the gas-phase momentum equations. The analytical method for the gas-phase pressure drop calculation with f$_i$ $\neq$ f$_G$ has also been described using the liquid level prediction model developed in the present study.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.4
• /
• pp.377-383
• /
• 2005

This paper deals with the flow characteristics of air-water two-phase flow in a vertical tube of 10mm I.D. and 600mm in length at an adiabatic condition. The obtained experimental data were covered with the liquid superficial velocity ranging from 0.095m/s to 2.56m/s. and the gas superficial velocity ranging from 0.032m/s to 21.08m/s. The effects of the gas and liquid superficial velocity on the flow pattern transitions, frictional pressure drop, and film thickness and gas-liquid interface roughness were also examined. It was found that the film thickness increased and the liquid film wave length was more longer with the liquid superficial velocity $j_L$ increasing at $j_G$ constant. It was also showed that the frictional pressure drops were experienced in three regions. namely increasing region(bubbly flow), decreasing region (Taylor bubble and slug flows) and re-increasing region (annular flow).