No abstract

The macroscopic properties of polymeric fluids are inherited from the material properties of the fibers embedded in the solvent. The behavior of such passive fibers in flow has been of interest in a wide range of systems, including cellular mechanics, nutrient aquisition by diatom chains in the ocean, and industrial applications such as paper manufacturing. The rotational dynamics and shape evolution of fibers in shear depends upon the slenderness of the fiber and the non-dimensional "elasto-viscous" number that measures the ratio of the fluid's viscous forces to the fiber's elastic forces. For a small elasto-viscous number, the nearly-rigid fiber rotates in the shear, but when the elasto-viscous number reaches a threshhold, buckling occurs. For even larger elasto-viscous numbers, there is a transition to a "snaking behavior" where the fiber remains aligned with the shear axis, but its ends curl in, in opposite directions. These experimentally-observed behaviors have recently been characterized computationally using slender-body theory and immersed boundary computations. However, classical experiments with nylon fibers and recent experiments with actin filaments have demonstrated that for even larger elasto-viscous numbers, multiple buckling sites and coiling can occur. Using a regularized Stokeslet framework coupled with a kernel independent fast multipole method, we present simulations that capture these complex fiber dynamics.

Motivated by bacterial transport through porous media, here we study the swimming of an actuated, flexible helical filament in both three-dimensional free space and within a cylindrical tube whose diameter is much smaller than the length of the helix. The filament, at rest, has a native helical shape modeled after the geometry of a typical bacterial flagellar bundle. The finite length filament is a free swimmer, and is driven by an applied torque as well as a counter-torque (of equal strength and opposite direction) that represents a virtual cell body. We use a regularized Stokeslet framework to examine the shape changes of the flexible filament in response to the actuation as well as the swimming performance as a function of the nondimensional Sperm number that characterizes the elastohydrodynamic system. We also show that a modified Sperm number may be defined to characterize the swimming progression within a tube. Finally, we demonstrate that a helical filament whose axis is not aligned with the tube axis can exhibit centering behavior in the narrowest tubes.

We describe the implementation of optimal local radiation boundary condition sequences for second order finite difference approximations to Maxwell's equations and the scalar wave equation using the double absorbing boundary formulation. Numerical experiments are presented which demonstrate that the design accuracy of the boundary conditions is achieved and, for comparable effort, exceeds that of a convolutional perfectly matched layer with reasonably chosen parameters. An advantage of the proposed approach is that parameters can be chosen using an accurate a priori error bound.

No abstract