High order finite volume WENO schemes for the shallow water flows through channels with irregular geometry

The shallow water equations are widely used to model flows in rivers and coastal areas. In this paper, we consider the shallow water flows in open channels with irregular geometry and a non-flat bottom topography, and design high order finite volume weighted essentially non-oscillatory (WENO) methods. A special source term approximation is introduced so that the proposed methods can preserve the still water steady state exactly. We also employ a simple positivity-preserving limiter to provide efficient and robust simulations near the wetting and drying front. The proposed methods are well-balanced for the still water steady state solutions, preserve the non-negativity of the wet cross section, and are genuinely high order accurate in smooth regions for general solutions and essentially non-oscillatory for general solutions with discontinuities. Numerical examples are performed at the end to verify these properties.

[1]  Manuel Jesús Castro Díaz,et al.  High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems , 2006, Math. Comput..

[2]  Javier Murillo,et al.  Accurate numerical modeling of 1D flow in channels with arbitrary shape. Application of the energy balanced property , 2014, J. Comput. Phys..

[3]  Pilar García-Navarro,et al.  1‐D Open‐Channel Flow Simulation Using TVD‐McCormack Scheme , 1992 .

[4]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[5]  Chi-Wang Shu.  Total-variation-diminishing time discretizations , 1988 .

[6]  Rémi Abgrall,et al.  A comment on the computation of non-conservative products , 2010, J. Comput. Phys..

[7]  Yulong Xing,et al.  High-order well-balanced finite volume WENO schemes for shallow water equation with moving water , 2007, J. Comput. Phys..

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

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

[10]  J LeVequeRandall.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

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

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

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

[14]  Chi-Wang Shu,et al.  High-order finite volume WENO schemes for the shallow water equations with dry states , 2011 .

[15]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

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

[17]  B. Perthame,et al.  A kinetic scheme for the Saint-Venant system¶with a source term , 2001 .

[18]  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 Terms† , 2005 .

[19]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

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

[21]  Yulong Xing,et al.  Positivity-Preserving Well-Balanced Discontinuous Galerkin Methods for the Shallow Water Equations on Unstructured Triangular Meshes , 2013, J. Sci. Comput..

[22]  R. LeVeque,et al.  Balancing Source Terms and Flux Gradientsin High-Resolution Godunov Methods : The Quasi-Steady Wave-Propogation AlgorithmRandall , 1998 .

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

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

[25]  Doron Levy,et al.  CENTRAL-UPWIND SCHEMES FOR THE SAINT-VENANT SYSTEM , 2002 .

[26]  Smadar Karni,et al.  A CENTRAL SCHEME FOR SHALLOW WATER FLOWS ALONG CHANNELS WITH IRREGULAR GEOMETRY , 2009 .

[27]  Yulong Xing,et al.  Exactly well-balanced discontinuous Galerkin methods for the shallow water equations with moving water equilibrium , 2014, J. Comput. Phys..

[28]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[29]  P. García-Navarro,et al.  On numerical treatment of the source terms in the shallow water equations , 2000 .

[30]  Smadar Karni,et al.  Shallow Water Flows in Channels , 2011, J. Sci. Comput..

[31]  Yulong Xing,et al.  High order finite difference WENO schemes with the exact conservation property for the shallow water equations , 2005 .

[32]  Xiangxiong Zhang,et al.  On maximum-principle-satisfying high order schemes for scalar conservation laws , 2010, J. Comput. Phys..

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

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