A HLLC scheme for nonconservative hyperbolic problems. Application to turbidity currents with sediment transport

¹ Dpto. de Análisis Matemático, Facultad de Ciencias, Universidad de Málga, Campus de Teatinos, s/n, 29071 Málaga, Spain
castor@anamat.cie.uma.es; pares@anamat.cie.uma.es
² Dpto. Matemática Aplicada I, ETS Arquitectura, Universidad de Sevilla, Avda. Reina Mercedes No. 2, 41012 Sevilla, Spain
edofer@us.es; gnarbona@us.es
³ Dpto. de Matemáticas, Universidad de Córdoba, Campus de Rabanales, 14071 Córdoba, Spain
tomas.morales@uco.es

Received: 24 June 2011
Revised: 4 April 2012

Abstract

The goal of this paper is to obtain a well-balanced, stable, fast, and robust HLLC-type approximate Riemann solver for a hyperbolic nonconservative PDE system arising in a turbidity current model. The main difficulties come from the nonconservative nature of the system. A general strategy to derive simple approximate Riemann solvers for nonconservative systems is introduced, which is applied to the turbidity current model to obtain two different HLLC solvers. Some results concerning the non-negativity preserving property of the corresponding numerical methods are presented. The numerical results provided by the two HLLC solvers are compared between them and also with those obtained with a Roe-type method in a number of 1d and 2d test problems. This comparison shows that, while the quality of the numerical solutions is comparable, the computational cost of the HLLC solvers is lower, as only some partial information of the eigenstructure of the matrix system is needed.