Balancing source terms and flux gradients in high-resolution Godunov methods: The quasi-steady wave-propagation algorithm (Q1275177)

We develop a variant of the wave-propagation algorithm which addresses the title problem by introducing a Riemann problem at the center of each grid cell whose flux difference exactly cancels the source term. This leads to modified Riemann problems at the cell edges in which the jump now corresponds to perturbations from the steady state. Computing waves and limiters based on the solution to these Riemann problems give high-resolution results. The one- and two-dimensional shallow water equations for flow over arbitrary bottom topography are used as an example, through the ideas apply to many other systems. The method is easily implemented in the software package CLAWPACK. \(\copyright\) Academic Press.