Hydrodynamic cascade drives tumbling in sheared colloidal rod suspensions

本研究揭示，在先前被认为可忽略的半稀浓度区间内，流体动力学相互作用驱动剪切胶体棒悬浮液发生集体翻滚级联，该过程破坏了流动取向并显著增加了粘度，从而需要对现有本构模型进行修正。

要点

想象一个拥挤的舞池，每个人都握着一根长而僵硬的棍子。如果音乐缓慢且人群稀疏，每位舞者可以自由地旋转手中的棍子，主要受自身随机运动的引导。但当音乐加速、人群变得更加密集时会发生什么？

本文正是研究这一情境，只不过研究对象不是舞者，而是悬浮在液体中的微观棒状颗粒（胶体棒）；而“音乐”则对应液体被搅拌或“剪切”的过程。

以下是研究人员发现的故事，以简明的方式阐述：

旧有观念：“液体太稀薄，无关紧要”

长期以来，科学家认为，当这些棒状颗粒处于半密集状态（既不太拥挤也不太稀疏）时，它们之间的液体就像一位沉默的旁观者。他们相信，如果你推动液体，棒状颗粒会像溪流中的树叶一样顺着流动方向排列，而液体自身的运动并不会真正改变棒状颗粒的行为。他们认为棒状颗粒主要是独立的，只有发生物理接触时才会相互碰撞。

新发现：“多米诺效应”

研究人员利用强大的计算机模拟来观察这些棒状颗粒的运动。他们发现，液体并非一位沉默的旁观者。事实上，它就像一支混乱乐队的指挥。

他们发现的机制如下：

1. 翻滚：当液体快速流动时，一根棒状颗粒试图顺着流动方向排列。但就在它即将完美对齐时，它会被推离轨道，不得不“翻滚”（翻转）以重新开始这一过程。
2. 涟漪：当一根棒状颗粒翻滚时，它会搅动周围的液体，形成微小的漩涡或涟漪。
3. 级联：这股涟漪击中邻近的另一根棒状颗粒，迫使它也发生翻滚。第二根棒状颗粒随后搅动液体，导致第三根棒状颗粒也发生翻滚。
4. 连锁反应：这就形成了一种级联。一次翻滚会触发邻近颗粒之间的一系列翻滚连锁反应。

作者将这种现象称为“流体动力学级联”。它就像一场多米诺骨牌游戏，液体是那只看不见的手将它们全部推倒，而不是它们自己倒下。

令人惊讶的结果

由于这种多米诺效应，棒状颗粒的行为与科学家的预测大相径庭：

• 它们无法对齐：棒状颗粒并没有全部指向同一方向（这会使液体更容易流动），而是不断被邻居的翻滚动作推离对齐状态。最终它们指向各个方向，包括垂直于流动方向的侧面。
• 液体变稠：由于棒状颗粒不断翻滚并努力保持对齐，液体变得难以搅拌。“粘度”（稠度）急剧上升。
• 应力发生变化：液体施加的力以特定方式发生变化，这与近期关于病毒状棒状颗粒的真实实验结果相符，而之前的理论无法解释这一现象。

类比：交通堵塞

将棒状颗粒想象成高速公路上的汽车。

• 旧理论：如果汽车高速行驶，它们都会保持在各自的车道内平稳行驶。它们之间的空气无关紧要。
• 新发现：当一辆汽车为了避开碰撞而突然变道（翻滚）时，它会产生一股气流，将旁边的汽车也推得变道。那辆车又推动下一辆车。突然间，整条高速公路变成了一片混乱，汽车左右摇摆。交通急剧减速（粘度增加），汽车不再沿直线行驶。

为何这很重要

该论文声称，长期以来，科学家们忽略了这些棒状颗粒之间的“风”（流体动力学相互作用），因为他们认为这种作用太弱而不重要。这项研究证明，在高速和特定密度下，这种“风”实际上是导致混乱的主要驱动力。

这一发现解释了为什么一些真实世界的实验（例如涉及病毒颗粒的实验）显示出粘稠且混乱的行为，而旧有的数学模型无法预测这一点。作者得出结论，我们需要重写描述这些材料的规则（本构模型），承认液体本身会产生连锁反应，从而决定整个群体的运动方式。

简而言之：液体不仅仅是背景；它是积极的驱动者，将一群独立的棒状颗粒转变为一个混乱、翻滚的群体，使流体变得比我们想象的更加粘稠和复杂。

技术摘要

技术摘要：流体动力学级联驱动剪切胶体棒悬浮液中的翻滚运动

