Multidimensional HLLC Riemann solver for unstructured meshes - With application to Euler and MHD flows

The goal of this paper is to formulate genuinely multidimensional HLL and HLLC Riemann solvers for unstructured meshes by extending our prior papers on the same topic for logically rectangular meshes Balsara (2010, 2012) [4,5]. Such Riemann solvers operate at each vertex of a mesh and accept as an input the set of states that come together at that vertex. The mesh geometry around that vertex is also one of the inputs of the Riemann solver. The outputs are the resolved state and multidimensionally upwinded fluxes in both directions. A formulation which respects the detailed geometry of the unstructured mesh is presented. Closed-form expressions are provided for all the integrals, making it particularly easy to implement the present multidimensional Riemann solvers in existing numerical codes. While it is visually demonstrated for three states coming together at a vertex, our formulation is general enough to treat multiple states (or zones with arbitrary geometry) coming together at a vertex. The present formulation is very useful for two-dimensional and three-dimensional unstructured mesh calculations of conservation laws. It has been demonstrated to work with second to fourth order finite volume schemes on two-dimensional unstructured meshes. On general triangular grids an arbitrary number of states might come together at a vertex of the primal mesh, while for calculations on the dual mesh usually three states come together at a grid vertex. We apply the multidimensional Riemann solvers to hydrodynamics and magnetohydrodynamics (MHD) on unstructured meshes. The Riemann solver is shown to operate well for traditional second order accurate total variation diminishing (TVD) schemes as well as for weighted essentially non-oscillatory (WENO) schemes with ADER (Arbitrary DERivatives in space and time) time-stepping. Several stringent applications for compressible gasdynamics and magnetohydrodynamics are presented, showing that the method performs very well and reaches high order of accuracy in both space and time. The present multidimensional Riemann solver is cost-competitive with traditional, one-dimensional Riemann solvers. It offers the twin advantages of isotropic propagation of flow features and a larger CFL number. Please see http://www.nd.edu/~dbalsara/Numerical-PDE-Course for a video introduction to multidimensional Riemann solvers.

[1]  Rémi Abgrall,et al.  On essentially non-oscillatory schemes on unstructured meshes: analysis and implementation , 1994 .

[2]  Katharine Gurski,et al.  An HLLC-Type Approximate Riemann Solver for Ideal Magnetohydrodynamics , 2001, SIAM J. Sci. Comput..

[3]  Michael Dumbser,et al.  Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems , 2007, J. Comput. Phys..

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

[5]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[6]  Dinshaw S. Balsara.  A two-dimensional HLLC Riemann solver for conservation laws: Application to Euler and magnetohydrodynamic flows , 2012, J. Comput. Phys..

[7]  R. LeVeque.  Wave Propagation Algorithms for Multidimensional Hyperbolic Systems , 1997 .

[8]  Dinshaw Balsara,et al.  Divergence-free adaptive mesh refinement for Magnetohydrodynamics , 2001 .

[9]  Gerald Warnecke,et al.  Finite volume evolution Galerkin methods for Euler equations of gas dynamics , 2002 .

[10]  E. Toro,et al.  Solution of the generalized Riemann problem for advection–reaction equations , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[11]  A. Scott,et al.  Solitons and the Inverse Scattering Transform (Mark J. Ablowitz and Harvey Segur) , 1983 .

[12]  James P. Collins,et al.  Numerical Solution of the Riemann Problem for Two-Dimensional Gas Dynamics , 1993, SIAM J. Sci. Comput..

[13]  Michael Fey,et al.  Multidimensional Upwinding. Part II. Decomposition of the Euler Equations into Advection Equations , 1998 .

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

[15]  Rémi Abgrall,et al.  A Cell-Centered Lagrangian Scheme for Two-Dimensional Compressible Flow Problems , 2007, SIAM J. Sci. Comput..

[16]  Shengtai Li.  An HLLC Riemann solver for magneto-hydrodynamics , 2005 .

