GPU driven finite difference WENO scheme for real time solution of the shallow water equations

Research output: Contribution to journalArticle

2 Citations (Scopus)
33 Downloads (Pure)

Abstract

The shallow water equations are applicable to many common engineering problems involving modelling of waves dominated by motions in the horizontal directions (e.g. tsunami propagation, dam breaks). As such events pose substantial economic costs, as well as potential loss of life, accurate real-time simulation and visualization methods are of great importance. For this purpose, we propose a new finite difference scheme for the 2D shallow water equations that is specifically formulated to take advantage of modern GPUs. The new scheme is based on the so-called Picard integral formulation of conservation laws combined with Weighted Essentially Non-Oscillatory reconstruction. The emphasis of the work is on third order in space and second order in time solutions (in both single and double precision). Further, the scheme is well-balanced for bathymetry functions that are not surface piercing and can handle wetting and drying in a GPU-friendly manner without resorting to long and specific case-by-case procedures. We also present a fast single kernel GPU implementation with a novel boundary condition application technique that allows for simultaneous real-time visualization and single precision simulations even on large ( > 2000 × 2000) grids on consumer-level hardware - the full kernel source codes are also provided online at https://github.com/pparna/swe_pifweno3.
Original languageEnglish
Pages (from-to)107-120
Number of pages14
JournalComputers and Fluids
Volume161
Early online date21 Nov 2017
DOIs
Publication statusPublished - 15 Jan 2018

Fingerprint

Visualization
Water
Bathymetry
Piercing
Tsunamis
Dams
Wetting
Conservation
Drying
Boundary conditions
Hardware
Economics
Graphics processing unit
Costs

Cite this

@article{a7a87c59c09040d08d319828db794551,
title = "GPU driven finite difference WENO scheme for real time solution of the shallow water equations",
abstract = "The shallow water equations are applicable to many common engineering problems involving modelling of waves dominated by motions in the horizontal directions (e.g. tsunami propagation, dam breaks). As such events pose substantial economic costs, as well as potential loss of life, accurate real-time simulation and visualization methods are of great importance. For this purpose, we propose a new finite difference scheme for the 2D shallow water equations that is specifically formulated to take advantage of modern GPUs. The new scheme is based on the so-called Picard integral formulation of conservation laws combined with Weighted Essentially Non-Oscillatory reconstruction. The emphasis of the work is on third order in space and second order in time solutions (in both single and double precision). Further, the scheme is well-balanced for bathymetry functions that are not surface piercing and can handle wetting and drying in a GPU-friendly manner without resorting to long and specific case-by-case procedures. We also present a fast single kernel GPU implementation with a novel boundary condition application technique that allows for simultaneous real-time visualization and single precision simulations even on large ( > 2000 × 2000) grids on consumer-level hardware - the full kernel source codes are also provided online at https://github.com/pparna/swe_pifweno3.",
author = "P. Parna and K. Meyer and R. Falconer",
year = "2018",
month = "1",
day = "15",
doi = "10.1016/j.compfluid.2017.11.012",
language = "English",
volume = "161",
pages = "107--120",
journal = "Computers and Fluids",
issn = "0045-7930",
publisher = "Elsevier Limited",

}

GPU driven finite difference WENO scheme for real time solution of the shallow water equations. / Parna, P.; Meyer, K.; Falconer, R.

In: Computers and Fluids, Vol. 161, 15.01.2018, p. 107-120.

Research output: Contribution to journalArticle

TY - JOUR

T1 - GPU driven finite difference WENO scheme for real time solution of the shallow water equations

AU - Parna, P.

AU - Meyer, K.

AU - Falconer, R.

PY - 2018/1/15

Y1 - 2018/1/15

N2 - The shallow water equations are applicable to many common engineering problems involving modelling of waves dominated by motions in the horizontal directions (e.g. tsunami propagation, dam breaks). As such events pose substantial economic costs, as well as potential loss of life, accurate real-time simulation and visualization methods are of great importance. For this purpose, we propose a new finite difference scheme for the 2D shallow water equations that is specifically formulated to take advantage of modern GPUs. The new scheme is based on the so-called Picard integral formulation of conservation laws combined with Weighted Essentially Non-Oscillatory reconstruction. The emphasis of the work is on third order in space and second order in time solutions (in both single and double precision). Further, the scheme is well-balanced for bathymetry functions that are not surface piercing and can handle wetting and drying in a GPU-friendly manner without resorting to long and specific case-by-case procedures. We also present a fast single kernel GPU implementation with a novel boundary condition application technique that allows for simultaneous real-time visualization and single precision simulations even on large ( > 2000 × 2000) grids on consumer-level hardware - the full kernel source codes are also provided online at https://github.com/pparna/swe_pifweno3.

AB - The shallow water equations are applicable to many common engineering problems involving modelling of waves dominated by motions in the horizontal directions (e.g. tsunami propagation, dam breaks). As such events pose substantial economic costs, as well as potential loss of life, accurate real-time simulation and visualization methods are of great importance. For this purpose, we propose a new finite difference scheme for the 2D shallow water equations that is specifically formulated to take advantage of modern GPUs. The new scheme is based on the so-called Picard integral formulation of conservation laws combined with Weighted Essentially Non-Oscillatory reconstruction. The emphasis of the work is on third order in space and second order in time solutions (in both single and double precision). Further, the scheme is well-balanced for bathymetry functions that are not surface piercing and can handle wetting and drying in a GPU-friendly manner without resorting to long and specific case-by-case procedures. We also present a fast single kernel GPU implementation with a novel boundary condition application technique that allows for simultaneous real-time visualization and single precision simulations even on large ( > 2000 × 2000) grids on consumer-level hardware - the full kernel source codes are also provided online at https://github.com/pparna/swe_pifweno3.

U2 - 10.1016/j.compfluid.2017.11.012

DO - 10.1016/j.compfluid.2017.11.012

M3 - Article

VL - 161

SP - 107

EP - 120

JO - Computers and Fluids

JF - Computers and Fluids

SN - 0045-7930

ER -