LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs

<p>LISFLOOD-FP 8.0 includes second-order discontinuous Galerkin (DG2) and first-order finite-volume (FV1) solvers of the two-dimensional shallow-water equations for modelling a wide range of flows, including rapidly propagating, supercritical flows, shock waves or flows over very smooth surfac...

Full description

Bibliographic Details
Main Authors: J. Shaw, G. Kesserwani, J. Neal, P. Bates, M. K. Sharifian
Format: Article
Language:English
Published: Copernicus Publications 2021-06-01
Series:Geoscientific Model Development
Online Access:https://gmd.copernicus.org/articles/14/3577/2021/gmd-14-3577-2021.pdf
_version_ 1831544479966822400
author J. Shaw
G. Kesserwani
J. Neal
P. Bates
M. K. Sharifian
author_facet J. Shaw
G. Kesserwani
J. Neal
P. Bates
M. K. Sharifian
author_sort J. Shaw
collection DOAJ
description <p>LISFLOOD-FP 8.0 includes second-order discontinuous Galerkin (DG2) and first-order finite-volume (FV1) solvers of the two-dimensional shallow-water equations for modelling a wide range of flows, including rapidly propagating, supercritical flows, shock waves or flows over very smooth surfaces. The solvers are parallelised on multi-core CPU and Nvidia GPU architectures and run existing LISFLOOD-FP modelling scenarios without modification. These new, fully two-dimensional solvers are available alongside the existing local inertia solver (called ACC), which is optimised for multi-core CPUs and integrates with the LISFLOOD-FP sub-grid channel model. The predictive capabilities and computational scalability of the new DG2 and FV1 solvers are studied for two Environment Agency benchmark tests and a real-world fluvial flood simulation driven by rainfall across a 2500 <span class="inline-formula">km<sup>2</sup></span> catchment. DG2's second-order-accurate, piecewise-planar representation of topography and flow variables enables predictions on coarse grids that are competitive with FV1 and ACC predictions on 2–4 times finer grids, particularly where river channels are wider than half the grid spacing. Despite the simplified formulation of the local inertia solver, ACC is shown to be spatially second-order-accurate and yields predictions that are close to DG2. The DG2-CPU and FV1-CPU solvers achieve near-optimal scalability up to 16 CPU cores and achieve greater efficiency on grids with fewer than 0.1 million elements. The DG2-GPU and FV1-GPU solvers are most efficient on grids with more than 1 million elements, where the GPU solvers are 2.5–4 times faster than the corresponding 16-core CPU solvers. LISFLOOD-FP 8.0 therefore marks a new step towards operational DG2 flood inundation modelling at the catchment scale. LISFLOOD-FP 8.0 is freely available under the GPL v3 license, with additional documentation and case studies at <span class="uri">https://www.seamlesswave.com/LISFLOOD8.0</span> (last access: 2 June 2021).</p>
first_indexed 2024-12-17T01:09:57Z
format Article
id doaj.art-8c5193523bac4d9c8e31398ee9fa3cfb
institution Directory Open Access Journal
issn 1991-959X
1991-9603
language English
last_indexed 2024-12-17T01:09:57Z
publishDate 2021-06-01
publisher Copernicus Publications
record_format Article
series Geoscientific Model Development
spelling doaj.art-8c5193523bac4d9c8e31398ee9fa3cfb2022-12-21T22:09:10ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032021-06-01143577360210.5194/gmd-14-3577-2021LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUsJ. Shaw0G. Kesserwani1J. Neal2P. Bates3M. K. Sharifian4Department of Civil and Structural Engineering, The University of Sheffield, Western Bank, Sheffield, UKDepartment of Civil and Structural Engineering, The University of Sheffield, Western Bank, Sheffield, UKSchool of Geographical Sciences, University of Bristol, Bristol, UKSchool of Geographical Sciences, University of Bristol, Bristol, UKDepartment of Civil and Structural Engineering, The University of Sheffield, Western Bank, Sheffield, UK<p>LISFLOOD-FP 8.0 includes second-order discontinuous Galerkin (DG2) and first-order finite-volume (FV1) solvers of the two-dimensional shallow-water equations for modelling a wide range of flows, including rapidly propagating, supercritical flows, shock waves or flows over very smooth surfaces. The solvers are parallelised on multi-core CPU and Nvidia GPU architectures and run existing LISFLOOD-FP modelling scenarios without modification. These new, fully two-dimensional solvers are available alongside the existing local inertia solver (called ACC), which is optimised for multi-core CPUs and integrates with the LISFLOOD-FP sub-grid channel model. The predictive capabilities and computational scalability of the new DG2 and FV1 solvers are studied for two Environment Agency benchmark tests and a real-world fluvial flood simulation driven by rainfall across a 2500 <span class="inline-formula">km<sup>2</sup></span> catchment. DG2's second-order-accurate, piecewise-planar representation of topography and flow variables enables predictions on coarse grids that are competitive with FV1 and ACC predictions on 2–4 times finer grids, particularly where river channels are wider than half the grid spacing. Despite the simplified formulation of the local inertia solver, ACC is shown to be spatially second-order-accurate and yields predictions that are close to DG2. The DG2-CPU and FV1-CPU solvers achieve near-optimal scalability up to 16 CPU cores and achieve greater efficiency on grids with fewer than 0.1 million elements. The DG2-GPU and FV1-GPU solvers are most efficient on grids with more than 1 million elements, where the GPU solvers are 2.5–4 times faster than the corresponding 16-core CPU solvers. LISFLOOD-FP 8.0 therefore marks a new step towards operational DG2 flood inundation modelling at the catchment scale. LISFLOOD-FP 8.0 is freely available under the GPL v3 license, with additional documentation and case studies at <span class="uri">https://www.seamlesswave.com/LISFLOOD8.0</span> (last access: 2 June 2021).</p>https://gmd.copernicus.org/articles/14/3577/2021/gmd-14-3577-2021.pdf
spellingShingle J. Shaw
G. Kesserwani
J. Neal
P. Bates
M. K. Sharifian
LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs
Geoscientific Model Development
title LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs
title_full LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs
title_fullStr LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs
title_full_unstemmed LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs
title_short LISFLOOD-FP 8.0: the new discontinuous Galerkin shallow-water solver for multi-core CPUs and GPUs
title_sort lisflood fp 8 0 the new discontinuous galerkin shallow water solver for multi core cpus and gpus
url https://gmd.copernicus.org/articles/14/3577/2021/gmd-14-3577-2021.pdf
work_keys_str_mv AT jshaw lisfloodfp80thenewdiscontinuousgalerkinshallowwatersolverformulticorecpusandgpus
AT gkesserwani lisfloodfp80thenewdiscontinuousgalerkinshallowwatersolverformulticorecpusandgpus
AT jneal lisfloodfp80thenewdiscontinuousgalerkinshallowwatersolverformulticorecpusandgpus
AT pbates lisfloodfp80thenewdiscontinuousgalerkinshallowwatersolverformulticorecpusandgpus
AT mksharifian lisfloodfp80thenewdiscontinuousgalerkinshallowwatersolverformulticorecpusandgpus