A Residual-Based Shock Capturing Scheme for the Continuous/Discontinuous Spectral Element Solution of the 2D Shallow Water Equations

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.

[1]  Francis X. Giraldo,et al.  Analysis of adaptive mesh refinement for IMEX discontinuous Galerkin solutions of the compressible Euler equations with application to atmospheric simulations , 2014, J. Comput. Phys..

[2]  Timothy C. Warburton,et al.  A GPU-accelerated continuous and discontinuous Galerkin non-hydrostatic atmospheric model , 2019, Int. J. High Perform. Comput. Appl..

[3]  Mark Ainsworth,et al.  Dispersive and Dissipative Properties of Discontinuous Galerkin Finite Element Methods for the Second-Order Wave Equation , 2006, J. Sci. Comput..

[4]  Clint Dawson,et al.  A Discontinuous Galerkin Method for Three-Dimensional Shallow Water Equations , 2005, J. Sci. Comput..

[5]  Francis X. Giraldo,et al.  High‐order semi‐implicit time‐integrators for a triangular discontinuous Galerkin oceanic shallow water model , 2009 .

[6]  Jean-Luc Guermond,et al.  Stabilized Spectral Element Approximation of the Saint Venant System Using the Entropy Viscosity Technique , 2015 .

[7]  Ethan J. Kubatko,et al.  hp Discontinuous Galerkin methods for advection dominated problems in shallow water flow , 2006 .

[8]  Francis X. Giraldo,et al.  Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES , 2015, J. Comput. Phys..

[9]  Claes Eskilsson,et al.  An hp‐adaptive discontinuous Galerkin method for shallow water flows , 2011 .

[10]  Johan Hoffman,et al.  Residual‐based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods , 2013 .

[11]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[12]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[13]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[14]  Manuel Jesús Castro Díaz,et al.  On a well-balanced high-order finite volume scheme for shallow water equations with topography and dry areas , 2007, J. Comput. Phys..

[15]  Francis X. Giraldo,et al.  Simulation of Shallow Water Jets with a Unified Element-based Continuous/Discontinuous Galerkin Model with Grid Flexibility on the Sphere , 2015 .

[16]  Natalja Rakowsky,et al.  Operational tsunami modelling with TsunAWI – recent developments and applications , 2013 .

[17]  Emil M. Constantinescu,et al.  Acceleration of the IMplicit–EXplicit nonhydrostatic unified model of the atmosphere on manycore processors , 2017, Int. J. High Perform. Comput. Appl..

[18]  E. Hendricks,et al.  Evaluation of the Utility of Static and Adaptive Mesh Refinement for Idealized Tropical Cyclone Problems in a Spectral Element Shallow-Water Model , 2015 .

[19]  Qiuhua Liang,et al.  Dynamically adaptive grid based discontinuous Galerkin shallow water model , 2012 .

[20]  J. Kelly,et al.  Linear and nonlinear ultrasound simulations using the discontinuous Galerkin method. , 2018, The Journal of the Acoustical Society of America.

[21]  S. Soares-Frazão,et al.  Simulation of flow in compound open-channel using a discontinuous Galerkin finite-element method with Smagorinsky turbulence closure , 2014 .

[22]  M. Carpenter,et al.  Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations , 2003 .

[23]  Olivier Delestre,et al.  SWASHES: a compilation of shallow water analytic solutions for hydraulic and environmental studies , 2011, 1110.0288.

[24]  Tuomas Kärnä,et al.  A flux-limiting wetting–drying method for finite-element shallow-water models, with application to the Scheldt Estuary , 2009 .

[25]  V. Legat,et al.  A fully implicit wetting-drying method for DG-FEM shallow water models, with an application to the Scheldt Estuary , 2011 .

[26]  Emil M. Constantinescu,et al.  Implicit-Explicit Formulations of a Three-Dimensional Nonhydrostatic Unified Model of the Atmosphere (NUMA) , 2013, SIAM J. Sci. Comput..

[27]  Jan S. Hesthaven,et al.  Nodal high-order discontinuous Galerkin methods for the spherical shallow water equations , 2002 .

[28]  M. Taylor.  The Spectral Element Method for the Shallow Water Equations on the Sphere , 1997 .

