A robust 2D shallow water model for solving flow over complex topography using homogenous flux method

A robust Godunov-type numerical scheme solver is proposed for solving 2D SWEs and is applied to simulate flow over complex topography with wetting and drying. In reality, the topography is usually complex and irregular; therefore, to avoid the numerical errors generated by such features, a Homogenous Flux Method is used to handle the bed slope term in the SWEs. The method treats the bed slope term as a flux to be incorporated into the flux gradient and so maintains the balance between the two in a Godunov-type shock-capturing scheme. The main advantages of the method are: first, it is simple and easy to implement; second, numerical experiments demonstrate that it can handle discontinuous or vertical bed topography without any special treatment and third, it is applicable to both steady and unsteady flows. It is demonstrated how the approach set out here can be applied to the nonlinear hyperbolic system of the SWEs. The two-dimensional hyperbolic system is then solved by use of a second-order total-variation-diminishing version of the weighted average flux method in conjunction with a Harten-Lax-van Leer-Contract approximate Riemann solver incorporating the new flux gradient term. Several benchmark tests are presented to validate the model and the approach is verified against experimental measurements from the European Union Concerted Action on Dam Break Modelling project. These show very good agreement. Finally, the method is applied to a volcano-induced outburst flood over an initially dry channel with complex irregular topography to demonstrate the technique's capability in simulating a real flood.

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

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

[3]  Yi Liu,et al.  A robust well-balanced finite volume model for shallow water flows with wetting and drying over irregular terrain , 2011 .

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

[5]  János Józsa,et al.  Solution-adaptivity in modelling complex shallow flows , 2007 .

[6]  Ó. Knudsen,et al.  An unusual jökulhlaup resulting from subglacial volcanism, Sólheimajökull, Iceland , 2010 .

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

[8]  G. Stelling,et al.  On the construction of computational methods for shallow water flow problems , 1983 .

[9]  M. Berzins,et al.  An unstructured finite-volume algorithm for predicting flow in rivers and estuaries , 1998 .

[10]  Valerio Caleffi,et al.  Finite volume method for simulating extreme flood events in natural channels , 2003 .

[11]  E. F. Toro,et al.  Riemann problems and the WAF method for solving the two-dimensional shallow water equations , 1992, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

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

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

[14]  Brett F. Sanders,et al.  Integration of a shallow water model with a local time step , 2008 .

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

[16]  Roger Alexander Falconer,et al.  Mathematical modelling of jet-forced circulation in reservoirs and harbours , 1977 .

[17]  Sandra Soares Frazao,et al.  Dam break in channels with 90 degrees bend , 2002 .

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

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

[20]  Qiuhua Liang,et al.  Well-balanced RKDG2 solutions to the shallow water equations over irregular domains with wetting and drying , 2010 .

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

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

[23]  Pilar García-Navarro,et al.  Unsteady free surface flow simulation over complex topography with a multidimensional upwind technique , 2003 .

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

[25]  Derek M. Causon,et al.  Numerical solutions of the shallow water equations with discontinuous bed topography , 2002 .

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

[27]  Ó. Knudsen,et al.  Flash flood at Sólheimajökull heralds the reawakening of an Icelandic subglacial volcano , 2000 .

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

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

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

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

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

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

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

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

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

[37]  Yves Zech,et al.  Dam Break in Channels with 90° Bend , 2002 .

[38]  A. J. Crossley,et al.  Time accurate local time stepping for the unsteady shallow water equations , 2005 .

[39]  Enrique Domingo Fernández-Nieto,et al.  Extension of WAF Type Methods to Non-Homogeneous Shallow Water Equations with Pollutant , 2008, J. Sci. Comput..

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

[41]  Derek M. Causon,et al.  Numerical prediction of dam-break flows in general geometries with complex bed topography , 2004 .

[42]  Eleuterio F. Toro,et al.  A weighted average flux method for hyperbolic conservation laws , 1989, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[43]  Sandra Soares Frazao,et al.  Dam-break Flow through Sharp Bends - Physical Model and 2D Boltzmann Model Validation , 1999 .