Artificial Viscosity Technique: A Riemann-Solver-Free Method for 2D Urban Flood Modelling on Complex Topography

This study deals with a simulation of two-dimensional urban flood problems on complex topography based on a cell-centred finite volume (CCFV) scheme. Unlike many numerical models that use a Riemann solver to deal with discontinuities due to the rapid change of a flow regime (caused by very shallow water) or due to wet–dry problems, an artificial viscosity technique is used in this study to tackle numerical instabilities caused by such discontinuities. This technique is a Riemann-solver-free method for solving the shallow water equations and is constructed from a combination of a Laplacian and a biharmonic operator, in which the variable scaling factor is devised by using the spectral radius of the Jacobian matrix. For a time discretisation, the Runge–Kutta fourth-order scheme is then used to achieve high-order accuracy. In order to avoid a computational overhead, this Runge–Kutta scheme is applied in its hybrid formulation, in which the artificial viscosity is only computed once per time step. Another advantage of our technique is a simple computation of the convective flux which is performed only by averaging the left and right states of every edge instead of evaluating complex if-then-else statements as required in the Riemann solver such as Harten-Lax-van Leer-contact (HLLC) scheme. Other improvement aspects address both the proper treatment of the friction source term when dealing with very shallow water on very rough beds and an advanced wet–dry technique which is solely applied in an edge-based fashion. Our results show that this artificial viscosity technique is highly accurate for solving the shallow water equations. Also, we show that this technique is cheaper than the HLLC scheme and entails a much less computational complexity.

[1]  Application of Finite Volume Cell Center Method with Wet and Dry Treatment in Hydrodynamic Flow Modeling , 2012 .

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

[3]  Ioannis K. Nikolos,et al.  Performance and Comparison of Cell-Centered and Node-Centered Unstructured Finite Volume Discretizations for Shallow Water Free Surface Flows , 2011 .

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

[5]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

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

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

[8]  V. A. Prokof'ev,et al.  State-Of-The-Art Numerical Schemes Based on the Control Volume Method for Modeling Turbulent Flows and Dam-Break Waves , 2002 .

[9]  Abdolmajid Mohammadian,et al.  Simulation of shallow flows over variable topographies using unstructured grids , 2006 .

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

[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]  Masayuki Fujihara,et al.  Godunov-Type Solution of Curvilinear Shallow-Water Equations , 2000 .

[13]  Javier Murillo,et al.  Time step restrictions for well‐balanced shallow water solutions in non‐zero velocity steady states , 2009 .

[14]  Ioannis K. Nikolos,et al.  A novel multidimensional solution reconstruction and edge‐based limiting procedure for unstructured cell‐centered finite volumes with application to shallow water dynamics , 2013 .

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

[16]  J. Kuerten,et al.  Improved shock‐capturing of Jameson's scheme for the Euler equations , 1991 .

[17]  Doron Levy,et al.  CENTRAL-UPWIND SCHEMES FOR THE SAINT-VENANT SYSTEM , 2002 .

[18]  D. Mavriplis.  Multigrid solution of the two-dimensional Euler equations on unstructured triangular meshes , 1987 .

[19]  Bobby Minola Ginting,et al.  A two-dimensional artificial viscosity technique for modelling discontinuity in shallow water flows , 2017 .

[20]  R. C. Swanson.  On Some Numerical Dissipation Schemes , 1998 .

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

[22]  Guohua Liu,et al.  Development of a Cell-Centered Godunov-Type Finite Volume Model for Shallow Water Flow Based on Unstructured Mesh , 2014 .

[23]  A. Jameson,et al.  Finite volume solution of the two-dimensional Euler equations on a regular triangular mesh , 1985 .

[24]  Jingming Hou,et al.  A robust well-balanced model on unstructured grids for shallow water flows with wetting and drying over complex topography , 2013 .

[25]  G. Petrova,et al.  A SECOND-ORDER WELL-BALANCED POSITIVITY PRESERVING CENTRAL-UPWIND SCHEME FOR THE SAINT-VENANT SYSTEM ∗ , 2007 .