[29]  Michael J. Briggs,et al.  Laboratory experiments of tsunami runup on a circular island , 1995 .

[30]  Yulong Xing,et al.  Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations , 2010 .

[31]  Qiuhua Liang,et al.  A conservative high‐order discontinuous Galerkin method for the shallow water equations with arbitrary topography , 2011 .

[32]  Georges Kesserwani,et al.  Multiwavelet-based grid adaptation with discontinuous Galerkin schemes for shallow water equations , 2015, J. Comput. Phys..

[33]  Claes Eskilsson,et al.  Method of moving frames to solve the shallow water equations on arbitrary rotating curved surfaces , 2017, J. Comput. Phys..

[34]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[35]  Francis X. Giraldo,et al.  A spectral element shallow water model on spherical geodesic grids , 2001 .

[36]  Francis X. Giraldo,et al.  Continuous and discontinuous Galerkin methods for a scalable three-dimensional nonhydrostatic atmospheric model: Limited-area mode , 2012, J. Comput. Phys..

[37]  J. J. Stoker.  Water Waves: The Mathematical Theory with Applications , 1957 .

[38]  Qiuhua Liang,et al.  Locally Limited and Fully Conserved RKDG2 Shallow Water Solutions with Wetting and Drying , 2012, J. Sci. Comput..

[39]  Stephen J. Thomas,et al.  A Discontinuous Galerkin Global Shallow Water Model , 2005, Monthly Weather Review.

[40]  Francis X. Giraldo,et al.  Mass conservation of the unified continuous and discontinuous element-based Galerkin methods on dynamically adaptive grids with application to atmospheric simulations , 2015, J. Comput. Phys..

[41]  Francis X. Giraldo,et al.  Efficient construction of unified continuous and discontinuous Galerkin formulations for the 3D Euler equations , 2016, J. Comput. Phys..

[42]  John C. Butcher,et al.  The effective order of singly-implicit Runge-Kutta methods , 1999, Numerical Algorithms.

[43]  Hervé Vandeven,et al.  Family of spectral filters for discontinuous problems , 1991 .

[44]  Georges Kesserwani,et al.  Benchmarking a multiresolution discontinuous Galerkin shallow water model: Implications for computational hydraulics , 2015 .

[45]  W. Thacker.  Some exact solutions to the nonlinear shallow-water wave equations , 1981, Journal of Fluid Mechanics.

[46]  Paul A. Ullrich,et al.  The spectral element method (SEM) on variable-resolution grids: evaluating grid sensitivity and resolution-aware numerical viscosity , 2014 .

[47]  D. Lilly.  On the numerical simulation of buoyant convection , 1962 .

[48]  Guido Boffetta,et al.  Two-Dimensional Turbulence , 2012 .

[49]  Maojun Li,et al.  A Positivity-Preserving Well-Balanced Central Discontinuous Galerkin Method for the Nonlinear Shallow Water Equations , 2017, J. Sci. Comput..

[50]  J. Boyd,et al.  A staggered spectral element model with application to the oceanic shallow , 1995 .

[51]  E. Toro.  Shock-Capturing Methods for Free-Surface Shallow Flows , 2001 .

[52]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[53]  Costas E. Synolakis,et al.  The runup of solitary waves , 1987, Journal of Fluid Mechanics.

[54]  Francis X. Giraldo,et al.  Strong scaling for numerical weather prediction at petascale with the atmospheric model NUMA , 2015, Int. J. High Perform. Comput. Appl..

[55]  Hong Ma,et al.  A spectral element basin model for the shallow water equations , 1993 .

[56]  S. R. Brus,et al.  A Comparison of Artificial Viscosity, Limiters, and Filters, for High Order Discontinuous Galerkin Solutions in Nonlinear Settings , 2016, J. Sci. Comput..

[57]  Jean-Luc Guermond,et al.  Viscous Regularization of the Euler Equations and Entropy Principles , 2012, SIAM J. Appl. Math..

[58]  R Gandham,et al.  GPU Accelerated Discontinuous Galerkin Methods for Shallow Water Equations , 2014, 1403.1661.

[59]  M. Vázquez-Cendón.  Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .