Three dimensional HLL Riemann solver for conservation laws on structured meshes; Application to Euler and magnetohydrodynamic flows

In this paper we build on our prior work on multidimensional Riemann solvers by detailing the construction of a three-dimensional HLL Riemann solver. As with the two-dimensional Riemann solver, this is accomplished by introducing a constant resolved state between the states being considered, which introduces sufficient dissipation for systems of conservation laws. Closed form expressions for the resolved fluxes are provided to facilitate numerical implementation. This is accomplished by introducing a novel derivation of the multidimensional Riemann solver. The novelty consists of integrating Lagrangian fluxes across moving surfaces. This makes the problem easier to visualize in three dimensions. (A video introduction to multidimensional Riemann solvers is available on http://www.nd.edu/~dbalsara/Numerical-PDE-Course.)A robust and efficient second order accurate numerical scheme for three dimensional Euler and MHD flows is presented. The scheme is built on the current three-dimensional Riemann solver and has been implemented in the author's RIEMANN code. We demonstrate that schemes that are based on the three-dimensional Riemann solver permit multidimensional discontinuities to propagate more isotropically on resolution-starved meshes. The number of zones updated per second by this scheme on a modern processor is shown to be cost competitive with schemes that are based on a one-dimensional Riemann solver. However, the present scheme permits larger timesteps in three dimensions because of its inclusion of genuinely three-dimensional effects in the flow. For MHD problems it is not necessary to double the dissipation when evaluating the edge-centered electric fields.Accuracy analysis for three-dimensional Euler and MHD problems shows that the scheme meets its design accuracy. Several stringent test problems involving Euler and MHD flows are also presented and the scheme is shown to perform robustly on all of them. For the very first time, we present the formulation and solution of three-dimensional Riemann problems.

[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]  Gerald Warnecke,et al.  Finite volume evolution Galerkin methods for Euler equations of gas dynamics , 2002 .

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

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

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

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

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

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

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

[11]  Dinshaw S. Balsara,et al.  Comparing Coarray Fortran (CAF) with MPI for several structured mesh PDE applications , 2015, J. Comput. Phys..

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

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

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

[15]  Edouard Audit,et al.  A simple two-dimensional extension of the HLL Riemann solver for hyperbolic systems of conservation laws , 2015, J. Comput. Phys..

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

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

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

[19]  Michael Dumbser,et al.  A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws , 2014, J. Comput. Phys..

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

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

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

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

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

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

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

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

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

[29]  Michael Dumbser,et al.  Divergence-free MHD on unstructured meshes using high order finite volume schemes based on multidimensional Riemann solvers , 2015, J. Comput. Phys..

[30]  Michael Dumbser,et al.  Lagrangian ADER-WENO finite volume schemes on unstructured triangular meshes based on genuinely multidimensional HLL Riemann solvers , 2013, J. Comput. Phys..

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

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

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

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

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

[36]  J. Stone,et al.  An unsplit Godunov method for ideal MHD via constrained transport , 2005, astro-ph/0501557.

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

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

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

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

[41]  Michael Dumbser,et al.  Multidimensional Riemann problem with self-similar internal structure. Part II - Application to hyperbolic conservation laws on unstructured meshes , 2015, J. Comput. Phys..

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

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

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

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

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

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

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

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

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

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

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

[53]  Edouard Audit,et al.  A Simple Two-Dimensional Extension of the HLLE Riemann Solver for Gas Dynamics , 2014 .

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

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

[56]  Rémi Abgrall,et al.  Multidimensional HLLC Riemann solver for unstructured meshes - With application to Euler and MHD flows , 2014, J. Comput. Phys..

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

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

[59]  James M. Stone,et al.  An unsplit Godunov method for ideal MHD via constrained transport in three dimensions , 2007, J. Comput. Phys..

[60]  Zhiliang Xu,et al.  Divergence-Free WENO Reconstruction-Based Finite Volume Scheme for Solving Ideal MHD Equations on Triangular Meshes , 2011, 1110.0860.

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

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

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

[64]  Michael Dumbser,et al.  High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics , 2013, 1310.7256.

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