Abstract

We develop a method for computing a nearly singular integral, such as a double layer potential due to sources on a curve in the plane, evaluated at a point near the curve. The approach is to regularize the singularity and obtain a preliminary value from a standard quadrature rule. Then we add corrections for the errors due to smoothing and discretization, which are found by asymptotic analysis. We prove an error estimate for the corrected value, uniform with respect to thepoint of evaluation. One application is a simple method for solving the Dirichlet problem for Laplace's equation on a grid covering an irregular region in the plane, similar to an earlier method of A. Mayo [SIAM J. Sci. Statist. Comput., 6 (1985), pp. 144--157]. This approach could also be used to compute the pressure gradient due to a force on a moving boundary in an incompressible fluid. Computational examples are given for the double layer potential and for the Dirichlet problem.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call