Characteristic-based and interface-sharpening algorithm for high-order simulations of immiscible compressible multi-material flows

The present work focuses on the simulation of immiscible compressible multi-material flows with the MieGrneisen-type equation of state governed by the non-conservative five-equation model [1]. Although low-order single fluid schemes have already been adopted to provide some feasible results, the application of high-order schemes (introducing relatively small numerical dissipation) to these flows may lead to results with severe numerical oscillations. Consequently, attempts to apply any interface-sharpening techniques to stop the progressively more severe smearing interfaces for a longer simulation time may result in an overshoot increase and in some cases convergence to a non-physical solution occurs. This study proposes a characteristic-based interface-sharpening algorithm for performing high-order simulations of such flows by deriving a pressure-equilibrium-consistent intermediate state (augmented with approximations of pressure derivatives) for local characteristic variable reconstruction and constructing a general framework for interface sharpening. First, by imposing a weak form of the jump condition for the non-conservative five-equation model, we analytically derive an intermediate state with pressure derivatives treated as additional parameters of the linearization procedure. Based on this intermediate state, any well-established high-order reconstruction technique can be employed to provide the state at each cell edge. Second, by designing another state with only different reconstructed values of the interface function at each cell edge, the advection term in the equation of the interface function is discretized twice using any common algorithm. The difference between the two discretizations is employed consistently for interface compression, yielding a general framework for interface sharpening. Coupled with the fifth-order improved accurate monotonicity-preserving scheme [2] for local characteristic variable reconstruction and the tangent of hyperbola for the interface capturing scheme [3] for designing other reconstructed values of the interface function, the present algorithm is examined using some typical tests, with the MieGrneisen-type equation of state used for characterizing the materials of interest in both one- and two-dimensional spaces. The results of these tests verify the effectiveness of the present algorithm: essentially non-oscillatory and interface-sharpened results are obtained.

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

[2]  D. Fu,et al.  A 5th order monotonicity-preserving upwind compact difference scheme , 2011 .

[3]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[4]  Huijing Yuan,et al.  Experimental studies of surface waves inside a cylindrical container , 2011, Journal of Fluid Mechanics.

[5]  Feng Xiao,et al.  An Eulerian interface sharpening algorithm for compressible two-phase flow: The algebraic THINC approach , 2014, J. Comput. Phys..

[6]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

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

[8]  G. M. Ward,et al.  A hybrid, center-difference, limiter method for simulations of compressible multicomponent flows with Mie-Grüneisen equation of state , 2010, J. Comput. Phys..

[9]  Grégoire Allaire,et al.  A five-equation model for the simulation of interfaces between compressible fluids , 2002 .

[10]  Feng Xiao,et al.  A simple algebraic interface capturing scheme using hyperbolic tangent function , 2005 .

[11]  Jianxian Qiu,et al.  On the construction, comparison, and local characteristic decomposition for high-Order central WENO schemes , 2002 .

[12]  Luigi Vigevano,et al.  An Evaluation of Roe's Scheme Generalizations for Equilibrium Real Gas Flows , 1997 .

[13]  Zhiwei He,et al.  An improved accurate monotonicity-preserving scheme for the Euler equations , 2016 .

[14]  Carlos Pantano,et al.  A diffuse interface model with immiscibility preservation , 2013, 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]  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..

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

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

[19]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[20]  Maolin Tang A Hybrid , 2010 .

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

[22]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[23]  Huanan Yang,et al.  An artificial compression method for ENO schemes - The slope modification method. [essentially nonoscillatory , 1990 .

[24]  Bruno Després,et al.  Contact Discontinuity Capturing Schemes for Linear Advection and Compressible Gas Dynamics , 2002, J. Sci. Comput..

[25]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[26]  A. Harten,et al.  The artificial compression method for computation of shocks and contact discontinuities. I - Single conservation laws , 1977 .

[27]  肖锋,et al.  Revisit To the Thinc Scheme: A Simple Algebraic VOF Algorithm , 2011 .

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

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

[30]  Jack R. Edwards,et al.  An investigation of interface-sharpening schemes for multi-phase mixture flows , 2009, J. Comput. Phys..

[31]  S. Zalesak Introduction to “Flux-Corrected Transport. I. SHASTA, A Fluid Transport Algorithm That Works” , 1997 .

[32]  Xian Liang,et al.  Direct numerical simulation of compressible turbulent flows , 2010 .

[33]  A. Harten,et al.  The artificial compression method for computation of shocks and contact discontinuities: III. Self , 1978 .

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

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

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

[37]  H. Huynh,et al.  Accurate Monotonicity-Preserving Schemes with Runge-Kutta Time Stepping , 1997 .

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

[39]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

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

[41]  Rémi Abgrall,et al.  A Simple Method for Compressible Multifluid Flows , 1999, SIAM J. Sci. Comput..

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

[43]  Shiyi Chen,et al.  Experimental investigation of freely falling thin disks. Part 1. The flow structures and Reynolds number effects on the zigzag motion , 2013, Journal of Fluid Mechanics.

[44]  Nikolaus A. Adams,et al.  Anti-diffusion interface sharpening technique for two-phase compressible flow simulations , 2012, J. Comput. Phys..

[45]  Gunilla Kreiss,et al.  A conservative level set method for two phase flow II , 2005, J. Comput. Phys..

[46]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

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

[48]  Ratnesh K. Shukla,et al.  Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows , 2014, J. Comput. Phys..

[49]  Paul Glaister,et al.  An approximate linearised Riemann solver for the three-dimensional Euler equations for real gases using operator splitting , 1988 .