Rusanov-type Schemes for Hyperbolic Equations: Wave-Speed Estimates, Monotonicity and Stability

TL;DR

This paper investigates how errors in wave speed estimates affect the monotonicity and stability of Rusanov-type schemes for hyperbolic equations, highlighting that overestimation is generally safer than underestimation.

Contribution

It provides a theoretical analysis of wave speed estimate errors on Rusanov schemes, revealing the impacts on stability and monotonicity, and offers guidance for more reliable wave speed estimation.

Findings

Underestimating wave speeds can cause loss of monotonicity and stability.

Overestimating wave speeds preserves monotonicity and enhances stability.

Overestimation is preferable to underestimation for scheme robustness.

Abstract

HLL-type schemes constitute a large hierarchy of numerical methods, in the finite volume and discontinuous Galerkin finite element frameworks, for solving hyperbolic equations. The hierarchy of fluxes includes Rusanov schemes, HLL schemes, HLLC schemes, and other variations. All of these schemes rely on wave speed estimates. Recent work has shown that most wave speed estimates in current use underestimate the true wave speeds. In the present paper we carry out a theoretical study of the consequences arising from errors in the wave speed estimates, on the monotonicity and stability properties of the derived schemes. For the simplest case of the hierarchy, that is Rusanov-type schemes, we carry out a detailed analysis in terms of the linear advection equation in one and two space dimensions. It is found that errors from underestimates of the wave speed could cause loss of monotonicity, severe reduction of the stability limit, and even loss of stability. Errors from overestimates, though preserving monotonicity, will cause a reduction of the stability limit. We find that overestimation is preferable to underestimation, for two reasons. First, schemes from overestimation are monotone, and second, their stability regions are larger than those from underestimation, for equivalent displacements from the exact speed. The findings of this paper may prove useful in raising awareness of the potential pitfalls of a seemingly simple practical computational task, that of providing wave speed estimates. Our reported findings may also motivate subsequent studies for complex non-linear hyperbolic systems, requiring estimates for two or more waves, such as in HLL and HLLC schemes.