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.

Classification: 74S10, 35L60, 74G15

Keywords: two-layer shallow water, nonconservative system, complex eigenvalues, time-splitting, entropy inequality, well-balanced scheme, nonnegativity

@article{M2AN_2008__42_4_683_0, author = {Bouchut, Fran\c cois and Tom\'as Morales de Luna}, title = {An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique}, publisher = {EDP-Sciences}, volume = {42}, number = {4}, year = {2008}, pages = {683-698}, doi = {10.1051/m2an:2008019}, zbl = {pre05318506}, mrnumber = {2437779}, language = {en}, url = {http://www.numdam.org/item/M2AN_2008__42_4_683_0} }

Bouchut, François; Tomás Morales de Luna. An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Volume 42 (2008) no. 4, pp. 683-698. doi : 10.1051/m2an:2008019. http://www.numdam.org/item/M2AN_2008__42_4_683_0/

[1] Computations of compressible multifluids. J. Comput. Phys. 169 (2001) 594-623. | MR 1836526 | Zbl 1033.76029

and ,[2] A relaxation scheme for the two-layer shallow water system, in Proceedings of the 11th International Conference on Hyperbolic Problems (Lyon, 2006), Hyperbolic problems: theory, numerics, applications, S. Benzoni-Gavage and D. Serre Eds., Springer (2007) 135-144. | Zbl pre05258327

and ,[3] A well-balanced positivity preserving “second-order” scheme for shallow water flows on unstructured meshes. J. Comput. Phys. 206 (2005) 311-333. | MR 2135839 | Zbl 1087.76072

and ,[4] Kinetic schemes for Saint-Venant equations with source terms on unstructured grids. INRIA Report, RR-3989 (2000).

, and ,[5] A fast and stable well-balanced scheme with hydrostatic reconstruction for shallow water flows. SIAM J. Sci. Comput. 25 (2004) 2050-2065 (electronic). | MR 2086830 | Zbl 1133.65308

, , , and ,[6] A wave propagation method for conservation laws and balance laws with spatially varying flux functions. SIAM J. Sci. Comput. 24 (2002) 955-978 (electronic) | MR 1950520 | Zbl 1034.65068

, , and ,[7] A relaxation method for two-phase flow models with hydrodynamic closure law. Numer. Math. 99 (2005) 411-440. | MR 2117734 | Zbl pre02146841

, , , and ,[8] Travelling wave solutions of a convective diffusive system with first and second order terms in nonconservation form, in Hyperbolic problems: theory, numerics, applications, Vol. I (Zürich, 1998), Internat. Ser. Numer. Math. 129, Birkhäuser, Basel (1999) 74-54. | MR 1715732 | Zbl 0934.35030

and ,[9] Nonlinear stability of finite volume methods for hyperbolic conservation laws and well-balanced schemes for sources, Frontiers in Mathematics. Birkhäuser Verlag, Basel (2004). | MR 2128209 | Zbl 1086.65091

,[10] Nonlinear dynamics of rotating shallow water: methods and advances, Edited Series on Advances in Nonlinear Science and Complexity. Elsevier (2007).

, , , and ,[11] A $Q$-scheme for a class of systems of coupled conservation laws with source term. Application to a two-layer 1-D shallow water system. ESAIM: M2AN 35 (2001) 107-127. | Numdam | MR 1811983 | Zbl 1094.76046

, and ,[12] Ideal shocks in a 2-layer flow. II: Under a passive layer. Tellus 53A (2001) 146-167.

and ,[13] On the propagation of internal bores. J. Fluid Mech. 331 (1997) 81-106. | Zbl 0952.76078

, and ,[14] A note on hydraulic theory of internal bores. Dyn. Atm. Oceans 28 (1998) 1-7.

and ,[15] On the well-balance property of Roe's method for nonconservative hyperbolic systems. Applications to shallow-water systems. ESAIM: M2AN 38 (2004) 821-852. | Numdam | MR 2104431 | Zbl 1130.76325

and ,[16] Numerical modeling of two-phase gravitational granular flows with bottom topography, in Proc. of HYP06, Lyon, France (2007). | Zbl pre05258395

, , and ,[17] A kinetic scheme for the Saint-Venant system with a source term. Calcolo 38 (2001) 201-231. | MR 1890353 | Zbl 1008.65066

and ,[18] Theoretical considerations on the motion of salt and fresh water, in Proc. of the Minn. Int. Hydraulics Conv., Joint meeting IAHR and Hyd. Div. ASCE (1953) 321-333.

and ,