Numerical Computation of Doubly-Periodic Stokes Flow Bounded by a Plane with Applications to Nodal Cilia. (6th July 2017)

• Record Type: Journal Article
• Authors: Hoffmann, Franz; Cortez, Ricardo
• Is Part Of: Communications in computational physics. Volume 22:Number 3(2017:Sep.)
• Page Start: 620
• Page End: 642
• Publication Date: 2017-07-06
• Subjects: 76D07, -- 76Z99, -- 92-08; Stokes flow, -- Ewald summation, -- FFT; Mathematical physics -- Data processing -- Periodicals; Physics -- Data processing -- Periodicals; 530.150285
• DOI: 10.4208/cicp.OA-2016-0151
• Languages: English
• ISSNs: 1815-2406

Abstract

Abstract: A numerical method is presented for the computation of externally forced Stokes flows bounded by the plane z =0 and satisfying periodic boundary conditions in the x and y directions. The motivation for this work is the simulation of flows generated by cilia, which are hair-like structures attached to the surface of cells that generate flows through coordinated beating. Large collections of cilia on a surface can be modeled using a doubly-periodic domain. The approach presented here is to derive a regularized version of the fundamental solution of the incompressible Stokes equations in Fourier space for the periodic directions and physical space for the z direction. This analytical expression forû ( k, m ; z ) can then be used to compute the fluid velocityu ( x, y, z ) via a two-dimensional inverse fast Fourier transform for any fixed value of z . Repeating the computation for multiple values of z leads to the fluid velocity on a uniform grid in physical space. The zero-flow condition at the plane z =0 is enforced through the use of images. The performance of the method is illustrated by numerical examples of particle transport by nodal cilia, which verify optimal particle transport for parameters consistent with previous studies. The results also show that for two cilia in the periodic box, out-of-phase beating produces considerablemore particle transport than in-phase beating.