On general and complete multidimensional Riemann solvers for nonlinear systems of hyperbolic conservation laws

In this work, we introduce a framework to design multidimensional Riemann solvers for nonlinear systems of hyperbolic conservation laws on general unstructured polygonal Voronoi-like tessellations. In this framework we propose two simple but complete solvers. The first method is a direct extension of the Osher-Solomon Riemann solver to multiple space dimensions. Here, the multidimensional numerical dissipation is obtained by integrating the absolute value of the flux Jacobians over a dual triangular mesh around each node of the primal polygonal grid. The required nodal gradient is then evaluated on a local computational simplex involving the d + 1 neighbors meeting at each corner. The second method is a genuinely multidimensional upwind flux. By introducing a fluctuation form of finite volume methods with corner fluxes, we show an equivalence with residual distribution schemes (RD). This naturally allows to construct genuinely multidimensional upwind corner fluxes starting from existing and well-known RD fluctuations. In particular, we explore the use of corner fluxes obtained from the so-called N scheme, i.e. the Roe's original optimal multidimensional upwind advection scheme. Both methods use the full eigenstructure of the underlying hyperbolic system and are therefore complete by construction. A simple higher order extension up to fourth order in space and time is then introduced at the aid of a CWENO reconstruction in space and an ADER approach in time, leading to a family of high order accurate fully-discrete one-step schemes based on genuinely multidimensional Riemann solvers. We present applications of our new numerical schemes to several test problems for the compressible Euler equations of gas-dyanamics. In particular, we show that the proposed schemes are at the same time carbuncle-free and preserve certain stationary shear waves exactly.