• 대한기계학회논문집B
• /
• 제32권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.

• 천문학회지
• /
• 제34권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.

• 한국가시화정보학회지
• /
• 제18권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.

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

• 자원환경지질
• /
• 제41권6호
• /
• pp.731-746
• /
• 2008

대부분의 분을 흐름 접근 방식을 따르는 다상흐름모델들은 주로 유체나 매질의 비 압축성을 가정하거나 완전한 삼상 오염물을 다루지 않고 이상 유체 (물 또는 공기, 물 또는 불용융성유체)의 거동만을 다룬다. 그러나 본 연구에서는 유체 및 매질의 압축성과 완전한 삼상 오염물을 가정하는, 분율흐름접근방식을 따르는 수학적 지배방정식을 개발하고, 이를 토대로 압축성 유체 및 매질을 고려한 삼상흐름 수치모의프로그램을 개발하였다 개발된 삼상흐름 수치모의프로그램 (Compressible Multiphase Flow Simulator, CMPS)을 검증하기 위해서 기존에 개발된 비압축성 유체 및 매질을 고려한 삼상흐름수치모의 프로그램인 MPS (Suk and Yeh, 2007; Suk and Yeh, 2008) 및 해석해를 통해서 간접적으로 비교 검증하였다. 비교결과 CMPS와 MPS의 결과와 해석해들은 서로 잘 일치하였다. 따라서 CMPS는 압축성 유체 및 매질을 고려한 삼상흐름 수치모의를 구현 할 수 있는 가능성을 가진다.

• Ocean Systems Engineering
• /
• 제1권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년 춘계학술대회논문집
• /
• 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.

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

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

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1996년도 추계학술대회 논문집
• /
• 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.