Simple and efficient solution of the shallow water equations with source terms

A simple and efficient method to solve the one-dimensional shallow water equations with source terms is presented. To avoid a fractional step method for the discretization of the source terms, a homogeneous form of the shallow water equations is proposed and well-known conservative numerical schemes are modified to solve the new form of the equations. The modification to the homogeneous form equations combines the source terms with the flux term and solves them by the same solution structure of the numerical scheme. As a result, the source terms are automatically discretized to achieve perfect balance with flux terms without any special treatment and the method does not introduce numerical errors. The proposed method is verified against several benchmark tests and shows good agreement with analytical solutions. Copyright © 2009 John Wiley & Sons, Ltd.

[1]  A. Jha,et al.  High resolution flux‐difference‐splitting scheme on adaptive grid for open‐channel flows , 2001 .

[2]  José M. Corberán,et al.  Construction of second-order TVD schemes for nonhomogeneous hyperbolic conservation laws , 2001 .

[3]  Y. Zech,et al.  Treatment of Natural Geometry in Finite Volume River Flow Computations , 2003 .

[4]  Bernd Einfeld On Godunov-type methods for gas dynamics , 1988 .

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

[6]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

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

[8]  Robert W. MacCormack,et al.  Numerical solution of the interaction of a shock wave with a laminar boundary layer , 1971 .

[9]  P. García-Navarro,et al.  McCormack's method for the numerical simulation of one-dimensional discontinuous unsteady open channel flow , 1992 .

[10]  P. Lax,et al.  Systems of conservation laws , 1960 .

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

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

[13]  Pilar García-Navarro,et al.  Dam-break flow simulation : some results for one-dimensional models of real cases , 1999 .

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

[15]  Pilar García-Navarro,et al.  Improving simple explicit methods for unsteady open channel and river flow , 2004 .

[16]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

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

[18]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

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

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

[21]  Reinaldo García,et al.  Numerical solution of the St. Venant equations with the MacCormack finite‐difference scheme , 1986 .