WENO Scheme with Subcell Resolution for Computing Nonconservative Euler Equations with Applications to One-Dimensional Compressible Two-Medium Flows

High order path-conservative schemes have been developed for solving nonconservative hyperbolic systems in (Parés, SIAM J. Numer. Anal. 44:300–321, 2006; Castro et al., Math. Comput. 75:1103–1134, 2006; J. Sci. Comput. 39:67–114, 2009). Recently, it has been observed in (Abgrall and Karni, J. Comput. Phys. 229:2759–2763, 2010) that this approach may have some computational issues and shortcomings. In this paper, a modification to the high order path-conservative scheme in (Castro et al., Math. Comput. 75:1103–1134, 2006) is proposed to improve its computational performance and to overcome some of the shortcomings. This modification is based on the high order finite volume WENO scheme with subcell resolution and it uses an exact Riemann solver to catch the right paths at the discontinuities. An application to one-dimensional compressible two-medium flows of nonconservative or primitive Euler equations is carried out to show the effectiveness of this new approach.

[1]  Zhi J. Wang,et al.  Spectral (finite) volume method for conservation laws on unstructured grids IV: extension to two-dimensional systems , 2004 .

[2]  Richard Saurel,et al.  A Riemann Problem Based Method for the Resolution of Compressible Multimaterial Flows , 1997 .

[3]  Charlie H. Cooke,et al.  On the Riemann problem for liquid or gas‐liquid media , 1994 .

[4]  I. Toumi A weak formulation of roe's approximate riemann solver , 1992 .

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

[6]  B. Larrouturou How to preserve the mass fractions positivity when computing compressible multi-component flows , 1991 .

[7]  Jianxian Qiu,et al.  Simulations of compressible two-medium flow by Runge-Kutta discontinuous Galerkin methods with the ghost fluid method , 2008 .

[8]  Manuel Jesús Castro Díaz,et al.  High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems , 2006, Math. Comput..

[9]  Carlos Parés Madroñal,et al.  Numerical methods for nonconservative hyperbolic systems: a theoretical framework , 2006, SIAM J. Numer. Anal..

[10]  Z. Wang,et al.  Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids , 2002 .

[11]  Smadar Karni Viscous shock profiles and primitive formulations , 1992 .

[12]  Chi-Wang Shu,et al.  A technique of treating negative weights in WENO schemes , 2000 .

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

[14]  Eleuterio F. Toro,et al.  HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow , 2010, J. Comput. Phys..

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

[16]  E. Toro,et al.  A path-conservative method for a five-equation model of two-phase flow with an HLLC-type Riemann solver , 2011 .

[17]  Jiaquan Gao,et al.  How to prevent pressure oscillations in multicomponent flow calculations , 2000, Proceedings Fourth International Conference/Exhibition on High Performance Computing in the Asia-Pacific Region.

[18]  Andrew B Wardlaw,et al.  Underwater Explosion Test Cases , 1998 .

[19]  J. Smoller Shock Waves and Reaction-Diffusion Equations , 1983 .

[20]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

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

[22]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[23]  Chi-Wang Shu,et al.  High order finite difference methods with subcell resolution for advection equations with stiff source terms , 2012, J. Comput. Phys..

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

[25]  Chi-Wang Shu,et al.  High Order Strong Stability Preserving Time Discretizations , 2009, J. Sci. Comput..

[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]  Khoon Seng Yeo,et al.  The simulation of compressible multi-medium flow. I. A new methodology with test applications to 1D gas–gas and gas–water cases , 2001 .

[29]  Michael Dumbser,et al.  A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems , 2011, J. Sci. Comput..

[30]  Tiegang Liu,et al.  Runge-Kutta discontinuous Galerkin methods for compressible two-medium flow simulations: One-dimensional case , 2007, J. Comput. Phys..

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

[32]  Boo Cheong Khoo,et al.  The ghost fluid method for compressible gas-water simulation , 2005 .

[33]  Manuel Jesús Castro Díaz,et al.  Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes , 2008, J. Comput. Phys..

[34]  G. D. Maso,et al.  Definition and weak stability of nonconservative products , 1995 .

[35]  Ronald Fedkiw,et al.  An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows , 1999 .

[36]  Boo Cheong Khoo,et al.  The simulation of compressible multi-medium flow: II. Applications to 2D underwater shock refraction , 2001 .

[37]  Chi-Wang Shu,et al.  Fast Sweeping Fifth Order WENO Scheme for Static Hamilton-Jacobi Equations with Accurate Boundary Treatment , 2010, J. Sci. Comput..

[38]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[39]  R. LeVeque Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[40]  Barry Koren,et al.  Modelling, Analysis and Simulation Mas Modelling, Analysis and Simulation a Pressure-invariant Conservative Godunov-type Method for Barotropic Two-fluid Flows , 2022 .

[41]  Chi-Wang Shu,et al.  An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method , 2010, J. Comput. Phys..

[42]  Boo Cheong Khoo,et al.  Ghost fluid method for strong shock impacting on material interface , 2003 .

[43]  Boo Cheong Khoo,et al.  A Real Ghost Fluid Method for the Simulation of Multimedium Compressible Flow , 2006, SIAM J. Sci. Comput..

[44]  Manuel Jesús Castro Díaz,et al.  High Order Extensions of Roe Schemes for Two-Dimensional Nonconservative Hyperbolic Systems , 2009, J. Sci. Comput..

[45]  S. Osher,et al.  Computing interface motion in compressible gas dynamics , 1992 .

[46]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[47]  Smadar Karni,et al.  Multicomponent Flow Calculations by a Consistent Primitive Algorithm , 1994 .

[48]  Jianxian Qiu,et al.  RKDG methods with WENO type limiters and conservative interfacial procedure for one-dimensional compressible multi-medium flow simulations , 2011 .

[49]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[50]  A. Harten ENO schemes with subcell resolution , 1989 .

[51]  T. Hou,et al.  Why nonconservative schemes converge to wrong solutions: error analysis , 1994 .

[52]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

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

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

[55]  Rémi Abgrall,et al.  A comment on the computation of non-conservative products , 2010, J. Comput. Phys..

[56]  Chi-Wang Shu,et al.  A Conservative Modification to the Ghost Fluid Method for Compressible Multiphase Flows , 2011 .

[57]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .