In this paper, we propose a unified framework, the Hessian discretisation method (HDM), which is based on four discrete elements (called altogether a Hessian discretisation) and a few intrinsic indicators of accuracy, independent of the considered model. An error estimate is obtained, using only these intrinsic indicators, when the HDM framework is applied to linear fourth order problems. It is shown that HDM encompasses a large number of numerical methods for fourth order elliptic problems: finite element methods (conforming and non-conforming) as well as finite volume methods. We also use the HDM to design a novel method, based on conforming ℙ₁ finite element space and gradient recovery operators. Results of numerical experiments are presented for this novel scheme and for a finite volume scheme.
History
Journal title
Journal of Scientific Computing
Volume
78
Issue
3
Pagination
1405-1437
Publisher
Springer
Language
en, English
College/Research Centre
Faculty of Science
School
School of Mathematical and Physical Sciences
Rights statement
This is a post-peer-review, pre-copyedit version of an article published in Journal of Scientific Computing. The final authenticated version is available online at: http://dx.doi.org/10.1007/s10915-018-0814-7.