🔬 condensed matter

Energy non-equipartition in vibrofluidized particles

本研究利用三维离散元模拟表明，在振动流化颗粒系统中，真实的切向弹簧刚度参数会导致平动与转动模式间的能量均分失效，尤其对于高摩擦系数的颗粒而言。

原作者： Alok Tiwari, Manaswita Bose, V. Kumaran

发布于 2026-05-18

📖 1 分钟阅读☕ 轻松阅读

原作者： Alok Tiwari, Manaswita Bose, V. Kumaran

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

想象一个装满成千上万颗微小弹球的盒子。现在，想象你非常快速地上下摇晃这个盒子的底部。弹球四处跳跃，相互碰撞，并从盒壁上弹开。在物理学中，这被称为“振动流化”系统。

通常，科学家预期如果你足够长时间地摇晃这些弹球，能量将会均匀分布。这被称为能量均分定理。可以把它想象成一群朋友分享一个披萨：最终，每个人都会得到相等的一块。在这种情况下，“披萨”就是运动能量。其中一部分能量使弹球向前、向后、向上和向下移动（称为平动能量），另一部分则使它们旋转（称为转动能量）。

标准规则指出：如果弹球完全相同且系统稳定，那么用于移动的能量应等于用于旋转的能量。

核心问题
本文的作者问道：“如果弹球并非完全光滑会怎样？如果它们像砂纸一样粗糙呢？”他们想要探究“粗糙度”（摩擦）和弹球的“刚度”是否会打破能量均分的规则。

为了找出答案，他们并没有使用真实的弹球盒。相反，他们使用名为 LAMMPS 的程序在计算机内部构建了一个巨大的虚拟三维世界。他们模拟了数百万次球形粒子之间的碰撞。

两个关键要素
研究人员在模拟中改变了两个主要因素：

摩擦（粗糙度）： 他们测试了几乎完全光滑（像冰一样）到极其粗糙（像砂纸一样）的弹球。
刚度比（ $\kappa$ ）： 这有点像弹球的“弹性”。当两个弹球相撞时，它们会轻微变形。研究人员观察了弹球在受到侧面撞击（切向）时的刚度与受到正面撞击（法向）时的刚度之比。他们测试了两种特定的刚度“配方”：
- 配方 A： 特定的 2/7 比例。
- 配方 B： 特定的 3/4 比例。

他们的发现

1. 光滑情况（“冰面”场景）
当弹球非常光滑且摩擦系数很低时，能量确实保持了分离状态。弹球四处移动，但几乎不旋转。这就像一群人在光滑的冰场上：他们可以轻易滑行，但无法获得足够的抓地力来旋转。移动能量和旋转能量截然不同。

2. 粗糙情况配合配方 A（“砂纸”成功）
当他们使用2/7 刚度比并使弹球非常粗糙（高摩擦）时，神奇的事情发生了。弹球终于实现了能量的均等分配！

类比： 想象弹球穿着厚重的橡胶靴。当它们相互碰撞时，摩擦力如此之大，以至于它们会“粘”在一起一瞬间。这种粘附使得向前的运动能够轻易转化为旋转。
结果： 移动能量和旋转能量变得几乎相等。“披萨”终于被公平地分享了。

3. 粗糙情况配合配方 B（“失效”的规则）
这里情况变得奇怪了。当他们使用3/4 刚度比并使弹球非常粗糙时，能量并没有均等分配，尽管弹球是粗糙的。

类比： 想象弹球穿着既粘又带有奇怪弹簧机制的靴子。当它们碰撞时，它们不仅仅是粘在一起；它们会以一种浪费能量的方式弹跳和滑动。
结果： 弹球保持着带有大量能量的向前运动，但拒绝将旋转提升到相同的水平。移动能量与旋转能量的比率稳定在约5 比 1。“披萨”没有被分享；一群朋友得到了五块，而另一群只得到了一块。

为什么会发生这种情况？
作者仔细研究了碰撞过程，以理解为何配方 B 未能实现能量共享。

在成功的案例（配方 A）中，弹球主要进入“粘附”模式，即它们短暂地相互保持接触。这是一种将移动能量高效转化为旋转能量的方式。
在失败的案例（配方 B）中，弹球主要进入“滑动”模式。它们相互摩擦但没有粘附。这产生了大量的摩擦热（能量损失），却未能有效地将向前运动转化为旋转。能量在能够被共享之前就被摩擦消耗掉了。

结论
该论文得出结论：粒子的刚度（特别是侧向刚度与正向刚度的比率）与其粗糙度同样重要。

你不能假设粗糙的粒子会自动均分它们的能量。根据材料具体的“弹性”不同，粗糙粒子要么完美地共享能量（如 2/7 的情况），要么囤积能量，导致旋转运动微弱（如 3/4 的情况）。

简而言之：仅凭粗糙度并不能保证公平。材料内部的“弹簧”特性同样重要。

