平面泊肃叶流中固体颗粒的聚焦、定向和去聚焦

《AIChE Journal?AIChE》：Solid particles focusing, orienting, and de-focusing in planar Poiseuille flow

　　

研究人员进行了平面泊肃叶流（planar Poiseuille flow）中中性浮力颗粒（neutrally buoyant particles）的三维、时间依赖的格子玻尔兹曼（Lattice-Boltzmann, LB）模拟。颗粒为球形或圆柱形（后者具有长径比（length-to-diameter ratio, L/D）在1/3至2的范围内），体积分数（volume fraction, φ）从几乎为零（单个颗粒）变化到0.15。通道的雷诺数（Reynolds number, Re）在15至150范围内。研究人员首先展示了与先前研究在单颗粒行为（横向迁移（lateral migration）和取向（orientation））方面的一致性。然后，研究人员展示了在较高固体载荷下，颗粒-颗粒相互作用（particle-particle interactions）在多大程度上影响颗粒的聚焦，以及对于圆柱形颗粒，其取向如何分布。观察到的、通常强烈不均匀的颗粒分布以及显著的取向分布对悬浮液的流动阻力（flow resistance）有明显影响。研究人员借助通道高度上的颗粒和液体剪切应力（shear stress）分布来解释后者。

**论文解读：平面泊肃叶流中固体颗粒的聚焦、定向与去聚焦行为的数值模拟研究**

**研究背景、问题与意义**

在层流剪切流中，固体颗粒会因惯性升力（inertial lift）而横跨流线迁移，产生优先集中（“聚焦”）和优先取向现象，这对悬浮液的流变行为（rheological behavior）具有潜在影响。早期实验和计算研究主要关注单个孤立颗粒在惯性剪切流中的行为，但实际应用中（如微流控芯片中的颗粒分离与排序）往往涉及多颗粒系统。随着固体载荷的增加，颗粒间相互作用可能显著改变单颗粒聚焦行为，导致散焦（de-focusing）甚至根本性的位置分布变化，进而影响分离效率或形状选择性。此外，Saffman曾推测非球形颗粒的优先取向和优先浓度会对悬浮液的有效粘度（effective viscosity）产生影响。已有实验（如Han & Kim在管流中的磁共振成像研究）表明，当总体固体体积分数低于0.1时，颗粒大致迁移到单颗粒平衡位置；而更高载荷下颗粒行为完全不同，倾向于集中在管道中心。然而，对于平面泊肃叶流（planar Poiseuille flow）中多颗粒系统的定量研究，尤其是颗粒-颗粒相互作用如何影响聚焦和取向分布，以及这些分布如何影响流动阻力和有效粘度，仍缺乏系统的数值模拟结果。为此，研究人员开展了这项工作，旨在通过高分辨率的三维时间依赖格子玻尔兹曼（LB）模拟，量化球形和圆柱形颗粒在平面泊肃叶流中的迁移、取向及其对宏观流变特性的影响。论文发表在《AIChE Journal》。

**研究开展与结论**

研究人员首先对单个颗粒（球形和圆柱形）在平面泊肃叶流中的横向迁移和取向演化进行了模拟，复现了文献中的经典结果，验证了方法的可靠性。随后，将颗粒数量从单个增加到多颗粒系统，总体固体体积分数φ从0.00625变化至0.15，雷诺数Re从15变化至150。研究系统地考察了颗粒-颗粒相互作用对聚焦位置分布、取向分布的影响，并分析了这些分布对悬浮液流动阻力及有效粘度（通过粘度系数μ_e）的影响。主要结论包括：（1）在低固体载荷下，球体迁移到单颗粒稳定位置，形成清晰的聚焦峰；随着φ增加，峰变宽，当φ足够高时（尤其在低Re下），浓度分布转为中心峰值（core-peaking）。（2）圆柱形颗粒的取向分布强烈依赖于长径比和φ：短盘（L/D≤2/3）倾向与剪切方向（y方向）对齐并自旋，长杆（L/D≥3/2）倾向垂直于剪切方向并翻滚；中间长径比颗粒的最终取向取决于初始条件。随着φ增加，取向分布变宽，优先取向减弱。（3）所有颗粒加载均导致体积流量降低，即流动阻力增加。球体悬浮液的有效粘度系数μ_e随Re降低而显著减小（从Re=150到15下降近3倍），这一趋势与颗粒-颗粒及颗粒-壁面短程相互作用产生的法向动量传递（应力）有关；圆柱体悬浮液在Re=150时μ_e随长径比增大而单调递减，反映了优先取向的影响。该研究阐明了在（微）流控装置中，其他颗粒的存在如何导致聚焦质量下降（散焦）并最终使颗粒（无论尺寸或形状）趋向于通道中心，为未来多分散系统分离装置的设计提供了基础。

