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

Summary, in English

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.

• Author: Johan Helsing; Rikard Ojala
• Department/s: Mathematics (Faculty of Engineering); Harmonic Analysis and Applications
• Publishing year: 2008
• Language: English
• Pages: 8820-8840
• Publication/Series: Journal of Computational Physics
• Volume: 227
• Issue: 20
• Full text: - 524 kB; Download statistics
• Links: http://www.maths.lth.se/na/staff/helsing/JCP08b.pdf; http://dx.doi.org/10.1016/j.jcp.2008.06.022
• Document type: Journal article
• Publisher: Elsevier
• Topic: Mathematics
• Keywords: Corner singularity Multiphase material Triple-junction Integral equation Mesh grading Conductivity
• Status: Published
• Research group: Harmonic Analysis and Applications; Harmonic Analysis and Applications
• ISBN/ISSN/Other: ISSN: 0021-9991