涡旋晶体中的颗粒捕获

《Journal of Fluid Mechanics》：Particle trapping in vortex crystals

　　

研究人员研究了二维无黏与黏性涡旋晶体中惯性颗粒的运动，其中涡（vortex）位于正多边形顶点处，并确定了颗粒发生长期捕获的条件。在含中心涡的晶体中，研究人员发现除先前分析已得到的捕获点外还存在其他捕获点，并进一步发现颗粒轨迹可逼近一个可用渐近方法描述的极限环（limit cycle）。研究人员表明这些新发现的固定点在黏度存在时暂时持续，并研究了黏性涡合并后的不同转变过程。研究发现，在中等雷诺数（Reynolds number）下，涡晶体的环状合并形成环状涡层，可将颗粒困在原点为中心的圆盘内；在较大雷诺数下，成对涡合并导致涡数更少的晶体，颗粒聚集于新形成涡晶体的固定点处。

研究背景方面，2016年朱诺（Juno）航天器对木星极区的观测重新引发流体物理界对涡旋晶体的关注。涡旋晶体是无黏二维流动中出现的点涡规则排布结构，可绕确定角速度刚体式旋转，曾在超流氦、纯电子柱、浮于液面的旋转圆盘以及大雷诺数平面流和受迫二维湍流中被观察到。已有研究多关注无中心涡、N个卫星涡位于正多边形顶点的构型，并发现小斯托克斯数（Stokes number, St）重颗粒会在反气旋回流胞（anticyclonic recirculation cell）中心附近出现平衡点与“笼效应（cage effect）”。但含中心涡的晶体中，反气旋胞数量、颗粒捕获点结构、极限环及黏性演化尚不清楚，因此研究人员开展本研究。

主要关键技术方法如下：在旋转参考系中建立相对流函数与相对流线，利用小St下流体驻点经隐函数定理延拓得到惯性颗粒平衡点；通过流体速度梯度矩阵特征值推导驻点类型与粒子平衡稳定性；沿θ=π/N轴将速度分量化为多项式P_N(r,γ_c)=0求解驻点；在黏性情形采用黏性涡演化与涡合并数值模拟研究中雷诺数依赖的转捩；样本为理论点涡构型与数值流场，无生物样本队列。

研究结果部分，小标题2 涡旋晶体中的反气旋胞与颗粒捕获中，研究人员通过分析相对流线及P_N(r,γ_c)=0，得出无中心涡时有N+1个反气旋胞，含中心涡γ_c>0时沿θ=π/N轴出现新驻点并可共有2N个反气旋胞。小标题2.2 颗粒动力学与平衡点中，由方程u(X)+St X=0定义惯性颗粒平衡，借流体驻点特征值μ²=I(X_eq⁰)判断：I₀>0为双曲不稳定，μ²<0且μ²>?1时小St平衡点渐近稳定，从而快速甄别捕获点。小标题2.3 反气旋胞与捕获中，展示N=5时γ_c=0有6个胞、Γ_c=Γ/4有10个胞，说明中心涡改变胞数与颗粒晶格。小标题4 旋转系内一维吸引流线中，研究人员发现晶体内存在吸引流线这一one-dimensional attractor，可用渐近方法刻画。小标题5 黏度影响中，中等Re下环状合并成环状涡层并把颗粒保持在原点圆盘内；大Re下成对合并使颗粒聚到新晶体固定点；新固定点在黏度下暂存。

讨论与结论部分，研究人员总结：含中心涡显著改变反气旋胞结构与颗粒“质量晶体（mass crystal）”构型；除已知外部双曲捕获点外还存在新的稳定或极限环型捕获行为；黏性会暂时保留无黏捕获点，但Re决定合并路径与最终聚集模式。论文结论可译为：研究人员研究了位于正多边形顶点上的二维无黏与黏性涡旋晶体中惯性颗粒的运动，并确定长期颗粒捕获发生的条件。在含中心涡的晶体中，除先前分析发现的捕获点外还存在其他捕获点，且颗粒轨迹可逼近可用渐近方法描述的极限环。这些新固定点在黏度存在时暂时持续；中等雷诺数下环状合并形成环状涡层，将颗粒捕获于以原点为中心的圆盘内；较大雷诺数下成对涡合并产生涡数更少的晶体，颗粒聚集于新晶体的固定点。该研究深化了对行星极区涡晶、二维湍流颗粒输运与黏性涡合并后聚集规律的认识，论文发表于《Journal of Fluid Mechanics》。