Comparison of Shallow Water Solvers: Applications for Dam-Break and Tsunami Cases with Reordering Strategy for Efficient Vectorization on Modern Hardware

We investigate in this paper the behaviors of the Riemann solvers (Roe and Harten-Lax-van Leer-Contact (HLLC) schemes) and the Riemann-solver-free method (central-upwind scheme) regarding their accuracy and efficiency for solving the 2D shallow water equations. Our model was devised to be spatially second-order accurate with the Monotonic Upwind Scheme for Conservation Laws (MUSCL) reconstruction for a cell-centered finite volume scheme—and be temporally fourth-order accurate using the Runge–Kutta fourth-order method. Four benchmark cases of dam-break and tsunami events dealing with highly-discontinuous flows and wet–dry problems were simulated. To this end, we applied a reordering strategy for the data structures in our code supporting efficient vectorization and memory access alignment for boosting the performance. Two main features are pointed out here. Firstly, the reordering strategy employed has enabled highly-efficient vectorization for the three solvers investigated on three modern hardware (AVX, AVX2, and AVX-512), where speed-ups of 4.5–6.5× were obtained on the AVX/AVX2 machines for eight data per vector while on the AVX-512 machine we achieved a speed-up of up to 16.7× for 16 data per vector, all with singe-core computation; with parallel simulations, speed-ups of up to 75.7–121.8× and 928.9× were obtained on AVX/AVX2 and AVX-512 machines, respectively. Secondly, we observed that the central-upwind scheme was able to outperform the HLLC and Roe schemes 1.4× and 1.25×, respectively, by exhibiting similar accuracies. This study would be useful for modelers who are interested in developing shallow water codes.

[1]  Ralf-Peter Mundani,et al.  Artificial Viscosity Technique: A Riemann-Solver-Free Method for 2D Urban Flood Modelling on Complex Topography , 2018 .

[2]  Michael Bader,et al.  Vectorization of Riemann solvers for the single- and multi-layer Shallow Water Equations , 2018, 2018 International Conference on High Performance Computing & Simulation (HPCS).

[3]  Ayal Zaks,et al.  Auto-vectorization of interleaved data for SIMD , 2006, PLDI '06.

[4]  T. Gallouët,et al.  Some approximate Godunov schemes to compute shallow-water equations with topography , 2003 .

[5]  Ralf-Peter Mundani,et al.  Parallel Flood Simulations for Wet–Dry Problems Using Dynamic Load Balancing Concept , 2019, Journal of Computing in Civil Engineering.

[6]  I. Nikolos,et al.  An unstructured node-centered finite volume scheme for shallow water flows with wet/dry fronts over complex topography , 2009 .

[7]  Dongfang Liang,et al.  Urban flood modeling using shallow water equations with depth-dependent anisotropic porosity , 2016 .

[8]  Qiuhua Liang,et al.  An efficient unstructured MUSCL scheme for solving the 2D shallow water equations , 2015, Environ. Model. Softw..

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

[10]  Michael J. Briggs,et al.  Laboratory experiments of tsunami runup on a circular island , 1995 .

[11]  Ralf-Peter Mundani,et al.  Parallel Simulations of Shallow Water Solvers for Modelling Overland Flows , 2018 .

[12]  Alexander Kurganov,et al.  Well-balanced positivity preserving cell-vertex central-upwind scheme for shallow water flows , 2016 .

[13]  Ji Yueh Liu,et al.  Hybrid OpenMP/AVX acceleration of a Split HLL Finite Volume Method for the Shallow Water and Euler Equations , 2015 .

[14]  Bobby Minola Ginting,et al.  Central-upwind scheme for 2D turbulent shallow flows using high-resolution meshes with scalable wall functions , 2019, Computers & Fluids.

[15]  J. Hou,et al.  An efficient and stable hydrodynamic model with novel source term discretization schemes for overland flow and flood simulations , 2017 .

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

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

[18]  Günter Blöschl,et al.  A two‐dimensional numerical scheme of dry/wet fronts for the Saint‐Venant system of shallow water equations , 2015 .

[19]  Carlos Parés,et al.  NUMERICAL TREATMENT OF WET/DRY FRONTS IN SHALLOW FLOWS WITH A MODIFIED ROE SCHEME , 2006 .

[20]  Alexander Kurganov,et al.  Semidiscrete Central-Upwind Schemes for Hyperbolic Conservation Laws and Hamilton-Jacobi Equations , 2001, SIAM J. Sci. Comput..

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

[22]  Randall J. LeVeque,et al.  Validating Velocities in the GeoClaw Tsunami Model Using Observations near Hawaii from the 2011 Tohoku Tsunami , 2013, Pure and Applied Geophysics.

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

[24]  Haijun Yu,et al.  Efficient Finite-Volume Model for Shallow-Water Flows Using an Implicit Dual Time-Stepping Method , 2015 .

[25]  Sebastian Rettenberger,et al.  Vectorization of an augmented Riemann solver for the shallow water equations , 2014, 2014 International Conference on High Performance Computing & Simulation (HPCS).

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

[27]  Aart J. C. Bik,et al.  Automatic Intra-Register Vectorization for the Intel® Architecture , 2002, International Journal of Parallel Programming.

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

[29]  Y. Zech,et al.  Experimental study of dam-break flow against an isolated obstacle , 2007 .

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

[31]  A. Duran,et al.  A robust and well‐balanced scheme for the 2D Saint‐Venant system on unstructured meshes with friction source term , 2015 .

[32]  Luis Cea,et al.  A simple and efficient unstructured finite volume scheme for solving the shallow water equations in overland flow applications , 2015 .

[33]  S. F. Davis.  Simplified second-order Godunov-type methods , 1988 .

[34]  J. Duan,et al.  Two-dimensional depth-averaged finite volume model for unsteady turbulent flow , 2012 .