An hp‐adaptive discontinuous Galerkin method for shallow water flows

An adaptive spectral/hp discontinuous Galerkin method for the two-dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non-conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p−1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h-type refinement, the parent element is subdivided into four similar sibling elements. The time-stepping is performed using a third-order Runge–Kutta scheme. The performance of the hp-adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p-adaptivity is more efficient than h-adaptivity with respect to degrees of freedom and computational time.

[1]  Francesca Rapetti,et al.  Basics and some applications of the mortar element method , 2005 .

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

[3]  Jean-François Remacle,et al.  High-order h-adaptive discontinuous Galerkin methods for ocean modelling , 2007 .

[4]  Brett F. Sanders,et al.  Non-reflecting boundary flux function for finite volume shallow-water models , 2002 .

[5]  Clint Dawson,et al.  A Performance Comparison of Continuous and Discontinuous Finite Element Shallow Water Models , 2009, J. Sci. Comput..

[6]  Francis X. Giraldo,et al.  High-order triangle-based discontinuous Galerkin methods for hyperbolic equations on a rotating sphere , 2006, J. Comput. Phys..

[7]  John P. Boyd,et al.  Equatorial Solitary Waves. Part I: Rossby Solitons , 1980 .

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

[9]  David Jon Furbish,et al.  Numerical Solution of the Dam-Break Problem with a Discontinuous Galerkin Method , 2004 .

[10]  Hong Li,et al.  The discontinuous Galerkin finite element method for the 2D shallow water equations , 2001 .

[11]  M. Berzins,et al.  An unstructured finite-volume algorithm for predicting flow in rivers and estuaries , 1998 .

[12]  Onno Bokhove,et al.  Space discontinuous Galerkin method for shallow water flows : kinetic and HLLC flux, and potential vorticity generation , 2007 .

[13]  Ethan J. Kubatko,et al.  Dynamic p-adaptive Runge–Kutta discontinuous Galerkin methods for the shallow water equations , 2009 .

[14]  Alain Vande Wouwer,et al.  Upwinding in the method of lines , 2001 .

[15]  Francis Loth,et al.  Spectral Element Methods for Transitional Flows in Complex Geometries , 2002, J. Sci. Comput..

[16]  Qiuhua Liang,et al.  Shallow flow simulation on dynamically adaptive cut cell quadtree grids , 2007 .

[17]  Francis X. Giraldo,et al.  A high‐order triangular discontinuous Galerkin oceanic shallow water model , 2008 .

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

[19]  Spencer J. Sherwin,et al.  On discontinuous Galerkin methods , 2003 .

[20]  Jean-François Remacle,et al.  An adaptive discretization of shallow‐water equations based on discontinuous Galerkin methods , 2006 .

[21]  David A. Kopriva,et al.  A Staggered-Grid Multidomain Spectral Method for the Compressible Navier-Stokes Equations , 1998 .

[22]  Spencer J. Sherwin,et al.  A triangular spectral/hp discontinuous Galerkin method for modelling 2D shallow water equations , 2004 .

[23]  J. Oden,et al.  A Posteriori Error Estimation in Finite Element Analysis: Oden/A Posteriori , 2000 .

[24]  Richard L. Stockstill,et al.  Refinement Indicator for Mesh Adaption in Shallow-Water Modeling , 2006 .

[25]  H. Stommel,et al.  The westward intensification of wind‐driven ocean currents , 1948 .

[26]  J. Remacle,et al.  Shock detection and limiting with discontinuous Galerkin methods for hyperbolic conservation laws , 2004 .