An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment. Bouchut, F. & Morales de Luna, T. ESAIM: Mathematical Modelling and Numerical Analysis (ESAIM: M2AN), 2008.
An entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment [link]Paper  abstract   bibtex   
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.
@Article{bouchut08entropy,
  author   = {Bouchut, Fran{\c{c}}ois and Morales de Luna, Tom{\'a}s},
  title    = {{A}n entropy satisfying scheme for two-layer shallow water equations with uncoupled treatment},
  journal  = {ESAIM: Mathematical Modelling and Numerical Analysis (ESAIM: M2AN)},
  year     = {2008},
  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.},
  keywords = {Two-layer shallow water, nonconservative system, complex eigenvalues, time-splitting, entropy inequality, well-balanced scheme, nonnegativity},
  url      = {http://dx.doi.org/10.1051/m2an:2008019},
}

Downloads: 0