BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Jens Markus Melenk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
The inverse of the stiffness matrix of the time-harmonic Maxwell equation with perfectly conducting boundary conditions is approximated in the blockwise low-rank format of ${\mathcal H}$-matrices. We prove that root exponential convergence in the block rank can be achieved if the block structure conforms to a standard admissibility criterion.
On the reference tetrahedron $\widehat K$, we define three projection-based interpolation operators on $H^2(\widehat K)$, ${\mathbf H}^1(\widehat K,\operatorname{\mathbf{curl}})$, and ${\mathbf H}^1(\widehat K,\operatorname{div})$. These operators are projections onto space of polynomials, they have the commuting diagram property and feature the optimal convergence rate as the polynomial degree increases in $H^{1-s}(\widehat K)$, ${\mathbf H}^{-s}(\widehat K,\operatorname{\mathbf{curl}})$, ${\mathbf H}^{-s}(\widehat K,\operatorname{div})$ for $0 \leq s \leq 1$.
The time-harmonic Maxwell equations at high wavenumber $k$ are discretized by edge elements of degree $p$ on a mesh of width $h$. For the case of a ball and exact, transparent boundary conditions, we show quasi-optimality of the Galerkin method under the $k$-explicit scale resolution condition that a) $kh/p$ is sufficient small and b) $p/\log k$ is bounded from below.
We develop a regularization operator based on smoothing on a locally defined length scale. This operator is defined on $L_1$ and has approximation properties that are given by the local regularity of the function it is applied to and the local length scale. Additionally, the regularized function satisfies inverse estimates commensurate with the approximation orders. By combining this operator with a classical hp-interpolation operator, we obtain an hp-Clément type quasi-interpolation operator, i.e., an operator that requires minimal smoothness of the function to be approximated but has the expected approximation properties in terms of the local mesh size and polynomial degree. As a second application, we consider residual error estimates in hp-boundary element methods that are explicit in the local mesh size and the local approximation order.
We propose and analyze an overlapping Schwarz preconditioner for the $p$ and $hp$ boundary element method for the hypersingular integral equation in 3D. We consider surface triangulations consisting of triangles. The condition number is bounded uniformly in the mesh size $h$ and the polynomial order $p$. The preconditioner handles adaptively refined meshes and is based on a local multilevel preconditioner for the lowest order space. Numerical experiments on different geometries illustrate its robustness.
The $hp$-version of the finite element method is applied to a singularly perturbed reaction-diffusion equation posed in one- and two-dimensional domains with analytic boundary. On suitably designed \emph{Spectral Boundary Layer meshes}, robust exponential convergence in a balanced norm is shown. This balanced norm is stronger than the energy norm in that the boundary layers are $O(1)$ uniformly in the singular perturbation parameter. Robust exponential convergence in the maximum norm is also established. The theoretical findings are illustrated with two numerical experiments.
We consider the question of approximating the inverse $\mathbf W = \mathbf V^{-1}$ of the Galerkin stiffness matrix $\mathbf V$ obtained by discretizing the simple-layer operator $V$ with piecewise constant functions. The block partitioning of $\mathbf W$ is assumed to satisfy any of the standard admissibility criteria that are employed in connection with clustering algorithms to approximate the discrete BEM operator $\mathbf V$. We show that $\mathbf W$ can be approximated by blockwise low-rank matrices such that the error decays exponentially in the block rank employed. Similar exponential approximability results are shown for the Cholesky factorization of $\mathbf V$.
For the reference triangle or tetrahedron $T$, we study the stability properties of the $L^2(T)$-projection $Π_N$ onto the space of polynomials of degree $N$. We show $\|Π_N u\|_{L^2(\partial T)}^2 \leq C \|u\|_{L^2(T)} \|u\|_{H^1(T)}$. This implies optimal convergence rates for the approximation error $\|u - Π_N u\|_{L^2(\partial T)}$ for all $u \in H^k(T)$, $k > 1/2$.
We study the question of approximability for the inverse of the FEM stiffness matrix for (scalar) second order elliptic boundary value problems by blockwise low rank matrices such as those given by the H-matrix format. We show that exponential convergence in the local block rank can be achieved. We also show that exponentially accurate LU-decompositions in the H-matrix format are possible for the stiffness matrices arising in the FEM. Unlike prior works, our analysis avoids any coupling of the block rank r and the mesh width h and also covers mixed Dirichlet-Neumann-Robin boundary conditions.
We consider a coupled system of two singularly perturbed reaction-diffusion equations, with two small parameters $0< ε\le μ\le 1$, each multiplying the highest derivative in the equations. The presence of these parameters causes the solution(s) to have \emph{boundary layers} which overlap and interact, based on the relative size of $ε$ and $% μ$. We construct full asymptotic expansions together with error bounds that cover the complete range $0 < ε\leq μ\leq 1$. For the present case of analytic input data, we derive derivative growth estimates for the terms of the asymptotic expansion that are explicit in the perturbation parameters and the expansion order.
We consider a (possibly) nonlinear interface problem in 2D and 3D, which is solved by use of various adaptive FEM-BEM coupling strategies, namely the Johnson-Nédélec coupling, the Bielak-MacCamy coupling, and Costabel's symmetric coupling. We provide a framework to prove that the continuous as well as the discrete Galerkin solutions of these coupling methods additionally solve an appropriate operator equation with a Lipschitz continuous and strongly monotone operator. Therefore, the coupling formulations are well-defined, and the Galerkin solutions are quasi-optimal in the sense of a Céa-type lemma. For the respective Galerkin discretizations with lowest-order polynomials, we provide reliable residual-based error estimators. Together with an estimator reduction property, we prove convergence of the adaptive FEM-BEM coupling methods. A key point for the proof of the estimator reduction are novel inverse-type estimates for the involved boundary integral operators which are advertized. Numerical experiments conclude the work and compare performance and effectivity of the three adaptive coupling procedures in the presence of generic singularities.
We prove inverse-type estimates for the four classical boundary integral operators associated with the Laplace operator. These estimates are used to show convergence of an h-adaptive algorithm for the coupling of a finite element method with a boundary element method which is driven by a weighted residual error estimator.
We review the stability properties of several discretizations of the Helmholtz equation at large wavenumbers. For a model problem in a polygon, a complete $k$-explicit stability (including $k$-explicit stability of the continuous problem) and convergence theory for high order finite element methods is developed. In particular, quasi-optimality is shown for a fixed number of degrees of freedom per wavelength if the mesh size $h$ and the approximation order $p$ are selected such that $kh/p$ is sufficiently small and $p = O(\log k)$, and, additionally, appropriate mesh refinement is used near the vertices. We also review the stability properties of two classes of numerical schemes that use piecewise solutions of the homogeneous Helmholtz equation, namely, Least Squares methods and Discontinuous Galerkin (DG) methods. The latter includes the Ultra Weak Variational Formulation.