TY - JOUR
T1 - Remarks on regularized stokeslets in slender body theory
AU - Ohm, Laurel
N1 - Funding Information:
Funding: This research was funded by NSF postdoctoral fellowship DMS-2001959.
Publisher Copyright:
© 2021 by the author. Licensee MDPI, Basel, Switzerland.
PY - 2021/8
Y1 - 2021/8
N2 - We remark on the use of regularized Stokeslets in the slender body theory (SBT) approximation of Stokes flow about a thin fiber of radius ɛ > 0. Denoting the regularization parameter by δ, we consider regularized SBT based on the most common regularized Stokeslet plus a regularized doublet correction. Given sufficiently smooth force data along the filament, we derive L∞ bounds for the difference between regularized SBT and its classical counterpart in terms of δ, ɛ, and the force data. We show that the regularized and classical expressions for the velocity of the filament itself differ by a term proportional to log(δ/ɛ); in particular, δ = ɛ is necessary to avoid an O(1) discrepancy between the theories. However, the flow at the surface of the fiber differs by an expression proportional to log(1 + δ2 /ɛ2), and any choice of δ ∝ ɛ will result in an O(1) discrepancy as ɛ → 0. Consequently, the flow around a slender fiber due to regularized SBT does not converge to the solution of the well-posed slender body PDE which classical SBT is known to approximate. Numerics verify this O(1) discrepancy but also indicate that the difference may have little impact in practice.
AB - We remark on the use of regularized Stokeslets in the slender body theory (SBT) approximation of Stokes flow about a thin fiber of radius ɛ > 0. Denoting the regularization parameter by δ, we consider regularized SBT based on the most common regularized Stokeslet plus a regularized doublet correction. Given sufficiently smooth force data along the filament, we derive L∞ bounds for the difference between regularized SBT and its classical counterpart in terms of δ, ɛ, and the force data. We show that the regularized and classical expressions for the velocity of the filament itself differ by a term proportional to log(δ/ɛ); in particular, δ = ɛ is necessary to avoid an O(1) discrepancy between the theories. However, the flow at the surface of the fiber differs by an expression proportional to log(1 + δ2 /ɛ2), and any choice of δ ∝ ɛ will result in an O(1) discrepancy as ɛ → 0. Consequently, the flow around a slender fiber due to regularized SBT does not converge to the solution of the well-posed slender body PDE which classical SBT is known to approximate. Numerics verify this O(1) discrepancy but also indicate that the difference may have little impact in practice.
KW - Error analysis
KW - Regularized Stokeslets
KW - Slender body theory
UR - http://www.scopus.com/inward/record.url?scp=85112743484&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85112743484&partnerID=8YFLogxK
U2 - 10.3390/fluids6080283
DO - 10.3390/fluids6080283
M3 - Article
AN - SCOPUS:85112743484
SN - 2311-5521
VL - 6
JO - Fluids
JF - Fluids
IS - 8
M1 - 283
ER -