Gradient correction for polygonal and polyhedral finite elements

Cameron Talischi, Anderson Pereira, Ivan F.M. Menezes, Glaucio H. Paulino

Research output: Contribution to journalArticlepeer-review

36 Scopus citations

Abstract

Previous studies have shown that the commonly used quadrature schemes for polygonal and polyhedral finite elements lead to consistency errors that persist under mesh refinement and subsequently render the approximations non-convergent. In this work, we consider minimal perturbations to the gradient field at the element level in order to restore polynomial consistency and recover optimal convergence rates when the weak form integrals are evaluated using quadrature. For finite elements of arbitrary order, we state the accuracy requirements on the underlying volumetric and boundary quadrature rules and discuss the properties of the resulting corrected gradient operator. We compare the proposed approach with the pseudo-derivative method developed by Belytschko and co-workers and, for linear elliptic problems, with our previous remedy that involves splitting of polynomial and non-polynomial of elemental energy bilinear form. We present several numerical results for linear and nonlinear elliptic problems in two and three dimensions that not only confirm the recovery of optimal convergence rates but also suggest that the global error levels are close to those of approximations obtained from exact evaluation of the weak form integrals.

Original languageEnglish (US)
Pages (from-to)728-747
Number of pages20
JournalInternational Journal for Numerical Methods in Engineering
Volume102
Issue number3-4
DOIs
StatePublished - Apr 20 2015
Externally publishedYes

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Engineering(all)
  • Applied Mathematics

Keywords

  • Generalized barycentric coordinates
  • Gradient correction
  • Polygonal and polyhedral finite elements
  • Pseudo-derivatives
  • Quadrature error

Fingerprint

Dive into the research topics of 'Gradient correction for polygonal and polyhedral finite elements'. Together they form a unique fingerprint.

Cite this