A single-layer based numerical method for the slender body boundary value problem

William H. Mitchell, Henry G. Bell, Yoichiro Mori, Laurel Ohm, Daniel Spirn

Research output: Contribution to journalArticlepeer-review

1 Scopus citations

Abstract

Fluid flows containing dilute or dense suspensions of thin fibers are widespread in biological and industrial processes. To describe the motion of a thin immersed fiber, or to describe the forces acting on it, it is convenient to work with one-dimensional fiber centerlines and force densities rather than two-dimensional surfaces and surface tractions. Slender body theories offer ways to model and simulate the motion of immersed fibers using only one-dimensional data. However, standard formulations can break down when the fiber surface comes close to intersecting itself or other fibers. In this paper we introduce a numerical method for a recently derived three-dimensional slender body boundary value problem that can be stated entirely in terms of a one-dimensional distribution of forces on the centerline. The method is based on a new completed single-layer potential formulation of fluid velocity which removes the nullspace associated with the unmodified single layer potential. We discretize the model and present numerical results demonstrating the good conditioning and improved performance of the method in the presence of near-intersections. To avoid the modeling and numerical choices involved with free ends, we consider closed fibers.

Original languageEnglish (US)
Article number110865
JournalJournal of Computational Physics
Volume450
DOIs
StatePublished - Feb 1 2022

All Science Journal Classification (ASJC) codes

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Keywords

  • Integral equations
  • Numerical methods
  • Slender body theory
  • Stokes flows

Fingerprint

Dive into the research topics of 'A single-layer based numerical method for the slender body boundary value problem'. Together they form a unique fingerprint.

Cite this