371 research outputs found
Space Decompositions and Solvers for Discontinuous Galerkin Methods
We present a brief overview of the different domain and space decomposition
techniques that enter in developing and analyzing solvers for discontinuous
Galerkin methods. Emphasis is given to the novel and distinct features that
arise when considering DG discretizations over conforming methods. Connections
and differences with the conforming approaches are emphasized.Comment: 2 pages 2 figures no table
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
Multilevel Preconditioning of Discontinuous-Galerkin Spectral Element Methods, Part I: Geometrically Conforming Meshes
This paper is concerned with the design, analysis and implementation of
preconditioning concepts for spectral Discontinuous Galerkin discretizations of
elliptic boundary value problems. While presently known techniques realize a
growth of the condition numbers that is logarithmic in the polynomial degrees
when all degrees are equal and quadratic otherwise, our main objective is to
realize full robustness with respect to arbitrarily large locally varying
polynomial degrees degrees, i.e., under mild grading constraints condition
numbers stay uniformly bounded with respect to the mesh size and variable
degrees. The conceptual foundation of the envisaged preconditioners is the
auxiliary space method. The main conceptual ingredients that will be shown in
this framework to yield "optimal" preconditioners in the above sense are
Legendre-Gauss-Lobatto grids in connection with certain associated anisotropic
nested dyadic grids as well as specially adapted wavelet preconditioners for
the resulting low order auxiliary problems. Moreover, the preconditioners have
a modular form that facilitates somewhat simplified partial realizations. One
of the components can, for instance, be conveniently combined with domain
decomposition, at the expense though of a logarithmic growth of condition
numbers. Our analysis is complemented by quantitative experimental studies of
the main components.Comment: 41 pages, 11 figures; Major revision: rearrangement of the contents
for better readability, part on wavelet preconditioner adde
- …