Dietmar Gallistl

Donald L. Brown, Dietmar Gallistl, Daniel Peterseim

This paper presents a multiscale Petrov-Galerkin finite element method for time-harmonic acoustic scattering problems with heterogeneous coefficients in the high-frequency regime. We show that the method is pollution- free also in the case of heterogeneous media provided that the stability bound of the continuous problem grows at most polynomially with the wave number k. By generalizing classical estimates of [Melenk, Ph.D. Thesis 1995] and [Hetmaniuk, Commun. Math. Sci. 5 (2007)] for homogeneous medium, we show that this assumption of polynomially wave number growth holds true for a particular class of smooth heterogeneous material coefficients. Further, we present numerical examples to verify our stability estimates and implement an example in the wider class of discontinuous coefficients to show computational applicability beyond our limited class of coefficients.

Dietmar Gallistl, Daniel Peterseim

This paper aims at bridging existing theories in numerical and analytical homogenization. For this purpose the multiscale method of Målqvist and Peterseim [Math. Comp. 2014], which is based on orthogonal subspace decomposition, is reinterpreted by means of a discrete integral operator acting on standard finite element spaces. The exponential decay of the involved integral kernel motivates the use of a diagonal approximation and, hence, a localized piecewise constant coefficient. In a periodic setting, the computed localized coefficient is proved to coincide with the classical homogenization limit. An a priori error analysis shows that the local numerical model is appropriate beyond the periodic setting when the localized coefficient satisfies a certain homogenization criterion, which can be verified a posteriori. The results are illustrated in numerical experiments.

Dietmar Gallistl

This paper studies formulations of second-order elliptic partial differential equations in nondivergence form on convex domains as equivalent variational problems. The first formulation is that of Smears \& Süli [SIAM J.\ Numer.\ Anal.\ 51(2013), pp.\ 2088--2106.], and the second one is a new symmetric formulation based on a least-squares functional. These formulations enable the use of standard finite element techniques for variational problems in subspaces of $H^2$ as well as mixed finite element methods from the context of fluid computations. Besides the immediate quasi-optimal a~priori error bounds, the variational setting allows for a~posteriori error control with explicit constants and adaptive mesh-refinement. The convergence of an adaptive algorithm is proved. Numerical results on uniform and adaptive meshes are included.

Dietmar Gallistl, Daniel Peterseim

This paper proposes a numerical upscaling procedure for elliptic boundary value problems with diffusion tensors that vary randomly on small scales. The resulting effective deterministic model is given through a quasilocal discrete integral operator, which can be further compressed to an effective partial differential operator. Error estimates consisting of a priori and a posteriori terms are provided that allow one to quantify the impact of uncertainty in the diffusion coefficient on the expected effective response of the process.

Dietmar Gallistl, Patrick Henning, Barbara Verfürth

If an elliptic differential operator associated with an $\mathbf{H}(\mathrm{curl})$-problem involves rough (rapidly varying) coefficients, then solutions to the corresponding $\mathbf{H}(\mathrm{curl})$-problem admit typically very low regularity, which leads to arbitrarily bad convergence rates for conventional numerical schemes. The goal of this paper is to show that the missing regularity can be compensated through a corrector operator. More precisely, we consider the lowest order Nédélec finite element space and show the existence of a linear corrector operator with four central properties: it is computable, $\mathbf{H}(\mathrm{curl})$-stable, quasi-local and allows for a correction of coarse finite element functions so that first-order estimates (in terms of the coarse mesh-size) in the $\mathbf{H}(\mathrm{curl})$ norm are obtained provided the right-hand side belongs to $\mathbf{H}(\mathrm{div})$. With these four properties, a practical application is to construct generalized finite element spaces which can be straightforwardly used in a Galerkin method. In particular, this characterizes a homogenized solution and a first order corrector, including corresponding quantitative error estimates without the requirement of scale separation.

Dietmar Gallistl, Daniel Peterseim

We present and analyze a pollution-free Petrov-Galerkin multiscale finite element method for the Helmholtz problem with large wave number $κ$ as a variant of [Peterseim, ArXiv:1411.1944, 2014]. We use standard continuous $Q_1$ finite elements at a coarse discretization scale $H$ as trial functions, whereas the test functions are computed as the solutions of local problems at a finer scale $h$. The diameter of the support of the test functions behaves like $mH$ for some oversampling parameter $m$. Provided $m$ is of the order of $\log(κ)$ and $h$ is sufficiently small, the resulting method is stable and quasi-optimal in the regime where $H$ is proportional to $κ^{-1}$. In homogeneous (or more general periodic) media, the fine scale test functions depend only on local mesh-configurations. Therefore, the seemingly high cost for the computation of the test functions can be drastically reduced on structured meshes. We present numerical experiments in two and three space dimensions.

Daniele Boffi, Dietmar Gallistl, Francesca Gardini et al.

It is shown that the h-adaptive mixed finite element method for the discretization of eigenvalue clusters of the Laplace operator produces optimal convergence rates in terms of nonlinear approximation classes. The results are valid for the typical mixed spaces of Raviart-Thomas or Brezzi-Douglas-Marini type with arbitrary fixed polynomial degree in two and three space dimensions.

Dietmar Gallistl

This paper studies the nonconforming Morley finite element approximation of the eigenvalues of the biharmonic operator. A new $C^1$ conforming companion operator leads to an $L^2$ error estimate for the Morley finite element method which directly compares the $L^2$ error with the error in the energy norm and, hence, can dispense with any additional regularity assumptions. Furthermore, the paper presents new eigenvalue error estimates for nonconforming finite elements that bound the error of (possibly multiple or clustered) eigenvalues by the approximation error of the computed invariant subspace. An application is the proof of optimal convergence rates for the adaptive Morley finite element method for eigenvalue clusters.