Non-hydrostatic pressure shallow flows: GPU implementation using finite volume and finite difference scheme

We consider the depth-integrated non-hydrostatic system derived by Yamazaki et al. An efficient formally second-order well-balanced hybrid finite volume finite difference numerical scheme is proposed. The scheme consists of a two-step algorithm based on a projection-correction type scheme initially...

Full description

Saved in:
Bibliographic Details
Published inApplied mathematics and computation Vol. 338; pp. 631 - 659
Main Authors Escalante, C., Morales de Luna, T., Castro, M.J.
Format Journal Article
LanguageEnglish
Published Elsevier Inc 01.12.2018
Subjects
Finite difference
Finite volume
GPU
Non-hydrostatic Shallow-Water
Tsunami simulation
Wave breaking
Finite difference
Non-hydrostatic Shallow-Water
Finite volume
Tsunami simulation
GPU
Wave breaking
Online AccessGet full text

Cover

More Information
Summary:We consider the depth-integrated non-hydrostatic system derived by Yamazaki et al. An efficient formally second-order well-balanced hybrid finite volume finite difference numerical scheme is proposed. The scheme consists of a two-step algorithm based on a projection-correction type scheme initially introduced by Chorin–Temam [15]. First, the hyperbolic part of the system is discretized using a Polynomial Viscosity Matrix path-conservative finite volume method. Second, the dispersive terms are solved by means of compact finite differences. A new methodology is also presented to handle wave breaking over complex bathymetries. This adapts well to GPU-architectures and guidelines about its GPU implementation are introduced. The method has been applied to idealized and challenging experimental test cases, which shows the efficiency and accuracy of the method.
ISSN:0096-3003
1873-5649
DOI:10.1016/j.amc.2018.06.035