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

DOI10.1137/18M1191816MaRDI QIDQ4628392zbMATH OpenWikidataFDO

Authors Meurig T. Gallagher, Debajyoti Choudhuri, David J. Smith

Publication date 13 March 2019

Published in SIAM Journal on Scientific Computing (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/1806.01560

zbMATH Keywords

boundary integral equationquadrature error boundsregularized stokesletsStokes flow

Mathematics Subject Classification ID

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

Abstract: The method of regularized stokeslets is a powerful numerical method 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 method contains three sources of numerical error, the regularization error associated from 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 discretisation, 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: with disjoint sets ($delta>0$) being close to linear in $h_q$ and insensitive to $varepsilon$, and contained sets ($delta=0$) being quadratic in $h_q$ with inverse dependence on $varepsilon$. The practical implications of these error bounds are discussed with 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 ($deltageq0$) are revealed to be proportional to the sum of the errors for each case.

Recommendations

• A nearest-neighbour discretisation of the regularized stokeslet boundary integral equation
• Reduction of the regularization error of the method of regularized stokeslets for a rigid object immersed in a three-dimensional Stokes flow
• Regularized single and double layer integrals in 3D Stokes flow
• A novel regularization for higher accuracy in the solution of the 3-dimensional Stokes flow
• Analysis of the immersed boundary method for a finite element Stokes problem

Cites work

• A boundary element regularized Stokeslet method applied to cilia- and flagella-driven flow
• A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D
• A fast algorithm for simulating vesicle flows in three dimensions
• A fast integral equation method for solid particles in viscous flow using quadrature by expansion
• A fast numerical method for computing doubly-periodic regularized Stokes flow in 3D
• A fast, high-order algorithm for the solution of surface scattering problems: Basic implementation, tests, and applications
• A general system of images for regularized stokeslets and other elements near a plane wall
• A high-order 3D boundary integral equation solver for elliptic PDEs in smooth domains
• A highly accurate boundary integral equation method for surfactant-laden drops in 3D
• A nearest-neighbour discretisation of the regularized stokeslet boundary integral equation
• An efficient algorithm for hydrodynamical interaction of many deformable drops
• Asymptotic analysis for close evaluation of layer potentials
• Boundary Integral and Singularity Methods for Linearized Viscous Flow
• Gershgorin's Theorem and the Zeros of Polynomials
• scientific article; zbMATH DE number 1852568 (Why is no real title available?)
• Slender body theory for Stokes flows with regularized forces
• The method of images for regularized Stokeslets
• The method of regularized stokeslets
• The method of regularized Stokeslets in three dimensions:  Analysis, validation, and application to helical swimming
• Weakening accuracy dependence with the regularization parameter in the method of regularized Stokeslets

Cited in

(3)

• A nearest-neighbour discretisation of the regularized stokeslet boundary integral equation
• Passively parallel regularized stokeslets
• Single-particle dynamics in a low-Reynolds-number fluid under spherical confinement

Describes a project that uses

Uses Software

• BEMLIB

This page was built for publication: Sharp quadrature error bounds for the nearest-neighbor discretization of the regularized stokeslet boundary integral equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4628392)