885 research outputs found
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
Multilevel quadrature for elliptic problems on random domains by the coupling of FEM and BEM
Elliptic boundary value problems which are posed on a random domain can be
mapped to a fixed, nominal domain. The randomness is thus transferred to the
diffusion matrix and the loading. While this domain mapping method is quite
efficient for theory and practice, since only a single domain discretisation is
needed, it also requires the knowledge of the domain mapping.
However, in certain applications, the random domain is only described by its
random boundary, while the quantity of interest is defined on a fixed,
deterministic subdomain. In this setting, it thus becomes necessary to compute
a random domain mapping on the whole domain, such that the domain mapping is
the identity on the fixed subdomain and maps the boundary of the chosen fixed,
nominal domain on to the random boundary.
To overcome the necessity of computing such a mapping, we therefore couple
the finite element method on the fixed subdomain with the boundary element
method on the random boundary. We verify the required regularity of the
solution with respect to the random domain mapping for the use of multilevel
quadrature, derive the coupling formulation, and show by numerical results that
the approach is feasible
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
On local Fourier analysis of multigrid methods for PDEs with jumping and random coefficients
In this paper, we propose a novel non-standard Local Fourier Analysis (LFA)
variant for accurately predicting the multigrid convergence of problems with
random and jumping coefficients. This LFA method is based on a specific basis
of the Fourier space rather than the commonly used Fourier modes. To show the
utility of this analysis, we consider, as an example, a simple cell-centered
multigrid method for solving a steady-state single phase flow problem in a
random porous medium. We successfully demonstrate the prediction capability of
the proposed LFA using a number of challenging benchmark problems. The
information provided by this analysis helps us to estimate a-priori the time
needed for solving certain uncertainty quantification problems by means of a
multigrid multilevel Monte Carlo method
Schnelle Löser für Partielle Differentialgleichungen
The workshop Schnelle Löser für partielle Differentialgleichungen, organised by Randolph E. Bank (La Jolla), Wolfgang Hackbusch (Leipzig), and Gabriel Wittum (Frankfurt am Main), was held May 22nd–May 28th, 2011. This meeting was well attended by 54 participants with broad geographic representation from 7 countries and 3 continents. This workshop was a nice blend of researchers with various backgrounds
- …