Skip to main content
    • Aa
    • Aa

Numerical Computation of Doubly-Periodic Stokes Flow Bounded by a Plane with Applications to Nodal Cilia

  • Franz Hoffmann (a1) and Ricardo Cortez (a1)

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 considerablemore particle transport than in-phase beating.

Corresponding author
*Corresponding author. Email addresses: (F. Hoffmann), (R. Cortez)
Hide All

Communicated by Boo-Cheong Khoo

Linked references
Hide All

This list contains references from the content that can be linked to their source. For a full set of references and notes please see the PDF or HTML where available.

[1] J. Ainley , S. Durkin , R. Embid , P. Boindala , and R. Cortez . The method of images for regularized Stokeslets. J. Comp. Phy, 227:46004616, 2008.

[2] C. W. J. Beenakker . Ewald sum of the Rotne-Prager tensor. J. Chem. Phys, 85(3), 1986.

[3] J. R. Blake . A note on the image system for a Stokeslet in a no-slip boundary. Proc. Cambridge Phil. Soc, 70:303310, 1971.

[4] J. R. Blake . A spherical envelope approach to ciliary propulsion. J. FluidMech, 46(1):199208, 1971.

[5] J. R. Blake . A model for the micro-structure in ciliated organisms. J. Fluid Mech, 55(1):123, 1972.

[6] J. J. Breunig , J. I. Arellano , and P. Rakic . Cilia in the brain: going with the flow. Nat. Neurosci, 13(6):654655, 2010.

[7] C-Y. Chen , C-Y. Chen , C-Y. Lin , and Y-T. Hu . Magnetically actuated artificial cilia for optimum mixing performance in microfluidics. Lab Chip, 13:28342839, 2013.

[8] R. Cortez . The method of regularized Stokeslets. SIAM J. Sci. Comp, 23(4):12041225, 2001.

[9] R. Cortez , L. Fauci , and A. Medovikov . The method of regularized Stokeslets in Three Dimensions: Analysis,Validation and Application to Helical Swimming. Phys. Fluids, 17:114, 2005.

[10] R. Cortez and F. Hoffmann . A fast numerical method for computing doubly-priodic regularized Stokes flow in 3D. J. Comp. Phy, 258, 2014.

[11] R. Cortez and M. Nicholas . Slender body theory for stokes flows with regularized forces. Comm. App.Math. Com. Sc, 7(1):3362, 2012.

[12] R. Cortez and D. Varela . A general system of images for regularized stokeslets and other elements near a plane wall. J. Comp. Phy, 285:4154, 2015.

[14] Y. Ding , J. C. Nawroth , M. J. McFall-Ngai , and E. Kanso . Mixing and transport by ciliary carpets: a numerical study. J. Fluid Mech, 743:124140, 2014.

[15] M. T. Downton and H. Stark . Beating kinematics of magneticallly actuated cilia. EPL, 85(4):44002–p1–p6, 2009.

[16] J. Elgeti and G. Gompper . Emergence of metachronal waves in cilia arrays. Proc. Nat. A. Sci, 110(12):44704475, 2013.

[17] G. R. Fulford and J. R. Blake . Muco-ciliary transport in the lung. Theor. Biol, 121(4):381402, 1986.

[18] E. M. Gauger , M. T. Downton , and H. Stark . Fluid transport at low Reynolds number with magnetically actuated artificial cilia. Eur. Phys. J. E, 28:231242, 2009.

[19] S. Gueron and K. Levit-Gurevich . Computation of the internal forces in cilia: Application to ciliary motion, the effects of viscosity, and cilia interactions. Biophy. J, 74:16581676, 1998.

[20] S. Gueron and K. Levit-Gurevich . Energetic considerations of ciliary beating and the advantage of metachronal coordination. Proc. Natl. Acad. Sci, 96(22):1224012245, 1999.

[22] J. D. Hunter . Matplotlib: A 2D graphics environment. Comput. Sci. Eng, 9(3):9095, 2007.

[23] J. Hussong , N. Schorr , J. Belardi , O. Prucker , J. Rheb , and J. Westerweel . Experimental investigation of the flow induced by artificial cilia. Lab Chip, 11:20172022, 2011.

[25] A. Keißner and C. Brücker . Directional fluid transport along artificial ciliary surfaces with base-layer actuation of counter-rotating orbital beating patterns. Soft Matter, 8:53425349, 2012.

[26] S. N. Khaderi , J. M. J. den Toonder , and P. R. Onck . Microfluidic propulsion by the metachronal beating of magnetic artificial cilia: a numerical analysis. J. Fluid Mech, 688:4465, 2011.

[27] G. Kokot , M. Vilfan , N. Osterman , A. Vilfan , B. Kavčič , I. Poberaj , and D. Babič . Measurement of fluid flow generated by artificial cilia. Biomicrofluidics, 5:034103, 2011.

[28] P. Lenz and A. Ryskin . Collective effects in ciliar arrays. Phys. Bio, 3:285294, 2006.

[29] N. Liron and S. Mochon . The discrete-cilia approach to propulsion of ciliated micro-organisms. J. Fluid Mech, 75(3):593607, 1975.

[30] W. F. Marshall and S. Nonaka . Cilia: Tuning in to the cell's antenna. Curr. Bio, 16(15):604614, 2006.

[31] S. M. Mitran . Metachronal wave formation in a model of pulmonary cilia. Comput. Struct, 85:763774, 2007.

[32] H-N. Nguyen and K. Leiderman . Computation of the singular and regularized image systems for doubly-periodic stokes flow in the presence of a wall. J. of Comp. Phy, 297:442461, 2015.

[33] N. Osterman and A. Vilfan . Finding the ciliary beating pattern with optimal efficency. Proc. Nat. A. Sci, 108(38):1572715732, 2011.

[34] A. R. Shields , B. L. Fiser , B. A. Evans , M. R. Falco , S. Washburn , and R. Superfine . Biometric cilia arrays generate simultaneous pumping and mixing regimes. Proc. Nat. A. Sci, 107(36):1567015675, 2010.

[35] D.J. Smith , J. R. Blake , and E. A. Gaffney . Fluid mechanics of nodal flow due to embryonic primary cilia. J. R. Soc. Interface, 5:567573, 2008.

[36] D.J. Smith , E. A. Gaffney , and J. R. Blake . Discrete cilia modelling with singularity distributions: application to the embryonic node and the airway surface liquid. Bull. Math. Bio, 69:14771510, 2007.

[37] M. Villfan , A. Potočnik , B. Kavačič , N. Osterman , I. Poberaj , A. Vilfan , and D. Babič . Self-assembled artificial cilia. Proc. Nat. A. Sci, 107(5):18441847, 2010.

[38] Y. Wang , Y. Gao , H. Wyss , P. Anderson , and J. den Tonder . Out of the cleanroom, self-assembled magnetic artificial cilia. Lab Chip, 13:33603366, 2013.

Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Communications in Computational Physics
  • ISSN: 1815-2406
  • EISSN: 1991-7120
  • URL: /core/journals/communications-in-computational-physics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 9 *
Loading metrics...

Abstract views

Total abstract views: 36 *
Loading metrics...

* Views captured on Cambridge Core between 6th July 2017 - 27th July 2017. This data will be updated every 24 hours.