**主要关键技术方法**

研究人员采用基于格子玻尔兹曼（LB）方法（Somers & Eggels变体）的三维时间依赖数值模拟，在均匀立方网格上求解纳维-斯托克斯（Navier-Stokes）方程。固体颗粒表面由密集标记点表示，利用浸入边界（immersed boundary, IB）方法在标记点处施加局部力来实现颗粒表面无滑移（no-slip）条件。通过积分表面力分布得到流体对颗粒的总水动力和力矩，进而用牛顿-欧拉方程更新颗粒的平动和转动速度，并用四元数（quaternions）处理取向更新。为补偿有限网格分辨率，显式引入短程颗粒-颗粒和颗粒-壁面相互作用：润滑力（lubrication force）采用低雷诺数解析表达式，接触力（接触弹簧及干摩擦，摩擦系数0.25）防止体积重叠。模拟在平面泊肃叶流通道（高度H，宽度与长度方向采用周期性边界）中进行，颗粒（球形或圆柱形）随机无重叠释放，体积分数从0.00625到0.15。未考虑重力以限制参数空间维度。

**研究结果**

**4.1 单颗粒迁移**
通过追踪单个球体在不同初始z位置（z₀=0.2和0.4）下的迁移轨迹，研究人员发现球体最终迁移至与Re相关的稳定位置（Re=150时为0.18H，Re=15时为0.22H），与文献（如Feng等）结果一致。低Re时迁移时间更长。对于圆柱体（以L/D=1，Re=150为例），其最终取向对初始取向敏感：若初始中心线与剪切方向（y轴）对齐，则保持自旋并稳定在z≈0.226H；若初始中心线与x轴对齐，则发生翻滚（绕y轴旋转），稳定在z≈0.216H。通过分析取向向量分量时间序列，研究人员揭示了从初始取向到最终稳定的演化过程。表1总结了不同长径比和Re下圆柱体的最终z位置和取向状态：盘状（L/D≤2/3）最终自旋（spinning），杆状（L/D≥3/2）最终翻滚（tumbling）；中间长径比（L/D≈0.86和1）行为依赖于初始取向。

**4.2 多颗粒系统**
在该部分，研究人员首先考察了球形颗粒系统。通过随机释放240个球体（φ=0.15）至10个球体（φ=0.00625），监测稳定后的z方向固体体积分数分布φ_v(z)。低φ（0.0125和0.025）时球体聚焦在单颗粒稳定位置附近（峰位置约为0.215H），随φ增大峰变宽；当φ达到0.15时，浓度分布转为中心峰值（core-peaking）（见图4）。进一步改变Re和φ的矩阵（图5）显示，中心峰值现象在低Re低φ时即可出现，且随φ增加提前发生，与Han & Kim实验（轴对称管流）一致。对于圆柱形颗粒，以L/D=1为例，在低φ时颗粒强烈趋向于中心线与剪切垂直（翻滚模式），φ升高后取向分布变宽。L/D=0.5（盘）与L/D=2（杆）在φ=0.025时分别表现出自旋和翻滚的优先取向，随φ增加取向偏好减弱（图7）。研究人员通过引入取向角θ（中心线与y轴夹角）及其分布（对比各向同性分布sinθ）量化了取向变化。

