Construction of a low Mach finite volume scheme for the isentropic Euler system with porosity

Classical finite volume schemes for the Euler system are not accurate at low Mach number and some fixes have to be used and were developed in a vast literature over the last two decades. The question we are interested in in this article is: What about if the porosity is no longer uniform? We first show that this problem may be understood on the linear wave equation taking into account porosity. We explain the influence of the cell geometry on the accuracy property at low Mach number. In the triangular case, the stationary space of the Godunov scheme approaches well enough the continuous space of constant pressure and divergence-free velocity, while this is not the case in the Cartesian case. On Cartesian meshes, a fix is proposed and accuracy at low Mach number is proved to be recovered. Based on the linear study, a numerical scheme and a low Mach fix for the non-linear system, with a non-conservative source term due to the porosity variations, is proposed and tested.

[1]  Mai Duc Thanh,et al.  A Godunov-type scheme for the isentropic model of a fluid flow in a nozzle with variable cross-section , 2015, Appl. Math. Comput..

[2]  S. Schochet.  Fast Singular Limits of Hyperbolic PDEs , 1994 .

[3]  Laurent Gosse,et al.  A Well-Balanced Scheme Using Non-Conservative Products Designed for Hyperbolic Systems of Conservati , 2001 .

[4]  David Rochette,et al.  First- and second-order finite volume methods for the one-dimensional nonconservative Euler system , 2009, J. Comput. Phys..

[5]  M. Pelanti,et al.  A Roe-type scheme with low Mach number preconditioning for a two-phase compressible flow model with pressure relaxation , 2016 .

[6]  Pascal Galon,et al.  A Mach-Sensitive Splitting Approach for Euler-like Systems , 2018 .

[7]  D. Rochette,et al.  Numerical Investigations on the Pressure Wave Absorption and the Gas Cooling Interacting in a Porous Filter, During an Internal arc Fault in a Medium-Voltage Cell , 2008, IEEE Transactions on Power Delivery.

[8]  A. Majda,et al.  Singular limits of quasilinear hyperbolic systems with large parameters and the incompressible limit , 1981 .

[9]  P. Colella,et al.  A Projection Method for Low Speed Flows , 1999 .

[10]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[11]  P. Degond,et al.  All speed scheme for the low Mach number limit of the isentropic Euler equations , 2009, 0908.1929.

[12]  Philipp Birken,et al.  L2Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers , 2016 .

[13]  R. Nicolaides.  Analysis and convergence of the MAC scheme. I : The linear problem , 1992 .

[14]  Stéphane Clain,et al.  Numerical scheme to complete a compressible gas flow in variable porosity media , 2005 .

[15]  J. Bear,et al.  Introduction to Modeling of Transport Phenomena in Porous Media , 1990 .

[16]  James R. DeBonis,et al.  A Survey of the Isentropic Euler Vortex Problem Using High-Order Methods , 2015 .

[17]  Wang Chi-Shu,et al.  Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws , 1997 .

[18]  Jean-Marc Hérard,et al.  A sequel to a rough Godunov scheme: application to real gases , 2000 .

[19]  T. Gallouët,et al.  Some approximate Godunov schemes to compute shallow-water equations with topography , 2003 .

[20]  Christian Klingenberg,et al.  An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity , 2019, J. Comput. Phys..

[21]  Thierry Gallouët,et al.  Un schéma de Godunov approché , 1996 .

[22]  E. Turkel,et al.  PRECONDITIONING TECHNIQUES IN COMPUTATIONAL FLUID DYNAMICS , 1999 .

[23]  R. Nicolaides.  Direct discretization of planar div-curl problems , 1992 .

[24]  A. Majda.  Compressible fluid flow and systems of conservation laws in several space variables , 1984 .

[25]  Felix Rieper,et al.  On the dissipation mechanism of upwind-schemes in the low Mach number regime: A comparison between Roe and HLL , 2010, J. Comput. Phys..

[26]  Samuel Kokh,et al.  An all-regime Lagrange-Projection like scheme for the gas dynamics equations on unstructured meshes , 2016 .

[27]  Marica Pelanti,et al.  Low Mach number preconditioning techniques for Roe-type and HLLC-type methods for a two-phase compressible flow model , 2017, Appl. Math. Comput..

[28]  H. Guillard,et al.  On the behaviour of upwind schemes in the low Mach number limit , 1999 .

[29]  Mai Duc Thanh,et al.  Numerical Solutions to Compressible Flows in a Nozzle with Variable Cross-section , 2005, SIAM J. Numer. Anal..

[30]  Stéphane Dellacherie,et al.  Preliminary results for the study of the godunov scheme applied to the linear wave equation with porosity at low mach number , 2015 .

[31]  D. Arnold,et al.  A uniformly accurate finite element method for the Reissner-Mindlin plate , 1989 .

[32]  Stéphane Dellacherie,et al.  Construction of modified Godunov-type schemes accurate at any Mach number for the compressible Euler system , 2016 .

[33]  Simon Delmas,et al.  A low Mach correction able to deal with low Mach acoustics , 2019, J. Comput. Phys..

[34]  Felix Rieper,et al.  A low-Mach number fix for Roe's approximate Riemann solver , 2011, J. Comput. Phys..

[35]  Stéphane Dellacherie,et al.  Analysis of Godunov type schemes applied to the compressible Euler system at low Mach number , 2010, J. Comput. Phys..

[36]  Stéphane Clain,et al.  Two-dimensional computation of gas flow in a porous bed characterized by a porosity jump , 2006, J. Comput. Phys..

[37]  Xue-song Li,et al.  An All-Speed Roe-type scheme and its asymptotic analysis of low Mach number behaviour , 2008, J. Comput. Phys..

[38]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[39]  Wasilij Barsukow,et al.  A Numerical Scheme for the Compressible Low-Mach Number Regime of Ideal Fluid Dynamics , 2016, Journal of Scientific Computing.

[40]  Georg Bader,et al.  The influence of cell geometry on the accuracy of upwind schemes in the low mach number regime , 2009, J. Comput. Phys..

[41]  Chun-wei Gu,et al.  Development of Roe-type scheme for all-speed flows based on preconditioning method , 2009 .

[42]  M. Thanh,et al.  The Riemann Problem for Fluid Flows in a Nozzle with Discontinuous Cross-Section , 2003 .

[43]  Hervé Guillard.  On the behavior of upwind schemes in the low Mach number limit. IV: P0 approximation on triangular and tetrahedral cells , 2009 .

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

[45]  Wasilij Barsukow,et al.  A low-Mach Roe-type solver for the Euler equations allowing for gravity source terms , 2017 .

[46]  Stéphane Dellacherie,et al.  The influence of cell geometry on the Godunov scheme applied to the linear wave equation , 2010, J. Comput. Phys..