The influence of source terms on stability, accuracy and conservation in two‐dimensional shallow flow simulation using triangular finite volumes

The two-dimensional shallow water model is a hyperbolic system of equations considered well suited to simulate unsteady phenomena related to some surface wave propagation. The development of numerical schemes to correctly solve that system of equations finds naturally an initial step in two-dimensional scalar equation, homogeneous or with source terms. We shall first provide a complete formulation of the second-order finite volume scheme for this equation, paying special attention to the reduction of the method to first order as a particular case. The explicit first and second order in space upwind finite volume schemes are analysed to provide an understanding of the stability constraints, making emphasis in the numerical conservation and in the preservation of the positivity property of the solution when necessary in the presence of source terms. The time step requirements for stability are defined at the cell edges, related with the traditional Courant–Friedrichs–Lewy (CFL) condition. Copyright © 2007 John Wiley & Sons, Ltd.

[1]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

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

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

[4]  Mutsuto Kawahara,et al.  Finite element method for moving boundary problems in river flow , 1986 .

[5]  J. Cunge,et al.  Practical aspects of computational river hydraulics , 1980 .

[6]  Pilar García-Navarro,et al.  A HIGH-RESOLUTION GODUNOV-TYPE SCHEME IN FINITE VOLUMES FOR THE 2D SHALLOW-WATER EQUATIONS , 1993 .

[7]  Derek M. Causon,et al.  POSITIVELY CONSERVATIVE HIGH‐RESOLUTION CONVECTION SCHEMES FOR UNSTRUCTURED ELEMENTS , 1996 .

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

[9]  A. A. Khan.  Modeling flow over an initially dry bed , 2000 .

[10]  Javier Murillo,et al.  Extension of an explicit finite volume method to large time steps (CFL>1): application to shallow water flows , 2006 .

[11]  Pilar García-Navarro,et al.  Flux difference splitting and the balancing of source terms and flux gradients , 2000 .

[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]  Scott F. Bradford,et al.  Finite-Volume Model for Shallow-Water Flooding of Arbitrary Topography , 2002 .

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

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

[16]  B. Sanders,et al.  Unstructured grid finite-volume algorithm for shallow-water flow and scalar transport with wetting and drying , 2006 .

[17]  Pilar García-Navarro,et al.  Efficient construction of high‐resolution TVD conservative schemes for equations with source terms: application to shallow water flows , 2001 .

[18]  S. Osher,et al.  Triangle based adaptive stencils for the solution of hyperbolic conservation laws , 1992 .

[19]  Bram van Leer,et al.  On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe , 1984 .

[20]  Emmanuel Audusse,et al.  A well-balanced positivity preserving second-order scheme for shallow water flows on unstructured meshes , 2005 .

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

[22]  Monika Wierse,et al.  A new theoretically motivated higher order upwind scheme on unstructured grids of simplices , 1997, Adv. Comput. Math..