A dual singular complement method for the numerical solution of the Poisson equation with boundary data in non-convex domains
arXiv:1505.00414
Abstract
The very weak solution of the Poisson equation with boundary data is defined by the method of transposition. The finite element solution with regularized boundary data converges with order 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.