CEMay 1, 2018Code
Python Framework for HP Adaptive Discontinuous Galerkin Method for Two Phase Flow in Porous MediaAndreas 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.
NAJan 23, 2017
Constrained Reconstruction in MUSCL-type Finite Volume SchemesChristoph Gersbacher, Martin Nolte
In this paper we are concerned with the stabilization of MUSCL-type finite volume schemes in arbitrary space dimensions. We consider a number of limited reconstruction techniques which are defined in terms inequality-constrained linear or quadratic programming problems on individual grid elements. No restrictions to the conformity of the grid or the shape of its elements are made. In the special case of Cartesian meshes a novel QP reconstruction is shown to coincide with the widely used Minmod reconstruction. The accuracy and overall efficiency of the stabilized second-order finite volume schemes is supported by numerical experiments.