Frame-invariant directional vector limiters for discontinuous Galerkin methods

Second and higher order numerical approximations of conservation laws for vector fields call for the use of limiting techniques based on generalized monotonicity criteria. In this paper, we introduce a family of directional vertexbased slope limiters for tensor-valued gradients of formally second-order accurate piecewise-linear discontinuous Galerkin (DG) discretizations. The proposed methodology enforces local maximum principles for scalar products corresponding to projections of a vector field onto the unit vectors of a frame-invariant orthogonal basis. In particular, we consider anisotropic limiters based on singular value decompositions and the Gram-Schmidt orthogonalization procedure. The proposed extension to hyperbolic systems features a sequential limiting strategy and a global invariant domain fix. The pros and cons of different approaches to vector limiting are illustrated by the results of numerical studies for the two-dimensional shallow water equations and for the Euler equations of gas dynamics

Embedded discontinuous Galerkin transport schemes with localised limiters

Motivated by finite element spaces used for representation of temperature in the compatible finite element approach for numerical weather prediction, we introduce locally bounded transport schemes for (partially-)continuous finite element spaces. The underlying high-order transport scheme is constructed by injecting the partially-continuous field into an embedding discontinuous finite element space, applying a stable upwind discontinuous Galerkin (DG) scheme, and projecting back into the partially-continuous space; we call this an embedded DG scheme. We prove that this scheme is stable in L2 provided that the underlying upwind DG scheme is. We then provide a framework for applying limiters for embedded DG transport schemes. Standard DG limiters are applied during the underlying DG scheme. We introduce a new localised form of element-based flux-correction which we apply to limiting the projection back into the partially-continuous space, so that the whole transport scheme is bounded. We provide details in the specific case of tensor-product finite element spaces on wedge elements that are discontinuous P1/Q1 in the horizontal and continuous P2 in the vertical. The framework is illustrated with numerical tests

Das unstetige Galerkinverfahren für Strömungen mit freier Oberfläche und im Grundwasserbereich in geophysikalischen Anwendungen

Free surface flows and subsurface flows appear in a broad range of geophysical applications and in many environmental settings situations arise which even require the coupling of free surface and subsurface flows. Many of these application scenarios are characterized by large domain sizes and long simulation times. Hence, they need considerable amounts of computational work to achieve accurate solutions and the use of efficient algorithms and high performance computing resources to obtain results within a reasonable time frame is mandatory. Discontinuous Galerkin methods are a class of numerical methods for solving differential equations that share characteristics with methods from the finite volume and finite element frameworks. They feature high approximation orders, offer a large degree of flexibility, and are well-suited for parallel computing. This thesis consists of eight articles and an extended summary that describe the application of discontinuous Galerkin methods to mathematical models including free surface and subsurface flow scenarios with a strong focus on computational aspects. It covers discretization and implementation aspects, the parallelization of the method, and discrete stability analysis of the coupled model.Für viele geophysikalische Anwendungen spielen Strömungen mit freier Oberfläche und im Grundwasserbereich oder sogar die Kopplung dieser beiden eine zentrale Rolle. Oftmals charakteristisch für diese Anwendungsszenarien sind große Rechengebiete und lange Simulationszeiten. Folglich ist das Berechnen akkurater Lösungen mit beträchtlichem Rechenaufwand verbunden und der Einsatz effizienter Lösungsverfahren sowie von Techniken des Hochleistungsrechnens obligatorisch, um Ergebnisse innerhalb eines annehmbaren Zeitrahmens zu erhalten. Unstetige Galerkinverfahren stellen eine Gruppe numerischer Verfahren zum Lösen von Differentialgleichungen dar, und kombinieren Eigenschaften von Methoden der Finiten Volumen- und Finiten Elementeverfahren. Sie ermöglichen hohe Approximationsordnungen, bieten einen hohen Grad an Flexibilität und sind für paralleles Rechnen gut geeignet. Diese Dissertation besteht aus acht Artikeln und einer erweiterten Zusammenfassung, in diesen die Anwendung unstetiger Galerkinverfahren auf mathematische Modelle inklusive solcher für Strömungen mit freier Oberfläche und im Grundwasserbereich beschrieben wird. Die behandelten Themen umfassen Diskretisierungs- und Implementierungsaspekte, die Parallelisierung der Methode sowie eine diskrete Stabilitätsanalyse des gekoppelten Modells

HIGH ORDER BOUND-PRESERVING DISCONTINUOUS GALERKIN METHODS AND THEIR APPLICATIONS IN PETROLEUM ENGINEERING

This report contains researches in the theory of high-order bound-preserving (BP) discontinuous Galerkin (DG) method and their applications in petroleum engineering. It contains both theoretical analysis and numerical experiments. The compressible miscible displacements and wormhole propagation problem, arising in petroleum engineering, is used to describe the evolution of the pressure and concentrations of different components of fluid in porous media. The important physical features of concentration and porosity include their boundedness between 0 and 1, as well as the monotone increasing for porosity in wormhole propagation model. How to keep these properties in the simulation is crucial to the robustness of the numerical algorithm. In the first project, we develop high-order bound-preserving discontinuous Galerkin methods for the coupled system of compressible miscible displacements on triangular meshes. We consider the problem with multi-component fluid mixture and the (volumetric) concentration of the jth component,cj, should be between 0 and 1. The main idea is stated as follows. First, we apply the second-order positivity-preserving techniques to all concentrations c′ js and enforce P jcj= 1 simultaneously to obtain physically relevant boundedness for every components. Then, based on the second-order BP schemes, we use the second-order numerical fluxes as the lower order one to combine with high-order numerical fluxes to achieve the high-order accuracy. Finally, since the classical slope limiter cannot be applied to polynomial upper bounds, we introduce a new limiter to our algorithm. Numerical experiments are given to demonstrate the high-order accuracy and good performance of the numerical technique. In our second project, we propose high-order bound-preserving discontinuous Galerkin methods to keep the boundedness for the porosity and concentration of acid, as well as the monotone increasing for porosity. The main technique is to introduce a new variable r to replace the original acid concentration and use a consistent flux pair to deduce a ghost equation such that the positive-preserving technique can be applied on both original and deduced equations. A high-order slope limiter is used to keep a polynomial upper bound which changes over time for r. Moreover, the high-order accuracy is attained by the flux limiter. Numerical examples are given to demonstrate the high-order accuracy and bound-preserving property of the numerical technique

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