Skip to main content
    • Aa
    • Aa

A Comparison Study of Numerical Methods for Compressible Two-Phase Flows

  • Jianyu Lin (a1), Hang Ding (a1), Xiyun Lu (a1) and Peng Wang (a2)

In this article a comparison study of the numerical methods for compressible two-phase flows is presented. Although many numerical methods have been developed in recent years to deal with the jump conditions at the fluid-fluid interfaces in compressible multiphase flows, there is a lack of a detailed comparison of these methods. With this regard, the transport five equation model, the modified ghost fluid method and the cut-cell method are investigated here as the typical methods in this field. A variety of numerical experiments are conducted to examine their performance in simulating inviscid compressible two-phase flows. Numerical experiments include Richtmyer-Meshkov instability, interaction between a shock and a rectangle SF 6 bubble, Rayleigh collapse of a cylindrical gas bubble in water and shock-induced bubble collapse, involving fluids with small or large density difference. Based on the numerical results, the performance of the method is assessed by the convergence order of the method with respect to interface position, mass conservation, interface resolution and computational efficiency.

Corresponding author
*Corresponding author. Email: (J. Y. Lin), (H. Ding), (X. Y. Lu), (P. Wang)
Hide All
[1] AbgrallR., How to prevent pressure oscillations in multicomponent flow calculations: A quasi conservative approach, J. Comput. Phys., 125 (1996), pp. 150160.
[2] ShyueK.-M., An efficient shock-capturing algorithm for compressible multicomponent problems, J. Comput. Phys., 142(1) (1998), pp. 208242.
[3] SaurelR. and AbgrallR., A multiphase Godunov method for compressible multifluid and multiphase flows, J. Comput. Phys., 150(2) (1999), pp. 425467.
[4] AllaireG., ClercS. and KokhS., A five-equation model for the simulation of interfaces between compressible fluids, J. Comput. Phys., 181(2) (2002), pp. 577616.
[5] SaurelR., GavrilyukS. and RenaudF., A multiphase model with internal degrees of freedom: Application to shock-bubble interaction, J. Fluid Mech., 495 (2003), pp. 283321.
[6] MarquinaA. and MuletP., A flux-split algorithm applied to conservative models for multicomponent compressible flows, J. Comput. Phys., 185(1) (2003), pp. 120138.
[7] ChangC.-H. and LiouM.-S., Robust and accurate approach to computing compressible multiphase flow: Stratified flow model and AUSM+-up scheme, J. Comput. Phys., 225(1) (2007), pp. 840873.
[8] JohnsenT. and ColoniusE., Implementation of WENO schemes in compressible multicomponent flow problems, J. Comput. Phys., 219(4) (2006), pp. 715732.
[9] KokhS. and LagoutièreF., An anti-diffusive numerical scheme for the simulation of interfaces between compressible fluids by means of a five-equation model, J. Comput. Phys., 229(8) (2010), pp. 27732809.
[10] ShuklaR. K., PantanoC. and FreundJ. B., An interface capturing method for the simulation of multi-phase compressible flows, J. Comput. Phys., 229(19) (2010), pp. 74117439.
[11] SoK. K., HuX. Y. and AdamsN. A., Anti-diffusion interface sharpening technique for two-phase compressible flow simulations, J. Comput. Phys., 231(11) (2012), pp. 43044323.
[12] TiwariA., FreundJ. B. and PantanoC., A diffuse interface model with immiscibility preservation, J. Comput. Phys., 252 (2013), pp. 290309.
[13] ShuklaR. K., Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows, J. Comput. Phys., 276 (2014), pp. 508540.
[14] DeligantM., SpecklinM. and KhelladiS., A naturally anti-diffusive compressible two phases Kapila model with boundedness preservation coupled to a high order finite volume solver, Comput. Fluids, 114 (2015), pp. 265273.
[15] FedkiwR. P., AslamT., MerrimanB. and OsherS., A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method), J. Comput. Phys., 152(2) (1999), pp. 457492.
[16] LiuT. G., KhooB. C. and YeoK. S., Ghost fluid method for strong shock impacting on material interface, J. Comput. Phys., 190(2) (2003), pp. 651681.
[17] WangC. W., LiuT. G. and KhooB. C., A real ghost fluid method for the simulation of multimedium compressible flow, SIAM J. Sci. Comput., 28(1) (2006), pp. 278302.
[18] SambasivanS. K. and UdaykumarH. S., Ghost fluid method for strong shock interactions part 1: Fluid-fluid interfaces, AIAA J., 47(12) (2009), pp. 29072922.
[19] TerashimaH. and TryggvasonG., A front-tracking/ghost-fluid method for fluid interfaces in compressible flows, J. Comput. Phys., 228(11) (2009), pp. 40124037.
[20] HuX. Y., KhooB. C., AdamsN. A. and HuangF. L., A conservative interface method for compressible flows, J. Comput. Phys., 219(2) (2006), pp. 553578.
[21] ChangC.-H., DengX. and TheofanousT. G., Direct numerical simulation of interfacial instabilities: A consistent, conservative, all-speed, sharp-interface method, J. Comput. Phys., 242 (2013), pp. 946990.
[22] NourgalievR. R., LiouM.-S. and TheofanousT. G., Numerical prediction of interfacial instabilities: Sharp interface method (SIM), J. Comput. Phys., 227(8) (2008), pp. 39403970.
[23] KimH. and LiouM.-S., Adaptive Cartesian sharp interface method for three-dimensional multiphase flows, AIAA Paper, (2009), pp. 20094153.
[24] BoW. and GroveJ. W., A volume of fluid method based ghost fluid method for compressible multi-fluid flows, Comput. Fluids, 90 (2014), pp. 113122.
[25] LinJ., ShenY., DingH., LiuN. and LuX., Simulation of compressible two-phase flows with topology change of fluid-fluid interface by a robust cut-cell method, preparing.
[26] RoeP. L., Approximate Riemann solvers, parameter vectors and difference schemes, J. Comput. Phys., 43 (1981), pp. 357372.
[27] HartmannD., MeinkeM. and SchröderW., The constrained reinitialization equation for level set methods, J. Comput. Phys., 229(5) (2010), pp. 15141535.
[28] OsherS. and FedkiwR., Level Set Methods and Dynamic Implicit Surfaces, Springer, 2003.
[29] HartenA., LaxP. D. and Van LeerB., On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Rev., 25(1) (1983), pp. 3561.
[30] ShuC.-W. and OsherS., Efficient implementation of essentially non-oscillatory shock-capturing schemes, J. Comput. Phys., 77(2) (1988), pp. 439471.
[31] HolmesR. L., GroveJ.W. and SharpD. H., Numerical investigation of Richtmyer–Meshkov instability using front tracking, J. Fluid Mech., 301 (1995), pp. 5164.
[32] UllahM. A., MaoD.-K and GaoW.-B., Numerical simulations of Richtmyer–Meshkov instabilities using conservative front-tracking method, Appl. Math. Mech., 32(1) (2011), pp. 119132.
[33] BatesK. R., NikiforakisN. and HolderD., Richtmyer–Meshkov instability induced by the interaction of a shock wave with a rectangular block of SF6 , Phys. Fluids, 19 (2007), 036101.
[34] JohnsenT. and ColoniusE., Shock-induced collapse of a gas bubble in shockwave lithotripsy, J. Acoust. Soc. Am., 124(4) (2008), pp. 20112020.
[35] BallG. J., HowellB. P., LeightonT. G. and SchofieldM. J., Shock-induced collapse of a cylindrical air cavity in water: A free-Lagrange simulation, Shock Waves, 10(4) (2000), pp. 265276.
[36] HuX. Y. and KhooB. C., An interface interaction method for compressible multifluids, J. Comput. Phys., 198(1) (2004), pp. 3564.
[37] NourgalievR. R. and DinhT. N., Adaptive characteristics-based matching for compressible multifluid dynamics, J. Comput. Phys., 213(2) (2006), pp. 500529.
[38] HawkerN. A. and VentikosY., Interaction of a strong shockwave with a gas bubble in a liquid medium: A numerical study, J. Fluid Mech., 701 (2012), pp. 5997.
Recommend this journal

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

Advances in Applied Mathematics and Mechanics
  • ISSN: 2070-0733
  • EISSN: 2075-1354
  • URL: /core/journals/advances-in-applied-mathematics-and-mechanics
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: 57 *
Loading metrics...

Abstract views

Total abstract views: 171 *
Loading metrics...

* Views captured on Cambridge Core between 11th July 2017 - 24th October 2017. This data will be updated every 24 hours.