Adaptive quadtree simulation of shallow flows with wet-dry fronts over complex topography

Many natural terrains have complicated surface topography. The simulation of steep-fronted flows that occur after heavy rainfall flash floods or as inundation from dyke breaches is usually based on the non-linear shallow water equations in hyperbolic conservation form. Particular challenges to numerical modellers are posed by the need to balance correctly the flux gradient and source terms in Godunov-type finite volume shock-capturing schemes and by the moving wet–dry boundary as the flood rises or falls. This paper presents a Godunov-type shallow flow solver on adaptive quadtree grids aimed at simulating flood flows as they travel over natural terrain. By choosing the stage and discharge as dependent variables in the hyperbolic non-linear shallow water equations, a new deviatoric formulation is derived that mathematically balances the flux gradient and source terms in cases where there are wet–dry fronts. The new formulation is more general in application than previous a priori approaches. Three benchmark tests are used to validate the solver, and include steady flow over a submerged hump, flow disturbances propagating over an elliptical-shaped hump, and free surface sloshing motions in a vessel with a parabolic bed. The model is also used to simulate the propagation of a flood due to a dam break over an initially dry floodplain containing three humps.

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

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

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

[4]  Qiuhua Liang,et al.  Simulation of dam‐ and dyke‐break hydrodynamics on dynamically adaptive quadtree grids , 2004 .

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

[6]  Masayuki Fujihara,et al.  Adaptive Q-tree Godunov-type scheme for shallow water equations , 2001 .

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

[8]  C. Hirsch.  Computational methods for inviscid and viscous flows , 1990 .

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

[10]  P. Bates,et al.  A simple raster-based model for flood inundation simulation , 2000 .

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

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

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

[14]  Siegfried Müller,et al.  Solution of shallow water equations using fully adaptive multiscale schemes , 2005 .

[15]  Bharat K. Soni,et al.  Handbook of Grid Generation , 1998 .

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

[17]  C. Hirsch.  Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows , 1990 .

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

[19]  Stuart N. Lane,et al.  Urban fluvial flood modelling using a two‐dimensional diffusion‐wave treatment, part 2: development of a sub‐grid‐scale treatment , 2006 .

[20]  K. Anastasiou,et al.  SOLUTION OF THE 2D SHALLOW WATER EQUATIONS USING THE FINITE VOLUME METHOD ON UNSTRUCTURED TRIANGULAR MESHES , 1997 .

[21]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[22]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[23]  S. Lane,et al.  Urban fluvial flood modelling using a two‐dimensional diffusion‐wave treatment, part 1: mesh resolution effects , 2006 .

[24]  Alessandro Valiani,et al.  Divergence Form for Bed Slope Source Term in Shallow Water Equations , 2006 .

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

[26]  Michael B Abbott,et al.  Coastal, estuarial and harbour engineers' reference book , 1993 .

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

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

[29]  Manmohan Singh,et al.  Moving boundary shallow water flow above parabolic bottom topography , 2006 .

[30]  P. L. George,et al.  Automatic Mesh Generation: Application to Finite Element Methods , 1992 .

[31]  C. I. Moutzouris,et al.  Godunov-type solution of the shallow water equations on adaptive unstructured triangular grids , 2006 .

[32]  Matthew E. Hubbard,et al.  A 2D numerical model of wave run-up and overtopping , 2002 .

[33]  Eleuterio F. Toro,et al.  Experimental and numerical assessment of the shallow water model for two-dimensional dam-break type problems , 1995 .

[34]  Galina V. Muratova,et al.  Adaptive grid shallow water modeling , 2000 .

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

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

[37]  Qiuhua Liang,et al.  Shallow flow simulation on dynamically adaptive cut cell quadtree grids , 2007 .

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

[39]  Derek M. Causon,et al.  Calculation of shallow water flows using a Cartesian cut cell approach , 2000 .

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

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

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