Jens M. Melenk

Mario Amrein, Jens M. Melenk, Thomas P. Wihler

In this paper we develop an $hp$-adaptive procedure for the numerical solution of general, semilinear elliptic boundary value problems in 1d, with possible singular perturbations. Our approach combines both a prediction-type adaptive Newton method and an $hp$-version adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully $hp$-adaptive Newton-Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for various examples.

Jens M. Melenk, Stefan A. Sauter, Céline Torres

We present a stability and convergence theory for the lossy Helmholtz equation and its Galerkin discretization. The boundary conditions are of Robin type. All estimates are explicit with respect to the real and imaginary part of the complex wave number $ζ\in\mathbb{C}$, $\operatorname{Re}ζ\geq0$, $\left\vert ζ\right\vert \geq1$. For the extreme cases $ζ\in\operatorname*{i}\mathbb{R}$ and $ζ\in\mathbb{R}_{\geq0}$, the estimates coincide with the existing estimates in the literature and exhibit a seamless transition between these cases in the right complex half plane.

Lehel Banjai, Jens M. Melenk, Ricardo H. Nochetto et al.

We design and analyze several Finite Element Methods (FEMs) applied to the Caffarelli-Silvestre extension that localizes the fractional powers of symmetric, coercive, linear elliptic operators in bounded domains with Dirichlet boundary conditions. We consider open, bounded, polytopal but not necessarily convex domains $Ω\subset \mathbb{R}^d$ with $d=1,2$. For the solution to the extension problem, we establish analytic regularity with respect to the extended variable $y\in (0,\infty)$. We prove that the solution belongs to countably normed, power-exponentially weighted Bochner spaces of analytic functions with respect to $y$, taking values in corner-weighted Kondat'ev type Sobolev spaces in $Ω$. In $Ω\subset \mathbb{R}^d$, we discretize with continuous, piecewise linear, Lagrangian FEM ($P_1$-FEM) with mesh refinement near corners, and prove that first order convergence rate is attained for compatible data $f\in \mathbb{H}^{1-s}(Ω)$. We also prove that tensorization of a $P_1$-FEM in $Ω$ with a suitable $hp$-FEM in the extended variable achieves log-linear complexity with respect to $\mathcal{N}_Ω$, the number of degrees of freedom in the domain $Ω$. In addition, we propose a novel, sparse tensor product FEM based on a multilevel $P_1$-FEM in $Ω$ and on a $P_1$-FEM on radical-geometric meshes in the extended variable. We prove that this approach also achieves log-linear complexity with respect to $\mathcal{N}_Ω$. Finally, under the stronger assumption that the data is analytic in $\overlineΩ$, and without compatibility at $\partial Ω$, we establish exponential rates of convergence of $hp$-FEM for spectral, fractional diffusion operators. We also report numerical experiments for model problems which confirm the theoretical results. We indicate several extensions and generalizations of the proposed methods.

Jens M. Melenk, Thomas P. Wihler

We consider the approximation of singularly perturbed linear second-order boundary value problems by $hp$-finite element methods. In particular, we include the case where the associated differential operator may not be coercive. Within this setting we derive an a posteriori error estimate for a natural residual norm. The error bound is robust with respect to the perturbation parameter and fully explicit with respect to both the local mesh size $h$ and the polynomial degree $p$.