Second-order discontinuous Galerkin flood model: comparison with industry-standard finite volume models

Finite volume (FV) numerical solvers of the two-dimensional shallow water equations are core to industry-standard flood models. The second-order Discontinuous Galerkin (DG) alternative is well-known to perform better than first- and second-order FV to capture sharp flow fronts and converge faster at coarser resolutions, but DG2 models typically rely on local slope limiting to selectively damp numerical oscillations in the vicinity of shock waves. Yet flood inundation events are smooth and gradually-varying, and shock waves play only a minor role in flood inundation modelling. Therefore, this paper investigates two DG2 variants - with and without local slope limiting - to identify the simplest and most efficient DG2 configuration suitable for flood inundation modelling. The predictive capabilities of the DG2 variants are analysed for a synthetic test case involving advancing and receding waves representative of flood-like flow. The DG2 variants are then benchmarked against industry-standard FV models over six UK Environment Agency scenarios. Results indicate that the DG2 variant without local slope limiting closely reproduces solutions of the commercial models at twice as coarse a spatial resolution, and removing the slope limiter can halve model runtime. Results also indicate that DG2 can capture more accurate hydrographs incorporating small-scale transients over long-range simulations, even when hydrographs are measured far away from the flooding source. Accompanying details of software and data accessibility are provided

Well-resolved velocity fields using discontinuous Galerkin shallow water solutions

Computational models based on the depth-averaged shallow water equations (SWE) offer an efficient choice to analyse velocity fields around hydraulic structures. Second-order finite volume (FV2) solvers have often been used for this purpose subject to adding an eddy viscosity term at sub-meter resolution, but have been shown to fall short of capturing small-scale field transients emerging from wave-structure interactions. The second-order discontinuous Galerkin (DG2) alternative is significantly more resistant to the growth of numerical diffusion and leads to faster convergence rates. These properties make the DG2 solver a promising modelling tool for detailed velocity field predictions. This paper focuses on exploring this DG2 capability with reference to an FV2 counterpart for a selection of test cases that require well-resolved velocity field predictions. The findings of this work lead to identifying a particular setting for the DG2 solver that allows for obtaining more accurate and efficient depth-averaged velocity fields incorporating small-scale transients

Shallow-flow velocity predictions using discontinuous Galerkin solutions

Numerical solvers of the two-dimensional (2D) shallow water equations (2D-SWE) can be an efficient option to predict spatial distribution of velocity fields in quasi-steady flows past or throughout hydraulic engineering structures. A second-order finite-volume (FV2) solver spuriously elongates small-scale recirculating eddies within its predictions, unless sustained by an artificial eddy viscosity, while a third-order finite-volume (FV3) solver can distort the eddies within its predictions. The extra complexity in a second-order discontinuous Galerkin (DG2) solver leads to significantly reduced error dissipation and improved predictions at a coarser resolution, making it a viable contender to acquire velocity predictions in shallow flows. This paper analyses this predictive capability for a grid-based, open source DG2 solver with reference to FV2 or FV3 solvers for simulating velocity magnitude and direction at the submeter scale. The simulated predictions are assessed against measured velocity data for four experimental test cases. The results consistently indicate that the DG2 solver is a competitive choice to efficiently produce more accurate velocity distributions for the simulations dominated by smooth flow regions