导流管流化床包衣机中颗粒循环与高分位接触载荷的CFD-DEM评估：以白菜种子为对象的干态筛选研究

《Agriculture》：CFD-DEM Evaluation of Particle Circulation and High-Percentile Contact Loading in a Draft-Tube Fluidized-Bed Seed Coater for Chinese Cabbage Seeds

　　

稳定循环与有限机械载荷是小型蔬菜种子流化床包衣的关键要求。对于白菜种子而言，其质量小、几何形状不规则且机械敏感性高，使得仅依靠平均颗粒速度和床层膨胀高度等全局指标难以评估操作条件。本研究针对导流管流化床种子包衣机建立了双向耦合的计算流体力学-离散元法（CFD-DEM）模型。白菜种子采用七球团簇（seven-sphere clumps）表示，并采用L9(33)正交表作为筛选设计，考察入口空气速度、初始床层高度和底部循环入口间隙三个因素。评估综合了循环时间分布（CTD）、全局低速颗粒分数Rs、95分位法向接触力F95,n以及归一化高分位接触载荷比ηc,95。由于L9正交表无法解析交互作用或支持完整的二次模型，因子效应仅在所研究的范围内进行描述性解释，而非作为验证性的全局优化。初始床层高度对循环周期、Rs和F95,n产生了最大的描述性贡献。增大入口空气速度缩短了循环周期并降低了Rs，但增加了高分位接触载荷。在九个工况中，F95,n的范围为7.50至14.50 mN，ηc,95的范围为0.18%至0.35%；ηc,95仅作为归一化载荷比使用，而非经过验证的种子损伤概率。在等权重和接触载荷优先权重下，工况4排名第一，而在循环优先权重下工况7排名第一。因此，工况4被描述为在所测试参数范围内依赖权重的均衡候选方案。原型实验复现了三种干燥操作条件下循环周期的排序，支持了模拟与观测循环行为之间的定性一致性。可见破损率保持在0.5%以下，但该观察仅提供了与模拟接触载荷趋势的初步定性对应关系。所提出的框架旨在用于干态阶段筛选，不直接预测湿法包衣质量、粘附性、团聚或发芽性能。

论文解读

一、研究背景与意义

白菜作为重要的叶菜类蔬菜作物，其种子质量直接影响机械化播种、幼苗整齐度和田间成苗率。种子包衣技术被广泛用于改善种子处理、播种适应性、种子保护及功能性材料的递送。在机械处理和包衣过程中，种子损伤可能源于准静态压缩、短时冲击、重复加载、切向摩擦及这些机制的组合作用。损伤敏感性取决于品种、含水量、加载方向、加载速率、接触材料及重复接触次数。因此，从一种压缩构型获得的力值不能自动解释为另一种种子批次或操作过程的通用动态破损阈值。

流化床包衣是颗粒包衣的重要技术，Wurster型和导流管流化床提供了重复的颗粒循环和强烈的气固相互作用。在导流管流化床中，颗粒通过中央导流管被向上提升，在上部区域分散，沿环形下降区下降，再通过底部循环入口重新进入，形成中央向上输送、上部扩散、环形下降和底部再入的循环路径。然而，传统流化指标主要描述宏观床层行为，床层膨胀高度和平均颗粒速度虽有助于识别整体流化状态，但无法显示单个颗粒是否完成了有效的循环周期、有多少颗粒处于低速状态，以及颗粒接触载荷的上尾分布如何变化。因此，需要多指标联合评估框架。

已有CFD-DEM研究多集中于药丸或相对规则的模型颗粒，对于形状不规则且机械敏感的蔬菜种子，将循环描述符与谨慎解释的接触载荷指标相结合，用于早期干态条件筛选，构成了本研究的增量贡献。

二、研究内容概述

研究人员针对导流管流化床种子包衣机建立了双向耦合的CFD-DEM模型，以商业“优绿3号”白菜种子为代表材料，采用七球多球团簇表示种子几何形状。研究考察了入口空气速度（A）、初始床层高度（B）和底部循环入口间隙（C）三个操作因素，采用L9(3³)正交表进行筛选设计。评估框架联合比较了循环时间分布（CTD）、全局低速颗粒分数R_s、95分位法向接触力F_95,n和归一化高分位接触载荷比η_c,95。此外，研究还制作了原型机，对三种代表性干燥工况进行了循环行为对比实验验证。

三、关键技术方法

研究采用的主要技术方法包括：双向耦合的CFD-DEM数值模拟方法，其中气相采用三维不可压缩连续相模型，颗粒相采用经典DEM理论逐团簇解析；Hertz-Mindlin无滑移接触模型用于接触力计算；Gidaspow曳力模型用于气固耦合；ANSYS Fluent 2024 R1与Altair EDEM 2024.1作为计算平台；L9(3³)正交试验设计用于因子筛选；Wen-Yu最小流化速度关联式作为选速范围的量级参考；多指标范围归一化评分与权重敏感性分析用于综合评估；原型实验采用高速视频追踪示踪种子进行循环周期测量。

