Well-balanced RKDG2 solutions to the shallow water equations over irregular domains with wetting and drying

Abstract This paper presents a new one-dimensional (1D) second-order Runge–Kutta discontinuous Galerkin (RKDG2) scheme for shallow flow simulations involving wetting and drying over complex domain topography. The shallow water equations that adopt water level (instead of water depth) as a flow variable are solved by an RKDG2 scheme to give piecewise linear approximate solutions, which are locally defined by an average coefficient and a slope coefficient. A wetting and drying technique proposed originally for a finite volume MUSCL scheme is revised and implemented in the RKDG2 solver. Extra numerical enhancements are proposed to amend the local coefficients associated with water level and bed elevation in order to maintain the well-balanced property of the RKDG2 scheme for applications with wetting and drying. Friction source terms are included and evaluated using splitting implicit discretization, implemented with a physical stopping condition to ensure stability. Several steady and unsteady benchmark tests with/without friction effects are considered to demonstrate the performance of the present model.

[1]  G. Kesserwani,et al.  Application of a second‐order Runge–Kutta discontinuous Galerkin scheme for the shallow water equations with source terms , 2008 .

[2]  Mengping Zhang,et al.  An analysis of and a comparison between the discontinuous Galerkin and the spectral finite volume methods , 2005 .

[3]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

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

[5]  Onno Bokhove,et al.  Flooding and Drying in Discontinuous Galerkin Finite-Element Discretizations of Shallow-Water Equations. Part 1: One Dimension , 2005, J. Sci. Comput..

[6]  T. Gallouët,et al.  Some approximate Godunov schemes to compute shallow-water equations with topography , 2003 .

[7]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

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

[9]  Javier Murillo,et al.  Time step restrictions for well‐balanced shallow water solutions in non‐zero velocity steady states , 2009 .

[10]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

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

[12]  Qiuhua Liang,et al.  Adaptive quadtree simulation of shallow flows with wet-dry fronts over complex topography , 2009 .

[13]  G. Kesserwani,et al.  Well‐balancing issues related to the RKDG2 scheme for the shallow water equations , 2009 .

[14]  P. Sweby.  High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[15]  Dirk Schwanenberg,et al.  Die Runge-Kutta-Discontinuous-Galerkin-Methode zur Lösung konvektionsdominierter tiefengemittelter Flachwasserprobleme , 2003 .

[16]  Javier Murillo,et al.  Friction term discretization and limitation to preserve stability and conservation in the 1D shallow‐water model: Application to unsteady irrigation and river flow , 2008 .

[17]  Chi-Wang Shu,et al.  Numerical Comparison of WENO Finite Volume and Runge–Kutta Discontinuous Galerkin Methods , 2001, J. Sci. Comput..

[18]  Ethan J. Kubatko,et al.  A wetting and drying treatment for the Runge-Kutta discontinuous Galerkin solution to the shallow water equations , 2009 .

[19]  Mourad Heniche,et al.  A two-dimensional finite element drying-wetting shallow water model for rivers and estuaries , 2000 .

[20]  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..

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

[22]  Bram van Leer,et al.  Upwind and High-Resolution Methods for Compressible Flow: From Donor Cell to Residual-Distribution Schemes , 2003 .

[23]  D. Causon,et al.  The surface gradient method for the treatment of source terms in the shallow-water equations , 2001 .

[24]  Manmohan Singh,et al.  Moving boundary shallow water flow above parabolic bottom topography , 2006 .

[25]  Pilar García-Navarro,et al.  A numerical model for the flooding and drying of irregular domains , 2002 .

[26]  Vincenzo Casulli,et al.  A high‐resolution wetting and drying algorithm for free‐surface hydrodynamics , 2009 .

[27]  Nicolas G. Wright,et al.  Simple and efficient solution of the shallow water equations with source terms , 2010 .

[28]  I. Nikolos,et al.  An unstructured node-centered finite volume scheme for shallow water flows with wet/dry fronts over complex topography , 2009 .

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

[30]  Yulong Xing,et al.  A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a Class of Hyperbolic Systems with Source , 2006 .

[31]  Qiuhua Liang,et al.  Flood Simulation Using a Well-Balanced Shallow Flow Model , 2010 .

[32]  J. Murillo,et al.  Analysis of the Friction Term in the One-Dimensional Shallow-Water Model , 2007 .

[33]  Brett F. Sanders,et al.  Conservative Wetting and Drying Methodology for Quadrilateral Grid Finite-Volume Models , 2007 .

[34]  Benedict D. Rogers,et al.  Mathematical balancing of flux gradient and source terms prior to using Roe's approximate Riemann solver , 2003 .

[35]  Q. Liang,et al.  Numerical resolution of well-balanced shallow water equations with complex source terms , 2009 .

[36]  Vincent Guinot,et al.  Godunov-type Schemes: An Introduction for Engineers , 2003 .

[37]  P. García-Navarro,et al.  Godunov-type methods for free-surface shallow flows: A review , 2007 .

[38]  Chi-Wang Shu,et al.  The Runge-Kutta local projection $P^1$-discontinuous-Galerkin finite element method for scalar conservation laws , 1988, ESAIM: Mathematical Modelling and Numerical Analysis.

[39]  Bernardo Cockburn.  Discontinuous Galerkin methods , 2003 .

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

[41]  Jostein R. Natvig,et al.  Well-balanced finite volume schemes of arbitrary order of accuracy for shallow water flows , 2006, J. Comput. Phys..

[42]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

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

[44]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[45]  Tae Hoon Yoon,et al.  Finite volume model for two-dimensional shallow water flows on unstructured grids , 2004 .

[46]  Vijaya R. Ambati,et al.  Space-time discontinuous Galerkin discretization of rotating shallow water equations on moving grids , 2006 .

[47]  Clinton N Dawson,et al.  A discontinuous Galerkin method for two-dimensional flow and transport in shallow water , 2002 .

[48]  Yulong Xing,et al.  High order well-balanced finite volume WENO schemes and discontinuous Galerkin methods for a class of hyperbolic systems with source terms , 2006, J. Comput. Phys..

[49]  Pierre Fabrie,et al.  Evaluation of well‐balanced bore‐capturing schemes for 2D wetting and drying processes , 2007 .

[50]  Fabien Marche,et al.  A Positive Preserving High Order VFRoe Scheme for Shallow Water Equations: A Class of Relaxation Schemes , 2008, SIAM J. Sci. Comput..

[51]  Jean-François Remacle,et al.  Boundary discretization for high‐order discontinuous Galerkin computations of tidal flows around shallow water islands , 2009 .

[52]  D. Zhao,et al.  Finite‐Volume Two‐Dimensional Unsteady‐Flow Model for River Basins , 1994 .

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

[54]  Alexandre Ern,et al.  A well‐balanced Runge–Kutta discontinuous Galerkin method for the shallow‐water equations with flooding and drying , 2008 .

[55]  Pilar García-Navarro,et al.  Zero mass error using unsteady wetting–drying conditions in shallow flows over dry irregular topography , 2004 .