🔬 condensed matter

Rheology of dense suspensions of granular spherocylinders by particle-based simulation

本文提出了一种基于粒子的模拟模型，该模型成功预测了剪切流下致密颗粒长圆柱体悬浮液的粘度和微观结构，揭示了这些特性如何随体积分数和长径比而演变，并证实了有限的实验数据。

原作者： Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

发布于 2026-07-31

📖 1 分钟阅读☕ 轻松阅读

原作者： Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下这样一个世界：你所倾倒、泵送或混合的东西不仅仅是水或沙子，而是一团由漂浮在液体中的微小、刚性细棒组成的混沌汤。这就是流变学（rheology）的领域——研究物质如何流动的科学。你可能已经在现实中见过这种现象：比如一条被倒塌树木堵塞、形成巨大“木头堆积”的河流，或者是用于造纸的浓稠粘稠浆料。在这些情况下，颗粒由于过于拥挤而不断相互碰撞，它们的形状与数量一样重要。如果你尝试搅拌一桶弹珠，它们会轻易地滚过彼此；但如果你尝试搅拌一桶长而细的杆状物，它们会缠绕在一起、锁死在一起，突然间整个混合物就像固体一样。科学家们早已了解圆球在这些拥挤环境中的行为，但对于长棒状颗粒的规律一直以来都有些神秘，尤其是在它们紧密堆积时。理解这一点至关重要，因为这些“黏性棒状物”混合物无处不在，从工业制造到火山内部熔岩的流动皆是如此。

在这篇论文中，一个研究小组构建了一个超级智能的计算机模拟，作为这些拥挤棒状悬浮液的虚拟实验室。他们没有在实验室里混合真实的化学品，而是创造了一个充满数千个“球柱体”（spherocylinders）——即看起来像胶囊或药丸（两端为半球形的圆柱体）的粒子——的数字世界。他们将这些胶囊置于虚拟流体中漂浮，然后开始进行“剪切”（shearing），这是一种专业说法，指的是在保持底层静止的同时，将流体顶层向侧面拖动，就像在吐司上抹黄油一样。其目标是观察随着他们向空间内填充更多棒状物以及增加棒状物长度，混合物的粘度（viscosity）会如何变化。

模拟揭示了一些迷人的行为。首先，当他们开始剪切时，混合物会在一瞬间变得极其粘稠——出现了一个“粘度峰值”——然后才稳定下来。这是因为棒状物最初是朝向随机方向杂乱无章的，像是在进行一场“冲撞舞”（mosh pit）一样互相碰撞。一旦流动开始，棒状物就开始排列，就像士兵列队行进一样，混合物也变得更容易搅拌了。研究人员发现，随着他们填充的棒状物越多（体积份额， $\phi$ ），混合物就变得越厚，最终达到一个使其发生“阻塞”（jamming）并停止流动的点。有趣的是，棒状物越长（长径比 $A$ ，范围从代表球体的 1 到代表极长细棒的 20），这种阻塞发生得就越早。例如，对于长度为 5 的棒状物，阻塞点大约在体积份额 0.57 处；但对于长度为 20 的棒状物，混合物在约 0.41 的较低密度下就会发生阻塞。

团队还发现，这些棒状物不仅仅是排列，它们还会专门沿着流向进行排列。他们测量了一个“序参数”（order parameter, $S$ ）来观察它们的对齐程度。在低密度下，棒状物会摇摆并进行环形旋转（称为杰弗里轨道，Jeffery orbits）；但随着人群变得更加密集，它们会锁定在沿流向的直线路径上。然而，如果人群变得太密集，棒状物实际上会感到“挫败”，无法保持完美的对齐，从而导致有序度略微下降。研究人员还分解了其中的作用力，表明当混合物较稀时，液体自身的阻力起主要作用；但随着变得拥挤，棒状物之间的物理碰撞和摩擦成为了主导力量，接管了使混合物变稠的任务。

通过创建一个平衡了接触、摩擦和流体阻力的模型，作者提供了一个预测这些复杂混合物行为的新工具。他们并非仅仅靠猜测，而是利用不同的粒子数量、不同的棒状物长度和不同的填充密度进行了数千次模拟。他们的结果表明，虽然“物质越多，流体越稠”的基本理念成立，但颗粒的形状会显著改变混合物锁死的临界点。这项工作并不声称解决了流体力学的每一个谜题，但它为理解如何在现实世界中管理这些棘手的、充满棒状物的流体（从防止河流中的木头堆积到设计更好的工业浆料）提供了坚实的模拟基础。

技术摘要：基于粒子模拟的致密颗粒长柱体悬浮液流变学研究

问题陈述
棒状颗粒组成的致密悬浮液在自然和工业领域广泛存在，从河流中的木材阻塞到岩浆中的晶体悬浮液，再到造纸业中的纤维网络。虽然球形悬浮液的流变学特性已得到充分表征（特别是通过 Boyer 等人对非布朗悬浮液与干颗粒流的统一描述），但对于拉长形颗粒的行为理解仍不完整。与颗粒取向相关的额外自由度引入了复杂的微观结构现象，如杰弗里轨道（Jeffery orbits）和剪切诱导排列，这些现象显著改变了有效粘度和应力响应。由于制备和测量困难，高浓度区域（高固体体积分数 $\phi$ ）的现有实验数据非常有限。此外，以往的模拟方法面临权衡：斯托克斯动力学（Stokesian Dynamics, SD）由于需要进行处理多体流体力学的矩阵求逆，在处理致密系统时计算成本极高；而离散元法（Discrete Element Method, DEM）模型往往过度简化了颗粒几何形状（例如使用球链模型）或忽略了必要的流体力学作用力。因此，需要一种计算上可行且能忠实捕捉致密非布朗长柱体悬浮液在剪切下主导微观物理特性的粒子基模拟模型。

