A Very Easy High-Order Well-Balanced Reconstruction for Hyperbolic Systems with Source Terms

Abstract

When adopting high-order finite volume schemes based on MUSCL reconstruction techniques to approximate the weak solutions of hyperbolic systems with source terms, the preservation of the steady states turns out to be very challenging. Indeed, the designed reconstruction must preserve the steady states under consideration in order to get the required well-balancedness property. A priori, to capture such a steady state, one needs to solve some strongly nonlinear equations. Here, we design a very easy correction to high-order finite volume methods. This correction can be applied to any scheme of order greater than or equal to 2, such as a MUSCL-type scheme, and ensures that this scheme exactly preserves the steady solutions. In contrast to usual techniques, our scheme avoids the inversion of the nonlinear function that governs the steady solutions. Moreover, for underdetermined steady solutions, several nonlinear functions must be considered simultaneously. Since the derived correction only considers the evaluation of the governing nonlinear functions, we are able to deal with underdetermined stationary systems. Several numerical experiments illustrate the relevance of the proposed well-balanced correction, as well as its main limitation, namely the fact that it may fail at being both well-balanced and more than second-order accurate for a specific class of initial conditions.