Uniform shift estimates for transmission problems and optimal rates of convergence for the parametric Finite Element Method

TL;DR

This paper establishes uniform regularity and convergence rates for parametric elliptic PDEs with discontinuous coefficients, enabling efficient polynomial chaos approximations in uncertainty quantification.

Contribution

It provides new regularity results in weighted Sobolev spaces for elliptic transmission problems with discontinuities, leading to optimal convergence rates for polynomial chaos methods.

Findings

Uniform shift theorem for parametric elliptic PDEs.

Optimal algebraic convergence rates for Galerkin approximations.

Regularity results in broken weighted Sobolev spaces.

Abstract

Let $\Omega \subset \RR^d$, $d \geqslant 1$, be a bounded domain with piecewise smooth boundary $\partial \Omega $ and let $U$ be an open subset of a Banach space $Y$. Motivated by questions in "Uncertainty Quantification," we consider a parametric family $P = (P_y)_{y \in U}$ of uniformly strongly elliptic, second order partial differential operators $P_y$ on $\Omega$. We allow jump discontinuities in the coefficients. We establish a regularity result for the solution $u: \Omega \times U \to \RR$ of the parametric, elliptic boundary value/transmission problem $P_y u_y = f_y$, $y \in U$, with mixed Dirichlet-Neumann boundary conditions in the case when the boundary and the interface are smooth and in the general case for $d=2$. Our regularity and well-posedness results are formulated in a scale of broken weighted Sobolev spaces $\hat\maK^{m+1}_{a+1}(\Omega)$ of Babu\v{s}ka-Kondrat'ev type in $\Omega$, possibly augmented by some locally constant functions. This implies that the parametric, elliptic PDEs $(P_y)_{y \in U}$ admit a shift theorem that is uniform in the parameter $y\in U$. In turn, this then leads to $h^m$-quasi-optimal rates of convergence (i.e. algebraic orders of convergence) for the Galerkin approximations of the solution $u$, where the approximation spaces are defined using the "polynomial chaos expansion" of $u$ with respect to a suitable family of tensorized Lagrange polynomials, following the method developed by Cohen, Devore, and Schwab (2010).