• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.6
• /
• pp.429-445
• /
• 2008

Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of each phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe‘s approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.237-241
• /
• 2001

Our goal is to present a simple volume-of-fluid type interface-tracking algorithm to compressible two-phase flow in two space dimensions. The algorithm uses a uniform underlying Cartesian grid with some cells cut by the tracked interfaces into two subcells. A volume-moving procedure that consists of two basic steps: (1) the update of volume fractions in each grid cell at the end of the time step, and (2) the reconstruction of interfaces from discrete set of volume fractions, is employed to follow the dynamical behavior of the interface motion. As in the previous work with a surface-tracking procedure for general front tracking (LeVeque & Shyue 1995, 1996), a high resolution finite volume method is then applied on the resulting slightly nonuniform grid to update all the cell values, while the stability of the method is maintained by using a large time step wave propagation approach even in the presence of small cells and the use of a time step with respect to the uniform grid cells. A sample preliminary numerical result for an underwater explosion problem is shown to demonstrate the feasibility of the algorithm for practical problems.

• Journal of the Korean Society of Visualization
• /
• v.18 no.3
• /
• pp.35-41
• /
• 2020

We conduct numerical simulations of the interaction of a deformable structure with two-phase compressible flow. The finite volume method (FVM) is used to simulate fluid phenomena including a shock wave, a gas bubble, and the deformation of free surface. The deformation of a floating structure is computed with the finite element method (FEM). The compressible two-phase volume of fluid (VOF) method is used for the generation and development of a cavitation bubble, and the immersed boundary method (IBM) is used to impose the effect of the structure on the fluid domain. The result of the simulation shows the generation of a shock wave, and the expansion of the bubble. Also, the deformation of the structure due to the hydrodynamic loading by the explosion is identified.

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

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

• Ocean Systems Engineering
• /
• v.1 no.1
• /
• pp.31-57
• /
• 2011

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, a new numerical method was developed for accurate resolution of violent sloshing flow inside a three-dimensional LNG tank including wave breaking, jet formation, gas entrapping and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.13-17
• /
• 2006

RoeM and AUSMPW+ schemes are two of the most accurate and efficient schemes which are recently developed for the analysis of single phase gas dynamics. In this paper, we developed two-phase versions of these schemes for the analysis of gas-liquid large density ratio two-phase flow. We adopt homogeneous equilibrium model (HEM) using mass fraction to describe different two phases. In the Eulerian-Eulerian framework, HEM assumes dynamic and thermal equilibrium of the two phases in the same computational mesh. From the mixture equation of state (EOS), we derived new shock-discontinuity sensing term (SDST), which is commonly used in RoeM and AUSMPW+ for the stable numerical flux calculation. The proposed two-phase versions of RoeM and AUSMPW+ schemes are applied on several air-water two-phase test problems. In spite of the large discrepancy of material properties such as density, enthalpy, and speed of sound, the numerical results show that both schemes provide very satisfactory solutions.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.219-222
• /
• 2005

A compressible two-fluid two-phase flow computation model using the stiffened-gas equation of state is formulated. Since the conservation equation system is of mixed type, it gives complex eigenvalues. The sonic speeds obtained from the individual single phase have been simply used in the literature for the fastest wave speeds necessary in the HLL scheme. This method has worked fine but proved to be quite diffusive according to our test. To improve the accuracy, we here propose to utilize the analytic eigenvalues evaluated from an approximate Jacobian matrix lot the fastest wave speeds. The interfacial transfer terms were dropped in constituting the Jacobian matrix for this purpose. The present scheme proved efficient, robust and accurate in comparison with other existing methods. We solved the cavitating flow problem using the present scheme. The result shows more detailed wave structure in the cavitating process caused by the strong expansion waves.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1996.11a
• /
• pp.610-614
• /
• 1996

Sealing of lubricat-air mixture in the high performance machining conte is one of most the important characteristics to carry out enhanced lubrication. High speed spindle requires non-contact type of sealing mechanism. Evaluating an optimum seal design to minimize leakage is concerned in the aspect of flow control. Effect of geometry and leakage path are evaluated according to variation of sealing geometry, Velocity, pressure, turbulence intensity of profile is calculated to fina more efficient geometry and variables. This offers a methodological way of enhancement seal design for high speed spindle. The working fluid is regarded as two phases that are mixed flow of oil phase and air phase. It is more reasonable to simulate an oil jet or oil mist type high speed spindle lubrication. Turbulence and compressible flow model are used to evaluate a flow characteristic, This paper arranges a geometry of mostly used non-contact type seal and analyzes leakage characteristics to minimize a leakage on the same sealing area.