Discontinuous galerkin flood model formulation: Luxury or necessity?

The finite volume Godunov-type flood model formulation is the most comprehensive amongst those currently employed for flood risk modeling. The local Discontinuous Galerkin method constitutes a more complex, rigorous, and extended local Godunov-type formulation. However, the practical merit associated with such an increase in the level of complexity of the formulation is yet to be decided. This work makes the case for a second-order Runge-Kutta Discontinuous Galerkin (RKDG2) formulation and contrasts it with the equivalently accurate finite volume (MUSCL) formulation, both of which solve the Shallow Water Equations (SWE) in two space dimensions. The numerical complexity of both formulations are presented and their capabilities are explored for wide-ranging diagnostic and real-scale tests, incorporating all challenging features relevant to flood inundation modeling. Our findings reveal that the extra complexity associated with the RKDG2 model pays off by providing higher-quality solution behavior on very coarse meshes and improved velocity predictions. The practical implication of this is that improved accuracy for flood modeling simulations will result when terrain data are limited or of a low resolution. © 2014. American Geophysical Union

RKDG2 shallow-water solver on non-uniform grids with local time steps: Application to 1D and 2D hydrodynamics

This paper investigates local time stepping (LTS) with the RKDG2 (second-order Runge–Kutta Discontinuous Galerkin) non-uniform solutions of the inhomogeneous SWEs (shallow water equations) with source terms. A LTS algorithm – recently designed for homogenous hyperbolic PDE(s) – is herein reconsidered and improved in combination with the RKDG2 shallow-flow solver (LTS-RKDG2) including topography and friction source terms as well as wetting and drying. Two LTS-RKDG2 schemes that adapt 3 and 4 levels of LTSs are configured on 1D and/or 2D (quadrilateral) non-uniform meshes that, respectively, adopt 3 and 4 scales of spatial discretization. Selected shallow water benchmark tests are used to verify, assess and compare the LTS-RKDG2 schemes relative to their conventional Global Time Step RKDG2 alternatives (GTS-RKDG2) considering several issues of practical relevance to hydraulic modelling. Results show that the LTS-RKDG2 models could offer (depending on both the mesh setting and the features of the flow) comparable accuracy to the associated GTS-RKDG2 models with a savings in runtime of up to a factor of 2.5 in 1D simulations and 1.6 in 2D simulations

Benchmarking a multiresolution discontinuous Galerkin shallow water model: Implications for Computational hydraulics

Numerical modelling of wide ranges of different physical scales, which are involved in Shallow Water (SW) problems, has been a key challenge in computational hydraulics. Adaptive meshing techniques have been commonly coupled with numerical methods in an attempt to address this challenge. The combination of MultiWavelets (MW) with the Runge-Kutta Discontinuous Galerkin (RKDG) method offers a new philosophy to readily achieve mesh adaptivity driven by the local variability of the numerical solution, and without requiring more than one threshold value set by the user. However, the practical merits and implications of the MWRKDG, in terms of how far it contributes to address the key challenge above, are yet to be explored. This work systematically explores this, through the verification and validation of the MWRKDG for selected steady and transient benchmark tests, which involves the features of real SW problems. Our findings reveal a practical promise of the SW-MWRKDG solver, in terms of efficient and accurate mesh-adaptivity, but also suggest further improvement in the SWRKDG reference scheme to better intertwine with, and harness the prowess of, the MW-based adaptivity

Haar wavelet-based adaptive finite volume shallow water solver

This paper presents the formulation of an adaptive finite volume (FV) model for the shallow water equations. A Godunov-type reformulation combining the Haar wavelet is achieved to enable solutiondriven resolution adaptivity (both coarsening and refinement) by depending on the wavelet’s threshold value. The ability to properly model irregular topographies and wetting/drying are transferred from the (baseline) FV uniform mesh model, with no extra notable efforts. Selected hydraulic tests are employed to analyse the performance of the Haar wavelet FV shallow water solver considering adaptivity and practical issues including choice for the threshold value driving the adaptivity, mesh convergence study, shock and wet/dry front capturing abilities. Our findings show that Haar wavelet-based adaptive FV solutions offer great potential to improve the reliability of multiscale shallow water models

A discontinuous Galerkin approach for conservative modelling of fully nonlinear and weakly dispersive wave transformations

This work extends a robust second-order Runge-Kutta Discontinuous Galerkin (RKDG2) method to solve the fully nonlinear and weakly dispersive flows, within a scope to simultaneously address accuracy, conservativeness, cost-efficiency and practical needs. The mathematical model governing such flows is based on a variant form of the Green-Naghdi (GN) equations decomposed as a hyperbolic shallow water system with an elliptic source term. Practical features of relevance (i.e. conservative modelling over irregular terrain with wetting and drying and local slope limiting) have been restored from an RKDG2 solver to the Nonlinear Shallow Water (NSW) equations, alongside new considerations to integrate elliptic source terms (i.e. via a fourth-order local discretization of the topography) and to enable local capturing of breaking waves (i.e. via adding a detector for switching off the dispersive terms). Numerical results are presented, demonstrating the overall capability of the proposed approach in achieving realistic prediction of nearshore wave processes involving both nonlinearity and dispersion effects within a single model

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

