[论文解读] Caustics and Intermittency in Turbulent Suspensions of Heavy Particles

Jérémie Bec, Luca Biferale+3 more|ArXiv.org|May 8, 2009

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

一句话总结

本研究探讨了重惯性粒子湍流悬浮体中的速度统计特性，揭示了在粘性子区中由于焦散现象导致的速度差呈现准双分形标度，而在惯性区中则表现出与尺度相关的霍尔德指数。高分辨率模拟表明，焦散现象主导了高阶结构函数，而粒子动力学则受与尺度相关的斯托克斯数支配，这对理解云中早期雨滴形成具有重要意义。

ABSTRACT

The statistics of velocity differences between very heavy inertial particles suspended in an incompressible turbulent flow is found to be extremely intermittent. When particles are separated by distances within the viscous subrange, the competition between quiet regular regions and multi-valued caustics leads to a quasi bi-fractal behavior of the particle velocity structure functions, with high-order moments bringing the statistical signature of caustics. Contrastingly, for particles separated by inertial-range distances, the velocity-difference statistics is characterized in terms of a local Hölder exponent, which is a function of the scale-dependent particle Stokes number only. Results are supported by high-resolution direct numerical simulations. It is argued that these findings might have implications in the early stage of rain droplets formation in warm clouds.

研究动机与目标

• 量化焦散现象与偏倚集中度在重惯性粒子湍流中速度差统计特性中的作用。
• 确定间歇性与多尺度动力学如何影响稀疏悬浮体中碰撞率，这些悬浮体与暖云物理密切相关。
• 探究由于焦散现象与粒子惯性相互作用，速度结构函数是否表现出双分形或多重标度行为。
• 建立与尺度相关的斯托克斯数和惯性区中局部标度指数之间的联系。
• 评估这些发现对雨滴形成早期阶段的相关性，此时粒子碰撞至关重要。

提出的方法

• 在 Rλ ≈ 400 和 Rλ ≈ 185 条件下对不可压缩纳维-斯托克斯湍流进行直接数值模拟（DNS），以解析全部尺度范围。
• 通过运动方程 τ d²X/dt² = −dX/dt + u(X,t) 计算粒子轨迹，其中 τ 为斯托克斯时间，St = τ/τη 为斯托克斯数。
• 计算粒子对间距为 r 时的速度结构函数 Sₚ(r;St) = ⟨|Ṙ|ᵖ | R = r⟩，其中 |Ṙ| 表示纵向速度差。
• 采用鬼粒子方法计算碰撞核，其中 κ(r;St) ∝ r^{D₂−1} S₁(r;St)，将结构函数与碰撞率关联。
• 分析局部霍尔德指数 ξₚ(r;τ) = d(ln Sₚ)/d(ln r) 作为与尺度相关的斯托克斯数 St(r) = ε^{1/3}τ / r^{2/3} 的函数。
• 对不同 r 和 St 下的速度差概率密度函数（PDF）进行统计分析，揭示幂律核心与拉伸指数尾部，表明焦散现象的贡献。

实验结果

研究问题

• RQ1 在粘性子区（r ≪ η）中，焦散现象如何影响重惯性粒子间速度差统计特性？
• RQ2 由于间歇性焦散现象的影响，速度结构函数 Sₚ(r;St) 在多大程度上表现出双分形或多重标度行为？
• RQ3 局部标度指数 ξₚ(r;τ) 是否仅依赖于与尺度相关的斯托克斯数 St(r)，如自相似流模型所预测？
• RQ4 偏倚集中度与速度差之间的相互作用如何影响稀疏悬浮体中的碰撞核？
• RQ5 这些发现对暖云中雨滴形成早期阶段的相关性如何，此时粒子碰撞至关重要？

主要发现

• 在粘性子区（r ≪ η）中，速度结构函数 Sₚ(r;St) 由于规则区域与多值焦散现象之间的竞争，表现出准双分形行为，高阶矩主要由焦散现象贡献。
• 归一化速度差 σ = Ṙ/R 的概率密度函数（PDF）具有重尾核心，幂律行为为 ∝ |σ|^{-(1+ξ∞)}（当 |σ| ≲ σ*(r) 时），表明当 p > ξ∞ 时矩的发散。
• 当 |σ| ≳ σ*(r) 时，PDF 的拉伸指数尾部与焦散现象尺寸统计相关，可能导致多重标度，但其对结构函数的贡献尚不明确。
• 在惯性区（η ≪ r ≪ L）中，一阶局部霍尔德指数 ξ₁(r;τ) 在以与尺度相关的斯托克斯数 St(r) 为横坐标绘制时发生收敛，证实其对 r 和 τ 的主导作用。
• 当 St(r) 较小时，ξ₁(St(r)) 趋近于 K41 预测值（ξ₁ = 1/3），表明具有示踪粒子行为；而当 St(r) 较大时，其显著偏离，反映增强的间歇性。
• 结果表明，在接近单分散的雨滴的暖云早期阶段，即使在重力沉降占主导之前，焦散现象也可能显著增强粒子碰撞率。