Fast Computation of Electrostatic Potentials for Piecewise Constant Conductivities

05/30/2022
by   Kyle Bower, et al.
0

We present a novel numerical method for solving the elliptic partial differential equation problem for the electrostatic potential with piecewise constant conductivity. We employ an integral equation approach for which we derive a system of well-conditioned integral equations by representing the solution as a sum of single layer potentials. The kernel of the resulting integral operator is smooth provided that the layers are well-separated. The fast multiple method is used to accelerate the generalized minimal residual method solution of the integral equations. For efficiency, we adapt the grid of the Nyström method based on the spectral resolution of the layer charge density. Additionally, we present a method for evaluating the solution that is efficient and accurate throughout the domain, circumventing the close-evaluation problem. To support the design choices of the numerical method, we derive regularity estimates with bounds explicitly in terms of the conductivities and the geometries of the boundaries between their regions. The resulting method is fast and accurate for solving for the electrostatic potential in media with piecewise constant conductivities.

READ FULL TEXT

page 12

page 17

research
05/26/2023

A boundary integral equation method for the complete electrode model in electrical impedance tomography with tests on real-world data

We develop a boundary integral equation-based numerical method to solve ...
research
08/07/2019

An explicit numerical algorithm to the solution of Volterra integral equation of the second kind

This paper considers a numeric algorithm to solve the equation y(t)...
research
07/15/2020

Layer potential identities and subtraction techniques

When using boundary integral equation methods, we represent solutions of...
research
05/29/2023

Automatic Synthesis of Low-Complexity Translation Operators for the Fast Multipole Method

We demonstrate a new, hybrid symbolic-numerical method for the automatic...
research
09/14/2021

A fast, high-order numerical method for the simulation of single-excitation states in quantum optics

We consider the numerical solution of a nonlocal partial differential eq...
research
04/24/2023

Singularity swap quadrature for nearly singular line integrals on closed curves in two dimensions

This paper presents a quadrature method for evaluating layer potentials ...

Please sign up or login with your details

Forgot password? Click here to reset