Discontinuous Galerkin formulation for 2D hydrodynamic modelling: trade-offs between theoretical complexity and practical convenience

In the modelling of hydrodynamics, the Discontinuous Galerkin (DG) approach constitutes a more complex and modern alternative to the well-established finite volume method. The latter retains some desired practical features for modelling hydrodynamics, such as well-balancing between spatial flux and steep topography gradients, ability to incorporate wetting and drying processes, and computational affordability. In this context, DG methods were originally devised to solve the two-dimensional (2D) Shallow Water Equations (SWE) with irregular topographies and wetting and drying, albeit at reduction in the formulation’s complexity to often being second-order accurate (DG2). The aims of this paper are: (a) to outline a so-called “slope-decoupled” formulation of a standard 2D-DG2-SWE simulator in which theoretical complexity is deliberately reduced; (b) to highlight the capabilities of the proposed slopedecoupled simulator in providing a setting where the simplifying assumptions are verified within the formulation. Both the standard and the slope-decoupled 2D-DG2-SWE models adopt 2D modal basis functions for shaping local planar DG2 solutions on quadrilateral elements, by using an average coefficient and two slope coefficients along the Cartesian coordinates. Over a quadrilateral element, the stencil of the slope-decoupled 2D-DG2 formulation is simplified to remove the interdependence of slope-coefficients for both flow and topography approximations. The fully well-balanced character the slope-decoupled 2D-DG2-SWE planar solutions is theoretically studied. The performance of the latter is compared with the standard 2D-DG2 formulation in classical simulation tests. Other tests are conducted to diagnostically verify the conservative properties of the 2D-DG2-SWE method in scenarios involving sharp topography gradients and wet and/or dry zones. The analyses conducted offer strong evidence that the proposed slope-decoupled 2D-DG2-SWE simulator is very attractive for the development of robust flood models

(Multi)wavelets increase both accuracy and efficiency of standard Godunov-type hydrodynamic models

This paper presents a scaled reformulation of a robust second-order Discontinuous Galerkin (DG2) solver for the Shallow Water Equations (SWE), with guiding principles on how it can be naturally extended to fit into the multiresolution analysis of multiwavelets (MW). Multiresolution analysis applied to the flow and topography data enables the creation of an adaptive MWDG2 solution on a non-uniform grid. The multiresolution analysis also permits control of the adaptive model error by a single user-prescribed parameter. This results in an adaptive MWDG2 solver that can fully exploit the local (de)compression of piecewise-linear modelled data, and from which a first-order finite volume version (FV1) is directly obtainable based on the Haar wavelet (HFV1) for local (de)compression of piecewise-constant modelled data. The behaviour of the adaptive HFV1 and MWDG2 solvers is systematically studied on a number of well-known hydraulic tests that cover all elementary aspects relevant to accurate, efficient and robust modelling. The adaptive solvers are run starting from a baseline mesh with a single element, and their accuracy and efficiency are measured referring to standard FV1 and DG2 simulations on the uniform grid involving the finest resolution accessible by the adaptive solvers. Our findings reveal that the MWDG2 solver can achieve the same accuracy as the DG2 solver but with a greater efficiency than the FV1 solver due to the smoothness of its piecewise-linear basis, which enables more aggressive coarsening than with the piecewise-constant basis in the HFV1 solver. This suggests a great potential for the MWDG2 solver to efficiently handle the depth and breadth in resolution variability, while also being a multiresolution mesh generator. Accompanying model software and simulation data are openly available online

Multiwavelet-based grid adaptation with discontinuous Galerkin schemes for shallow water equations

We provide an adaptive strategy for solving shallow water equations with dynamic grid adaptation including a sparse representation of the bottom topography. A challenge in computing approximate solutions to the shallow water equations including wetting and drying is to achieve the positivity of the water height and the well-balancing of the approximate solution. A key property of our adaptive strategy is that it guarantees that these properties are preserved during the refinement and coarsening steps in the adaptation process.The underlying idea of our adaptive strategy is to perform a multiresolution analysis using multiwavelets on a hierarchy of nested grids. This provides difference information between successive refinement levels that may become negligibly small in regions where the solution is locally smooth. Applying hard thresholding the data are highly compressed and local grid adaptation is triggered by the remaining significant coefficients. Furthermore we use the multiresolution analysis of the underlying data as an additional indicator of whether the limiter has to be applied on a cell or not. By this the number of cells where the limiter is applied is reduced without spoiling the accuracy of the solution.By means of well-known 1D and 2D benchmark problems, we verify that multiwavelet-based grid adaptation can significantly reduce the computational cost by sparsening the computational grids, while retaining accuracy and keeping well-balancing and positivity