问题陈述
剪切流动下胶体棒悬浮液的动力学仍是软物质物理中的核心挑战，这主要归因于粒子相互作用的各向异性和长程特性。尽管流体动力学相互作用（HI）已被确立为球形胶体和乳液流变学的关键因素，但其在细长棒悬浮液中的作用历史上一直存在争议。Shaqfeh 和 Fredrickson 的开创性理论工作表明，对于大长径比的棒，半稀浓度区间的应力中 HI 的贡献可忽略不计，导致许多本构模型将其排除。然而，近期针对高佩克莱特数（$Pe$）下 fd 病毒悬浮液的实验观察到了高剪切粘度和微弱的流动取向，这与忽略 HI 的理论相矛盾。在半稀、高-$Pe$ 区间内，流体动力学耦合影响翻滚和取向动力学的具体机制仍知之甚少。

方法论
为填补这一空白，作者采用了基于粒子的布朗动力学（BD）模拟，研究剪切流动下的半稀胶体棒悬浮液。

• 模型系统： 棒被建模为由 $N=10$ 个无扭矩珠子组成的链，通过强弯曲和拉伸约束连接，总长度 $L = 2aN$。系统采用 Lees–Edwards 周期性边界条件。
• 流体动力学相互作用： 为应对全 HI 计算（$O(N^2)$）带来的高昂计算成本，作者实施了一种基于棒间距的截断判据，并结合 Rotne–Prager–Yamakawa (RPY) 迁移率张量。HI 仅在同根棒或相邻棒（定义为质心间距 $r_c \leq L$）的珠子之间计算。
• 参数空间： 模拟探索了广泛的佩克莱特数范围（$1 \leq Pe \leq 10^6$）以及半稀浓度区间（$nL^3 \leq 8.9$）。
• 对比分析： 将包含 HI 的模拟结果与明确忽略 HI 的“自由排水”（FD）模拟结果进行系统比较，以隔离流体动力学耦合的影响。

关键结果
研究表明，流体动力学相互作用触发了一种集体性的“级联”翻滚事件，从根本上改变了悬浮液的动力学和流变学：

1. 翻滚的流体动力学级联： 在半稀条件下，单根棒的翻滚会诱导流动扰动，促使邻近的棒发生翻滚。这产生了一种级联效应，使得翻滚事件变得相互关联。
2. 流动取向的破坏： 与高 $Pe$ 下棒与流动对齐的 FD 模拟不同，HI 情况下的流动取向序参数（$S$）显著降低。作者观察到在特定浓度下 $S(Pe)$出现峰值（$Pe \approx 10^3 - 2 \times 10^3$），这一特征在 FD 模型中不存在，但与近期的实验观察一致。
3. 涡度偏置： HI 驱动棒取向朝向涡度方向，而非速度梯度方向。这体现在第二法向应力差（$N_2$）的显著增加以及双轴性参数（$B$）的非单调行为上。
4. 翻滚频率的标度律： 翻滚频率（$\rho$）表现出两个区间。在较低 $Pe$下，它遵循布朗标度$\rho/Dr \sim Pe^{2/3}$。超过临界$Pe_c$ 后，HI 占主导地位，频率与 $Pe$呈线性标度（$\rho/Dr \sim Pe$），趋近于非布朗棒的行为。作者推导出了一个标度律（公式 2），基于浓度参数$\Theta = nL^3 / \ln(L/4a)$ 统一了这两个区间。
5. 流变学后果： 级联机制导致高 $Pe$下剪切粘度（$\eta^*$）和法向应力差显著增加。不同浓度的粘度曲线在通过推导出的标度函数$F(\Theta, Pe)$ 归一化后发生坍塌，表明流变学偏差是由级联动力学而非直接的多体应力贡献驱动的。

意义与主张
本文主张，流体动力学耦合不仅仅是次要效应，而是控制高度各向异性、半稀悬浮液集体动力学的关键机制。作者证明：

• 先前针对静态构型计算得出的 HI 对应力“可忽略”的贡献，未能解释 HI 对系统演化的间接动力学效应。
• 近期实验中观察到的与经典稀溶液理论的偏差（例如 fd 病毒悬浮液中增强的粘度和取向峰值），可以通过这种流体动力学促进的级联效应来解释。
• 现有的胶体棒本构模型在很大程度上忽略了半稀条件下的 HI，需要修订以纳入这些集体流体动力学现象。

该工作建立了一个将微观级联动力学与宏观流变学联系起来的框架，表明类似的集体流体动力学机制可能在其他高度各向异性和半柔性丝状系统中具有相关性。