方法论
作者开发了一种用于模拟非弯曲、非布朗长柱体（具有半球形端帽的圆柱体）在简单剪切流下致密悬浮液的粒子模拟模型。该模型运行在粒子惯性存在但可忽略不计的机制下。

运动方程： 使用 Velocity-Verlet 算法更新粒子轨迹，对平移和旋转运动方程进行积分。这些方程考虑了成对接触力、成对润滑力和单体流体力学作用力（阻力、升力和转矩）。
流体力学： 单体作用力采用长球体（prolate spheroids）的解析解作为长柱体的近似。这包括与取向相关的阻力和剪切诱导升力。流体力学转矩能够重现杰弗里轨道。
接触与润滑：
- 接触： 建模为带有库仑摩擦的胡克定律模型，包括法向和切向的弹簧-阻尼分量。
- 润滑： 开发了一种新型方案来计算长柱体之间的润滑力。为了解决端部轴向连接处曲率不连续的问题，模型同时计算三种相互作用分量（轴-轴、端-轴和端-端）。这些分量根据粒子的相对位置进行加权，以确保不同相互作用类型之间的平滑过渡，从而解决了先前模型（如 Yamane 等人、Butler 和 Shaqfeh）中存在的间断问题。
动态时间步长： 该研究的一个关键方法论创新是动态时间步长计算。时间步长由接触相互作用、润滑力和润滑转矩中最小的特征时间尺度决定。这在确保数值稳定性的同时，避免了不必要的小时间步长，提高了计算效率。
模拟设置： 系统在流动方向和梯度方向使用周期性边界条件，并使用 Lees-Edwards 边界进行剪切。采用半径比例为 1:1.4 的双分散混合物以防止结晶。模拟涵盖了长径比（ $A$ ）从 1 到 20 以及体积分数（ $\phi$ ）从稀疏到接近阻塞（jamming）的全过程。

关键结果

粘度行为： 模型预测稳态粘度（ $\eta_s$ $η_{s}$ ）随体积分数（ $\phi$ $ϕ$ ）和长径比（ $A$ $A$ ）的增加而系统性上升。
- 瞬态响应： 在剪切启动时，观察到粘度峰值，这归因于初始的随机颗粒取向分布。随着颗粒排列，该峰值消退并达到稳态。相对于稳态，该峰值的幅度随 $\phi$ 的增加而增大。
- 稳态： 稳态粘度遵循幂律发散 $\eta_s \sim (1 - \phi/\phi_m)^{-\varsigma}$ 。阻塞体积分数 $\phi_m$ 随 $A$ 的增加而降低（从 $A=5$ 时的 $\approx 0.57$ 降至更高 $A$ 时更低的值）。模型证实了稳态粘度的速率无关性，这与颗粒悬浮液的框架一致。
- 力贡献： 在低 $\phi$ 时，粘度由流体应力和润滑力主导。随着 $\phi$ 增加，接触力（包括法向和切向）变得占主导地位，并在 $\phi \approx 0.5$ 附近最终超过润滑贡献。
微观结构与排列：
- 有序参数： 取向有序参数（ $S$ ）随 $\phi$ 增加直至一个临界点，之后由于几何挫折（geometric frustration）阻止了完美的排列，该参数会略有下降。
- 排列方向： 颗粒优先沿流动方向（ $x$ 轴）排列，相对于流-梯度平面的倾斜角（ $\phi_v$ ）较小。该倾斜角随 $\phi$ 增加至某一阈值后开始减小。
- 长径比效应： 更高的长径比会导致在较低的体积分数下产生更强的排列，但也导致更低的阻塞点。对于极高的 $A$ ，倾斜角可以忽略不计。
- 瞬态与稳态： 初始粘度峰值与系统从随机状态向剪切排列状态转变所需的时间相关。

意义与主张
本文声称提供了一个探索致密颗粒棒状悬浮液流体力学、微观结构和流变学的基本工具。通过成功模拟从稀疏到致密机制的转变，并捕捉流体力学作用、摩擦接触和颗粒几何形状之间的相互作用，该模型佐证了现有的有限的颗粒棒实验流变学数据。

作者指出，其工作将此前已建立的球形颗粒悬浮液流变学理解扩展到了更广泛的长柱体颗粒类别。该模型提供了一个平台，用于未来研究如何将阻塞点与微观结构细节联系起来，调查正应力，以及探索在非均匀条件下的行为。论文将其贡献谦逊地定位为迈向实现这些材料完整流体力学描述的一步，并将阻塞点与排列及摩擦关系的全面量化推迟到未来的工作。它并不声称直接解决特定的工业加工问题，而是提供了研究这些系统底层物理机制的模拟基础设施。

类似论文