835 research outputs found
A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure
In this paper we formulate and test numerically a fully-coupled discontinuous
Galerkin (DG) method for incompressible two-phase flow with discontinuous
capillary pressure. The spatial discretization uses the symmetric interior
penalty DG formulation with weighted averages and is based on a wetting-phase
potential / capillary potential formulation of the two-phase flow system. After
discretizing in time with diagonally implicit Runge-Kutta schemes the resulting
systems of nonlinear algebraic equations are solved with Newton's method and
the arising systems of linear equations are solved efficiently and in parallel
with an algebraic multigrid method. The new scheme is investigated for various
test problems from the literature and is also compared to a cell-centered
finite volume scheme in terms of accuracy and time to solution. We find that
the method is accurate, robust and efficient. In particular no post-processing
of the DG velocity field is necessary in contrast to results reported by
several authors for decoupled schemes. Moreover, the solver scales well in
parallel and three-dimensional problems with up to nearly 100 million degrees
of freedom per time step have been computed on 1000 processors
Python framework for HP adaptive discontinuous Galerkin methods for two phase flow in porous media
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 a fully coupled implicit scheme. The implementation of the discretization is very flexible allowing to test different formulations of the two-phase flow model and adaptation strategies
Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures
International audienceWe design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters
Python Framework for HP Adaptive Discontinuous Galerkin Method for Two Phase Flow in Porous Media
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.Comment: Keywords: DG, hp-adaptivity, Two-phase flow, IMPES, Fully implicit,
Dune, Python, Porous media. 28 pages, 9 figures, various code snippet
Mixed and discontinuous finite volume element schemes for the optimal control of immiscible flow in porous media
We introduce a family of hybrid discretisations for the numerical
approximation of optimal control problems governed by the equations of
immiscible displacement in porous media. The proposed schemes are based on
mixed and discontinuous finite volume element methods in combination with the
optimise-then-discretise approach for the approximation of the optimal control
problem, leading to nonsymmetric algebraic systems, and employing minimum
regularity requirements. Estimates for the error (between a local reference
solution of the infinite dimensional optimal control problem and its hybrid
approximation) measured in suitable norms are derived, showing optimal orders
of convergence
- …