[论文解读] The S shape of a granular pile in a rotating drum

Nicolas Taberlet, Patrick Richard+1 more|arXiv (Cornell University)|Jul 4, 2006

Granular flow and fluidized beds被引用 13

一句话总结

本文表明，旋转颗粒鼓中S形自由表面主要由摩擦性端板引起，而非离心力或颗粒惯性。通过离散元方法（DEM）模拟和理论模型，作者推导出一个无量纲数Λ，用于量化端板的影响，表明在长鼓或无摩擦端板中S形消失，且数值验证了标度律。

ABSTRACT

The shape of a granular pile in a rotating drum is investigated. Using Discrete Elements Method (DEM) simulations we show that the "S shape" obtained for high rotation speed can be accounted for by the friction on the end plates. A theoretical model which accounts for the effect of the end plates is presented and the equation of the shape of the free surface is derived. The model reveals a dimensionless number which quantifies the influence of the end plates on the shape of the pile. Finally, the scaling laws of the system are discussed and numerical results support our conclusions.

研究动机与目标

• 探究旋转颗粒鼓中S形自由表面的成因，尽管已有理论提出离心力或颗粒惯性作为解释，但该现象仍未被完全阐明。
• 检验此前被忽视的假设：端板摩擦是导致S形（尤其在短鼓中）的主要原因。
• 建立理论模型，利用无量纲数Λ量化端板摩擦对颗粒堆形状的影响。
• 推导并验证描述系统行为的标度律，预测鼓的几何形状、转速和颗粒尺寸对S形的影响。

提出的方法

• 采用软球分子动力学方法的离散元方法（DEM）模拟，研究端板摩擦特性可控的旋转鼓中颗粒流动。
• 采用弹簧-阻尼器模型和正则化库仑摩擦定律，分别控制颗粒间及颗粒与壁面的相互作用，实现对颗粒和端板摩擦的独立调控。
• 推导自由表面坡度的理论方程：∂y/∂x = μ_i + μ_w Λ √(1 - x² - y²)，其中Λ = (d Ω²/g)^{1/4} (R/L)。
• 使用无量纲变量（x̃ = x/R, ỹ = y_surf/R）对模拟数据进行归一化，验证不同R、L和Ω下的标度律。
• 系统性地改变鼓长L和端板摩擦，以分离其对S形的影响，比较有摩擦与无摩擦端板的差异。
• 通过保持Λ恒定对数据进行重标度，实现数据点汇聚于单一曲线，实现标度律的数值验证。

实验结果

研究问题

• RQ1 当弗劳德数过低、离心力不足以解释时，旋转颗粒鼓中S形自由表面的成因是什么？
• RQ2 在短旋转鼓中，摩擦性端板在多大程度上影响颗粒堆自由表面的曲率？
• RQ3 能否通过包含端板摩擦的理论模型再现观测到的S形，并预测其对系统参数的依赖性？
• RQ4 控制颗粒堆形状的无量纲数是什么？鼓的尺寸、转速和颗粒尺寸的标度律如何与该数相关？

主要发现

• S形在旋转鼓中主要由端板摩擦引起，而非离心力或颗粒惯性，该结论经不同端板摩擦条件下的DEM模拟验证。
• 当移除端板摩擦时，S形消失，即使在高转速和低弗劳德数（Fr ≈ 10⁻⁴）下，自由表面也趋于平坦。
• 增加鼓长L可减弱端板影响：当L = 300d时，系统趋于平坦自由表面；而当L = 10d时，S形明显出现。
• 无量纲数Λ = (d Ω²/g)^{1/4} (R/L) 完全表征了系统的形状，Λ决定了相对于斜率为μ_i的平坦表面的偏离程度。
• 标度律经数值验证：Ω ∝ 1/R²，Ω ∝ L²，Ω ∝ 1/√d；当Λ保持恒定时，数据可汇聚于单一曲线。
• 该模型解释了Felix（2002）的实验观测，特别是转速与鼓半径和颗粒尺寸的标度关系。