A finite volume WENO scheme for immiscible inviscid two-phase flows

Abstract This article presents a two-phase finite volume (FV) method based on the concept of equilibrium volume fraction for simulating inviscid two-phase flows. A reduced three-equation hyperbolic system is adopted as the governing equations, which is closed with a barotropic equation of state for both phases. A modified weighted essentially non-oscillatory (WENO) scheme is proposed for reconstructing the fluid state at cell-interfaces, which can avoid strictly pressure or velocity oscillations near material interfaces in the test-case of the advection of an isolated interface, while retaining the convergence order of the WENO scheme in 1D and 2D cases. Besides, with the generalized Riemann invariants (GRI), the eigenstructure of the hyperbolic system is analyzed and an approximate Riemann flux solver is proposed based on the linearization of the GRI. The semi-discrete system is explicitly integrated in time by means of the classical 4 th -order Runge-Kutta (RK4) scheme. The proposed scheme is validated in a series of two-phase test-cases, such as advection of a two-phase vortex, linear sloshing, long-time wave propagation and dam-break flows, for which good agreements are obtained with the references.

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

[2]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[3]  Boniface Nkonga,et al.  Towards the direct numerical simulation of nucleate boiling flows , 2014 .

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

[5]  Benjamin Bouscasse,et al.  CN-Stream: Open-source library for nonlinear regular waves using stream function theory , 2019 .

[6]  M. J. Castro,et al.  FORCE schemes on unstructured meshes II: Non-conservative hyperbolic systems , 2010 .

[7]  Chi-Wang Shu,et al.  High order WENO and DG methods for time-dependent convection-dominated PDEs: A brief survey of several recent developments , 2016, J. Comput. Phys..

[8]  Tim Colonius,et al.  Finite-volume WENO scheme for viscous compressible multicomponent flows , 2014, J. Comput. Phys..

[9]  M. Dumbser,et al.  High order space–time adaptive ADER-WENO finite volume schemes for non-conservative hyperbolic systems , 2013, 1304.5408.

[10]  Xiangxiong Zhang,et al.  Positivity-preserving high order finite difference WENO schemes for compressible Euler equations , 2012, J. Comput. Phys..

[11]  D. Drew.  Mathematical Modeling of Two-Phase Flow , 1983 .

[12]  Chi-Wang Shu,et al.  On positivity preserving finite volume schemes for Euler equations , 1996 .

[13]  Eleuterio F. Toro,et al.  Finite-volume WENO schemes for three-dimensional conservation laws , 2004 .

[14]  John D. Fenton,et al.  A Fourier approximation method for steady water waves , 1981, Journal of Fluid Mechanics.

[15]  M. Baer,et al.  Mathematical analysis of a two-phase continuum mixture theory , 1992 .

[16]  Vincent Perrier,et al.  Runge-Kutta discontinuous Galerkin method for the approximation of Baer and Nunziato type multiphase models , 2012, J. Comput. Phys..

[17]  Richard Saurel,et al.  A general formulation for cavitating, boiling and evaporating flows , 2016 .

[18]  Philippe Villedieu,et al.  An accurate low-Mach scheme for a compressible two-fluid model applied to free-surface flows , 2013, J. Comput. Phys..

[19]  Guillaume Oger,et al.  Comparisons of weakly-compressible and truly incompressible approaches for viscous flow into a high-order Cartesian-grid finite volume framework , 2019, J. Comput. Phys. X.

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

[21]  Eric Johnsen,et al.  A high-order accurate five-equations compressible multiphase approach for viscoelastic fluids and solids with relaxation and elasticity , 2019, J. Comput. Phys..

[22]  P. B. Whalley,et al.  Measurements of the speed of sound in air-water flows , 1997 .

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

[24]  Antonio Souto-Iglesias,et al.  Experimental investigation of dynamic pressure loads during dam break , 2013, 1308.0115.

[25]  Eleuterio F. Toro,et al.  A flux splitting method for the Baer-Nunziato equations of compressible two-phase flow , 2016, J. Comput. Phys..

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

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

[28]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[29]  David Le Touzé,et al.  A weakly-compressible Cartesian grid approach for hydrodynamic flows , 2017, Comput. Phys. Commun..

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

[31]  Nikolaus A. Adams,et al.  Adaptive multi-resolution method for compressible multi-phase flows with sharp interface model and pyramid data structure , 2014, J. Comput. Phys..

[32]  Xiangxiong Zhang,et al.  On maximum-principle-satisfying high order schemes for scalar conservation laws , 2010, J. Comput. Phys..

[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]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[36]  P. Lax.  Hyperbolic systems of conservation laws II , 1957 .

[37]  Derek M. Causon,et al.  On the Choice of Wavespeeds for the HLLC Riemann Solver , 1997, SIAM J. Sci. Comput..

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

[39]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[40]  S. Osher,et al.  Simplified Discretization of Systems of Hyperbolic Conservation Laws Containing Advection Equations , 2000, Journal of Computational Physics.

[41]  Bing Wang,et al.  An incremental-stencil WENO reconstruction for simulation of compressible two-phase flows , 2017, International Journal of Multiphase Flow.

[42]  H. Bruce Stewart,et al.  Two-phase flow: Models and methods , 1984 .

[43]  Victor H. Ransom,et al.  Hyperbolic two-pressure models for two-phase flow , 1984 .