A Smoluchowski equation for a sheared suspension of frictionally interacting rods

想象一个这样的世界：我们饮用、涂抹或打印的物质不仅仅是简单的液体，而是在流体中漂浮着的、由无数细小且僵硬的棒状物组成的混乱舞会。这就是悬浮液（suspensions）的世界，这是一个研究粒子在混合于液体中如何运动的物理学分支。一个多世纪以来，科学家们已经了解了当这些棒状物彼此远离时，它们是如何像随风飘荡的落叶一样慵懒地漂流的。他们也知道，当这些棒状物变得拥挤，但仍能在一层薄薄的、滑溜的流体上相互滑动时（就像冰上的滑冰运动员），它们会如何表现。但还有第三种更混乱的情景：当人群变得如此密集，以至于“冰面”破碎了，会发生什么？当这些棒状物相互碰撞、失去其润滑涂层并开始相互磨擦时，会发生什么？这就是摩擦力（friction）。这正是现代流变学（研究流动行为的学科）所驱动的核心问题。理解这一点至关重要，因为它解释了为什么某些材料在受到过快搅拌时，会突然从流动的液体变成坚硬的固体块——这种现象可能会堵塞工业管道或导致3D打印机卡死。

在本文中，Genevieve Quiñones和Peter D. Olmsted致力于解决这些摩擦性碰撞带来的复杂数学问题。他们建立了一套新的规则，称为斯莫卢霍夫斯基方程（Smoluchowski equation），用以预测受剪切（搅拌）作用下的密集摩擦棒状物群体的行为。可以将其想象为将交通模拟从“车辆仅在光滑、无摩擦的道路上行驶”的模型，升级为“车辆可以发生侧滑、磨擦并锁死”的模型。作者指出，当这些棒状物变得拥挤时，它们不仅仅是在滑动；它们有时会粘附并刮擦，产生一种“摩擦转矩（frictional torque）”，从而改变它们的旋转方式。他们的研究表明，这种摩擦力起到了“刹车”的作用，抑制了棒状物通常会进行的剧烈翻滚运动，并迫使它们更整齐地与流向对齐。他们还提出了一种计算棒状物与其邻居碰撞次数的新方法，发现人群越有序，发生的碰撞就越少。虽然他们尚未建立物理实验装置来证明这一点，但他们的计算机模拟显示，加入摩擦力改变了棒状物的“舞步”，防止它们失控旋转，并可能解释了为什么某些粘稠流体在受到更强推力时会突然变得更加粘稠。

磨损棒状物的寓言

想象你正在一个拥挤的舞池中。如果空间很宽敞，每个人都可以自由地旋转和旋转。这就像是棒状物的稀释悬浮液（dilute suspension），其中的粒子彼此远离，只是偶尔发生碰撞。物理学家早已知道如何预测这种舞蹈；他们使用方程（如杰弗里方程/Jeffery's equation）来描述单个棒状物在流动液体中的旋转。

现在，想象舞池变得拥挤不堪。舞者不再仅仅是旋转，他们还在互相碰撞。在半稀释（semi-dilute）的人群中，他们可能仍能相互滑动，因为有一层薄薄的流体（如汗水或油脂）充当润滑剂。但在浓缩（concentrated）的人群中，那层润滑层可能会破裂。突然间，舞者不再是滑动，而是在互相磨蹭鞋底。这就是固体摩擦（solid friction）。

长期以来，科学家们一直难以编写这种磨损过程的数学模型。之前的尝试假设即使在棒状物磨损时，它们仍遵循与稀疏人群相同的“舞步”（杰弗里方程），只是增加了一个额外的“推力”。但当本文作者试图计算这些推力的平均效应时，他们发现数学计算的结果抵消为零。这在他们看来是不合理的，因为我们从现实生活中知道，摩擦会减慢物体的速度并改变其运动方式。

新方法：“能量账单”

Quiñones和Olmsted决定尝试另一种策略。他们没有去猜测棒状物的运动方式然后再加入摩擦力，而是从能量入手。他们使用了一种称为**瑞利函数（Rayleighian）**的方法，这本质上是一种计算整个系统“能量账单”的方法。

可以将瑞利函数想象成一个总成本计算器。它累加了两部分：

1. 自由能（Free Energy）： 棒状物希望如何排列（例如磁铁的排列方式）。
2. 耗散（Dissipation）： 由于摩擦和阻力而损失到热能中的能量。

作者意识到，要找到棒状物的真实运动，必须使这个总“能量账单”最小化。当他们在账单中计入固体摩擦（磨损）和边界润滑摩擦（带薄膜的滑动）产生的能量损失时，神奇的事情发生了。数学自然地产生了一个新的“转矩”（扭转力），而这个力并不为零。

他们发现，摩擦力起到了旋转制动器的作用。当棒状物相互磨损时，它们旋转变得更加困难。这改变了斯莫卢霍夫斯基方程，即那个预测所有棒状物取向随时间变化的母方程。

他们的发现：舞步改变了

作者运行了计算机模拟，以观察搅拌这些摩擦棒状物时会发生什么。以下是他们的发现：

• “皮划艇运动”停止了： 在无摩擦的人群中，棒状物经常进行一种被称为“皮划艇运动（kayaking）”的动作，即进行复杂的、面外轨道的倾斜和旋转。作者发现，当开启摩擦力时，这种剧烈的旋转被抑制了。棒状物停止了皮划艇式的旋转，转而倾向于与流向保持一致，就像一群鱼在直线游动。
• 与“堵塞”的联系： 他们认为，这种摩擦制动是**不连续剪切增稠（Discontinuous Shear Thickening, DST）**的关键因素。DST是一种奇特的现象，即某种液体（如玉米淀粉和水混合物）在受到猛烈撞击时会突然变成固体。作者的模型表明，随着搅拌速度加快，棒状物之间的摩擦增加，“刹车”锁死，材料发生堵塞。
• 计数碰撞： 他们还提出了一个新公式，用于估算一个棒状物接触其邻居的次数。他们发现，接触次数取决于人群的有序程度。如果棒状物完美对齐（高有序度），它们接触的邻居就较少；如果它们杂乱无章（无序），则接触更多。这是对之前仅适用于混乱无序人群的理念的一种推广。

模型的局限性

需要注意的是，本文并未声称实现的内容。作者谨慎地指出，他们的模型假设棒状物始终在相互滑动。他们目前还无法模拟棒状物由于过度卡住而形成巨大的、静止集群（静态摩擦）的瞬间。在现实世界中，这些卡住的集群很可能是材料完全堵塞的主要原因。作者承认，由于静态摩擦不会产生热量（耗散），因此不会出现在他们的能量账单中，所以他们的模型遗漏了这种“卡住”的部分。

此外，他们的结果是基于模拟和数学推导，而非新的物理实验。他们暗示，其研究结果解释了为什么摩擦很重要，但尚未在实验室中实际测量。他们还假设悬浮液是完美混合且均匀的，忽略了摩擦可能导致棒状物在特定位置聚集的可能性。

总结

简单来说，Quiñones和Olmsted建立了一个新的数学工具箱，终于能够解释厚流体中粒子间的“磨损”。他们证明了，如果忽略摩擦，就会错过阻止棒状物剧烈旋转的“刹车”。通过引入这种摩擦，他们找到了一种自洽的方法，可以预测这些拥挤的棒状物如何运动，以及它们会对容器施加多少应力。虽然他们还没有解决流体为何会完全堵塞的谜团（这需要理解“卡住”的状态），但他们提供了一个关键的拼图碎片：即在堵塞发生前，通过摩擦转矩让“舞蹈”慢下来的机制。