Source author record

Thomas Apel

Thomas Apel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2020arXiv

A nonconforming pressure-robust finite element method for the Stokes equations on anisotropic meshes

Most classical finite element schemes for the (Navier-)Stokes equations are neither pressure-robust, nor are they inf-sup stable on general anisotropic triangulations. A lack of pressure-robustness may lead to large velocity errors, whenever the Stokes momentum balance is dominated by a strong and complicated pressure gradient. It is a consequence of a method, which does not exactly satisfy the divergence constraint. However, inf-sup stable schemes can often be made pressure-robust just by a recent, modified discretization of the exterior forcing term, using $\mathbf{H}(\operatorname{div})$-conforming velocity reconstruction operators. This approach has so far only been analyzed on shape-regular triangulations. The novelty of the present contribution is that the reconstruction approach for the Crouzeix-Raviart method, which has a stable Fortin operator on arbitrary meshes, is combined with results on the interpolation error on anisotropic elements for reconstruction operators of Raviart-Thomas and Brezzi-Douglas-Marini type, generalizing the method to a large class of anisotropic triangulations. Numerical examples confirm the theoretical results in a 2D and a 3D test case.

preprint2019arXiv

Brezzi-Douglas-Marini interpolation of any order on anisotropic triangles and tetrahedra

Recently, the $\vec{H}(\operatorname{div})$-conforming finite element families for second order elliptic problems have come more into focus, since due to hybridization and subsequent advances in computational efficiency their use is no longer mainly theoretical. Their property of yielding exactly divergence-free solutions for mixed problems makes them interesting for a variety of applications, including incompressible fluids. In this area, boundary and interior layers are present, which demand the use of anisotropic elements. While for the Raviart-Thomas interpolation of any order on anisotropic tetrahedra optimal error estimates are known, this contribution extends these results to the Brezzi-Douglas-Marini finite elements. Optimal interpolation error estimates are proved under two different regularity conditions on the elements, which both relax the standard minimal angle condition. Additionally a numerical application on the Stokes equations is presented to illustrate the findings.

preprint2016arXiv

Adapted numerical methods for the numerical solution of the Poisson equation with $L^2$ boundary data in non-convex domains

The very weak solution of the Poisson equation with $L^2$ boundary data is defined by the method of transposition. The finite element solution with regularized boundary data converges in the $L^2(Ω)$-norm with order $1/2$ in convex domains but has a reduced convergence order in non-convex domains although the solution remains to be contained in $H^{1/2}(Ω)$. The reason is a singularity in the dual problem. In this paper we propose and analyze, as a remedy, both a standard finite element method with mesh grading and a dual variant of the singular complement method. The error order 1/2 is retained in both cases also with non-convex domains. Numerical experiments confirm the theoretical results.

preprint2015arXiv

A dual singular complement method for the numerical solution of the Poisson equation with $L^2$ boundary data in non-convex domains

The very weak solution of the Poisson equation with $L^2$ boundary data is defined by the method of transposition. The finite element solution with regularized boundary data converges with order $1/2$ in convex domains but has a reduced convergence order in non-convex domains. As a remedy, a dual variant of the singular complement method is proposed. The error order of the convex case is retained. Numerical experiments confirm the theoretical results.

preprint2015arXiv

Discretization of the Poisson equation with non-smooth data and emphasis on non-convex domains

Several approaches are discussed how to understand the solution of the Dirichlet problem for the Poisson equation when the Dirichlet data are non-smooth such as if they are in $L^2$ only. For the method of transposition (sometimes called very weak formulation) three spaces for the test functions are considered, and a regularity result is proved. An approach of Berggren is recovered as the method of transposition with the second variant of test functions. A further concept is the regularization of the boundary data combined with the weak solution of the regularized problem. The effect of the regularization error is studied. The regularization approach is the simplest to discretize. The discretization error is estimated for a sequence of quasi-uniform meshes. Since this approach turns out to be equivalent to Berggren's discretization his error estimates are rendered more precisely. Numerical tests show that the error estimates are sharp, in particular that the order becomes arbitrarily small when the maximal interior angle of the domain tends to $2π$.