• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.109-114
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.187-195
• /
• 2012

This paper deals with numerical computations of cryogenic flows around turbopump inducer. Firstly, we introduce numerical methods to compute compressible/incompressible cryogenic two-phase flow. As a validation problem, computation results of 2 dimensional/axi-symmetric cryogenic flow will be presented. In this process, various cavitation model will be compared. Finally, numerical simulation of 3 dimensional turbopump inducer will be presented.

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

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Economic and Environmental Geology
• /
• v.41 no.6
• /
• pp.731-746
• /
• 2008

Most multiphase flow simulators following fractional flow approach assume incompressibility of fluid and matrix or consider only two phase flow (water and air, water and NAPL). However, in this study, mathematical governing equations were developed for fully compressible three-phase flow using fractional flow based approach. Also, fully compressible multiphase flow simulator (CMPS) considering compressibilities of matrix and fluid was developed using the mathematical governing equations. In order to verify CMPS, the CMPS were compared with analytical solution and the existing multiphase flow simulator, MPS, which had been developed for simulating incompressible multiphase flow (Suk and Yeh 2007; Suk and Yeh 2008). According to the results, solutions of CMPS and MPS and analytical solutions are well matched each other. Thus, it is found that CMPS has the capability of simulating compressible three phase flow phenomena assuming compressibilities of fluids and matrix.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.11a
• /
• pp.221-224
• /
• 2009

Since the coupling of cavitation modeling with turbulent flow is the difficulty topic, a numerical simulation for two phase flow remains as one of the challenging issues in the society. This research focuses on the development of numerical code to deal with incompressible two phase flow around conical body combined with cavitation model suggested by Kunz et al. with k-e turbulent model. The simulation results are compared to experimental data to verify the validity of the developed code. The calculation results show very good agreement with experimental observations. Also, the calculation of cavitation in cryogenic fluid is being done by implementing the temperature sensitivity in government equations and it is still in the progress. This code have been being further extended to 3D compressible two phase flow for the study on the fluid dynamics around inducers and impellers in turbo pump system.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.9-16
• /
• 2013

Cavitation on an axi-symmetric hemispherical head-form body was studied using an Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. To consider compressibility effects on the vapor phase and cavity interface, a pressure-based compressible flow CFD code was developed. To validate the developed CFD code, cavitating flow around the hemispherical head-form body was simulated using pressure-based incompressible and compressible CFD codes and validated against existing experimental data in the three-way comparison. The cavity shedding behavior, length of re-entrant jet, drag history, and Strouhal number of the hemispherical head-form body were compared between two CFD codes. The results, in this paper, suggested that the computations of cavitating flow with compressibility effects improve the description of cavity dynamics.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.653-662
• /
• 2007

Due to the difficulty raised from the coupling of cavitation modeling with turbulent flow, numerical simulation for two phase flow remains one of the challenging issues in the society. This research focuses on the development of numerical code to deal with incompressible two phase flow around 2D hydrofoil by combing the cavitation model suggested by Kunz et al. with $k-{\varepsilon}$ turbulent model. The simulation results are compared to experimental data to verify the validity of the developed code. Also, the comparison of the calculation results is made with LES results to evaluate the capability of $k-{\varepsilon}$ turbulence model. The calculation results show very good agreement with experimental observations even though this code can not grasp the small scaled bubbles in the calculation wheres LES can hold the real physics. This code will be extended to 3D compressible two phase flow for the study on the fluid dynamics in the inducers and impellers.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1821-1826
• /
• 2003

In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.1
• /
• pp.56-65
• /
• 2018

Pressure oscillation caused by the compressibility of entrapped air in dam break flow is analyzed using an open source code, which is a two-phase compressible code for non-isothermal immiscible fluids. Since compressible flows are computed based on a pressure-based method, the code can handle the equation of state of barotropic fluid, which is virtually incompressible. The computed time variation of pressure is compared with other experimental and computational results. The present result shows good agreements with other results until the air is entrapped. As the entrapped air bubbles pulsate, pressure oscillations are predicted and the pressure oscillations damp out quickly. Although the compressibility parameter of water has been varied for a wide range, it has no effects on the computed results, because the present equation of state for water is so close to that of incompressible fluid. Grid independency test for computed time variation of pressure shows that all results predict similar period of pressure oscillation and quick damping out of the oscillation, even though the amplitude of pressure oscillation is sensitive to the velocity field at the moment of the entrapping. It is observed that as pressure inside the entrapped air changes quickly, the pressure field in the neighboring water adjusts instantly, because the sound of speed is much higher in water. It is confirmed that the period of pressure oscillation is dominated by the added mass of neighboring water. It is found that the temperature oscillation of the entrapped air is critical to the quick damping out of the oscillations, due to the fact that the time averaged temperature inside the entrapped air is higher than that of surrounding water, which is almost constant.