Corner singularities for elliptic problems: Integral equations, graded meshes, quadrature, and compressed inverse preconditioning

Research output: Contribution to journalArticle

Abstract

We take a fairly comprehensive approach to the problem of solving elliptic partial differential equations numerically using integral equation methods on domains where the boundary has a large number of corners and branching points. Use of non-standard integral equations, graded meshes, interpolatory quadrature, and compressed inverse preconditioning are techniques that are explored, developed, mixed, and tested on some familiar problems in materials science. The recursive compressed inverse preconditioning, the major novelty of the paper, turns out to be particularly powerful and, when it applies, eliminates the need for mesh grading completely. In an electrostatic example for a multiphase granular material with about two thousand corners and triple junctions and a conductivity ratio between phases up to a million we compute a common functional of the solution with an estimated relative error of 10-12. In another example, five times as large but with a conductivity ratio of only a hundred, we achieve an estimated relative error of 10-14.

Details

Authors
Organisations
Research areas and keywords

Subject classification (UKÄ) – MANDATORY

  • Mathematics

Keywords

  • Corner singularity Multiphase material Triple-junction Integral equation Mesh grading Conductivity
Original languageEnglish
Pages (from-to)8820-8840
JournalJournal of Computational Physics
Volume227
Issue number20
Publication statusPublished - 2008
Publication categoryResearch
Peer-reviewedYes

Bibliographic note

The information about affiliations in this record was updated in December 2015. The record was previously connected to the following departments: Numerical Analysis (011015004)

Total downloads

No data available