arrow
Volume 14, Issue 4
A Fast Accurate Boundary Integral Method for Potentials on Closely Packed Cells

Wenjun Ying & J. Thomas Beale

Commun. Comput. Phys., 14 (2013), pp. 1073-1093.

Published online: 2013-10

Export citation
  • Abstract

Boundary integral methods are naturally suited for the computation of harmonic functions on a region having inclusions or cells with different material properties. However, accuracy deteriorates when the cell boundaries are close to each other. We present a boundary integral method in two dimensions which is specially designed to maintain second order accuracy even if boundaries are arbitrarily close. The method uses a regularization of the integral kernel which admits analytically determined corrections to maintain accuracy. For boundaries with many components we use the fast multipole method for efficient summation. We compute electric potentials on a domain with cells whose conductivity differs from that of the surrounding medium. We first solve an integral equation for a source term on the cell interfaces and then find values of the potential near the interfaces via integrals. Finally we use a Poisson solver to extend the potential to a regular grid covering the entire region. A number of examples are presented. We demonstrate that increased refinement is not needed to maintain accuracy as interfaces become very close.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-14-1073, author = {Wenjun Ying and J. Thomas Beale}, title = {A Fast Accurate Boundary Integral Method for Potentials on Closely Packed Cells}, journal = {Communications in Computational Physics}, year = {2013}, volume = {14}, number = {4}, pages = {1073--1093}, abstract = {

Boundary integral methods are naturally suited for the computation of harmonic functions on a region having inclusions or cells with different material properties. However, accuracy deteriorates when the cell boundaries are close to each other. We present a boundary integral method in two dimensions which is specially designed to maintain second order accuracy even if boundaries are arbitrarily close. The method uses a regularization of the integral kernel which admits analytically determined corrections to maintain accuracy. For boundaries with many components we use the fast multipole method for efficient summation. We compute electric potentials on a domain with cells whose conductivity differs from that of the surrounding medium. We first solve an integral equation for a source term on the cell interfaces and then find values of the potential near the interfaces via integrals. Finally we use a Poisson solver to extend the potential to a regular grid covering the entire region. A number of examples are presented. We demonstrate that increased refinement is not needed to maintain accuracy as interfaces become very close.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.210612.240113a}, url = {http://global-sci.org/intro/article_detail/cicp/7193.html} }
TY - JOUR T1 - A Fast Accurate Boundary Integral Method for Potentials on Closely Packed Cells AU - Wenjun Ying & J. Thomas Beale JO - Communications in Computational Physics VL - 4 SP - 1073 EP - 1093 PY - 2013 DA - 2013/10 SN - 14 DO - http://doi.org/10.4208/cicp.210612.240113a UR - https://global-sci.org/intro/article_detail/cicp/7193.html KW - AB -

Boundary integral methods are naturally suited for the computation of harmonic functions on a region having inclusions or cells with different material properties. However, accuracy deteriorates when the cell boundaries are close to each other. We present a boundary integral method in two dimensions which is specially designed to maintain second order accuracy even if boundaries are arbitrarily close. The method uses a regularization of the integral kernel which admits analytically determined corrections to maintain accuracy. For boundaries with many components we use the fast multipole method for efficient summation. We compute electric potentials on a domain with cells whose conductivity differs from that of the surrounding medium. We first solve an integral equation for a source term on the cell interfaces and then find values of the potential near the interfaces via integrals. Finally we use a Poisson solver to extend the potential to a regular grid covering the entire region. A number of examples are presented. We demonstrate that increased refinement is not needed to maintain accuracy as interfaces become very close.

Wenjun Ying and J. Thomas Beale. (2013). A Fast Accurate Boundary Integral Method for Potentials on Closely Packed Cells. Communications in Computational Physics. 14 (4). 1073-1093. doi:10.4208/cicp.210612.240113a
Copy to clipboard
The citation has been copied to your clipboard