Sharp quadrature error bounds for the nearest-neighbor discretization of the regularized stokeslet boundary integral equation

Meurig Thomas Gallagher, Debajyoti Choudhuri, David Smith

Research output: Contribution to journalArticlepeer-review

4 Citations (Scopus)
130 Downloads (Pure)


The method of regularized stokeslets is a powerful numerical approach to solve the Stokes flow equations for problems in biological fluid mechanics. A recent variation of this method incorporates a nearest-neighbor discretization to improve accuracy and efficiency while maintaining the ease of implementation of the original meshless method. This new method contains three sources of numerical error: the regularization error associated with using the regularized form of the boundary integral equations (with parameter ϵ), and two sources of discretization error associated with the force and quadrature discretizations (with lengthscales h f and h q ). A key issue to address is the quadrature error; initial work has not fully explained observed numerical convergence phenomena. In the present manuscript we construct sharp quadrature error bounds for the nearest-neighbor discretization, noting that the error for a single evaluation of the kernel depends on the smallest distance (δ ) between these discretization sets. The quadrature error bounds are described for two cases: disjoint sets ((δ > 0) that are close to linear in h q and insensitive to ϵ, and contained sets ((δ = 0) that are quadratic in h q with inverse dependence on ϵ. The practical implications of these error bounds are discussed in reference to the condition number of the matrix system for the nearestneighbor method, with the analysis revealing that the condition number is insensitive to ϵ for disjoint sets, and grows linearly with ϵ for contained sets. Error bounds for the general case (δ ≥ 0) are revealed to be proportional to the sum of the errors for each case.

Original languageEnglish
Pages (from-to)B139–B152
Number of pages14
JournalSIAM Journal on Scientific Computing
Issue number1
Publication statusPublished - 19 Feb 2019


  • quadrature error bounds
  • regularized stokeslets
  • stoeks flow
  • boundary integral


Dive into the research topics of 'Sharp quadrature error bounds for the nearest-neighbor discretization of the regularized stokeslet boundary integral equation'. Together they form a unique fingerprint.

Cite this