The method of images for regularized Stokeslets

The image system for the method of regularized Stokeslets is developed and implemented. The method uses smooth localized functions to approximate a delta distribution in the derivation of the fluid flow due to a concentrated force. In order to satisfy zero-flow boundary conditions at a plane wall, the method of images derived for a standard (singular) Stokeslet is extended to give exact cancellation of the regularized flow at the wall. As the regularization parameter vanishes, the expressions reduce to the known images for singular Stokeslets. The advantage of the regularized method is that it gives bounded velocity fields even for isolated forces or for distributions of forces along curves. These are useful in the simulation of ciliary beats, flagellar motion, and particle suspensions. The expression relating force and velocity can be inverted to find the forces that generate a given velocity boundary condition. The latter is exemplified by modeling a cilium as a filament moving in a three-dimensional flow. The cilium velocity at various times is constructed from known data and used to determine the force field along the filament. Those forces can then reproduce the flow everywhere. The validity of the method is evaluated by computing the drag on a sphere moving near a wall. Comparisons with known expressions for the drag show that the method gives accurate results for spheres even within a distance from the wall equal to the surface discretization size. (C) 2008 Elsevier Inc. All rights reserved.