The Influence of Higher Order FEM Discretisations on Multigrid Convergence

Abstract Quadratic and even higher order finite elements are interesting candidates for the numerical solution of partial differential equations (PDEs) due to their improved approximation properties in comparison to linear approaches. The systems of equations that arise from the discretisation of the underlying (elliptic) PDEs are often solved by iterative solvers like preconditioned Krylow-space methods, while multigrid solvers are still rarely used – which might be caused by the high effort that is associated with the realisation of the necessary data structures as well as smoothing and intergrid transfer operators. In this note, we discuss the numerical analysis of quadratic conforming finite elements in a multigrid solver. Using the “correct” grid transfer operators in conjunction with a quadratic finite element approximation allows to formulate an improved approximation property which enhances the (asymptotic) behaviour of multigrid: If m denotes the number of smoothing steps, the convergence rates behave asymptotically like O(1/m2) in contrast to O(1/m) for linear FEM.