Andreas Dedner

Andreas Dedner, Birane Kane, Robert Klöfkorn et al.

In this paper we present a framework for solving two phase flow problems in porous media. The discretization is based on a Discontinuous Galerkin method and includes local grid adaptivity and local choice of polynomial degree. The method is implemented using the new Python frontend Dune-FemPy to the open source framework Dune. The code used for the simulations is made available as Jupyter notebook and can be used through a Docker container. We present a number of time stepping approaches ranging from a classical IMPES method to fully coupled implicit scheme. The implementation of the discretization is very flexible allowing for test different formulations of the two phase flow model and adaptation strategies.

Luke Benfield, Andreas Dedner

The solution of partial differential equations (PDEs) on complex domains often presents a significant computational challenge by requiring the generation of fitted meshes. The Diffuse Domain Method (DDM) is an alternative which reformulates the problem on a larger, simple domain where the complex geometry is represented by a smooth phase-field function. This paper introduces and analyses several new DDM methods for solving problems with Dirichlet boundary conditions. We derive two new methods from the mixed formulation of the governing equations. This approach transforms the essential Dirichlet conditions into natural boundary conditions. Additionally, we develop coercive formulations based on Nitsche's method, and provide proofs of coercivity for all new and key existing approximations. Numerical experiments demonstrate the improved accuracy of the new methods, and reveal the balance between $L^2$ and $H^1$ errors. The practical effectiveness of this approach is demonstrated through the simulation of the incompressible Navier-Stokes equations on a benchmark fluid dynamics problems.

Andreas Dedner, Pravin Madhavan, Björn Stinner

We extend the discontinuous Galerkin (DG) framework to a linear second-order elliptic problem on a compact smooth connected and oriented surface. An interior penalty (IP) method is introduced on a discrete surface and we derive a-priori error estimates by relating the latter to the original surface via the lift introduced in Dziuk (1988). The estimates suggest that the geometric error terms arising from the surface discretisation do not affect the overall convergence rate of the IP method when using linear ansatz functions. This is then verified numerically for a number of test problems. An intricate issue is the approximation of the surface conormal required in the IP formulation, choices of which are investigated numerically. Furthermore, we present a generic implementation of test problems on surfaces.

Chinedu Nwaigwe, Andreas Dedner

One dimensional (1D) simulations of the flow and flooding of open channels are known to be inaccurate as the flow is multi-dimensional in nature, especially at the flooded regions. However, multi-dimensional simulations, even in two dimensions (2D), are computationally expensive, hence the problem of efficiently coupling 2D and 1D simulations for the flow and flooding of open channels has been the subject of much research and is investigated in this paper. We adopt a 1D model with coupling term for the channel flow and the 2D shallow water flow model for the floodplain. The 1D model with coupling term is derived by integrating the 3D Free Surface Euler equations but without imposing any restriction on the channel width variations. %leading to a new coupling term. Finite volume methods are formulated for both the 2D and 1D models including a discrete coupling term in closed form. Coupling is achieved through the discrete coupling term in the 1D model and the lateral numerical fluxes in the 2D model. Since the lateral discharge in the channel cannot be guaranteed to be zero during flooding, we aim to recover the lateral variation by computing two lateral discharges over each cross section and propose to use an ad-hoc model based on the $y$-discharge equation in the 2D model for this purpose. We then propose the numerical scheme for this ad-hoc model following the hydrostatic reconstruction philosophy. Then, we show that the resulting method, named Horizontal Coupling Method (HCM), is well-balanced; we introduce the no-numerical flooding property and also show that the method satisfies the property. Three numerical test cases are used to verify the performance of the method. The results show that the method performs well in both accuracy and efficiency and also approximates the channel lateral discharges with very good accuracy and little computational overhead.

Andreas Dedner, Chinedu Nwaigwe

In this paper, we propose a novel approach for coupling 2D/1D shallow water flow models. Efficiently coupling these models is vital for simulating the flow and flooding of open channels. Currently, existing methods couple the models either at the channel lateral boundaries (lateral methods) or at the location, along the channel flow direction, where the two sub-domains intersect (frontal methods). We classify these methods as horizontal methods since the coupling points are on the horizontal plane. The limitations of these methods include their inability to recover the 2D channel flow structure during flooding without losing efficiency, and their inability to switch between lateral and frontal methods, based on the problem. Here, we propose a new paradigm which is to think of the channel flow as a two layer flow. This leads to the vertical coupling method (VCM) which is able to recover 2D channel flow structure without losing efficiency, switch types depending on the flow and is a superset of some existing methods. The VCM is based on a user chosen elevation above which the user considers the channel to be full. From this elevation, all other flow quantities are defined including the two layers. The flows in the lower and upper layers are assumed to be 1D and 2D respectively, and the appropriate flow models with exchange terms are derived. The 2D shallow water flow model is retained for the floodplains. Finite volume methods (FVM) are formulated for the flood model; the FVM with the operator splitting approach are also formulated to solve and couple the two layer channel flow models. We show that the resulting method (i) is well-balanced (ii) preserves the "no-numerical" flooding property (iii) preserves conservation properties and (iv) adapts to the flow situation.

Andreas Dedner, Jan Giesselmann

We present reliable a posteriori estimators for some fully discrete schemes applied to nonlinear systems of hyperbolic conservation laws in one space dimension with strictly convex entropy. The schemes are based on a method of lines approach combining discontinuous Galerkin spatial discretization with single- or multi-step methods in time. The construction of the estimators requires a reconstruction in time for which we present a very general framework first for odes and then apply the approach to conservation laws. The reconstruction does not depend on the actual method used for evolving the solution in time. Most importantly it covers in addition to implicit methods also the wide range of explicit methods typically used to solve conservation laws. For the spatial discretization, we allow for standard choices of numerical fluxes. We use reconstructions of the discrete solution together with the relative entropy stability framework, which leads to error control in the case of smooth solutions. We study under which conditions on the numerical flux the estimate is of optimal order pre-shock. While the estimator we derive is computable and valid post-shock for fixed meshsize, it will blow up as the meshsize tends to zero. This is due to a breakdown of the relative entropy framework when discontinuities develop. We conclude with some numerical benchmarking to test the robustness of the derived estimator.

Andreas Dedner, Pravin Madhavan

We extend the discontinuous Galerkin (DG) framework to the analysis of first-order hyperbolic and advection-dominated problems posed on implicitly defined surfaces. The focus will be on the hyperbolic part, which is discretised using a "discrete surface" generalisation of the jump-stabilised upwind flux. A key issue arising in the analysis (which does not appear in the planar setting) is the treatment of the discrete velocity field, choices of which play an important role in the stability of the scheme. We then prove optimal error estimates in an appropriate norm given a number of assumptions on the discrete velocity field, which are then investigated and discussed in more detail. The theoretical results are verified numerically for a number of test problems exhibiting advection-dominated behaviour.

Andreas Dedner, Tristan Pryer

We extend the finite element method introduced by Lakkis and Pryer [2011] to approximate the solution of second order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin (DG) framework. This is done by viewing the NVFEM as a mixed method whereby the finite element Hessian is an auxiliary variable in the formulation. Representing the finite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems. Furthermore, the system matrix is very easy to assemble, Thus this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach. We conduct a stability and consistency analysis making use of the unified framework set out in Arnold et. al. [2001]. We also give an apriori analysis of the method. The analysis applies to any consistent representation of the finite element Hessian, thus is applicable to the previous works making use of continuous Galerkin approximations.