**4.3 流动阻力与有效粘度**
研究人员通过监测平均体积流量Q_s，计算表观有效粘度μ_e using Einstein型关系：μ_e=μ₀(1+αφ)，其中α为粘度系数。对于球体系统（图8），α显著依赖Re：从Re=150的约4.5降至Re=15的约1.5。α随φ变化呈现轻微最小值（出现在浓度分布由聚焦向中心峰值过渡的φ附近）。通过分解剪切应力贡献（图9），发现颗粒-颗粒和颗粒-壁面（ppw）相互作用产生的应力是主要增粘贡献，且该应力随Re降低而显著减小，解释了α的Re依赖性。对于圆柱体系统（Re=150，图10），α随L/D增大而单调递减（盘L/D=0.5的α最高，杆L/D=2的α最低），归因于优先取向效应：低φ时取向强烈，φ升高后取向分布变宽导致α趋近。

**讨论与结论翻译**

**结论（翻译自论文第5节）：**
本文报告了由体力/压力梯度驱动、位于两平行平板间的层流固-液流动的数值模拟。从单个颗粒开始，本文重现了文献中关于其在平板间优先位置以及（对于圆柱颗粒）相对于剪切方向的取向的结果。通过增加颗粒数量，观察到固体体积分数分布从靠近稳定单颗粒位置的优先固体浓度逐渐过渡到对通道中心的偏好。类似的效应在球形颗粒的实验研究中也有观察到，其过渡方式与模拟中总体固体体积分数φ和雷诺数Re的依赖关系相同。长径比L/D≤2/3的圆柱颗粒优先使其中心线与剪切（y）方向对齐，并围绕中心线旋转（自旋）。由于直接的颗粒-颗粒相互作用（碰撞）以及周围颗粒引起的液体流动扰动，取向分布随固体体积分数增加而变得不明显。L/D≥3/2的圆柱颗粒演化成其中心线垂直于剪切方向，并在通道中翻滚。与L/D≤2/3圆柱一样，L/D≥3/2圆柱的取向分布在固体体积分数增加时也会变宽。在L/D≈0.86至1范围内，单个圆柱的最终取向可以是翻滚或自旋，取决于其初始条件。无论颗粒类型如何，研究人员观察到增加颗粒数量时体积流量减小。这种增加的流动阻力被量化为固-液混合物有效粘度μ_e的增加。随后，受启发于已有文献，研究人员采用以下假设来标度从单相粘度μ₀到有效粘度的变化：μ_e=μ₀(1+αφ)，其中φ为总体固体体积分数。对于球体悬浮液，该标度较好地捕捉了φ的影响；在研究的任何雷诺数下，α在φ从0.01到0.15范围内变化不超过30%（该最大变化出现在Re=30时）。然而，研究人员确实观察到α随φ变化的系统趋势：在φ值略低于固体体积分数分布过渡到中心峰值的点处，α出现最小值。球体悬浮液的α强烈依赖于通道雷诺数：从Re=150到Re=15，α几乎减小了3倍。研究表明，这可归因于颗粒显著贡献了壁法向（z方向）的动量传递，特别是短程颗粒-颗粒和颗粒-壁面相互作用（润滑与碰撞）相关的力，而该力被证明是雷诺数的强函数。对于Re=150时的圆柱悬浮液，α随φ增加呈单调递减趋势，初步归因于优先取向效应：低φ时取向分布非常显著，φ增加时取向分布显著变宽。本研究展示了（微）流控装置中颗粒聚焦质量可能因其他颗粒的存在而受到影响，这些颗粒导致散焦，并最终使颗粒（无论尺寸或形状）倾向于向通道中心聚集。考虑到这一点，未来工作的一个方向是在泊肃叶流装置中引入多种颗粒尺寸和/或形状，并研究该装置在这些方面的区分能力。应注意，以上仅是计算研究。尽管基于以往经验（包括对所用模拟方法的实验验证）以及单颗粒迁移结果与文献的一致性，研究人员对数值方法有信心，但一项探究所观察到的转变（例如浓度分布的去聚焦和中心峰值；非球形颗粒取向角分布的展宽）的实验计划将是有意义且关键的。