技术摘要：振动流化颗粒中的能量非均分

问题陈述
能量均分定理认为动能在流体的所有自由度中均匀分布，但在颗粒系统中经常偏离。尽管先前的研究已记录了涉及光滑颗粒（导致涨落动能的各向异性）和粗糙颗粒的振动流化床中的非均分现象，但真实颗粒相互作用参数的具体影响仍是研究课题。具体而言，本工作探讨了切向与法向弹簧刚度系数之比（ $\kappa = k_t/k_n$ ）如何影响平动模式与转动模式之间的能量分配。本研究旨在确定真实的颗粒接触模型（特别是涉及高摩擦和特定刚度比的那些模型）是否会导致能量均分，还是导致持续的非均分。

方法论
作者利用开源软件 LAMMPS 采用离散元法（DEM）模拟了一个由球形颗粒组成的三维垂直振动颗粒系统。

相互作用模型：接触力采用线性弹簧 - 阻尼器（LSD）模型进行建模。法向弹簧刚度常数（ $k_n = 10^8 mg/d$ ）被选为足够大，以确保碰撞主要为二体碰撞（ $t_c/t_f \ll 1$ ）。
参数：研究将摩擦系数（ $\mu$ ）从近乎光滑（ $\mu \approx 0$ ）变化至近乎完全粗糙（ $\mu = 10$ ）。主要研究了两个刚度比 $\kappa$ 值： $2/7$ 和 $3/4$ 。此外，还针对 $0.67 \le \kappa < 1$ 进行了额外模拟。
模拟协议：模拟分两个阶段运行：首先是初始弛豫阶段（ $7 \times 10^7 t_c$ ）以达到稳态，随后使用较小的时间步长（ $\Delta t = t_c/100$ ）进行数据收集。
分析：分析瞬时线速度和角速度，以确定平均涨落平动动能（ $T$ ）和转动动能（ $T_R$ ）的分布。系综平均值在超过 $10^3$ 个构型上计算。比率 $K = K_{ET}/K_{ER}$ （平动动能与转动动能之比）被用作能量均分的主要度量指标。

主要结果

对刚度比（ $\kappa$ ）的依赖性：能量比率 $K$ 的行为 critically 依赖于 $\kappa$ 。
- 对于 $\kappa = 2/7$ ：随着摩擦系数 $\mu$ 的增加，比率 $K$ 单调递减。在高摩擦下（ $\mu = 10$ ）， $K$ 趋近于 1，表明平动模式与转动模式之间实现了能量均分。
- 对于 $\kappa = 3/4$ ：行为是非单调的。在随 $\mu$ 增加而初始减小后， $K$ 上升并稳定在约 5 的水平。这表明即使在非常高的摩擦系数下，能量均分的缺失依然持续存在。
- 针对 $0.67 \le \kappa < 1$ 的模拟证实，在此范围内具有高摩擦系数的颗粒会持续存在非均分现象。
碰撞动力学与能量交换：
- 对单次碰撞的分析表明，对于 $\kappa = 2/7$ ，处于“粘着”机制（其中转动恢复系数 $\beta \approx 1$ ）的接触比例随 $\mu$ 增加而增加。在 $\mu = 1$ 时，近 75% 的接触处于这种能量守恒机制中，从而促进了能量均分。
- 对于 $\kappa = 3/4$ ，处于粘 - 滑机制的接触比例在 $\mu > 1$ 时稳定在约 55%，且 $\beta < 1$ 。因此，能量交换比率 $\Theta = E_{exch}/\Delta E$ 保持小于 1，并在非常粗糙的颗粒情况下进一步减小。耗散超过了模式间的能量交换，从而阻止了能量均分。
理论验证：平衡颗粒温度（ $T_o$ ）的模拟结果与理论预测（基于主导阶矩展开）表现出极好的一致性，适用于完全光滑颗粒以及稀薄状态下的粗糙颗粒，前提是使用了适当的缩放参数（涉及 $\beta$ ）。

意义与主张
本文主张，刚度比 $\kappa$ 的选择是预测真实颗粒宏观流动行为的关键因素。研究表明，在振动流化粗糙颗粒系统中，能量均分并非必然结果；相反，它取决于摩擦与切向 - 法向刚度比之间的具体相互作用。

具体而言，作者得出结论：

在 $\kappa = 2/7$ 且高摩擦条件下观察到能量均分，这是由于能量守恒的粘着碰撞占主导地位。
在 $0.67 \le \kappa < 1$ （包括 $\kappa = 3/4$ ）条件下观察到非均分，因为这些机制中的粘 - 滑碰撞不守恒能量，导致平动动能与转动动能之间存在显著差异。

这些发现表明，真实的颗粒流变学建模需要仔细考虑刚度比，因为它从根本上改变了系统内的能量分布机制。该工作未提出新的实验装置或应用，而是系统分析了微观相互作用参数如何决定宏观能量分布。

类似论文