An accurate, robust, and easy-to-implement method for integration over arbitrary polyhedra: Application to embedded interface methods

Y. Sudhakar, J. P. Moitinho de Almeida, Wolfgang A. Wall

Research output: Contribution to journalArticlepeer-review

48 Scopus citations

Abstract

We present an accurate method for the numerical integration of polynomials over arbitrary polyhedra. Using the divergence theorem, the method transforms the domain integral into integrals evaluated over the facets of the polyhedra. The necessity of performing symbolic computation during such transformation is eliminated by using one dimensional Gauss quadrature rule. The facet integrals are computed with the help of quadratures available for triangles and quadrilaterals. Numerical examples, in which the proposed method is used to integrate the weak form of the Navier-Stokes equations in an embedded interface method (EIM), are presented. The results show that our method is as accurate and generalized as the most widely used volume decomposition based methods. Moreover, since the method involves neither volume decomposition nor symbolic computations, it is much easier for computer implementation. Also, the present method is more efficient than other available integration methods based on the divergence theorem. Efficiency of the method is also compared with the volume decomposition based methods and moment fitting methods. To our knowledge, this is the first article that compares both accuracy and computational efficiency of methods relying on volume decomposition and those based on the divergence theorem.

Original languageEnglish
Pages (from-to)393-415
Number of pages23
JournalJournal of Computational Physics
Volume273
DOIs
StatePublished - 15 Sep 2014

Keywords

  • Complex volumes
  • Divergence theorem
  • Embedded interface method
  • Enriched partition of unity method
  • Extended finite element method
  • Numerical integration
  • Polyhedra

Fingerprint

Dive into the research topics of 'An accurate, robust, and easy-to-implement method for integration over arbitrary polyhedra: Application to embedded interface methods'. Together they form a unique fingerprint.

Cite this