TY - JOUR T1 - Numerical Computation of Doubly-Periodic Stokes Flow Bounded by a Plane with Applications to Nodal Cilia JO - Communications in Computational Physics VL - 3 SP - 620 EP - 642 PY - 2017 DA - 2017/09 SN - 22 DO - http://doi.org/10.4208/cicp.OA-2016-0151 UR - https://global-sci.org/intro/article_detail/cicp/9975.html KW - AB -

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 velocity u(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 considerable more particle transport than in-phase beating.