Department of Mechanical Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
J R Soc Interface. 2012 Mar 7;9(68):571-85. doi: 10.1098/rsif.2011.0355. Epub 2011 Aug 24.
The emergence of coherent structures, large-scale flows and correlated dynamics in suspensions of motile particles such as swimming micro-organisms or artificial microswimmers is studied using direct particle simulations. A detailed model is proposed for a slender rod-like particle that propels itself in a viscous fluid by exerting a prescribed tangential stress on its surface, and a method is devised for the efficient calculation of hydrodynamic interactions in large-scale suspensions of such particles using slender-body theory and a smooth particle-mesh Ewald algorithm. Simulations are performed with periodic boundary conditions for various system sizes and suspension volume fractions, and demonstrate a transition to large-scale correlated motions in suspensions of rear-actuated swimmers, or Pushers, above a critical volume fraction or system size. This transition, which is not observed in suspensions of head-actuated swimmers, or Pullers, is seen most clearly in particle velocity and passive tracer statistics. These observations are consistent with predictions from our previous mean-field kinetic theory, one of which states that instabilities will arise in uniform isotropic suspensions of Pushers when the product of the linear system size with the suspension volume fraction exceeds a given threshold. We also find that the collective dynamics of Pushers result in giant number fluctuations, local alignment of swimmers and strongly mixing flows. Suspensions of Pullers, which evince no large-scale dynamics, nonetheless display interesting deviations from the random isotropic state.
采用直接粒子模拟研究了游动颗粒(如游泳微生物或人工微泳者)悬浮液中相干结构、大规模流动和相关动力学的出现。为一种细长的杆状颗粒提出了详细的模型,该颗粒通过在其表面施加规定的切向应力来推动自身在粘性流体中运动,并设计了一种使用细长体理论和光滑粒子网格 Ewald 算法在这种颗粒的大规模悬浮液中高效计算流体动力学相互作用的方法。模拟在各种系统尺寸和悬浮液体积分数下采用周期性边界条件进行,并在临界体积分数或系统尺寸以上的后驱游动者(或推进器)悬浮液中证明了向大规模相关运动的转变,而在前驱游动者(或拉拔器)悬浮液中则没有观察到这种转变。这种转变在颗粒速度和被动示踪剂统计中最为明显。这些观察结果与我们之前的平均场动理论的预测一致,其中之一表明,当线性系统尺寸与悬浮液体积分数的乘积超过给定阈值时,推进器的均匀各向同性悬浮液中将出现不稳定性。我们还发现,推进器的集体动力学导致了巨大的数量波动、泳者的局部对齐和强烈的混合流动。尽管没有大规模动力学的拉拔器悬浮液仍表现出与随机各向同性状态的有趣偏离。