A higher‐order unsplit 2D direct Eulerian finite volume method for two‐material compressible flows based on the MOOD paradigms

SUMMARY A higher-order unsplit multi-dimensional discretization of the diffuse interface model for two-material compressible flows proposed by R. Saurel, F. Petitpas and R. A. Berry in 2009 is developed. The proposed higher-order method is based on the concepts of the Multidimensional Optimal Order Detection (MOOD) method introduced in three recent papers for single-material flows. The first-order unsplit multi-dimensional Finite Volume discretization presented by SPB serves as foundation for the development of the higher-order unlimited schemes. Specific detection criteria along with a novel decrementing algorithm for the MOOD method are designed in order to deal with the complexity of multi-material flows. Numerically, we compare errors and computational times on several 1D problems (stringent shock tube and cavitation problems) computed on 2D meshes with the second- and fourth-order MOOD methods using a classical MUSCL method as reference. Several simulations of a 2D shocked R22 bubble in the air are also presented on Cartesian and unstructured meshes with the second- and fourth-order MOOD methods, and qualitative comparisons confirm the conclusions obtained with 1D problems. These numerical results demonstrate the robustness of the MOOD approach and the interest of using more than second-order methods even for locally singular solutions of complex physics models. Copyright © 2014 John Wiley & Sons, Ltd.

[1]  Yong-Tao Zhang,et al.  Resolution of high order WENO schemes for complicated flow structures , 2003 .

[2]  Richard Saurel,et al.  A multiphase model with internal degrees of freedom: application to shock–bubble interaction , 2003, Journal of Fluid Mechanics.

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

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

[5]  Stéphane Clain,et al.  Monoslope and multislope MUSCL methods for unstructured meshes , 2010, J. Comput. Phys..

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

[7]  Rémi Abgrall,et al.  Efficient numerical approximation of compressible multi-material flow for unstructured meshes , 2003 .

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

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

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

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

[12]  Michael Dumbser,et al.  A New Family of High Order Unstructured MOOD and ADER Finite Volume Schemes for Multidimensional Systems of Hyperbolic Conservation Laws , 2014 .

[13]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[14]  E. Toro.  The Equations of Fluid Dynamics , 1999 .

[15]  Stéphane Clain,et al.  Improved detection criteria for the multi-dimensional optimal order detection (MOOD) on unstructured meshes with very high-order polynomials , 2012 .

[16]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[17]  J. Grove.  Pressure-velocity equilibrium hydrodynamic models , 2010 .

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

[19]  John W. Grove,et al.  A volume of fluid method based ghost fluid method for compressible multi-fluid flows , 2014 .

[20]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

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

[22]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[23]  Rémi Abgrall,et al.  Discrete equations for physical and numerical compressible multiphase mixtures , 2003 .

[24]  Stéphane Clain,et al.  A high-order finite volume method for systems of conservation laws - Multi-dimensional Optimal Order Detection (MOOD) , 2011, J. Comput. Phys..

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

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

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

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

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

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

[31]  Stéphane Clain,et al.  The Multidimensional Optimal Order Detection method in the three‐dimensional case: very high‐order finite volume method for hyperbolic systems , 2013 .

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

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