Nodal cilia-driven flow: Development of a computational model of the nodal cilia axoneme

Cilia-driven nodal flow is important in the determination of left-right asymmetry in the body. Several theoretical and computational models have been proposed to explain the mechanics of ciliary motion, although the full mechanism remains unknown. Here, we developed a three-dimensional nodal cilia axoneme model using a finite element-boundary element coupling method, and investigated the mechanics of nodal ciliary motion. We found that the rotational orbit was strongly dependent on the dynein activation frequency. We also investigated flow field generated by the ciliary rotation, and the flow strength decayed as r^-3 at the far field from the cilium. Our numerical results also suggest that experimentally observed tilt angle θ=2π/9 is sufficiently large to make a leftward flow. These findings are helpful in better understanding ciliary motion and nodal flow.