Simulating flexible fiber suspensions using a scalable immersed boundary algorithm

Jeffrey K. Wiens,John M. Stockie
DOI: https://doi.org/10.1016/j.cma.2015.02.026
2014-07-29
Abstract:We present an approach for numerically simulating the dynamics of flexible fibers in a three-dimensional shear flow using a scalable immersed boundary (IB) algorithm based on Guermond and Minev's pseudo-compressible fluid solver. The fibers are treated as one-dimensional Kirchhoff rods that resist stretching, bending, and twisting, within the generalized IB framework. We perform a careful numerical comparison against experiments on single fibers performed by S. G. Mason and co-workers, who categorized the fiber dynamics into several distinct orbit classes. We show that the orbit class may be determined using a single dimensionless parameter for low Reynolds flows. Lastly, we simulate dilute suspensions containing up to hundreds of fibers using a distributed- memory computer cluster. These simulations serve as a stepping stone for studying more complex suspension dynamics including non-dilute suspensions and aggregation of fibers (also known as flocculation).
Fluid Dynamics,Computational Engineering, Finance, and Science
What problem does this paper attempt to address?