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

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

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.330-340
• /
• 2005

Compressible two-phase flow is analyzed based on the arbitrary Lagrangian-Eulerian (ALE) formulation. For water, Tamman type stiffened equation of state is used. Numerical fluxes are calculated using the ALE two-phase Godunov scheme which assumes only that the speed of sound and pressure can be provided whenever density and internal energy are given. Effects of the approximations of a material interface speed are Investigated h method Is suggested to assign a rigid body boundary condition effectively To validate the developed code, several well-known problems are calculated and the results are compared with analytic or other numerical solutions including a single material Sod shock tube problem and a gas/water shock tube problem The code is applied to analyze the refraction and transmission of shock waves which are impacting on a water-gas interface from gas or water medium.

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