Skip to main content Accessibility help

A Conservative Modification to the Ghost Fluid Method for Compressible Multiphase Flows

  • Wei Liu (a1), Li Yuan (a1) and Chi-Wang Shu (a2)


A conservative modification to the ghost fluid method (GFM) is developed for compressible multiphase flows. The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance. We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell. The modification procedure can be used on the GFM with any base schemes. In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization. The level set method is used to capture the interface. Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.


Corresponding author



Hide All
[1]Abgrall, R., How to prevent pressure oscillations in multicomponent flow calculations: a quasi-conservative approach, J. Comput. Phys., 125 (1996), 150160.
[2]Abgrall, R. and Karni, S., Computations of compressible multifluids, J. Comput. Phys., 169 (2001), 594623.
[3]Adalsteinsson, D. and Sethian, J.A., A fast level set method for propagating interfaces, J. Comput. Phys., 118 (1995), 269277.
[4]Bailey, D.A., Sweby, P.K. and Glaister, P., A ghost fluid, moving finite volume plus continuous remap method for compressible Euler flow, Int. J. Numer. Methods Fluids, 47 (2005), 833840.
[5]Fedkiw, R.P., Coupling an Eulerian fluid calculation to a Lagrangian solid calculation with the ghost fluid method, J. Comput. Phys., 175 (2002), 200224.
[6]Fedkiw, R.P., Aslam, T., Merriman, B. and Osher, S., A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the Ghost Fluid Method), J. Comput. Phys., 152 (1999), 457492.
[7]Fedkiw, R.P., Marquina, A. and Merriman, B., An isobaric fix for the overheating problem in multimaterial compressible flows, J. Comput. Phys., 148 (1999), 545578.
[8]Glimm, J., Grove, J.M., Li, X.L. and Zhao, N., Simple front tracking, in: Chen, G.Q., Di Bebe-detto, E. (Eds.), Contemporary Mathematics, vol. 238, AMS, Providence, RI, 1999.
[9]Glimm, J., Marchesin, D. and McBryan, O., A numerical method for two phase flow with an unstable interface, J. Comput. Phys., 39 (1981), 179200.
[10]Haas, J.-F. and Sturtevant, B., Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities, J. Fluid Mech., 181 (1987), 4176.
[11]Jiang, G. and Peng, D.-P., Weighted ENO schemes for Hamilton-Jacobi equations, SIAM J. Sci. Comput., 21 (2000), 21262143.
[12]Jiang, G. and Shu, C.-W., Efficient implementation of weighted ENO schemes, J. Comput. Phys., 126 (1996), 202228.
[13]Johnsen, E. and Colonius, T., Implementation of WENO schemes in compressible multi-component flow problems, J. Comput. Phys., 219 (2006), 715732.
[14]Liu, T.G., Khoo, B.C. and Wang, C.W., The ghost fluid method for compressible gas-water simulation, J. Comput. Phys., 204 (2005), 193221.
[15]Liu, T.G., Khoo, B.C. and Yeo, K.S., The simulation of compressible multi-medium flow. Part I: a new methodology with applications to 1D gas-gas and gas-water cases, Computers & Fluids, 30 (2001), 291314.
[16]Liu, T.G., Khoo, B.C. and Yeo, K.S., The simulation of compressible multi-medium flow. Part II: applications to 2D underwater shock refraction, Computers & Fluids, 30 (2001), 315337.
[17]Liu, T.G., Khoo, B.C. and Yeo, K.S., Ghost fluid method for strong shock impacting on material interface, J. Comput. Phys., 190 (2003), 651681.
[18]Marquina, A. and Mulet, P., A flux-split algorithm applied to conservative models for multi-component compressible flows, J. Comput. Phys., 185 (2003), 120138.
[19]Mulder, W., Osher, S. and Sethian, J.A., Computing interface motion in compressible gas dynamics, J. Comput. Phys., 100 (1992), 209228.
[20]Osher, S., Fedkiw, R.P., Level Set Methods and Dynamic Implicit Surfaces, Springer, 2003.
[21]Osher, S., Sethian, J.A., Fronts propagating with curvature-dependent speed: algorithms based on Hamilton Jacobi formulations, J. Comput. Phys., 79 (1988), 1249.
[22]Qiu, J., Liu, T. and Khoo, B.C., Simulations of compressible two-medium flow by Runge-Kutta discontinuous Galerkin methods with the Ghost Fluid Method, Commun. Comput. Phys., 3 (2008), 479504.
[23]Quirk, J. and Karni, S., On the dynamics of a shock-bubble interaction, J. Fluid Mech., 318 (1996), 129163.
[24]Sethian, J.A., Level Set Methods and Fast Marching Methods, Cambridge University Press, 1999.
[25]Shu, C.-W., Total-variation-diminishing time discretizations, SIAM J. Sci. Stat. Comput., 9 (1988), 10731084.
[26]Shu, C.-W. and Osher, S., Efficient implementation of essentially non-oscillatory shock-capturing schemes, J. Comput. Phys., 77 (1988), 439471.
[27]Sussman, M., Smereka, P. and Osher, S., A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114 (1994), 146159.
[28]Wang, C.W., Liu, T.G. and Khoo, B.C., A real-ghost fluid method for the simulation of multimedium compressible flow, SIAM J. Sci. Comput., 28 (2006), 278302.
Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Communications in Computational Physics
  • ISSN: 1815-2406
  • EISSN: 1991-7120
  • URL: /core/journals/communications-in-computational-physics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed