Mpemba effect in a sheared granular gas with velocity-dependent restitution

利用动理学理论，本研究证明，具有速度依赖恢复系数的稀薄剪切颗粒气体同时表现出温度和粘度姆佩姆巴效应，即初始温度较高的系统比温度较低的系统弛豫得更快，而速度依赖性引入的内在时间尺度使得多条弛豫曲线能够发生交叉。

要点

以下是用通俗语言和日常类比对这篇论文的解读。

核心理念：弹跳球盒子里的“热水”之谜

你可能听说过姆潘巴效应。这是一种反直觉的现象：热水有时比冷水结冰更快。这听起来不可能，但它之所以发生，是因为“热”水具有不同的内部结构或历史，一旦关上冷冻室的门，这种特性反而有助于它更快地冷却。

这篇论文探讨了同样的奇怪现象是否也会发生在颗粒气体中。想象一个装满成千上万颗微小、坚硬钢球的盒子，它们在四处弹跳。与真实的气体分子不同，这些钢球每次相互碰撞时都会损失能量（它们不会完美地弹跳）。为了让它们保持运动，科学家会对盒子施加“剪切”，这意味着他们向右滑动盒子的顶部，向左滑动盒子的底部，像搅拌机一样不断搅动这些钢球。

研究人员提出了这样一个问题：如果你有两个装有这些弹跳球的盒子，其中一个比另一个“更热”（运动速度更快），那么较热的那个盒子是否真的能比较冷的那个更快地进入平稳、稳定的节奏？

两个起点

为了验证这一点，他们设置了两种不同的场景（协议），这两种场景最终都达到了完全相同的“最终状态”（特定的搅拌速度）：

1. “已搅拌”起点（FS 协议）： 想象一个已经搅拌了很长时间的球盒。它们正以一种特定、有序但混乱的模式运动。然后，搅拌速度突然改变。
2. “静止”起点（FI 协议）： 想象一个刚刚静止不动（或正在自行冷却）且未受搅拌的球盒。在完全相同的时刻，搅拌以与第一个盒子相同的新速度启动。关键在于，这个盒子里的球初始具有更高的温度（运动速度更快），而第一个盒子里的球速度较慢。

结果：热者胜出

在正常世界中，你可能会预期较冷的盒子更快达到最终的稳定状态。但是，就像热水结冰的把戏一样，较热的盒子（“静止”起点）追上了并超过了较冷的盒子。

• 为什么？ “已搅拌”的盒子由于之前的搅拌，积累了大量的内部应力和“坏习惯”。当速度改变时，它必须解开这些旧模式，这拖慢了它的速度。
• “静止”的盒子虽然更热，但它是从一张白纸开始的（没有内部应力）。它能够更有效地吸收新的搅拌运动，并更快地进入节奏，尽管它开始时拥有更多的能量。

这就是温度姆潘巴效应：拥有更多能量的系统反而更快地弛豫到了稳定状态。

转折：“粘度”把戏

论文还发现了一个更奇怪的现象。不仅仅是球的温度（速度）表现出这种效应，粘度（气体对搅拌感觉有多“稠”或阻力多大）也是如此。

通常，当你改变搅拌流体的速度时，其厚度会平滑地变化。但在这里，研究人员观察到粘度曲线多次相互交叉。较热的系统不仅仅是一次性超过了较冷的系统；它之字形地超过它，然后可能又落后，接着再次超过，最后才安定下来。

秘密成分：“弹性”开关

为什么会发生这种情况？关键在于他们对钢球应用的一条特殊规则：弹性取决于它们撞击的猛烈程度。

• 轻柔撞击： 钢球非常有弹性（像超级球）。
• 猛烈撞击： 钢球弹性较差（像一团粘土）。

这在物理上创造了一个“开关”。因为钢球在不同速度下表现不同，这给系统引入了第二个时钟或时间尺度。

这就像一辆拥有两种不同挡位的汽车。如果你只有一档，汽车会平稳加速。但如果有一辆汽车会根据行驶速度突然换挡，加速过程就会变得生涩且复杂。正是钢球物理特性中的这种“换挡”，导致弛豫曲线多次交叉，从而产生了多重姆潘巴效应。

结论

这篇论文证明，在一个“弹性”取决于速度的弹跳球气体中：

1. 较热的系统可以比冷却的系统更快地弛豫到稳定状态（温度姆潘巴效应）。
2. 气体的“厚度”也可以表现出这种效应（粘度姆潘巴效应）。
3. 由于速度依赖性的弹性，这些系统在走向稳定的过程中可以多次交叉，这是更简单的模型中未曾见过的行为。