[17]  E. F. Toro,et al.  The development of a Riemann solver for the steady supersonic Euler equations , 1994, The Aeronautical Journal (1968).

[18]  Michael Dumbser,et al.  Efficient implementation of ADER schemes for Euler and magnetohydrodynamical flows on structured meshes - Speed comparisons with Runge-Kutta methods , 2013, J. Comput. Phys..

[19]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[20]  Dinshaw S. Balsara.  Multidimensional HLLE Riemann solver: Application to Euler and magnetohydrodynamic flows , 2010, J. Comput. Phys..

[21]  P. Roe,et al.  On Godunov-type methods near low densities , 1991 .

[22]  Burton Wendroff,et al.  A two-dimensional HLLE riemann solver and associated godunov-type difference scheme for gas dynamics☆ , 1999 .

[23]  Dinshaw Balsara,et al.  Second-Order-accurate Schemes for Magnetohydrodynamics with Divergence-free Reconstruction , 2003, astro-ph/0308249.

[24]  Dinshaw S. Balsara,et al.  Linearized Formulation of the Riemann Problem for Adiabatic and Isothermal Magnetohydrodynamics , 1998 .

[25]  Michael Fey,et al.  Multidimensional Upwinding. Part I. The Method of Transport for Solving the Euler Equations , 1998 .

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

[27]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[28]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[29]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[30]  Dinshaw S. Balsara.  Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics , 2009, J. Comput. Phys..

[31]  Bernd Einfeld.  On Godunov-type methods for gas dynamics , 1988 .

[32]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[33]  C. Munz,et al.  Hyperbolic divergence cleaning for the MHD equations , 2002 .

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

[35]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[36]  Rémi Abgrall.  APPROXIMATION DU PROBLEME DE RIEMANN VRAIMENT MULTIDIMENSIONNEL DES EQUATIONS D'EULER PAR UNE METHODE DE TYPE ROE (I) : LA LINEARISATION , 1994 .

[37]  J. Saltzman,et al.  An unsplit 3D upwind method for hyperbolic conservation laws , 1994 .

[38]  P. Colella.  A Direct Eulerian MUSCL Scheme for Gas Dynamics , 1985 .

[39]  Dinshaw S. Balsara,et al.  Notes on the Eigensystem of Magnetohydrodynamics , 1996, SIAM J. Appl. Math..

[40]  Philip L. Roe,et al.  Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics , 1986 .

[41]  D. Balsara,et al.  A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations , 1999 .

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

[43]  Aramais R. Zakharian,et al.  Two-dimensional Riemann solver for Euler equations of gas dynamics , 2001 .

[44]  S. Orszag,et al.  Small-scale structure of two-dimensional magnetohydrodynamic turbulence , 1979, Journal of Fluid Mechanics.

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

[46]  Gérard Gallice,et al.  Roe Matrices for Ideal MHD and Systematic Construction of Roe Matrices for Systems of Conservation Laws , 1997 .

[47]  Eleuterio F. Toro,et al.  AOn WAF-Type Schemes for Multidimensional Hyperbolic Conservation Laws , 1997 .

[48]  Michael Dumbser,et al.  Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics , 2008, Journal of Computational Physics.

[49]  Dinshaw S. Balsara,et al.  Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics , 2012, J. Comput. Phys..

[50]  Alexandre J. Chorin,et al.  Random choice solution of hyperbolic systems , 1976 .

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

[52]  Eleuterio F. Toro,et al.  ADER schemes for three-dimensional non-linear hyperbolic systems , 2005 .

[53]  Eleuterio F. Toro,et al.  ADER: Arbitrary High Order Godunov Approach , 2002, J. Sci. Comput..

[54]  K. Kusano,et al.  A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics , 2005 .

[55]  Philip L. Roe,et al.  A multidimensional flux function with applications to the Euler and Navier-Stokes equations , 1993 .