Error estimates of finite element method about elliptic problems with singular righthand side (Q2371495)

The theory of finite element methods for elliptic problems with singular righthand side: \[ \begin{cases} Lu= \Sigma_{i,j}\frac{\partial}{\partial x_i} (a^{ij} \frac{\partial G_{x_0}}{\partial x_j})= \delta(x-x_0) &\text{in } \Omega\\ G_{x_0}=0 &\text{on } \partial \Omega \end{cases}\tag{1} \] \(n=2, x_0 \in \Omega, \rho(x_0, \partial \Omega) \geq C, \delta(0)=+ \infty, \delta(x)=0, x\neq0, \int_\Omega \delta(x)\,dx=1\) is studied. It is assumed that \(G_{x_0}\) is the weak solution for (1) and \(G^h_{x_0}\) is a \(k\)th-order Lagrangian finite element solution for (1) on basis of quasi-uniform partition, \(\mathcal{T}_h\) of \(\Omega\) with mesh size \(h\). Main result: The author shows that there exists a positive constant \(C\) (which is independent of \(h\)) such that \[ \| G_{x_0}-G^h_{x_0}\| _{L^2(\Omega - \Omega_r)}\leq C h^{k+1}r^{-k}, \] \[ \| \nabla G_{x_0}-G^h_{x_0}\| _{L^2(\Omega - \Omega_r)}\leq C h^{k}r^{-k}, \] where \(r>C_1h| \ln h| , \Omega_r=\{x'\mid | x'-x_0| <r\}, \partial \Omega_r=\{x'\mid | x'-x_0| =r\}.\) Finally the results of some numerical examples to investigate these error estimates are also presented.