这是一个纯粹关于颗粒材料中能量与应力如何相互作用的数学和物理发现，表明“更热”并不总是意味着“更慢安定”。

技术摘要

技术摘要：具有速度依赖恢复系数的剪切颗粒气体中的姆潘巴效应

问题陈述
姆潘巴效应最初是在水结冰过程中被观察到的，它描述了一种反直觉的现象：初始温度较高的系统比初始温度较低的相同系统更快地达到平衡。尽管该效应已被公认为驱动经典和量子系统中非平衡弛豫的普遍特征，但流变学效应和速度依赖耗散在颗粒气体中产生这种行为的具体作用仍 largely 未被探索。具体而言，速度依赖的恢复系数（它引入了固有的速度尺度和额外的弛豫时间）如何影响剪切颗粒气体的弛豫路径，尚不明确。

方法论
作者研究了一种由具有速度依赖恢复系数的无摩擦硬球组成的稀薄剪切颗粒气体。该研究采用两种主要方法：

1. 动理学理论：基于 Grad 矩方法，作者推导出一组封闭的耦合演化方程，用于描述颗粒温度（$T$）、温度各向异性（$\Delta T$）和剪切应力（$P_{xy}$）。该模型引入了阶梯式速度依赖恢复系数 $e(v_n)$，其中包含一个特征阈值速度 $v_c$。这使得能够解析计算碰撞能量耗散率（$\zeta$）和有效碰撞频率（$\nu$）。
2. 直接模拟蒙特卡洛（DSMC）：为了验证理论预测，作者在均匀简单剪切流下，利用 SLLOD 动力学和 Lees–Edwards 边界条件，对 $N=10^3$ 个粒子进行了 DSMC 模拟。

该研究比较了在剪切率突然改变至目标值 $\dot{\gamma}_{tar}$ 后，两种不同的弛豫协议：

• 从剪切态（FS）协议：系统从具有初始剪切率 $\dot{\gamma}_{ini}$ 的稳态剪切状态开始。
• 从各向同性态（FI）协议：系统从零剪切应力的未剪切、各向同性状态开始。

主要贡献与结果
本文证明了该系统中存在温度和粘度姆潘巴效应，其特征是弛豫曲线的交叉。

• 温度姆潘巴效应：作者表明，尽管 FI 系统的初始温度较高（$T^{(FI)}_{ini} \gtrsim T^{(FS)}_{ini}$），但它能比 FS 系统更快地弛豫到最终的稳态温度。这种机制源于能量注入与耗散之间的初始平衡：FS 系统由于负剪切应力而最初经历受抑的弛豫，而 FI 系统从零剪切应力开始，经历快速的初始耗散。
• 粘度姆潘巴效应：研究识别出一种“粘度姆潘巴效应”，即两种协议的剪切粘度弛豫曲线发生交叉。这在 FI 协议中尤为显著，因为粘度从零开始演化，产生强烈的瞬态响应，从而促进了交叉。
• 多重姆潘巴效应：一个重要的发现是弛豫曲线中出现多重交叉。与通常仅表现出单次交叉的具有恒定恢复系数的系统不同，速度依赖的恢复系数引入了一个额外的固有时间尺度。这导致了非单调的弛豫动力学，使得轨迹可以交叉多次。相图揭示了特定的参数区域（取决于初始剪切率和温度），在该区域内效应会发生两次。
• 流变学联系：多重姆潘巴效应的出现与稳态流变学表现出对剪切率非单调（S 形）依赖的参数区域相关，具体而言，当恢复系数参数满足 $e_1 < e_2$ 时。

意义与主张
本文声称提供了一种最小化的、解析可处理的动理学机制，用于解释驱动颗粒气体中的多重姆潘巴效应。作者强调，这些发现建立了反常弛豫与非平衡流变学之间的直接联系。至关重要的是，该机制纯粹是动理学的；它不依赖于致密堆积、摩擦接触或阻塞效应。该研究表明，通过调节微观耗散机制，特别是通过恢复系数的速度依赖性，可以在颗粒流中实现受控的姆潘巴现象。与以往关于稳态流变学的研究不同，本研究聚焦于瞬态弛豫动力学，提供了一个框架来理解内部自由度（应力和各向异性）如何与速度依赖耗散耦合，从而支配弛豫路径。