Numerical simulation of compressible two-phase flow using a diffuse interface method

Abstract In this article, a high-resolution diffuse interface method is investigated for simulation of compressible two-phase gas–gas and gas–liquid flows, both in the presence of shock wave and in flows with strong rarefaction waves similar to cavitations. A Godunov method and HLLC Riemann solver is used for discretization of the Kapila five-equation model and a modified Schmidt equation of state (EOS) is used to simulate the cavitation regions. This method is applied successfully to some one- and two-dimensional compressible two-phase flows with interface conditions that contain shock wave and cavitations. The numerical results obtained in this attempt exhibit very good agreement with experimental results, as well as previous numerical results presented by other researchers based on other numerical methods. In particular, the algorithm can capture the complex flow features of transient shocks, such as the material discontinuities and interfacial instabilities, without any oscillation and additional diffusion. Numerical examples show that the results of the method presented here compare well with other sophisticated modeling methods like adaptive mesh refinement (AMR) and local mesh refinement (LMR) for one- and two-dimensional problems.

[1]  Richard Saurel,et al.  Modelling detonation waves in condensed energetic materials: multiphase CJ conditions and multidimensional computations , 2009 .

[2]  Hiroshi Terashima,et al.  A front-tracking/ghost-fluid method for fluid interfaces in compressible flows , 2009, J. Comput. Phys..

[3]  S. F. Davis.  Simplified second-order Godunov-type methods , 1988 .

[4]  Richard Saurel,et al.  A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation , 2001, Journal of Fluid Mechanics.

[5]  Barry Koren,et al.  A new formulation of Kapila's five-equation model for compressible two-fluid flow, and its numerical treatment , 2010, J. Comput. Phys..

[6]  H. Udaykumar,et al.  Sharp interface simulations with Local Mesh Refinement for multi-material dynamics in strongly shocked flows , 2010 .

[7]  Miltiadis V. Papalexandris,et al.  An exact Riemann solver for compressible two-phase flow models containing non-conservative products , 2007, J. Comput. Phys..

[8]  Svend Tollak Munkejord,et al.  A Numerical Study of Two-Fluid Models with Pressure and Velocity Relaxation , 2010 .

[9]  M. Baer,et al.  A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials , 1986 .

[10]  R. Abgrall,et al.  A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows , 1999 .

[11]  J. Haas,et al.  Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities , 1987, Journal of Fluid Mechanics.

[12]  D. Stewart,et al.  Two-phase modeling of deflagration-to-detonation transition in granular materials: Reduced equations , 2001 .

[13]  Gretar Tryggvason,et al.  A front-tracking method with projected interface conditions for compressible multi-fluid flows , 2010 .

[14]  Richard Saurel,et al.  A relaxation-projection method for compressible flows. Part II: Artificial heat exchanges for multiphase shocks , 2007, J. Comput. Phys..

[15]  J. Freund,et al.  An interface capturing method for the simulation of multi-phase compressible flows , 2010, J. Comput. Phys..

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

[17]  T. Etoh,et al.  Cavitation induced by low-speed underwater impact , 2009 .

[18]  Klaus Hannemann,et al.  Computations of shock wave propagation with local mesh adaptation , 2009 .

[19]  Keh-Ming Shyue,et al.  An Efficient Shock-Capturing Algorithm for Compressible Multicomponent Problems , 1998 .

[20]  Chi-Wang Shu,et al.  An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method , 2010, J. Comput. Phys..

[21]  J. Sethian.  Level set methods : evolving interfaces in geometry, fluid mechanics, computer vision, and materials science , 1996 .

[22]  Michael Dumbser,et al.  A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems , 2011, J. Sci. Comput..

[23]  Frédéric Lagoutière,et al.  An anti-diffusive numerical scheme for the simulation of interfaces between compressible fluids by means of a five-equation model , 2010, J. Comput. Phys..

[24]  Tiegang Liu,et al.  Application of a one-fluid model for large scale homogeneous unsteady cavitation: The modified Schmidt model , 2006 .

[25]  Keh-Ming Shyue,et al.  A high-resolution mapped grid algorithm for compressible multiphase flow problems , 2010, J. Comput. Phys..

[26]  Shamsul Qamar,et al.  A high order kinetic flux-vector splitting method for the reduced five-equation model of compressible two-fluid flows , 2009, J. Comput. Phys..

[27]  Meng-Sing Liou,et al.  A robust and accurate approach to computing compressible multiphase flow: Stratified flow model and AUSM+-up scheme , 2007, J. Comput. Phys..

[28]  A. Huerta,et al.  Arbitrary Lagrangian–Eulerian Methods , 2004 .

[29]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[30]  James J. Quirk,et al.  On the dynamics of a shock–bubble interaction , 1994, Journal of Fluid Mechanics.

[31]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[32]  A. B. Wood,et al.  A textbook of sound , 1930 .

[33]  Hervé Guillard,et al.  A five equation reduced model for compressible two phase flow problems , 2005 .

[34]  Eleuterio F. Toro,et al.  HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow , 2010, J. Comput. Phys..

[35]  Richard Saurel,et al.  Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures , 2009, J. Comput. Phys..

[36]  Richard Saurel,et al.  A relaxation-projection method for compressible flows. Part I: The numerical equation of state for the Euler equations , 2007, J. Comput. Phys..

[37]  Ivan B. Bazhlekov,et al.  Interaction of a deformable bubble with a rigid wall at moderate Reynolds numbers , 1990, Journal of Fluid Mechanics.

[38]  Jeffrey W. Banks,et al.  A high-resolution Godunov method for compressible multi-material flow on overlapping grids , 2006, J. Comput. Phys..

[39]  Nikolaus A. Adams,et al.  On the HLLC Riemann solver for interface interaction in compressible multi-fluid flow , 2009, J. Comput. Phys..

[40]  Kiumars Mazaheri,et al.  Moment of fluid interface reconstruction method in multi-material arbitrary Lagrangian Eulerian (MMALE) algorithms , 2009 .

[41]  Steven F. Son,et al.  Two-Phase Modeling of DDT in Granular Materials: Reduced Equations , 2000 .

[42]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[43]  Eric Johnsen,et al.  Implementation of WENO schemes in compressible multicomponent flow problems , 2005, J. Comput. Phys..

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

[45]  David P. Schmidt,et al.  A moving mesh interface tracking method for 3D incompressible two-phase flows , 2007, J. Comput. Phys..

[46]  Theo G. Theofanous,et al.  Adaptive characteristics-based matching for compressible multifluid dynamics , 2006, J. Comput. Phys..

[47]  E. Puckett,et al.  Second-Order Accurate Volume-of-Fluid Algorithms for Tracking Material Interfaces , 2013 .