A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows

TL;DR

This paper introduces a robust, genuinely two-dimensional HLLC Riemann solver for conservation laws, applicable to Euler and MHD flows, with a second-order Godunov scheme for 3D applications, enhancing stability and accuracy.

Contribution

It develops a new two-dimensional HLLC Riemann solver with sub-structure, improving stability and applicability to complex flows, and integrates it into a second-order 3D Godunov scheme.

Findings

Achieves stable multi-dimensional flux calculations

Enables accurate modeling of Euler and MHD flows

Cost-competitive with traditional schemes

Abstract

In this paper we present a genuinely two-dimensional HLLC Riemann solver. On logically rectangular meshes, it accepts four input states that come together at an edge and outputs the multi-dimensionally upwinded fluxes in both directions. This work builds on, and improves, our prior work on two-dimensional HLL Riemann solvers. The HLL Riemann solver presented here achieves its stabilization by introducing a constant state in the region of strong interaction, where four one-dimensional Riemann problems interact vigorously with one another. A robust version of the HLL Riemann solver is presented here along with a strategy for introducing sub-structure in the strongly-interacting state. Introducing sub-structure turns the two-dimensional HLL Riemann solver into a two-dimensional HLLC Riemann solver. The sub-structure that we introduce represents a contact discontinuity which can be oriented in any direction relative to the mesh. The Riemann solver presented here is general and can work with any system of conservation laws. We also present a second order accurate Godunov scheme that works in three dimensions and is entirely based on the present multidimensional HLLC Riemann solver technology. The methods presented are cost-competitive with traditional higher order Godunov schemes.