Automated parameters for troubled-cell indicators using outlier detection

In Vuik and Ryan (2014) we studied the use of troubled-cell indicators for discontinuity detection in nonlinear hyperbolic partial differential equations and introduced a new multiwavelet technique to detect troubled cells. We found that these methods perform well as long as a suitable, problem-dependent parameter is chosen. This parameter is used in a threshold which decides whether or not to detect an element as a troubled cell. Until now, these parameters could not be chosen automatically. The choice of the parameter has impact on the approximation: it determines the strictness of the troubled-cell indicator. An inappropriate choice of the parameter will result in detection (and limiting) of too few or too many elements. The optimal parameter is chosen such that the minimal number of troubled cells is detected and the resulting approximation is free of spurious oscillations. In this paper we will see that for each troubled-cell indicator the sudden increase or decrease of the indicator value with respect to the neighboring values is important for detection. Indication basically reduces to detecting the outliers of a vector (one dimension) or matrix (two dimensions). This is done using Tukey's boxplot approach to detect which coefficients in a vector are straying far beyond others (Tukey, 1977). We provide an algorithm that can be applied to various troubled-cell indication variables. Using this technique the problem-dependent parameter that the original indicator requires is no longer necessary as the parameter will be chosen automatically

Detecting Discontinuities Over Triangular Meshes Using Multiwavelets

It is well known that solutions to nonlinear hyperbolic PDEs develop discontinuities in time. The generation of spurious oscillations in such regions can be prevented by applying a limiter in the troubled zones. In earlier work, we constructed a multiwavelet troubled-cell indicator for one and (tensor-product) two dimensions (SIAM J. Sci. Comput. 38(1):A84–A104, 2016). In this paper, we investigate multiwavelet troubled-cell indicators on structured triangular meshes. One indicator uses a problem-dependent parameter; the other indicator is combined with outlier detection

Recent Developments in the Numerics of Nonlinear Hyperbolic Conservation Laws

The development of reliable numerical methods for the simulation of real life problems requires both a fundamental knowledge in the field of numerical analysis and a proper experience in practical applications as well as their mathematical modeling. Thus, the purpose of the workshop was to bring together experts not only from the field of applied mathematics but also from civil and mechanical engineering working in the area of modern high order methods for the solution of partial differential equations or even approximation theory necessary to improve the accuracy as well as robustness of numerical algorithms

Shock capturing techniques for hp-adaptive finite elements

The aim of this work is to propose an hp-adaptive algorithm for discontinuous Galerkin methods that is capable to detect the discontinuities and sharp layers and avoid the spurious oscillation of the solution around them. In order to control the spurious oscillations, artificial viscosity is used with the particularity that it is only applied around the layers where the solution changes abruptly. To do so, a novel troubled-cell detector has been developed in order to mark the elements around those layers and to impose linear order in them. The detector takes advantage of the evolution of the value of the gradient through the adaptive process.Peer ReviewedPostprint (published version

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

(Multi)wavelets increase both accuracy and efficiency of standard Godunov-type hydrodynamic models: robust 2D approaches

Multiwavelets (MW) enable the compression, analysis and assembly of model data on a multiresolution grid within Godunov-type solvers based on second-order discontinuous Galerkin (DG2) and first-order finite volume (FV1) methods. Multiwavelet adaptivity has been studied extensively with one-dimensional (1D) hydrodynamic models (Kesserwani et al., 2019), revealing that MWDG2 can be 20 times faster than uniform DG2 and 2 times faster than uniform FV1, while preserving the accuracy and robustness of the underlying formulation. The potential of the MWDG2 scheme has yet to be studied for two-dimensional (2D) modelling, but this requires a design that robustly and efficiently solves the 2D shallow water equations (SWE) with complex source terms and wetting and drying. This paper presents a two-dimensional MWDG2 scheme that: (1) adopts a slope-decoupled DG2 solver as a reference scheme, for its ability to deliver well-balanced piecewise-planar solutions shaped by a simplified 3-component basis; and, (2) adapts the multiresolution analysis of multiwavelets for compatibility with the slope-decoupled DG2 basis. A scaled reformulation of slope-decoupled DG2 is presented alongside two multiwavelet approaches that yield MWDG2 schemes with similar properties, and a Haar wavelet FV1 (HFV1) variant for adapting piecewise-constant model data. The performance of the adaptive HFV1 and MWDG2 solvers is explored alongside their uniform counterparts, while analysing their accuracy, efficiency, grid-coarsening ability, reliability in handling wet-dry fronts across steep bed-slopes, and ability to capture features relevant to practical hydraulic modelling. The results indicate a particular multiwavelet approach that allows the MWDG2 scheme to exploit its grid-coarsening ability for the widest range of flow types. Results also indicate that the proposed (multi)wavelet-based adaptive schemes are even more efficient for the 2D case. Accompanying model software is openly available online

Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2018

This open access book features a selection of high-quality papers from the presentations at the International Conference on Spectral and High-Order Methods 2018, offering an overview of the depth and breadth of the activities within this important research area. The carefully reviewed papers provide a snapshot of the state of the art, while the extensive bibliography helps initiate new research directions