Adaptive finite volume methods with well‐balanced Riemann solvers for modeling floods in rugged terrain: Application to the Malpasset dam‐break flood (France, 1959)

The simulation of advancing flood waves over rugged topography, by solving the shallow-water equations with well-balanced high-resolution finite volume methods and block-structured dynamic adaptive mesh refinement (AMR), is described and validated in this paper. The efficiency of block-structured AMR makes large-scale problems tractable, and allows the use of accurate and stable methods developed for solving general hyperbolic problems on quadrilateral grids. Features indicative of flooding in rugged terrain, such as advancing wet―dry fronts and non-stationary steady states due to balanced source terms from variable topography, present unique challenges and require modifications such as special Riemann solvers. A well-balanced Riemann solver for inundation and general (non-stationary) flow over topography is tested in this context. The difficulties of modeling floods in rugged terrain, and the rationale for and efficacy of using AMR and well-balanced methods, are presented. The algorithms are validated by simulating the Malpasset dam-break flood (France, 1959), which has served as a benchmark problem previously. Historical field data, laboratory model data and other numerical simulation results (computed on static fitted meshes) are shown for comparison. The methods are implemented in GEOCLAW, a subset of the open-source CLAWPACK software. All the software is freely available at www. clawpack. org.

[1]  J. Hervouet,et al.  Malpasset dam-break revisited with two-dimensional computations , 1999 .

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

[3]  S. Savage,et al.  The motion of a finite mass of granular material down a rough incline , 1989, Journal of Fluid Mechanics.

[4]  Alistair G.L. Borthwick,et al.  The shallow flow equations solved on adaptive quadtree grids , 2001 .

[5]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

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

[7]  Jacques Sainte-Marie,et al.  Derivation of a non-hydrostatic shallow water model; Comparison with Saint-Venant and Boussinesq systems , 2008 .

[8]  R. LeVeque,et al.  Adaptive Mesh Refinement Using Wave-Propagation Algorithms for Hyperbolic Systems , 1998 .

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

[10]  F. Bouchut.  Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws: and Well-Balanced Schemes for Sources , 2005 .

[11]  Javier Murillo,et al.  The influence of source terms on stability, accuracy and conservation in two‐dimensional shallow flow simulation using triangular finite volumes , 2007 .

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

[13]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[14]  David L. George,et al.  Finite volume methods and adaptive refinement for tsunami propagation and inundation , 2006 .

[15]  David L. George,et al.  Augmented Riemann solvers for the shallow water equations over variable topography with steady states and inundation , 2008, J. Comput. Phys..

[16]  R. LeVeque,et al.  HIGH-RESOLUTION FINITE VOLUME METHODS FOR THE SHALLOW WATER EQUATIONS WITH BATHYMETRY AND DRY STATES , 2008 .

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

[18]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[19]  Richard M. Iverson,et al.  Elements of an Improved Model of Debris-flow Motion , 2009 .

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

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

[22]  R. LeVeque,et al.  FINITE VOLUME METHODS AND ADAPTIVE REFINEMENT FOR GLOBAL TSUNAMI PROPAGATION AND LOCAL INUNDATION. , 2006 .

[23]  Randall J. LeVeque,et al.  A Wave Propagation Method for Conservation Laws and Balance Laws with Spatially Varying Flux Functions , 2002, SIAM J. Sci. Comput..

[24]  R. Iverson,et al.  U. S. Geological Survey , 1967, Radiocarbon.

[25]  R. LeVeque,et al.  High-Resolution Methods and Adaptive Refinement for Tsunami Propagation and Inundation , 2008 .

[26]  Valerio Caleffi,et al.  Case Study: Malpasset Dam-Break Simulation using a Two-Dimensional Finite Volume Method , 2002 .

[27]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[28]  Richard M. Iverson,et al.  Granular avalanches across irregular three-dimensional terrain: 1. Theory and computation , 2004 .

[29]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[30]  R. LeVeque.  Wave Propagation Algorithms for Multidimensional Hyperbolic Systems , 1997 .

[31]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[32]  Kwok Fai Cheung,et al.  Depth‐integrated, non‐hydrostatic model for wave breaking and run‐up , 2009 .