Extending linear finite elements to quadratic precision on arbitrary meshes

详细信息    查看全文

• 作者：Daniel Duque ; ^a ; ^{daniel.duque@upm.es} ; Pep Españ ; ol^b ; Jaime Arturo de la Torre^b
• 关键词：Finite element method ; Quadratic precision ; Linear precision ; Laplacian
• 刊名：Applied Mathematics and Computation
• 出版年：2017
• 出版时间：15 May 2017
• 年：2017
• 卷：301
• 期：Complete
• 页码：201-213
• 全文大小：899 K
• 卷排序：301

文摘

The method of finite elements can provide discrete versions of differential operators in general arrangements of nodes. Linear finite elements defined on a triangulation, like Delaunay’s, constitute a typical route towards this aim, as well as being a tool for interpolation. We discuss a procedure to build interpolating functions with quadratic precision from this functional set. The idea is to incorporate (extend) the products of the original functions in order to go from linear precision (the property that linear functions are reconstructed exactly) to quadratic precision (quadratic functions are reconstructed exactly). This procedure is applied to the standard Galerkin approach, in order to study the Poisson problem and the direct evaluation of the Laplacian. We discuss convergence of these problems in 1D and 2D. The extended method is shown to be superior in 2D, featuring convergence for the direct evaluation of the Laplacian in distorted lattices.