[论文解读] Clustering of heavy particles in the inertial range of turbulence

Jérémie Bec, Massimo Cencini+1 more|ArXiv.org|Jun 14, 2006

Particle Dynamics in Fluid Flows参考文献 1被引用 3

一句话总结

本文研究了惯性尺度湍流中重粒子的团聚行为，表明粒子分布受局部斯托克斯数控制，并因流体速度场的粗糙性而减弱。通过采用δ-相关、自相似的速度场模型，本文证明了速度差概率密度函数（PDF）中幂律尾部——由罕见的近距离接近事件驱动——是团聚现象的根源，且随着赫ölder指数h的减小，团聚效率降低。

ABSTRACT

A statistical description of heavy particles suspended in incompressible rough self-similar flows is developed. It is shown that, differently from smooth flows, particles do not form fractal clusters. They rather distribute inhomogeneously with a statistics that only depends on a local Stokes number, given by the ratio between the particles' response time and the turnover time associated to the observation scale. Particle clustering is reduced when increasing the fluid roughness. Heuristic arguments supported by numerics are used to explain this effect in terms of the algebraic tails of the probability density function of the velocity difference between two particles. These tails are a signature of events during which particle couples approach each other very closely.

研究动机与目标

• 理解湍流惯性尺度中重惯性粒子的统计分布，其中粒子响应时间破坏了自相似性。
• 分离并分析耗散动力学（粘性阻力）在粒子团聚中的作用，排除持久涡旋或离心喷射的影响。
• 研究载体流的粗糙性——通过赫ölder指数h量化——如何通过相对粒子运动的统计特性调节粒子团聚。
• 建立速度差概率密度函数（PDF）中幂律尾部与粒子团聚机制之间的联系。

提出的方法

• 将载体流速度场建模为平稳、均匀、各向同性、高斯分布、时间δ-相关的场，具有由赫ölder指数h表征的自相似结构。
• 通过朗之万型方程（公式1）定义粒子对的相对运动，其中粒子响应时间τ和流体速度差δu(𝐫, t)决定动力学行为。
• 引入尺度相关的局部斯托克斯数S(ℓ) = D₁τ / ℓ²⁽¹⁻ʰ⁾，该参数捕捉了在尺度ℓ上粒子响应时间与流体周转时间的比值。
• 通过变量替换，将粒子对的动力学简化为三维随机系统（公式6–8），从而能够分析相空间中的大环路轨迹。
• 通过估计粒子在极端区域停留的时间比例，分析纵向速度差X的PDF，推导出当x ≪ −1时，P< (x) ∝ |x|⁻²/ʰ的幂律尾部分布。
• 通过模拟在不同h值下数值验证预测的幂律尾部指数α = 1 + 2/h，并与理论预测进行比较。

实验结果

研究问题

• RQ1 如何量化由赫ölder指数h表征的湍流速度场粗糙性，影响惯性尺度中重惯性粒子的团聚行为？
• RQ2 在排除涡旋喷射影响的前提下，耗散动力学（粘性阻力）在粒子团聚中起什么作用？
• RQ3 相对速度差PDF中的幂律尾部与粒子团聚机制有何关联？
• RQ4 局部斯托克斯数S(ℓ)在多大程度上能完全表征尺度ℓ上粒子团聚的统计特性？
• RQ5 观测到的速度差幂律尾部分布是否在三维流动中依然存在，还是仅限于所用的二维模型？

主要发现

• 惯性尺度中的粒子团聚行为并非尺度不变，也不会形成分形结构，这与黏性尺度中的行为相反。
• 粒子的统计分布仅依赖于局部斯托克斯数S(ℓ)，该参数结合了粒子响应时间τ和尺度ℓ上的流体周转时间。
• 随着赫ölder指数h减小（即流体速度场变得越粗糙），粒子团聚显著减弱，原因在于极端速度差事件发生的概率降低。
• 纵向速度差X的PDF表现出幂律尾部，其指数为α = 1 + 2/h，该结果在多个h值下均通过数值模拟得到验证。
• 速度差PDF中的幂律尾部源于相空间中罕见的大环路轨迹，这些轨迹对应于粒子极为接近的运动过程。
• 在缺乏持久涡旋的条件下，团聚机制主要由耗散动力学主导，且由于近距离接近事件被抑制，粗糙流体中的团聚效率更低。