Simulation of shallow flows over variable topographies using unstructured grids

Simulation of shallow flows over variable topographies is a challenging case for most available shock-capturing schemes. This problem arises because the source terms and flux gradients are not balanced in the numerical computations. Treatments for this problem generally work well on structured grids, but they are usually too expensive, and most of them are not directly applicable to unstructured grids. In this paper we propose two efficient methods to treat the source terms without upwinding and to satisfy the compatibility condition on unstructured grids. In the first method, the calculation of the bed slope source term is performed by employing a compatible approximation of water depth at the cell interfaces. In the second one, different components of the bed slope term are considered separately and a compatible discretization of the components is proposed. The present treatments are applicable for most schemes including the Roe's method without changing the performance of the original scheme for smooth topographies. Copyright © 2006 John Wiley & Sons, Ltd.

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

[2]  Kun Xu,et al.  A well-balanced gas-kinetic scheme for the shallow-water equations with source terms , 2002 .

[3]  M. Hanif Chaudhry,et al.  Simulation of one-dimensional dam-break flows , 1987 .

[4]  Kenneth G. Powell,et al.  Sonic-point capturing , 1989 .

[5]  Paul Glaister.  Approximate Riemann solutions of the shallow water equations , 1988 .

[6]  Fayssal Benkhaldoun,et al.  Exact solutions to the Riemann problem of the shallow water equations with a bottom step , 2001 .

[7]  Derek M. Causon,et al.  HIGH-RESOLUTION FINITE-VOLUME METHOD FOR SHALLOW WATER FLOWS , 1998 .

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

[9]  Derek M. Causon,et al.  Numerical simulation of wave overtopping of coastal structures using the non-linear shallow water equations , 2000 .

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

[11]  Alberto Guadagnini,et al.  Travel time and trajectory moments of conservative solutes in three dimensional heterogeneous porous media under mean uniform flow , 2005 .

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

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

[14]  Pilar García-Navarro,et al.  Balancing Source Terms and Flux Gradients in Finite Volume Schemes , 2001 .

[15]  Shi Jin,et al.  A steady-state capturing method for hyperbolic systems with geometrical source terms , 2001 .

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

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

[18]  Luka Sopta,et al.  ENO and WENO Schemes with the Exact Conservation Property for One-Dimensional Shallow Water Equations , 2002 .

[19]  Scott F. Bradford,et al.  Finite-Volume Model for Shallow-Water Flooding of Arbitrary Topography , 2002 .

[20]  Marinko Nujić,et al.  Efficient implementation of non-oscillatory schemes for the computation of free-surface flows , 1995 .

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