Implicit finite volume simulation of 2D shallow water flows in flexible meshes

Abstract In this work, an implicit method for solving 2D hyperbolic systems of equations is presented, focusing on the application to the 2D shallow water equations. It is based on the first order Roe’s scheme, in the framework of finite volume methods. A conservative linearization is done for the flux terms, leading to a non-structured matrix for unstructured meshes thus requiring iterative methods for solving the system. The validation is done by comparing numerical and exact solutions in both unsteady and steady cases. In order to test the applicability of the implicit scheme to real world situations, a laboratory scale tsunami simulation is carried out and compared to the experimental data. The implicit schemes have the advantage of the unconditional stability, but a quality loss in the transient solution can appear for high CFL numbers. The properties of the scheme are well suited for the simulation of unsteady shallow water flows over irregular topography using all kind of meshes.

[1]  Pilar García-Navarro,et al.  Implicit schemes with large time step for non‐linear equations: application to river flow hydraulics , 2004 .

[2]  Pilar García-Navarro,et al.  2D Zero-Inertia Model for Solution of Overland Flow Problems in Flexible Meshes , 2016 .

[3]  Pilar García-Navarro,et al.  An implicit method for water flow modelling in channels and pipes , 1994 .

[4]  Stephen C. Medeiros,et al.  Review of wetting and drying algorithms for numerical tidal flow models , 2013 .

[5]  Samir Kumar Das,et al.  Modeling of Shallow-Water Equations by Using Implicit Higher-Order Compact Scheme with Application to Dam-Break Problem , 2013 .

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

[7]  Nancy Nichols,et al.  Analytic Benchmark Solutions for Open-Channel Flows , 1997 .

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

[9]  Yousef Saad,et al.  ILUT: A dual threshold incomplete LU factorization , 1994, Numer. Linear Algebra Appl..

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

[11]  Javier Murillo,et al.  Augmented versions of the HLL and HLLC Riemann solvers including source terms in one and two dimensions for shallow flow applications , 2012, J. Comput. Phys..

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

[13]  Javier Murillo,et al.  Friction term discretization and limitation to preserve stability and conservation in the 1D shallow‐water model: Application to unsteady irrigation and river flow , 2008 .

[14]  A. M. Wasantha Lal,et al.  Performance Comparison of Overland Flow Algorithms , 1998 .

[15]  A. M. Wasantha Lal,et al.  Weighted Implicit Finite-Volume Model for Overland Flow , 1998 .

[16]  Qiuhua Liang,et al.  A positivity-preserving zero-inertia model for flood simulation , 2011 .

[17]  Qiuhua Liang,et al.  RKDG2 shallow-water solver on non-uniform grids with local time steps: Application to 1D and 2D hydrodynamics , 2015 .

[18]  Shuangcai Li,et al.  Fully-Coupled Modeling of Shallow Water Flow and Pollutant Transport on Unstructured Grids , 2012 .

[19]  Pilar García-Navarro,et al.  Diffusive-Wave Based Hydrologic-Hydraulic Model with Sediment Transport. I: Model Development , 2012 .

[20]  Michael Dumbser,et al.  A high order semi-implicit discontinuous Galerkin method for the two dimensional shallow water equations on staggered unstructured meshes , 2014, Appl. Math. Comput..

[21]  Implicit TVDLF Methods for Diffusion and Kinematic Flows , 2013 .

[22]  H. C. Yee,et al.  Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations. [in gas dynamics , 1985 .

[23]  V. Legat,et al.  A fully implicit wetting-drying method for DG-FEM shallow water models, with an application to the Scheldt Estuary , 2011 .

[24]  Joannes J. Westerink,et al.  A solution for the vertical variation of stress, rather than velocity, in a three-dimensional circulation model , 1991 .

[25]  J. Fernández-Pato,et al.  Simulación numérica de inundación de valles fluviales mediante un modelo difusivo implícito , 2016 .

[26]  Javier Murillo,et al.  Weak solutions for partial differential equations with source terms: Application to the shallow water equations , 2010, J. Comput. Phys..

[27]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[28]  Javier Murillo,et al.  Wave Riemann description of friction terms in unsteady shallow flows: Application to water and mud/debris floods , 2012, J. Comput. Phys..

[29]  Pilar García-Navarro,et al.  A 2D extension of a Large Time Step explicit scheme (CFL > 1) for unsteady problems with wet/dry boundaries , 2014, J. Comput. Phys..

[30]  Vincenzo Casulli,et al.  Semi-implicit finite difference methods for the two-dimensional shallow water equation , 1990 .

[31]  Jerónimo Puertas,et al.  Experimental validation of two-dimensional depth-averaged models for forecasting rainfall-runoff from precipitation data in urban areas , 2010 .