A residual-based shock capturing scheme for the continuous/discontinuous spectral element solution of the 2D shallow water equations. (April 2018)

• Record Type: Journal Article
• Authors: Marras, Simone; Kopera, Michal A.; Constantinescu, Emil M.; Suckale, Jenny; Giraldo, Francis X.
• Is Part Of: Advances in water resources. Volume 114(2018)
• Page Start: 45
• Page End: 63
• Publication Date: 2018-04
• Subjects: Shallow water equations -- Element-based Galerkin methods -- Dynamic artificial diffusion -- De-aliasing -- High-order methods -- Unified continuous/discontinuous Galerkin -- Large Eddy Simulation; Hydrology -- Periodicals; Hydrodynamics -- Periodicals; Hydraulic engineering -- Periodicals; 551.48
• DOI: 10.1016/j.advwatres.2018.02.003
• Languages: English
• ISSNs: 0309-1708

Abstract

Highlights: Novel dynamic artificial viscosity for the de-aliasing of the high order unified CG/DG solutions of the shallow water equations. The proposed diffusion model improves the continuous and discontinuous Galerkin solutions alike in the proximity of sharp fronts propagating on wet surface. The proposed model is mass conservative. 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.