×

Sharp quadrature error bounds for the nearest-neighbor discretization of the regularized stokeslet boundary integral equation. (English) Zbl 1412.35252

Summary: 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 \(\varepsilon\)), 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 (\(\delta\)) between these discretization sets. The quadrature error bounds are described for two cases: disjoint sets (\(\delta>0\)) that are close to linear in \(h_q\) and insensitive to \(\varepsilon\), and contained sets (\(\delta=0\)) that are quadratic in \(h_q\) with inverse dependence on \(\varepsilon\). The practical implications of these error bounds are discussed in reference to the condition number of the matrix system for the nearest-neighbor method, with the analysis revealing that the condition number is insensitive to \(\varepsilon\) for disjoint sets, and grows linearly with \(\varepsilon\) for contained sets. Error bounds for the general case (\(\delta\geq 0\)) are revealed to be proportional to the sum of the errors for each case.

MSC:

35Q35 PDEs in connection with fluid mechanics
65G99 Error analysis and interval analysis
65M15 Error bounds for initial value and initial-boundary value problems involving PDEs
76D07 Stokes and related (Oseen, etc.) flows
76M25 Other numerical methods (fluid mechanics) (MSC2010)
76Z05 Physiological flows

Software:

BEMLIB
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] L. A. Klinteberg and A.-K. Tornberg, {\it A fast integral equation method for solid particles in viscous flow using quadrature by expansion}, J. Comput. Phys., 326 (2016), pp. 420-445. · Zbl 1373.76037
[2] J. Ainley, S. Durkin, R. Embid, P. Boindala, and R. Cortez, {\it The method of images for regularized Stokeslets}, J. Comput. Phys., 227 (2008), pp. 4600-4616. · Zbl 1388.76060
[3] A. Barrero-Gil, {\it Weakening accuracy dependence with the regularization parameter in the method of regularized Stokeslets}, J. Comput. Appl. Math., 237 (2013), pp. 672-679. · Zbl 1426.76519
[4] H. E. Bell, {\it Gershgorin’s theorem and the zeros of polynomials}, Amer. Math. Monthly, 72 (1965), pp. 292-295. · Zbl 0134.02003
[5] O. P. Bruno and L. A. Kunyansky, {\it A fast, high-order algorithm for the solution of surface scattering problems: Basic implementation, tests, and applications}, J. Comput. Phys., 169 (2001), pp. 80-110. · Zbl 1052.76052
[6] C. Carvalho, S. Khatri, and A. D. Kim, {\it Asymptotic analysis for close evaluation of layer potentials}, J. Comput. Phys., 355 (2018), pp. 327-341. · Zbl 1380.65396
[7] R. Cortez, {\it The method of regularized Stokeslets}, SIAM J. Sci. Comput., 23 (2001), pp. 1204-1225, . · Zbl 1064.76080
[8] R. Cortez, L. Fauci, and A. Medovikov, {\it The method of regularized Stokeslets in three dimensions: Analysis, validation, and application to helical swimming}, Phys. Fluids, 17 (2005), 031504. · Zbl 1187.76105
[9] R. Cortez and F. Hoffmann, {\it A fast numerical method for computing doubly-periodic regularized Stokes flow in 3D}, J. Comput. Phys., 258 (2014), pp. 1-14. · Zbl 1349.76596
[10] R. Cortez and M. Nicholas, {\it Slender body theory for Stokes flows with regularized forces}, Commun. Appl. Math. Comput. Sci., 7 (2012), pp. 33-62. · Zbl 1452.76055
[11] R. Cortez and D. Varela, {\it A general system of images for regularized Stokeslets and other elements near a plane wall}, J. Comput. Phys., 285 (2015), pp. 41-54. · Zbl 1351.76022
[12] M. T. Gallagher and D. J. Smith, {\it Meshfree and efficient modeling of swimming cells}, Phys. Rev. Fluids, 3 (2018), 053101.
[13] S. Kim and S. J. Karrila, {\it Microhydrodynamics: Principles and Selected Applications}, Dover, Mineola, NY, 2013.
[14] C. Pozrikidis, {\it Boundary Integral and Singularity Methods for Linearized Viscous Flow}, Cambridge University Press, Cambridge, 1992. · Zbl 0772.76005
[15] C. Pozrikidis, {\it A Practical Guide to Boundary Element Methods with the Software Library BEMLIB}, Chapman & Hall/CRC, Boca Raton, FL, 2002. · Zbl 1019.65097
[16] M. W. Rostami and S. D. Olson, {\it Kernel-independent fast multipole method within the framework of regularized Stokeslets}, J. Fluid. Struct., 67 (2016), pp. 60-84.
[17] D. Smith, {\it A boundary element regularized Stokeslet method applied to cilia-and flagella-driven flow}, Proc. R. Soc. Lond. A, 465 (2009), pp. 3605-3626. · Zbl 1195.76452
[18] D. Smith, {\it A nearest-neighbour discretisation of the regularized Stokeslet boundary integral equation}, J. Comput. Phys., 358 (2018), pp. 88-102. · Zbl 1381.76052
[19] C. Sorgentone and A.-K. Tornberg, {\it A highly accurate boundary integral equation method for surfactant-laden drops in 3d}, J. Comput. Phys., 360 (2018), pp. 167-191. · Zbl 1391.76456
[20] S. Veerapaneni, D. Gueyffier, D. Zorin, and G. Biros, {\it A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in \(2\) D}, J. Comput. Phys., 228 (2009), pp. 2334-2353. · Zbl 1275.76175
[21] S. Veerapaneni, A. Rahimian, G. Biros, and D. Zorin, {\it A fast algorithm for simulating vesicle flows in three dimensions}, J. Comput. Phys., 230 (2011), pp. 5610-5634. · Zbl 1419.76475
[22] L. Ying, G. Biros, and D. Zorin, {\it A high-order 3D boundary integral equation solver for elliptic pdes in smooth domains}, J. Comput. Phys., 219 (2006), pp. 247-275. · Zbl 1105.65115
[23] A. Zinchenko and R. Davis, {\it An efficient algorithm for hydrodynamical interaction of many deformable drops}, J. Comput. Phys., 157 (2000), pp. 539-587. · Zbl 0961.76055
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.