Percolation of a rod-like particle in a static bed of spheres: trapping and passing

本研究通过数值模拟证明，棒状颗粒在静态球体床层中的渗透受控于由棒长和孔隙几何结构决定的捕获与通过机制之间的转变，其中较短的棒由于受几何捕获的影响较小，其移动速度几乎是较长棒的两倍。

要点

想象一个由成千上万颗紧密排列的大而光滑的弹珠组成的巨大、隐形的筛子。现在，想象你向这个筛子中丢入一把不同的物体：其中一些是微小的弹珠，另一些是长而光滑的棍子（比如未煮熟的意大利面或牙签）。

这篇论文是一个计算机模拟，它观察了在重力的作用下，这些“棍子”试图穿过弹珠之间的缝隙时会发生什么。研究人员想要了解为什么有些物体能完全穿透，而另一些却会被卡住。

以下是他们发现的研究结果，通过简单的概念进行了拆解：

1. 两种结果：“通过”与“被困”

当棍子落下时，它们最终会进入两个阵营之一：

• 通过者： 这些棍子找到了路径，在缝隙中扭动，并以稳定的速度穿过整个弹珠层。
• 被困者： 这些棍子落下了一段时间，但最终被卡住了。它们停止移动，并留在弹珠堆中。

论文发现，一根棍子是被困还是通过，主要取决于它的长度与弹珠之间缝隙大小的比例。

2. “钥匙与锁”的问题

把弹珠之间的缝隙想象成微小且不规则的门口。

• 短棍子就像小钥匙。它们可以轻松地扭转和旋转，以适应几乎任何门口。它们掉落得很快，因为不会被卡住。
• 长棍子则像长而坚硬的管道。要穿过一个门口，管道必须保持完美笔直并与开口对齐。如果它侧向撞击门框，就会被卡住。由于堆积中的缝隙是随机且杂乱的，长棍子经常会以错误的夹角撞到“门框”并发生卡顿。

3. 形状的“速度限制”

研究人员发现了一个令人惊讶的规则：短棍子的下落速度几乎是长棍子的两倍。

为什么？

• 短棍子的行为几乎和那些大弹珠一样。它们可以轻易地翻滚并从孔洞中滑过，没有太多麻烦。
• 长棍子必须进行大量的“舞蹈”。在下落过程中，它们必须不断旋转，以寻找适合其长度的缝隙。这种不断的扭转和旋转减慢了它们的进度。这就像是在拥挤的房间里行走：一个小孩子可以轻松穿梭于人群，但一个拿着长梯的高个子则必须停下、转身，并等待一条清晰的路径，这显著减慢了他们的进度。

4. “卡住”的瞬间

当一根棍子最终被困住时，它并不是像汽车撞墙那样瞬间停止。它会在很短的距离内（大约半个大弹珠的宽度）减速，然后才完全固定下来。

论文还观察了它们是如何被卡住的：

• 短棍子通常是以垂直站立的状态被卡住，横在弹珠的侧面。
• 长棍子则会以各种奇怪的角度被卡住。它们经常通过同时接触三到四个弹珠而被卡住，从而形成一个固定住它们的复杂“结”。

5. “魔力数字”

研究人员发现了一个特定的“临界点”。如果一根棍子的长度超过了大弹珠宽度的约一半，它开始有很高的概率被困住。如果它比这个长度短，它几乎总能穿过去。

大局观

核心结论是，形状与大小同样重要。 在一个由圆形弹珠组成的世界里，大小是决定你是否能掉下去的唯一因素。但当你引入长而细的形状时，规则就改变了。变得“长”会让你变慢，并且更有可能被卡住，这并不是因为你变重了，而是因为你很难与堆积物中杂乱、随机的孔洞对齐。

这有助于解释为什么在自然界或工业领域中，长形物体（如纤维或谷物）在与圆形物体（如沙子或药丸）混合时，其行为方式会有所不同。

技术摘要

技术摘要：棒状颗粒在静态球体床层中的渗透行为

问题陈述
细颗粒通过静态大颗粒床层的渗透是颗粒物分选的一个基本机制，与地球物理现象（如滑坡、雪崩）以及工业过程（如药物混合、食品加工）密切相关。虽然球形颗粒的渗透已被广泛研究，但现实世界的颗粒材料通常含有非球形颗粒，如棒状、纤维或片状颗粒。这些各向异性的形状引入了球形系统中不存在的、依赖于方向的约束、相互锁合以及独特的接触力学特性。本文旨在解决理解棒状颗粒如何渗透通过无序静态大球体床层的知识空白，特别研究了连续渗透与捕获之间的转变，以及几何各向异性如何影响渗透速度和动力学过程。

