Path-conservative central-upwind schemes for nonconservative hyperbolic systems

¹ Departamento de Análisis Matemático, Universidad de Málaga, 29080 Málaga, Spain
² Department of Mathematics, Southern University of Science and Technology, 518055 Shenzhen, China
³ Mathematics Department, Tulane University, New Orleans, LA 70118, USA
⁴ Departamento de Matemáticas, Universidad de Córdoba, Campus de Rabanales, 14071 Córdoba, Spain

^* Corresponding author: kurganov@math.tulane.edu

Received: 10 April 2018
Accepted: 20 December 2018

Abstract

We develop path-conservative central-upwind schemes for nonconservative one-dimensional hyperbolic systems of nonlinear partial differential equations. Such systems arise in a variety of applications and the most challenging part of their numerical discretization is a robust treatment of nonconservative product terms. Godunov-type central-upwind schemes were developed as an efficient, highly accurate and robust ``black-box’’ solver for hyperbolic systems of conservation and balance laws. They were successfully applied to a large number of hyperbolic systems including several nonconservative ones. To overcome the difficulties related to the presence of nonconservative product terms, several special techniques were proposed. However, none of these techniques was sufficiently robust and thus the applicability of the original central-upwind schemes was rather limited. In this paper, we rewrite the central-upwind schemes in the form of path-conservative schemes. This helps us (i) to show that the main drawback of the original central-upwind approach was the fact that the jump of the nonconservative product terms across cell interfaces has never been taken into account and (ii) to understand how the nonconservative products should be discretized so that their influence on the numerical solution is accurately taken into account. The resulting path-conservative central-upwind scheme is a new robust tool for both conservative and nonconservative hyperbolic systems. We apply the new scheme to the Saint-Venant system with discontinuous bottom topography and two-layer shallow water system. Our numerical results illustrate the good performance of the new path-conservative central-upwind scheme, its robustness and ability to achieve very high resolution.