Behaviour of the numerical entropy production of the one-and-a-half-dimensional shallow water equations

Authors

  • Sudi Mungkasi Australian National University
  • Stephen Gwyn Roberts Australian National University

DOI:

https://doi.org/10.21914/anziamj.v54i0.6243

Keywords:

numerical entropy production, smoothness indicator, refinement indicator, finite volume methods, shallow water equations, passive tracer, transverse velocity

Abstract

This article reports the behaviour of the numerical entropy production of the one-and-a-half-dimensional shallow water equations. The one-and-a-half-dimensional shallow water equations are the one-dimensional shallow water equations with a passive tracer or transverse velocity. The studied behaviour is with respect to the choice of numerical fluxes to evolve the mass, momentum, tracer-mass (transverse momentum), and entropy. When solving the one-and-a-half-dimensional shallow water equations using a finite volume method, we recommend the use of a double sided stencil flux for the mass and momentum, and in addition, a single sided stencil (upwind) flux for the tracer-mass. Having this recommended combination of fluxes, we use a double sided stencil entropy flux to compute the numerical entropy production, but this flux generates positive overshoots of the numerical entropy production. Positive overshoots of the numerical entropy production are avoided by use of a modified entropy flux, which satisfies a discrete numerical entropy inequality. References

Author Biographies

Sudi Mungkasi, Australian National University

Mathematical Sciences Institute, PhD Candidate

Stephen Gwyn Roberts, Australian National University

Mathematical Sciences Institute, Associate Professor

Published

2013-05-05

Issue

Section

Proceedings Computational Techniques and Applications Conference