On an angle-averaged Neumann-to-Dirichlet map for thin filaments

We consider the Laplace equation in the exterior of a thin filament in $\mathbb{R}^3$ and perform a detailed decomposition of a notion of slender body Neumann-to-Dirichlet (NtD) and Dirichlet-to-Neumann (DtN) maps along the filament surface. The decomposition is motivated by a filament evolution equation in Stokes flow for which the Laplace setting serves as an important toy problem. Given a general curved, closed filament with constant radius $\epsilon>0$, we show that both the slender body DtN and NtD maps may be decomposed into the corresponding operator about a straight, periodic filament plus lower order remainders. For the straight, periodic filament, both the slender body NtD and DtN maps are given by explicit Fourier multipliers and it is straightforward to compute their mapping properties. The remainder terms are lower order in the sense that they are small with respect to $\epsilon$ or smoother. While the strategy here is meant to serve as a blueprint for the Stokes setting, the Laplace problem may be of independent interest.Comment: 56 pages, 1 figur

An integral model based on slender body theory, with applications to curved rigid fibers

We propose a novel integral model describing the motion of curved slender fibers in viscous flow, and develop a numerical method for simulating dynamics of rigid fibers. The model is derived from nonlocal slender body theory (SBT), which approximates flow near the fiber using singular solutions of the Stokes equations integrated along the fiber centerline. In contrast to other models based on (singular) SBT, our model yields a smooth integral kernel which incorporates the (possibly varying) fiber radius naturally. The integral operator is provably negative definite in a non-physical idealized geometry, as expected from PDE theory. This is numerically verified in physically relevant geometries. We propose a convergent numerical method for solving the integral equation and discuss its convergence and stability. The accuracy of the model and method is verified against known models for ellipsoids. Finally, a fast algorithm for computing dynamics of rigid fibers with complex geometries is developed