Svetlana Tlupova

Svetlana Tlupova, J. Thomas Beale

We present a numerical method for computing the single layer (Stokeslet) and double layer (stresslet) integrals in Stokes flow. The method applies to smooth, closed surfaces in three dimensions, and achieves high accuracy both on and near the surface. The singular Stokeslet and stresslet kernels are regularized and, for the nearly singular case, corrections are added to reduce the regularization error. These corrections are derived analytically for both the Stokeslet and the stresslet using local asymptotic analysis. For the case of evaluating the integrals on the surface, as needed when solving integral equations, we design high order regularizations for both kernels that do not require corrections. This approach is direct in that it does not require grid refinement or special quadrature near the singularity, and therefore does not increase the computational complexity of the overall algorithm. Numerical tests demonstrate the uniform convergence rates for several surfaces in both the singular and near singular cases, as well as the importance of corrections when two surfaces are close to each other.

Luiz M. Faria, Carlos Perez-Arancibia, Svetlana Tlupova

This paper extends and analyzes the high-order kernel regularization framework of Beale & Tlupova (arXiv:2510.13639) to all four boundary integral operators of the Helmholtz Calderon calculus in three dimensions: the single-layer, double-layer, adjoint double-layer, and hypersingular operators. To the best of our knowledge, this work provides the first high-order kernel regularization of the hypersingular operator for both the Helmholtz and Laplace equations in three dimensions. The regularization replaces each singular kernel with a smooth modification constructed from error functions together with a polynomial correction whose coefficients are determined through moment conditions. Alongside the derivation of the regularizing functions, the paper provides a unified error analysis of the combined regularization and quadrature discretization procedure. By coupling the regularization parameter to the mesh size, the two error contributions can be balanced, leading to explicit overall convergence rates that depend jointly on the order of the regularization and the degree of exactness of the surface quadrature rule. A key practical feature of the method is its implementation simplicity. Once the regularizing functions are determined, the numerical task reduces entirely to the evaluation of smooth surface integrals using standard quadrature, without the need for element-local solves, singularity-specific precomputations, or specialized quadrature rules. Although the modified kernel is generally incompatible with kernel-specific fast methods, this limitation is addressed through H-matrix acceleration, applicable in a black-box manner. Numerical examples -- including verification of the predicted convergence rates and solution of sound-soft and sound-hard scattering problems by smooth obstacles -- demonstrate the accuracy and practicality of the proposed methodology.