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

• 한국전산유체공학회지
• /
• 제12권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.

• 한국분무공학회지
• /
• 제27권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.

• 해양환경안전학회지
• /
• 제15권1호
• /
• pp.63-69
• /
• 2009

가스상의 체적분율과 압력강하는 기액이상류에 대한 이해와 예측에 있어서 매우 중요한 인자이다. 또한 그것들은 산업용 대용량의 열교환시스템 및 선박에 설치되는 보일러 및 냉동시스템의 설계에 있어서 필수적인 항목이다. 따라서 본 논문에서는 파이프의 모든 경사각도에서 기액이상류 가스상의 체적분율과 압력손실을 예측할 수 있는 이론적 해석 방법을 제시한다. 여기서의 이론적 해석은 2유체 층상류 모델을 기초로 하고 있다. 또한 이론적 해석결과와 기존의 실험결과와 비교한 결과에 대해서도 제시한다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제15권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• 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년도 추계 학술대회논문집
• /
• 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.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1996년도 추계학술발표회논문집(1)
• /
• 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
• /
• 제29권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).