[论文解读] Velocity Distributions in Homogeneously Cooling and Heated Granular Fluids

Twan van Noije, M. H. Ernst|arXiv (Cornell University)|Mar 3, 1998

Granular flow and fluidized beds参考文献 3被引用 25

一句话总结

该论文利用Enskog-Boltzmann方程推导非弹性颗粒流体中的速度分布，重点研究四阶累积量以量化在均匀冷却和均匀加热稳态下非高斯行为。结果表明，在冷却状态下，高斯分布的偏离较小，但在加热系统中，高能尾部呈现 $\exp(-A c^{3/2})$ 的标度形式，其中 $A \sim 1/\sqrt{\epsilon}$，$\epsilon$ 为非弹性度。

ABSTRACT

We study the single particle velocity distribution for a granular fluid of inelastic hard spheres or disks, using the Enskog-Boltzmann equation, both for the homogeneous cooling of a freely evolving system and for the stationary state of a uniformly heated system, and explicitly calculate the fourth cumulant of the distribution. For the undriven case, our result agrees well with computer simulations of Brey et al. \cite{brey}. Corrections due to non-Gaussian behavior on cooling rate and stationary temperature are found to be small at all inelasticities. The velocity distribution in the uniformly heated steady state exhibits a high energy tail $\sim \exp(-A c^{3/2})$, where $c$ is the velocity scaled by the thermal velocity and $A\sim 1/\sqrt{\eps}$ with $\eps$ the inelasticity.

研究动机与目标

• 理解由于非弹性碰撞导致的颗粒流体中非高斯速度分布，挑战普遍假设的麦克斯韦分布。
• 计算自由冷却和均匀加热颗粒系统中速度分布的四阶累积量。
• 评估非高斯行为对非弹性颗粒气体冷却速率和稳态温度的影响。
• 使用非微扰方法推导均匀加热稳态下速度分布的高能尾部。
• 纠正早期四阶矩计算中的错误，并为基于矩的颗粒气体分析提供一致的框架。

提出的方法

• 使用 $d$ 维非弹性硬球或硬盘的Enskog-Boltzmann方程来模拟单粒子速度分布。
• 应用标度假设，即在冷却情况下，分布仅通过衰减的温度 $T(t)$ 依赖于时间。
• 通过切比雪夫多项式展开分布函数以计算矩，重点关注四阶累积量 $a_2$。
• 通过变换至质心与相对速度坐标系，计算角向积分 $\beta_n$，并推导矩 $\mu_p$ 的方程。
• 使用高斯加权平均 $\langle \cdot \rangle_0$ 计算涉及 $c_{12}$ 和 $C$ 速度的积分。
• 通过将矩表达式代入矩层次方程，求解冷却和加热情况下的 $a_2$。

实验结果

研究问题

• RQ1 非高斯行为（以四阶累积量量化）如何影响自由演化颗粒流体的冷却速率？
• RQ2 均匀加热颗粒流体中速度分布的高能尾部形式为何？
• RQ3 非高斯性引起的修正如何影响非弹性颗粒系统中的稳态温度和冷却速率？
• RQ4 高能尾部对非弹性度参数 $\epsilon = 1 - \alpha^2$ 的依赖关系如何？
• RQ5 在 $\mathcal{O}(\epsilon)$ 范围内，非微扰方法与微扰结果相比如何？

主要发现

• 在均匀冷却状态下，四阶累积量 $a_2$ 通过非微扰方法得到修正，解决了Goldshtein与Shapiro早期工作中存在的代数错误。
• 在所有非弹性度下，由于非高斯行为导致的冷却速率和稳态温度的修正均较小。
• 在均匀加热稳态下，速度分布表现出高能尾部，其标度形式为 $\exp(-A c^{3/2})$，其中 $A \sim 1/\sqrt{\epsilon}$，$\epsilon$ 为非弹性度。
• 尾部衰减系数 $A$ 与非弹性度的平方根成反比，表明在恢复系数较低时偏离高斯分布更显著。
• 加热状态下四阶累积量由一个包含 $\beta_3$、$\langle c_{12}^3 \rangle_0$ 以及依赖于维度和 $\alpha$ 的 $T_1$、$T_2$ 项的闭合方程导出。
• 对冷却速率 $\mathcal{O}(\epsilon)$ 修正结果与Sela和Goldhirsch的微扰计算高度一致，验证了非微扰方法的有效性。