Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
8883331 | Advances in Water Resources | 2018 | 19 Pages |
Abstract
The high-order numerical solution of the non-linear shallow water equations is susceptible to Gibbs oscillations in the proximity of strong gradients. In this paper, we tackle this issue by presenting a shock capturing model based on the numerical residual of the solution. Via numerical tests, we demonstrate that the model removes the spurious oscillations in the proximity of strong wave fronts while preserving their strength. Furthermore, for coarse grids, it prevents energy from building up at small wave-numbers. When applied to the continuity equation to stabilize the water surface, the addition of the shock capturing scheme does not affect mass conservation. We found that our model improves the continuous and discontinuous Galerkin solutions alike in the proximity of sharp fronts propagating on wet surfaces. In the presence of wet/dry interfaces, however, the model needs to be enhanced with the addition of an inundation scheme which, however, we do not address in this paper.
Related Topics
Physical Sciences and Engineering
Earth and Planetary Sciences
Earth-Surface Processes
Authors
Simone Marras, Michal A. Kopera, Emil M. Constantinescu, Jenny Suckale, Francis X. Giraldo,