The ghost fluid method for compressible gas-water simulation

An analysis is carried out for the ghost fluid method (GFM) based algorithm as applied to the gas-water Riemann problems, which can be construed as two single-medium GFM Riemann problems. It is found that the inability to provide correct and consistent Riemann waves in the respective real fluids by these two GFM Riemann problems may lead to inaccurate numerical results. Based on this finding, two conditions are suggested and imposed for the ghost fluid status in order to ensure that correct and consistent Riemann waves are provided in the respective real fluids during the numerical decomposition of the singularity. Using these two conditions to analyse some of the existing GFM-based algorithms such as the original GFM [J. Comput. Phys. 152 (1999) 457], the new version GFM [J. Comput. Phys. 166 (2001) 1; J. Comput. Phys. 175 (2002) 200] and the modified GFM (MGFM) [J. Comput. Phys. 190 (2003) 651], it is found that there are ranges of conditions for each type of solution where either the original GFM or the new version GFM or both are unable to provide correct or consistent Riemann waves in one of the real fluids. Within these ranges, examples can be found such that either the original GFM or the new version GFM or both are unable to provide accurate results. The MGFM is also found to encounter difficulties when applied to nearly cavitating flow. Various examples are presented to demonstrate the conclusions obtained. The MGFM with proposed modification when applied to nearly cavitating flow is then found to be quite robust and can provide relatively reasonable results.

[1]  Barry Koren,et al.  Riemann-problem and level-set approaches for homentropic two-fluid flow computations , 2002 .

[2]  Hansong Tang,et al.  A Second-Order Accurate Capturing Scheme for 1D Inviscid Flows of Gas and Water with Vacuum Zones , 1996 .

[3]  Boo Cheong Khoo,et al.  Ghost fluid method for strong shock impacting on material interface , 2003 .

[4]  Ning Zhao,et al.  Conservative front tracking and level set algorithms , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Barry Koren,et al.  A pressure-invariant conservative Godunov-type method for barotropic two-fluid flows , 2003 .

[6]  Khoon Seng Yeo,et al.  The simulation of compressible multi-medium flow. I. A new methodology with test applications to 1D gas–gas and gas–water cases , 2001 .

[7]  P. Roe,et al.  On Godunov-type methods near low densities , 1991 .

[8]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[9]  Charlie H. Cooke,et al.  On the Riemann problem for liquid or gas‐liquid media , 1994 .

[10]  Boo Cheong Khoo,et al.  The simulation of compressible multi-medium flow: II. Applications to 2D underwater shock refraction , 2001 .

[11]  R. Fedkiw,et al.  Coupling an Eulerian fluid calculation to a Lagrangian solid calculation with the ghost fluid method , 2002 .

[12]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[13]  Ronald Fedkiw,et al.  An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows , 1999 .

[14]  R. Abgrall How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[15]  B. Larrouturou How to preserve the mass fractions positivity when computing compressible multi-component flows , 1991 .

[16]  Tariq D. Aslam,et al.  A level-set algorithm for tracking discontinuities in hyperbolic conservation laws , 2001 .

[17]  Richard Saurel,et al.  A Riemann Problem Based Method for the Resolution of Compressible Multimaterial Flows , 1997 .

[18]  Patrick Jenny,et al.  Correction of Conservative Euler Solvers for Gas Mixtures , 1997 .

[19]  Tariq D. Aslam,et al.  A Level Set Algorithm for Tracking Discontinuities in Hyperbolic Conservation Laws II: Systems of Equations , 2003, J. Sci. Comput..

[20]  Smadar Karni,et al.  Multicomponent Flow Calculations by a Consistent Primitive Algorithm , 1994 .

[21]  R. Fedkiw,et al.  A numerical method for two-phase flow consisting of separate compressible and incompressible regions , 2000 .

[22]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .