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

• Published in: Journal of computational physics Vol. 333; pp. 247 - 268
• Main Authors: He, Zhiwei, Tian, Baolin, Zhang, Yousheng, Gao, Fujie
• Format: Journal Article
• Language: English
• Published: Cambridge Elsevier Inc 15.03.2017; Elsevier Science Ltd
• Subjects: Algorithms; Compressibility; Computational physics; Computer simulation; Derivatives; Equations of state; Galling; High-order; Interface sharpening; Mathematical models; Mie–Grüneisen-type equation of state; Miscibility; Multi-material flow; Numerical analysis; Numerical dissipation; Reconstruction; Roe approximation; Sharpening; Studies; Roe approximation; Interface sharpening; High-order; Multi-material flow; Mie–Grüneisen-type equation of state
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2016.12.035

The present work focuses on the simulation of immiscible compressible multi-material flows with the Mie–Grüneisen-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 Mie–Grüneisen-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.