研究方法
本研究采用基于 MercuryDPM 代码的离散元法 (DEM) 模拟。系统由 3,000 个无摩擦球体（直径为 $d_L$）组成的静态无序床层组成，通过各向同性压缩和减压制备，以确保处于随机堆积状态。棒状颗粒被建模为由一系列共线的、无摩擦的“粘合”球体（直径为 $d_S$）组成的链，具有不同的长径比 ($A = l/d_S$) 和尺寸比 ($R = d_L/d_S$)。

关键模拟参数包括：

• 棒状几何形状： 棒状颗粒由组成球体的数量 ($n$)、长径比 ($A$) 和无量纲长度 $L^* = l/d_L$ 定义。研究探讨了 $1 \le A \le 10$ 且 $0.1 \le L^* \le 1$ 的情况，涵盖了多种 $R$ 值（包括临界几何捕获比 $R_t \approx 6.464$）。
• 相互作用模型： 法向力遵循线性弹簧-阻尼器模型。切向力设为零，以隔离几何效应。禁用棒-棒相互作用；棒仅与静态床层球体发生相互作用。
• 实验方案： 将 1,000 到 5,000 个互不干扰的棒状颗粒在重力作用下投入床层中。通过追踪轨迹来确定渗透深度、速度、取向和接触统计数据。根据阻尼参数计算内生恢复系数，使得球-球接触时的 $e \le 0.6$。

核心贡献与结果

1. 识别出两个截然不同的机制：
   研究确定了通过机制（棒状颗粒以恒定平均速度连续渗透）和捕获机制（棒状颗粒在渗透有限距离后停止运动）。这两个机制之间的转换受棒长度 ($L^*$) 和尺寸比 ($R$) 的控制。

   • 一个临界无量几种长度 $L_c^* \approx 0.52$（对于 $R=7$）标志着转换点。在该长度以下，棒状颗粒可以无限期渗透；在该长度以上，捕获概率随深度增加。
   • 被捕获棒状颗粒的特征渗透长度 $\lambda$ 遵循幂律发散，当 $L^*$ 从上方趋近于 $L_c^*$ 时：$\lambda/d_L \propto (L^* - L_c^*)^{-1}$。

2. 捕获机制与几何形状：

   • 接触数： 被捕获的棒状颗粒以平均接触数 $\langle Z \rangle = 5$ 被固定。由于在缺乏摩擦的情况下，棒状颗粒仍保留绕其长轴旋转的自由度，该数值低于球体的等静力极限 ($\langle Z \rangle = 6$)。捕获是由几何纠缠而非单纯的力平衡驱动的。
   • 取向： 短棒（$L^* \approx 0.57$）主要以垂直取向被阻截，其接触点位于床层球体的赤道附近。较长的棒状颗粒（$L^* = 1$）表现出更宽的取向分布和接触位置分布（从极点到赤道均有分布）。
   • 接触位置： 对于短棒，接触集中在北半球（$20^\circ \le \theta \le 60^\circ$）。随着棒长度的增加，接触点向赤道和南半球扩散，反映了增强的几何挫折感。

3. 渗透速度标度：

   • 无量纲渗透速度 $\tilde{v}_p = v_p / \sqrt{gd_L}$ 随着棒长度 ($L^*$) 和长径比 ($A$) 的增加而降低。由于几何约束较少，短棒的渗透速度几乎是长棒的两倍。
   • 当使用 $\sqrt{gd_L(R - R_p)}$ 进行缩放时（其中 $R_p \approx 4.5$ 代表特征孔径直径比），所有棒状几何形状（包括单球极限）的数据都坍缩到同一条曲线上。该标度统一了球体和棒状颗粒的动力学过程。
   • 速度对取向具有强烈的依赖性：垂直对齐的棒状颗粒比水平对齐的具有显著更高的速度。

4. 动能分布：
   通过对动能比 ($K_{rot}^{rms}/K_{tra}^{rms}$) 的分析发现，对于短棒，旋转能和平动能近似实现能量均分。随着棒长度的增加，由于孔隙网络限制了旋转自由度，旋转激发相对于平动受到抑制。

意义
本文证明了颗粒形状的各向异性通过引入不同于球形系统的动力学约束和阈值，从根本上改变了颗粒渗透行为。研究结果表明：

• 棒状颗粒的行为并不等同于球体；它们通过孔隙的能力严格取决于取向。
• 存在一个导致捕获的临界几何阈值，该阈值取决于棒长度和床层孔隙结构。
• 棒状颗粒的渗透速度可以通过基于球体渗透模型的标度律进行预测，前提是必须考虑几何约束。

这些结果为预测非球形颗粒在颗粒混合物中的输运和分选提供了定量框架，对于理解含有纤维或碎屑的自然系统以及优化涉及药丸、谷物或纤维添加剂的工业过程具有重要意义。这项工作为将颗粒分选模型扩展到更广泛的各向异性颗粒形状奠定了基础。