《Powder Technology》：Shape-dependent mixing kinetics of randomised particles in binary mixtures: An experimental study
编辑推荐：
本摘要报道了一项关于随机化不规则颗粒形状对二元混合物混合动力学影响的实验研究。研究人员采用球谐函数（spherical harmonics, SH）生成六种不规则颗粒几何形状，并系统研究了这六种几何形状的所有自混合与交叉混合组合。结果表明，凸度（convexi
本摘要报道了一项关于随机化不规则颗粒形状对二元混合物混合动力学影响的实验研究。研究人员采用球谐函数（spherical harmonics, SH）生成六种不规则颗粒几何形状，并系统研究了这六种几何形状的所有自混合与交叉混合组合。结果表明，凸度（convexity）与可见壁混合动力学（visible-wall mixing kinetics）表现出最强的关联性；曲率匹配（curvature matching）与拟合的混合速率参数相关联；而Corey形状因子（Corey shape factor）与混合行为的相关性有限。
## 论文解读：形状如何影响二元颗粒混合动力学？
### 1. 研究背景与问题
颗粒材料混合是众多工业过程中的核心单元操作，涵盖建筑中的混凝土搅拌、机械工程中的粉末冶金以及制药行业的粉末处理。高效混合对于确保产品均匀性、机械稳定性或活性成分的均匀分布至关重要。然而，尽管颗粒混合普遍存在，其性能预测仍具挑战性，原因在于颗粒流动的复杂动力学以及混合动力学对粒径、密度和形状等颗粒级特性的强敏感性。
经典的颗粒混合描述框架通常区分对流、扩散和剪切驱动机制。虽然颗粒形状影响所有这些机制，但其作用往往仅被隐式处理。多数研究优先关注粒径或密度差异，即便显式考虑形状，也多局限于球体、超二次曲面、棒状、立方体或椭球体等高理想化几何形状。因此，当前对混合动力学的理解主要局限于规则形态的子集。
随着三维成像和增材制造技术的普及，颗粒尺度几何不规则性对混合的影响重新受到关注。已有离散元模拟表明多面体颗粒的混合速度明显慢于球体；实验也显示不同三维打印形状的混合物会增强偏析。然而，也有研究认为重建砂粒的混合动力学差异可由球度单独捕捉。这些结果凸显了颗粒形状的重要性，同时也暴露出量化形状方面缺乏共识。特别是，现有研究未提供颗粒形态的连续描述。近年来已有方法能够连续变化地生成颗粒形状，同时引入几何不规则性，这为系统研究特定形状属性受控、非受控几何特征保持随机的混合实验提供了能力。
尽管已有众多研究考虑颗粒形状，其对颗粒混合影响的全面理解仍属空白。为填补这一知识缺口，研究人员探究了算法生成的不规则颗粒形状如何影响旋转鼓中的混合动力学。
### 2. 研究内容与结论
研究人员开发了一种自起始（ab initio）几何生成算法，能够生成具有任意受控形态属性的随机化三维颗粒几何族。这些形状随后以相近体积和材料密度进行三维打印，以减少这些属性在比较几何效应时的干扰。为在实验可控性与系统变化之间取得平衡，研究聚焦于六种定义形状的小规模二元交叉混合集合，仅采用一种填充水平和一种转速。由此能够将测量的混合轨迹与一组定量形状描述符相关联，识别在所探究参数范围内与混合速率和稳态组成最相关的颗粒几何因素。
研究得出以下核心结论：凸度与可见壁混合动力学表现出最强的关联性，凸度越高，混合越快且达到良好混合状态所需旋转次数越少；曲率匹配与拟合的混合速率参数相关联，但不影响总混合时间；Corey形状因子与宏观形状描述符在本研究窗口内与混合行为相关性有限。这些发现表明，颗粒形状可用于调控混合过程，例如通过提高颗粒凸度或改变表面曲率匹配来减少达到可比混合状态所需的旋转次数。
### 3. 关键技术方法
研究人员采用三阶段自起始几何生成方法：使用球谐函数（SH）表示三维颗粒表面，并通过控制缩放因子和可选择的百分位向量随机生成候选形状；随后以外层优化循环迭代调整形状描述符，使生成几何接近目标描述符。混合实验在定制旋转鼓中进行，鼓径与粒径比约为23，填充水平约40%，转速为30±1 rpm。红绿颗粒各210 g（约2400颗粒子），通过分层初始状态装载。图像采集使用全局快门相机，图像分割采用经人工标注数据训练的U-Net风格神经网络（MobileNetV2主干），随后进行闭运算与后处理。混合强度采用基于方差的Lacey混合指数（LMI），并利用加权像素频率计算。混合轨迹以双指数修正模型拟合，并通过高斯过程回归结合重复实验数据来估计参数及其不确定区间。共研究了六种形状的全部21种自混合和交叉混合组合，每种组合至少三次重复实验。
### 4. 研究结果
**4.1 轨迹概览与拟合质量**
通过对所有混合轨迹的拟合，模型能良好捕捉实验轨迹。参数`r0.8`（达到可见壁LMI=0.8所需旋转次数）范围从颗粒II的2.2转到颗粒VI的5.1转，差异显著。速率常数γ以颗粒VI最低（0.71 rot?1），最高出现在颗粒II与V的交叉混合中（γ=1.5 rot?1）。稳态区噪声随系统接近稳态而增大，这是红绿边界像素增多的结果。视觉检查验证了LMI值的合理性。总体上，颗粒I–IV的自混合与交叉混合动力学相似，而颗粒V和特别是VI混合更慢，且高度不规则颗粒需要更多旋转才能达到可比均匀度。
**4.2 回归模型关联**
研究人员对`r0.8`、γ、`M∞`和`M(0)`分别拟合普通最小二乘（OLS）回归模型。结果显示，形状描述符与`r0.8`、γ和`M∞`具有强的样本内一致性，而与`M(0)`的关联较弱。`M(0)`仅与实测体积`V_meas`相关，该体积作为制造偏差的控制变量，并非刻意改变参数；`M(0)`主要作为容许初始条件实验变化的干扰参数。
**4.3 各形状描述符的关联**
`M∞`（稳态LMI）与平均凸度`g_CO`、形状指数匹配`g_SI`和曲率匹配`g_CU`均呈显著正相关。`M∞`低于1是因为稳态时仍存在残余可见壁浓度方差。速率γ与`g_CO`正相关，与`g_CU`负相关，与`g_SI`无显著相关。`r0.8`与`g_CO`和`g_SI`负相关，与`g_CU`无明确关联。这说明`g_CU`与混合轨迹形状的变化相关，但不改变总混合时间。由于自混合中`g_CU`=1，观察到的关联意味着交叉混合具有更低的稳态值和更高的速率。
**4.4 标准化系数与机理讨论**
标准化系数显示，凸度是对混合动力学关联最强的描述符。这支持已有文献中近圆形光滑颗粒混合最快的结论。机制上，凸度可解释为颗粒迁移率的代理指标：不规则颗粒滚动和重定向阻力增加导致迁移率降低，从而减少对流输运、延长混合时间。形状指数匹配`g_SI`主要与`r0.8`相关而非γ，提示其可能涉及较慢的扩散混合组分；当颗粒表面特征分布互补时，微分散改善，可能缩短扩散混合尾段，减少达到LMI阈值所需旋转次数。这与多机制混合模型中扩散机制在对流层建立后成为限速步骤的概念一致。
**4.5 未观察到关联的描述符**
与常见文献报道不同，本研究未发现混合速度与平均Corey形状因子`g_CF`或`V_meas`之间存在显著相关。这可能是由于所探究的形状描述符窗口较窄，介观效应可能对混合行为有更强的关联。剩余未解释方差可能源于表面纹理差异，这与先前关于表面粗糙度效应的研究以及将堆积行为差异归因于宏观和介观尺度特征的研究一致。
### 5. 讨论与结论总结
讨论部分指出，研究发现可见壁混合行为与单颗粒和成对几何描述符均相关，其中凸度对混合速率和达到良好混合所需旋转次数关联最强，曲率匹配描述符则与混合轨迹的曲率相关，暗示局部颗粒重排的差异。通过结合自起始几何生成与受控旋转鼓实验，研究人员在所探究设计范围内建立了形状描述符与颗粒混合动力学之间的关联。
研究结论翻译如下：本研究表明，不规则颗粒材料的可见壁混合行为与单颗粒和成对几何描述符（包括凸度和基于曲率的匹配度量）相关。其中，凸度与混合速率及达到良好混合可见壁LMI所需旋转次数表现出最强关联。此外，基于曲率的匹配描述符与混合轨迹的曲率相关，这与局部颗粒重排的差异一致。通过将自起始几何生成与受控旋转鼓实验相结合，研究人员在所研究设计范围内建立了形状描述符与颗粒混合动力学间的关联。这些发现表明，颗粒形状原则上可用于调整混合过程，例如通过提高颗粒凸度或改变表面曲率匹配来减少达到可比可见壁混合状态所需的旋转次数。该应用设想需在更复杂的配方中直接验证。
