Local Error Estimates of the Finite Element Method for an Elliptic Problem with a Dirac Source Term
Provides rigorous error estimates for finite element methods applied to elliptic problems with Dirac sources, which are important for reduced models in physics, but the results are incremental as they extend existing L2-seminorm estimates to Hs-seminorms.
The paper proves quasi-optimal convergence in Hs-seminorm and optimal convergence in H1-seminorm for finite element solutions of elliptic problems with a Dirac source term on quasi-uniform meshes in 2D, using local error estimates and a weak Aubin-Nitsche lemma. Numerical experiments confirm the theory.
The solutions of elliptic problems with a Dirac measure in right-hand side are not H1 and therefore the convergence of the finite element solutions is suboptimal. Graded meshes are standard remedy to recover quasi-optimality, namely optimality up to a log-factor, for low order finite elements in L2-norm. Optimal (or quasi-optimal for the lowest order case) convergence has been shown in L2-seminorm, where the L2-seminorm is defined as the L2-norm on a subdomain which excludes the singularity. Here we show a quasi-optimal convergence for the Hs-seminorm, s \textgreater{} 0, and an optimal convergence in H1-seminorm for the lowest order case, on a family of quasi- uniform meshes in dimension 2. This question is motivated by the use of the Dirac measure as a reduced model in physical problems, and a high accuracy at the singularity of the finite element method is not required. Our results are obtained using local Nitsche and Schatz-type error estimates, a weak version of Aubin-Nitsche duality lemma and a discrete inf-sup condition. These theoretical results are confirmed by numerical illustrations.