A Nystr枚m method for weakly singular integral operators on surfaces

详细信息	查看全文 | 推荐本文 |

• 作者：James Bremer^a ; ^{bremer@math.ucdavis.edu} ; Zydrunas Gimbutas^b
• 关键词：Integral operators ; Elliptic boundary value problems ; Singular domains ; Scattering theory
• 刊名：Journal of Computational Physics
• 出版年：2012
• 期刊代码：136_00219991
• 类别：cp
• 出版时间：20 May, 2012
• 卷：231
• 期：14
• 页码：4885-4903
• 文件大小：628 K

摘要

We describe a modified Nystr枚m method for the discretization of the weakly singular boundary integral operators which arise from the formulation of linear elliptic boundary value problems as integral equations. Standard Nystr枚m and collocation schemes proceed by representing functions via their values at a collection of quadrature nodes. Our method uses appropriately scaled function values in lieu of such representations. This results in a scheme which is mathematically equivalent to Galerkin discretization in that the resulting matrices are related to those obtained by Galerkin methods via conjugation with well-conditioned matrices, but which avoids the evaluation of double integrals. Moreover, we incorporate a new mechanism for approximating the singular integrals which arise from the discretization of weakly singular integral operators which is considerably more efficient than standard methods. We illustrate the performance of our method with numerical experiments.