四、研究结果

循环时间分布特征： 通过CTD曲线和统计量分析，工况7具有最短的平均循环周期（1.182 s），工况4次之（1.276 s），而工况3最长（2.064 s）。低床层高度工况通常产生较短的众数和平均周期，较高床层高度产生更长的右尾。增大入口空气速度从4.7至5.5 m s^-1使平均周期从约1.692 s缩短至1.396 s。

全局低速颗粒分数演变： R_s在气流启动后迅速下降，约2.5 s后趋于准稳态平台。W₁窗口平均R_s范围为16.52%至27.69%，工况7最低，工况4次低，工况3最高。三个连续窗口的绝对范围为0.03至0.39个百分点，表明工况间差异远大于残余时间波动。阈值敏感性分析表明改变速度阈值不改变九种工况的排序。

耦合流场与低速区域可视化： 三种典型工况均形成中央上升射流，但射流强度和近壁低速区域存在差异。工况3的提升力较弱，底部和壁面附近低速区域更广泛；工况7提升力更强，颗粒相对速度更大；工况4在中央提升与环形下降之间呈现较连续的连接。

高分位法向接触载荷： 27次独立初始化运行中，每次运行的时间采样活动接触力记录数为32,347,918至83,806,152条。F_95,n范围为7.50至14.50 mN，F_99,n范围为11.20至22.63 mN，运行级观测最大值范围为29.58至71.38 mN。最高观测值71.38 mN仅为转移参考力F_cr = 4.134 N的1.73%。F_95,n和F_99,n均随入口空气速度和初始床层高度描述性增加，基于F_99,n的工况排序与基于F_95,n的排序一致。

描述性主效应与多准则评估： 初始床层高度对平均循环周期、R_s和F_95,n具有最大的描述性贡献。入口空气速度也有贡献，尤其是对F_95,n，而底部循环入口间隙在所测试范围内贡献较小。多指标评分显示，工况4在等权重（0.896）和接触载荷优先权重（0.899）下排名第一，工况7在循环优先权重下排名第一（0.925）。等权重下工况4与7的分数差异仅为0.022。

原型循环对比实验： 工况3、4、7的实验平均循环周期分别为2.31 s、1.39 s和1.27 s，与模拟值的偏差分别为11.9%、8.9%和7.4%，排序与模拟一致。可见破损率分别为0.47 ± 0.06%、0.23 ± 0.06%和0.37 ± 0.06%。精确运行级置换Kruskal-Wallis检验显示存在总体工况效应（H = 6.442，精确p = 0.0214），但Holm校正后的两两比较均不显著。

五、讨论与结论

讨论部分指出，研究结果揭示了循环效率与高分位接触载荷之间的权衡关系。较高的入口速度增强了中央提升力，缩短了平均循环时间并降低了全局低速分数，但同时增加了F_95,n和F_99,n，这与Hertz-Mindlin接触关系中更大的相对碰撞速度和重叠率物理一致。初始床层高度对所有三个指标产生最大的描述性贡献，较大的颗粒库存增加了底部再入区域的局部固相分数和接触频率。工况4在三种权重方案中的两种下排名第一，但等权重下与工况7的分数差异较小，因此被描述为依赖权重的均衡候选方案而非通用最优解。

研究结论翻译如下：建立了用于导流管流化床包衣机中不规则白菜种子团簇干态循环的双向耦合CFD-DEM模型。评估框架结合了CTD、全局低速颗粒分数R_s、F_95,n和归一化高分位接触载荷比η_c,95。描述性L9筛选表明，在所测试范围内，初始床层高度对平均循环周期、R_s和F_95,n的变异贡献最大。增大入口速度缩短了循环时间并降低了R_s，但增加了F_95,n。由于L9设计无法估计交互作用或完整的二次模型，不做出验证性p值或全局优化声明。多准则计算确定工况4在等权重和接触载荷优先权重下排名第一，而工况7在循环优先权重下排名第一。原型实验复现了工况3、4、7的平均循环周期排序，支持了模拟与观测循环行为之间的定性一致性。可见破损率保持在0.5%以下；精确运行级检验显示存在总体工况效应，但Holm校正后的两两比较不显著。当前证据限于干态循环，湿法包衣可能改变摩擦、恢复系数、质量、粘附性、液桥力和团聚行为。未来研究应采用扩展响应面设计、单因子数值收敛测试、控湿种子力学测试、事件级接触统计、空间驻留时间保持度量和湿法包衣验证。