3,120 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
Numerical analysis for the pure Neumann control problem using the gradient discretisation method
The article discusses the gradient discretisation method (GDM) for
distributed optimal control problems governed by diffusion equation with pure
Neumann boundary condition. Using the GDM framework enables to develop an
analysis that directly applies to a wide range of numerical schemes, from
conforming and non-conforming finite elements, to mixed finite elements, to
finite volumes and mimetic finite differences methods. Optimal order error
estimates for state, adjoint and control variables for low order schemes are
derived under standard regularity assumptions. A novel projection relation
between the optimal control and the adjoint variable allows the proof of a
super-convergence result for post-processed control. Numerical experiments
performed using a modified active set strategy algorithm for conforming,
nonconforming and mimetic finite difference methods confirm the theoretical
rates of convergence
Unified convergence analysis of numerical schemes for a miscible displacement problem
This article performs a unified convergence analysis of a variety of
numerical methods for a model of the miscible displacement of one
incompressible fluid by another through a porous medium. The unified analysis
is enabled through the framework of the gradient discretisation method for
diffusion operators on generic grids. We use it to establish a novel
convergence result in of the approximate
concentration using minimal regularity assumptions on the solution to the
continuous problem. The convection term in the concentration equation is
discretised using a centred scheme. We present a variety of numerical tests
from the literature, as well as a novel analytical test case. The performance
of two schemes are compared on these tests; both are poor in the case of
variable viscosity, small diffusion and medium to small time steps. We show
that upstreaming is not a good option to recover stable and accurate solutions,
and we propose a correction to recover stable and accurate schemes for all time
steps and all ranges of diffusion
- …