Numerical study of natural superconvergence in least-squares finite element methods for elliptic problems
Natural superconvergence of the least-squares finite element method is surveyed for the one- and two-dimensional Poisson equation. For two-dimensional problems, both the families of Lagrange elements and Raviart-Thomas elements have been considered on uniform triangular and rectangular meshes. Numerical experiments reveal that many superconvergence properties of the standard Galerkin method are preserved by the least-squares finite element method.