关于垂直振动细颗粒混合物中空气驱动的自发分离现象的研究报告
主要作者及发表信息
本研究报告基于N. Burtally、P. J. King及Michael R. Swift在英国诺丁汉大学（University of Nottingham）物理与天文学院（School of Physics and Astronomy）完成的一项原创性科学研究。该研究成果于2002年3月8日发表在《Science》期刊（Science, Vol 295, pp. 1877-1879）上，论文标题为《垂直振动细颗粒混合物中的自发空气驱动分离》（Spontaneous Air-Driven Separation in Vertically Vibrated Fine Granular Mixtures）。
学术背景与研究目的
颗粒材料（granular materials）在自然界和工业生产中普遍存在，其动力学行为是物理学和工程学的重要研究对象。对于较大颗粒，其运动主要由颗粒间及颗粒与容器壁的非弹性碰撞（inelastic collisions）控制，在振动条件下可表现出流动（flow）、对流（convection）、拱桥效应（arching）以及条纹、方形或六边形图案的形成（pattern formation）等现象。其中，广为人知的“巴西果效应”（Brazil nut effect）描述了较大颗粒在振动下上升至混合物顶部的现象。然而，自法拉第（Faraday）时代起，人们就已认识到空气会强烈影响细颗粒物（fine particulates）的运动，例如导致振动颗粒层中自发堆积（spontaneous formation of heaps）和细颗粒集合体的自发倾斜（spontaneous tilting）。尽管近期有研究报道空气对巴西果效应也存在影响，但对于这些空气驱动效应的基本机制，学界至今尚未达成普遍共识。
在颗粒复合材料（granular composites）的分离与分层（segregation and separation）方面，现有知识体系虽已相当丰富，但对许多涉及的物理过程仍缺乏清晰理解，且此前的注意力多集中于空气效应不占主导地位的大颗粒体系。因此，本研究旨在探索细颗粒混合物在空气中受垂直振动时的行为，特别是探究颗粒分离与空气影响的相互作用。研究者特别关注由性质极为相似、仅在密度上存在差异的组分构成的二元混合物所展现的分离行为，这在振动颗粒材料于周围流体存在的条件下代表了一种新颖的有序化机制（ordering mechanism），对矿物提取、制药和粉末加工等长期以来难以分离等粒径细颗粒物的行业具有重要意义。
详细研究流程
1. 实验材料与装置 研究采用的初始混合物由球形青铜（bronze）颗粒和球形玻璃（glass）颗粒组成，两者直径范围均为90至120微米（μm），但密度差异显著：青铜密度为8900千克/立方米（kg/m³），玻璃密度为2500千克/立方米（kg/m³），体积混合比例为25%的青铜对75%的玻璃。两种颗粒的动态休止角（dynamic angles of repose）非常接近，分别为玻璃23.9°±0.5°和青铜23.4°±0.5°；真空中测量显示，在低冲击速度下两者的法向恢复系数（coefficients of normal restitution）均接近1。混合物平均深度为20毫米（mm），盛放于高50毫米、水平内部尺寸为40毫米×10毫米的玻璃容器中。容器与一个能产生垂直正弦运动的机电换能器（electromechanical transducer）相连，其轴线和容器侧面与垂直方向的对准精度在1度以内。容器上附有加速度计（accelerometers）以监测运动状态。垂直振动可用无量纲参数Γ（gamma）来表征，定义为Γ = Aω²/g，其中A为振幅（amplitude），ω=2πf为角频率（angular frequency），g为重力加速度。因此，Γ代表容器最大加速度与重力加速度之比。为辅助摄影，实验中添加了微量抗静电剂（antistatic agent），已确认该操作不影响分离现象。
2. 实验参数范围与相图绘制 研究系统性地探索了Γ值从1到20、振动频率f从10赫兹（Hz）至170赫兹范围内的颗粒行为，并在频率-Γ（f-Γ）平面内绘制了包含多种现象的相图。图中标示的各个区域如下：区域A为对流和微弱分离区，在所有频率下当Γ略大于1时出现颗粒相对运动，但任何分离尝试均为对流所挫败；区域B为出现清晰分离边界（separation boundaries）的区域，其下限代表临界值Γc(f)，上限则代表几乎完全分离为富青铜相和富玻璃相的状态；区域C和D为分离边界和上表面发生振荡的区域，其中区域D特指简单的倾斜振荡（tilt oscillations）；区域E为上层颗粒被剧烈抛掷的区域；区域F为青铜夹层区，青铜在上下玻璃层之间形成中间层；区域G为高频反转区，青铜层移至容器底部附近。
3. 不同频率区间的分离演化观察 研究通过连续摄影记录了分离的演化过程，发现分离行为在约65赫兹以下的低频区和以上的高频区存在本质差异。
在低频区（f < 65 Hz），临界加速度Γc(f)近似随频率线性变化。初始均匀混合的颗粒在经历数次振动周期后迅速分离为多个具有渐趋清晰边界的区域，富青铜相优先出现在上表面附近。上表面随后隆起并发生对称性破缺倾斜。最终系统分为几乎全为玻璃的下层相和青铜比例极高的上层相，其分离边界极度清晰，宽度仅约一个颗粒直径。分离边界和上表面的形态并非静止，而是在区域C和D内发生振荡运动。在10至25赫兹、Γ略超临界值的区域D内，仅观察到周期性振荡，表现为系统在两个对称性破缺倾斜构型之间的周期性摇摆，运动完全局限于容器较大的高度和宽度平面内。当频率高于25赫兹时，需更高Γ值才能引发摇摆，且出现了涵盖两个水平平面的非周期性振荡。随Γ值进一步增至区域E，上层青铜颗粒在振动循环中被强力抛入空中，分离效果减弱。
在高频区（f > 65 Hz），Γc对频率的依赖性不强。当Γ超过Γc（区域B下限），清晰分离边界出现；至区域B上限，富青铜相与富玻璃相已完成分离。然而，富青铜相并非处于顶部，而是作为水平或略微倾斜的夹层存在于下层玻璃和上层玻璃之间（区域F）。随Γ增大，该富青铜层逐渐向表层靠近。进入区域G时，颗粒系统发生突变反转（abrupt inversion），青铜移向容器底部附近形成一薄层。在区域F上部和区域G，上层玻璃表面常显示出六边形图案（hexagonal patterning）。反转后，研究者常观察到小而活跃的青铜“水洼”（puddles）被困于玻璃表层上，其内青铜颗粒剧烈运动，与邻近玻璃颗粒的相对静止形成鲜明对比，小水洼有合并趋势。若合并后的青铜团块质量过大，颗粒玻璃表面似无法承载，团块便会形成小球穿过玻璃层沉降至主体青铜层。
在低频和高频区，分离边界均通过局部类扩散过程（diffusion-like processes）及随后的粗化（coarsening）而形成。上层或中间富青铜层一旦形成，富玻璃区和富青铜区内可各自观察到对流，但并无对流混合两区。在水平尺寸缩至10毫米×10毫米的容器中重复杂实验，青铜和玻璃仍表现出强烈的分离趋势，并可识别出低、高频区间，但简单的倾斜振荡受到抑制，且已形成的青铜层无法在Γ变化后通过反转过程移至新位置，而是通过玻璃颗粒穿过青铜层的扩散运动缓慢移动。
4. 空气压力与颗粒尺寸的对照实验 为验证空气作用，研究者缓慢排空容器内的空气。在低频区，随着气压降低，容器前后壁面的青铜向下对流日益显著。在20赫兹时气压降至10至15毫巴（mbar），或40赫兹时降至20至25毫巴，快速的全局对流即会阻止任何分离趋势；在稍低气压下，表面倾斜现象消失。在高频区（例如80赫兹时降至60毫巴，140赫兹时降至90毫巴），随气压降低，富青铜区逐渐被玻璃稀释并扩张，最终玻璃与青铜完全混合并由快速对流携带，整体分离现象消失。此外，在空气中增大颗粒尺寸，当颗粒直径增至400至500微米时，分离边界清晰度减弱，分离趋势基本消失。
5. 理论模型与机制分析 研究者将颗粒床视为多孔介质（porous medium），采用方程（1）对局部气压变化p’进行建模：φ(∂p’/∂t) = p₀∇·[(k/μ)∇p’] - p₀∇·u，其中φ为局部孔隙率（porosity），k为床层对空气的渗透率（permeability），μ为空气粘度（viscosity），u为颗粒速度，∇为梯度算符，p₀为环境气压。在高压强和低振动频率下，时间导数项可忽略，作用于颗粒的合力表现为参考系中的粘性阻尼（viscous damping）且与p₀无关。量级分析表明，当气压高于~10⁻⁶(νl²/a²)毫巴（其中a为颗粒半径，l为床层深度）时，粘性阻尼占主导，而实验观察到良好的分离恰发生于超过此值的阻尼主导区内。
由于青铜和玻璃球尺寸相近，在以相同速度运动时将经受同样的粘性阻尼力。但青铜密度更大，故该阻尼力对玻璃球运动的阻滞作用更强，导致青铜球在同等条件下更具移动性。基于此，研究者提出了如下分离机制：在振动下，初始均匀混合的床层因颗粒运动而膨胀。富青铜区因颗粒更活跃而比富玻璃区更为膨胀，组分不均导致颗粒密度空间变化。在多振动循环中，颗粒趋向于从致密区向膨胀区沿密度梯度方向扩散。然而，因青铜移动性更强，其净转移方向是从致密富玻璃区指向膨胀富青铜区，从而进一步增强效应。这种正反馈构成了空气存在下的动力学不稳定性（dynamical instability），驱动了青铜与玻璃的分离。分离后，青铜球极难重返致密的富玻璃区，而少量玻璃球则可进入更为膨胀的富青铜相。
分离后富青铜相的位置取决于两种竞争效应：较重青铜球因空气阻尼差异被抛掷得更高，而富青铜区一旦形成，需由下方碰撞提供支撑。在低频（因而振幅较大）时，前者似占主导，青铜富集于顶部；在高频时则需更详尽的动力学描述来预测粒子在分离状态下的分布。
主要研究结果
结论与研究价值
本研究确立了大密度差异、等粒径细颗粒混合物在空气中受垂直振动时会发生自发的、依赖于空气的分离这一新颖物理现象。其科学价值在于：深化了对颗粒介质与间隙流体耦合动力学的理解，揭示了有别于传统几何效应（如尺寸分离）的密度驱动分离新机制，并为非平衡态下的自组织（self-organization）和图案形成现象提供了新的模型体系。其应用价值广泛且直接：该发现为矿物浮选、制药工程中的药物组分分离、粉末冶金及化工粉体处理等长期受困于等粒径细颗粒分离难题的工业领域提供了全新的思路和潜在的技术路径。
研究亮点