How a nonconvergent recovered Hessian works in mesh adaptation
View/ Open
Issue Date
2014-01-01Author
Huang, Weizhang
Kamenski, Lennard
Publisher
Society for Industrial and Applied Mathematics
Type
Article
Article Version
Scholarly/refereed, publisher version
Metadata
Show full item recordAbstract
Hessian recovery has been commonly used in mesh adaptation for obtaining the required magnitude and direction information of the solution error. Unfortunately, a recovered Hessian from a linear finite element approximation is nonconvergent in general as the mesh is refined. It has been observed numerically that adaptive meshes based on such a nonconvergent recovered Hessian can nevertheless lead to an optimal error in the finite element approximation. This also explains why Hessian recovery is still widely used despite its nonconvergence. In this paper we develop an error bound for the linear finite element solution of a general boundary value problem under a mild assumption on the closeness of the recovered Hessian to the exact one. Numerical results show that this closeness assumption is satisfied by the recovered Hessian obtained with commonly used Hessian recovery methods. Moreover, it is shown that the finite element error changes gradually with the closeness of the recovered Hessian. This provides an explanation of how a nonconvergent recovered Hessian works in mesh adaptation.
Description
This is the published version, also available here: http://dx.doi.org/10.1137/120898796.
Collections
Citation
Kamenski, Lennard & Huang, Weizhang "How a nonconvergent recovered Hessian works in mesh adaptation." SIAM J. Numer. Anal., 52(4), 1692–1708. (17 pages). http://dx.doi.org/10.1137/120898796.
Items in KU ScholarWorks are protected by copyright, with all rights reserved, unless otherwise indicated.
We want to hear from you! Please share your stories about how Open Access to this item benefits YOU.