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.