Description

In this project we plan to develop and analyze new subspace correction (SC) methods for the numerical solution of coupled systems of partial differential equations (PDE). The focus is on nearly singular symmetric positive definite (SPD) and on indefinite problems.
We propose an integrated approach in which it is essential to use discretization techniques that preserve certain conservation laws, and to combine them with an adaptive solution process. In this way, one can design methods that perform optimally with respect to:

(i) accurate approximation of the unknown quantities;
(ii) obtaining the numerical solution in optimal time; and
(iii) scalability with respect to both, problem size and advances in computer hardware.

The present project has the following three interrelated Components (C1)–(C3) with a main emphasis on systems with highly oscillatory coefficients:

(C1):
SC methods for nearly incompressible elasticity and Stokes flow.
(C2):
SC methods for total variation minimization of discrete functionals arising in sparse data recovery.
(C3):
Auxiliary space and SC methods for elliptic problems with highly oscillatory coefficients.

The primary goal of the proposed research work is to contribute to extending the theory and applicability of subspace correction methods to the above-mentioned classes of problems.

Starting point of the research plan is the use and interplay of stable and accurate finite element schemes and of the efficient preconditioning of the related discrete problems.
In the present setting nonconforming and in particular discontinuous Galerkin (DG) finite element methods provide adequate discretization tools. Some of their most attractive properties and practical advantages over conforming methods are that

a) it is easy to extended DG methods to higher approximation order;
b) they are well suited to treat complex geometries in combination with unstructured and hybrid meshes;
c) they can be combined with any element type where the grids are also allowed to have hanging nodes;
d) they can easily handle adaptive strategies;
e) they have favorable properties in view of parallel computing.

A main disadvantage of DG discretizations is that they produce an excess of degrees of freedom (as compared to conforming methods of the same approximation order) which in general makes the solution of the arising linear systems more difficult and more time consuming. We therefore put strong efforts on devising new efficient and robust solution methods, covering wider classes of problems (see (C1)–(C3)) that arise from nonconforming and discontinuous Galerkin discretizations. The final aim is to adapt our methods to and to test them on industrial and multiphysics applications, e.g., in reservoir engineering, or in life science. Some of the problems in which we are particularly interested stem from micro-mechanics modeling of heterogeneous media, e.g., the modeling of fluid flow in porous media, the determination of the bio-mechanical properties of bones, or the reconstruction of (medical) images. Typically such problems involve parameters that lead to highly ill-conditioned systems of linear algebraic equations.