Higher order energy-corrected finite element methods
For computational scientists solving elliptic PDEs on non-smooth domains, this method improves accuracy without increasing polynomial degree, addressing a known bottleneck in finite element methods.
The paper generalizes the energy-correction method for higher order finite elements to handle reduced regularity in polygonal domains with re-entrant corners, eliminating the pollution effect and achieving optimal order estimates in weighted L2-norms, with a post-processing step recovering optimal order convergence in the standard L2-norm.
The regularity of the solution of elliptic partial differential equa- tions in a polygonal domain with re-entrant corners is, in general, reduced compared to the one on a smooth convex domain. This results in a best approximation property for standard norms which depend on the re-entrant corner but does not increase with the polynomial degree. Standard Galerkin approximations are moreover affected by a global pollution effect. Even in the far field no optimal error reduction can be observed. Here, we generalize the energy-correction method for higher order finite elements. It is based on a parameter-dependent local modification of the stiffness matrix. We will show firstly that for such modified finite element approximation the pollution effect does not occur and thus optimal order estimates in weighted L2-norms can be obtained. Two different modification techniques are introduced and illustrated numerically. Secondly we propose a simple post-processing step such that even with respect to the standard L2-norm optimal order convergence can be recovered.