An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment

¹ Département de Mathématiques et Applications, École Normale Supérieure & CNRS, 45 rue d'Ulm, 75230 Paris cedex 05, France. francois.bouchut@ens.fr
² Departamento de Anlálisis Matemático, Facultad de Ciencias, Universidad de Málaga, Campus de Teatinons s/n, 29071 Málaga, Spain. morales@anamat.cie.uma.es

Received: 5 March 2007
Revised: 11 December 2007

Abstract

We consider the system of partial differential equations governing the one-dimensional flow of two superposed immiscible layers of shallow water. The difficulty in this system comes from the coupling terms involving some derivatives of the unknowns that make the system nonconservative, and eventually nonhyperbolic. Due to these terms, a numerical scheme obtained by performing an arbitrary scheme to each layer, and using time-splitting or other similar techniques leads to instabilities in general. Here we use entropy inequalities in order to control the stability. We introduce a stable well-balanced time-splitting scheme for the two-layer shallow water system that satisfies a fully discrete entropy inequality. In contrast with Roe type solvers, it does not need the computation of eigenvalues, which is not simple for the two-layer shallow water system. The solver has the property to keep the water heights nonnegative, and to be able to treat vanishing values.