Anti-diffusion interface sharpening technique for two-phase compressible flow simulations

In this paper we propose an interface sharpening technique for two-phase compressible-flow simulations based on volume-of-fluid methods. The idea of sharpening the two-fluid interface is to provide a correction algorithm which can be applied as post-processing to the volume-fraction field after each time step. For this purpose an anti-diffusion equation, i.e. a diffusion equation with a positive diffusion coefficient, is solved to counter-act the numerical diffusion resulting from the underlying VOF discretization. The numerical stability and volume-fraction boundedness in solving the anti-diffusion equation are ensured by a specified discretization scheme. No interface reconstruction and interface normal calculation are required in this method. All flow variables are updated with the sharpened volume-fraction field for ensuring the consistency of the variables, and the update of the phase mass, momentum and energy is conservative. Numerical results for shock-tube and shock-bubble interactions based on the ideal-gas EOS and shock contact problems based on the Mie-Gruneisen EOS show an improved interface resolution. The large-scale interface structures are in good agreement with reference results, and finer small-scale interface structures are recovered in a consistent manner as the grid resolution increases. As compared with reference high grid-resolution numerical results based on AMR algorithms, the interface roll-up phenomena due to the Richtmyer-Meshkov instability and the Kelvin-Helmholtz instability are recovered reliably for shock-bubble interactions involving different ideal gases.

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

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

[3]  E. Puckett,et al.  A High-Order Godunov Method for Multiple Condensed Phases , 1996 .

[4]  Nikolaus A. Adams,et al.  Anti-diffusion method for interface steepening in two-phase incompressible flow , 2011, J. Comput. Phys..

[5]  Santhosh K. Shankar,et al.  Numerical Simulation of Multicomponent Shock Accelerated Flows and Mixing using Localized Artificial Diffusivity Method , 2010 .

[6]  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..

[7]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[8]  J. P. Boris,et al.  Vorticity generation by shock propagation through bubbles in a gas , 1988, Journal of Fluid Mechanics.

[9]  R. Bonazza,et al.  Shock-Bubble Interactions , 2011 .

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

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

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

[13]  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..

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

[15]  D. Drikakis,et al.  Numerical experiments using high-resolution schemes for unsteady, inviscid, compressible flows , 2004 .

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

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

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

[19]  Diego Rossinelli,et al.  High order finite volume methods on wavelet-adapted grids with local time-stepping on multicore architectures for the simulation of shock-bubble interactions , 2010, J. Comput. Phys..

[20]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

[21]  Kathy Prestridge,et al.  An experimental investigation of mixing mechanisms in shock-accelerated flow , 2008, Journal of Fluid Mechanics.

[22]  Phillip Colella,et al.  An adaptive multifluid interface-capturing method for compressible flow in complex geometries , 1995 .

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

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

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

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

[27]  Ravi Samtaney,et al.  On initial‐value and self‐similar solutions of the compressible Euler equations , 1996 .

[28]  L. Margolin.  Introduction to “An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds” , 1997 .

[29]  Alexander Kurganov,et al.  Interface tracking method for compressible multifluids , 2008 .

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

[31]  Song Jiang,et al.  A second-order γ-model BGK scheme for multimaterial compressible flows , 2007 .

[32]  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..

[33]  Soshi Kawai,et al.  A high‐resolution scheme for compressible multicomponent flows with shock waves , 2011 .

[34]  Ravi Samtaney,et al.  Circulation deposition on shock-accelerated planar and curved density-stratified interfaces: models and scaling laws , 1994, Journal of Fluid Mechanics.

[35]  Wenjun Sun,et al.  A γ‐DGBGK scheme for compressible multi‐fluids , 2011 .

[36]  M. Brouillette.  THE RICHTMYER-MESHKOV INSTABILITY , 2002 .

[37]  Edward E. Zukoski,et al.  A model for characterization of a vortex pair formed by shock passage over a light-gas inhomogeneity , 1994, Journal of Fluid Mechanics.

[38]  William J. Rider,et al.  An adaptive Riemann solver using a two-shock approximation 1 1 This work performed under the auspice , 1999 .

[39]  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..

[40]  Keh-Ming Shyue,et al.  A wave-propagation based volume tracking method for compressible multicomponent flow in two space dimensions , 2006, J. Comput. Phys..

[41]  Mark H. Anderson,et al.  A computational parameter study for the three-dimensional shock–bubble interaction , 2007, Journal of Fluid Mechanics.

[42]  Keh-Ming Shyue,et al.  A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Grüneisen equation of state , 2001 .

[43]  Ravi Samtaney,et al.  Baroclinic circulation generation on shock accelerated slow/fast gas interfaces , 1998 .