In this paper, we develop a postprocessing derivative recovery scheme for the finite element solution uh on general unstructured but shape regular triangulations. In the case of continuous piecewise polynomials of degree p ≥ 1, by applying the global L2 projection (Qh) and a smoothing operator (Sh), the recovered pth derivatives (S mhQh∂puh) superconverge to the exact derivatives (∂pu). Based on this technique we are able to derive a local error indicator depending only on the geometry of corresponding element and the (p + 1)st derivatives approximated by ∂SmhQh∂ puh. We provide several numerical examples illustrating the effectiveness of our schemes. We also observe that higher order elements are likely to require more conservative refinement strategies to create meshes corresponding to optimal orders of convergence. © 2007 Society for Industrial and Applied Mathematics.
|Original language||English (US)|
|Number of pages||15|
|Journal||SIAM Journal on Numerical Analysis|
|State||Published - Dec 1 2007|
Bibliographical noteGenerated from Scopus record by KAUST IRTS on 2023-02-15
ASJC Scopus subject areas
- Numerical Analysis