![[PDF]](/pb-assets/Icons/adobe-pdf-icon-1498649896243.png){kind=icon}
Abstract
We present a new graphics processing unit implementation of two second-order numerical schemes of the shallow water equations on Cartesian grids. Previous implementations are not fast enough to evaluate multiple scenarios for a robust, uncertainty-aware decision support. To tackle this, we exploit the capabilities of the NVIDIA Kepler architecture. We implement a scheme developed by Kurganov and Petrova (KP07) and a newer, improved version by Horváth et al. (HWP14). The KP07 scheme is simpler but suffers from incorrect high velocities along the wet/dry boundaries, resulting in small time steps and long simulation runtimes. The HWP14 scheme resolves this problem but comprises a more complex algorithm. Previous work has shown that HWP14 has the potential to outperform KP07, but no practical implementation has been provided. The novel shuffle-based implementation of HWP14 presented here takes advantage of its accuracy and performance capabilities for real-world usage. The correctness and performance are validated on real-world scenarios.
References
|
Acuña, M, Aoki, T (2009) Real-time tsunami simulation on multi-node GPU cluster. In: ACM/IEEE Conference on Supercomputing, 2009. Oregon Convention Center, Portland, Oregon, 14—20 November 2009. Google Scholar | |
|
Bollermann, A, Chen, G, Kurganov, A. (2013) A well-balanced reconstruction of wet/dry fronts for the shallow water equations. Journal of Scientific Computing56(2): 267–290. Google Scholar | Crossref | ISI | |
|
Brandvik, T, Pullan, G (2008) Acceleration of a 3D Euler solver using commodity graphics hardware. In: 46th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, USA, 7–10 January 2008, p. 607. Google Scholar | Crossref | |
|
Brodtkorb, AR, Sætra, ML (2012) Explicit shallow water simulations on GPUs: Guidelines and best practices. In: XIX International Conference on Water Resources, CMWR, Champaign, USA, 17–22 June 2012, pp. 17–22. Google Scholar | |
|
Brodtkorb, AR, Sætra, ML, Altinakar, M (2012) Efficient shallow water simulations on GPUs: Implementation, visualization, verification, and validation. Computers & Fluids55: 1–12. Google Scholar | Crossref | ISI | |
|
Brodtkorb, AR, Hagen, TR, Lie, K.-A. (2010) Simulation and visualization of the Saint-Venant system using GPUs. Computing and Visualization in Science13(7): 341–353. Google Scholar | Crossref | |
|
Causon, DM, Ingram, DM, Mingham, CG (2001) A cartesian cut cell method for shallow water flows with moving boundaries. Advances in Water Resources24(8): 899–911. Google Scholar | Crossref | ISI | |
|
Courant, R, Friedrichs, K, Lewy, H (1967) On the partial difference equations of mathematical physics. IBM Journal of Research and Development11(2): 215–234. Google Scholar | Crossref | ISI | |
|
Cook, S (2012) CUDA Programming: A Developer’s Guide to Parallel Computing with GPUs, Applications of GPU Computing Series. Waltham, MA: Elsevier Science, 2012. Google Scholar | |
|
de la Asunción, M, Mantas, JM, Castro, MJ (2010) Programming CUDA-based GPUs to simulate two-layer shallow water flows. In: Euro-Par 2010-Parallel Processing, 2010, Italy, August 3–September 3, pp. 353–364. Heidelberg, Germany: Springer. Google Scholar | Crossref | |
|
De La Asunción, M, Mantas, JM, Castro, MJ (2011) Simulation of one-layer shallow water systems on multicore and CUDA architectures. The Journal of Supercomputing58(2): 206–214. Google Scholar | Crossref | ISI | |
|
de la Asuncin, M, Castro, MJ, Fernndez-Nieto, E. (2013) Efficient GPU implementation of a two waves TVD-WAF method for the two-dimensional one layer shallow water system on structured meshes. Computers & Fluids80: 441–452. Google Scholar | Crossref | ISI | |
|
Duran, A, Liang, Q, Marche, F (2013) On the well-balanced numerical discretization of shallow water equations on unstructured meshes. Journal of Computational Physics235: 565–586. Google Scholar | Crossref | ISI | |
|
EU (2007) Directive 2007/60/EC of the European Parliament and of the Council of 23 October 2007 on the assessment and management of flood risks. Official Journal of the European UnionL288: 27–34. Google Scholar | |
|
Frazao, SS, Alcrudo, F, Goutal, N (1999) Dam-break test cases summary. In: The Proceedings of the 4th CADAM Meeting, Zaragoza, Spain, 18—19 November 1999. Google Scholar | |
|
George, D (2011) Adaptive finite volume methods with well-balanced Riemann solvers for modeling floods in rugged terrain: application to the Malpasset dam-break flood (France, 1959). International Journal for Numerical Methods in Fluids 66(8): 1000–1018. Google Scholar | Crossref | ISI | |
|
Gottlieb, S, Shu, C-W (1998) Total variation diminishing Runge-Kutta schemes. Mathematics of Computation of the American Mathematical Society67(221): 73–85. Google Scholar | Crossref | |
|
Gottlieb, D, Tadmor, E (1991) The CFL condition for spectral approximations to hyperbolic initial-boundary value problems. Mathematics of Computation56(194): 565–588. Google Scholar | Crossref | ISI | |
|
Hagen, TR, Lie, K-A, Natvig, JR (2006) Solving the Euler equations on graphics processing units. In: Alexandrov, VN, van Albada, GD, Sloot, PMA. (eds) Computational Science–ICCS 2006, Reading, UK, 28—31 May 2006, pp. 220–227. Berlin/Heidelberg: Springer. Google Scholar | Crossref | |
|
Hagen, TR, Hjelmervik, JM, Lie, K.-A (2005) Visual simulation of shallow-water waves. Simulation Modelling Practice and Theory 13(8): 716–726. Google Scholar | Crossref | ISI | |
|
Hervouet, J-M, Petitjean, A (1999) Malpasset dam-break revisited with two-dimensional computations. Journal of Hydraulic Research 37(6): 777–788. Google Scholar | Crossref | ISI | |
|
Horváth, Z, Waser, JW, Perdigão, RAP. (2014) A two-dimensional numerical scheme of dry/wet fronts for the Saint-Venant system of shallow water equations. International Journal for Numerical Methods in Fluids. doi:10.1002/fld.3983. Available at: http://dx.doi.org/10.1002/fld.3983 (accessed 2 February 2016). Google Scholar | |
|
Hou, J, Liang, Q, Zhang, H (2014) Multislope MUSCL method applied to solve shallow water equations. Computers & Mathematics with Applications 68(2): 2012–2027. Google Scholar | Crossref | ISI | |
|
Hou, J, Liang, Q, Simons, F. (2013) A 2d well-balanced shallow flow model for unstructured grids with novel slope source term treatment. Advances in Water Resources 52: 107–131. Google Scholar | Crossref | ISI | |
|
Klöckner, A, Warburton, T, Bridge, J. (2009) Nodal discontinuous Galerkin methods on graphics processors. Journal of Computational Physics 228(21): 7863–7882. Google Scholar | Crossref | ISI | |
|
Kurganov, A, Levy, D (2002) Central-upwind schemes for the saint-venant system. ESAIM: Mathematical Modelling and Numerical Analysis36(03): 397–425. Google Scholar | Crossref | ISI | |
|
Kurganov, A, Petrova, G (2007) A second-order well-balanced positivity preserving central-upwind scheme for the Saint-Venant system. Communications in Mathematical Sciences 5(1): 133–160. Google Scholar | Crossref | ISI | |
|
Lacasta, A, Morales-Hernández, M, Murillo, J. (2014) An optimized gpu implementation of a 2d free surface simulation model on unstructured meshes. Advances in Engineering Software 78: 1–15. Google Scholar | Crossref | ISI | |
|
Lacasta, A, Morales-Hernández, M, Murillo, J. (2015) GPU implementation of the 2D shallow water equations for the simulation of rainfall/runoff events. Environmental Earth Sciences 74(11): 7295–7305. Google Scholar | Crossref | ISI | |
|
Lastra, M, Mantas, JM, Ureña, C. (2009) Simulation of shallow-water systems using graphics processing units. Mathematics and Computers in Simulation80(3): 598–618. Google Scholar | Crossref | ISI | |
|
Li, G, Gao, J, Liang, Q (2013) A well-balanced weighted essentially non-oscillatory scheme for pollutant transport in shallow water. International Journal for Numerical Methods in Fluids71(12): 1566–1587. Google Scholar | Crossref | ISI | |
|
Liang, W-Y, Hsieh, T-J, Satria, MT. (2009) A GPU-based simulation of tsunami propagation and inundation. In: Arrems, H, Chang, S-L (eds) Algorithms and Architectures for Parallel Processing. Berlin/Heidelberg: Springer, pp. 593–603. Google Scholar | Crossref | |
|
Noelle, S, Pankratz, N, Puppo, G. (2006) Well-balanced finite volume schemes of arbitrary order of accuracy for shallow water flows. Journal of Computational Physics 213(2): 474–499. Google Scholar | Crossref | ISI | |
|
NVIDIA Corporation (2015) NVIDIA CUDA Compute Unified Device Architecture Programming Guide. NVIDIA Corporation, 2015. Google Scholar | |
|
Open TELEMAC-MASCARET , http://www.opentelemac.org/ (accessed 24 June 2015). Google Scholar | |
|
Perdigão, RA, Blöschl, G (2014) Spatiotemporal flood sensitivity to annual precipitation: evidence for landscape-climate coevolution. Water Resources Research50(7): 5492–5509. Google Scholar | Crossref | ISI | |
|
Russo, G (2005) Central schemes for conservation laws with application to shallow water equations. In: Rionero, S, Romano, G (eds) Trends and Applications of Mathematics to Mechanics. Milano, Italy: Springer, pp. 225–246. Google Scholar | Crossref | |
|
Sætra, M (2013) Shallow water simulation on GPUs for sparse domains. In: Cangiani, A, Davidchack, RL, Georgoulis, EH. (eds) Numerical Mathematics and Advanced Applications 2011. Heidelberg, Germany: Springer, pp. 673–680. Google Scholar | Crossref | |
|
Sætra, ML, Brodtkorb, AR, Lie, K-A (2014) Efficient GPU-implementation of adaptive mesh refinement for the shallow-water equations. Journal of Scientific Computing63(1): 23–48. Google Scholar | Crossref | ISI | |
|
Schindler, B, Waser, J, Ribicic, H. (2013) Multiverse data-flow control. Visualization and Computer Graphics, IEEE Transactions on19(6): 1005–1019. Google Scholar | Crossref | Medline | ISI | |
|
Shu, C-W (1988) Total-variation-diminishing time discretizations. SIAM Journal on Scientific and Statistical Computing9(6): 1073–1084. Google Scholar | Crossref | ISI | |
|
Vacondio, R, Dal Palù, A, Mignosa, P (2014) GPU-enhanced finite volume shallow water solver for fast flood simulations. Environmental Modelling & Software57: 60–75. Google Scholar | Crossref | ISI | |
|
Visdom—an integrated simulation and visualization system , http://visdom.at/ (accessed 24 June 2015). Google Scholar | |
|
Wang, P, Abel, T, Kaehler, R (2010) Adaptive mesh fluid simulations on GPU. New Astronomy15(7): 581–589. Google Scholar | Crossref | ISI | |
|
Waser, J, Konev, A, Sadransky, B. (2014) Many plans: multidimensional ensembles for visual decision support in flood management. Computer Graphics Forum33(3): 281–290. doi:10.1111/cgf.12384. Available at: http://dx.doi.org/10.1111/cgf.12384 (accessed 2 February 2016). Google Scholar | |
|
Wilt, N (2013) The CUDA handbook: A comprehensive guide to GPU programming. Boston, MA: Pearson Education, 2013. Google Scholar | |
|
Xing, Y, Shu, C-W (2005) High order finite difference WENO schemes with the exact conservation property for the shallow water equations. Journal of Computational Physics208(1): 206–227. Google Scholar | Crossref | ISI |
