[论文解读] Longitudinal and transverse velocity scaling exponents from merging of the Vortex filament and Multifractal models

K. P. Zybin, V. A. Sirota|arXiv (Cornell University)|Apr 6, 2012

Fluid Dynamics and Turbulent Flows参考文献 1被引用 3

一句话总结

本文通过将涡旋纤维模型与多重分形（MF）猜想相结合，为湍流中纵向与横向速度标度指数差异提供了一个统一的解释。基于物理论证和奇异性谱 D(h) 的抛物形式，该模型在无需拟合参数的情况下，推导出 ζₙ∥ 与 ζₙ⊥ 的解析表达式，与直接数值模拟（DNS）和实验数据高度一致，尤其在高阶矩以及纵向与横向结构函数的显著标度行为差异方面表现优异。

ABSTRACT

We suggest a simple explanation of the difference between transverse and longitudinal scaling exponents observed in experiments and simulations. Based on the Vortex filament model and Multifractal conjecture, we calculate both scaling exponents for any n without any fitting parameters and ESS anzatz. The results are in very good agreement with the data of simulations.

研究动机与目标

• 解决模拟与实验中长期存在的纵向与横向速度结构函数标度指数差异问题。
• 基于涡旋纤维的动力学及其几何特性，提供一个物理解释。
• 统一涡旋纤维模型与多重分形猜想，实现标度指数的一致推导。
• 在不依赖 ESS 假设或拟合参数的前提下，对任意 n 预测 ζₙ∥ 与 ζₙ⊥。
• 基于涡旋纤维的不同形态（轴对称与强弯曲）解释 ζₙ∥ 与 ζₙ⊥ 的差异，尤其是高阶 n 时的差异。

提出的方法

• 使用涡旋纤维模型对湍流速度结构函数进行建模，其中涡度集中在细长纤维上。
• 应用多重分形猜想，通过奇异性谱 D(h) 描述局部标度行为，满足 P ∝ l^{3−D(h)}。
• 利用大波动定理，将结构函数表示为关于 h 的积分形式：⟨Δvⁿ⟩ ∝ ∫ l^{nh} l^{3−D(h)} dμ(h)。
• 通过 ζₙ = min_h (nh + 3 − D(h)) 推导标度指数，其中假设 D(h) 为凹抛物线形式 D(h) = 3 − b(h − h₀)²。
• 通过基于纤维几何的差异，分别赋予横向与纵向贡献不同的 D(h) 函数：D⊥(h) 用于横向，D∥(h) 用于纵向。
• 利用边界条件 D⊥(0) = 3, D⊥(h₀) = 3, D∥(0) = 0, D∥(h₀) = 3 确定 D(h) 参数，从而为每种情况分别得出不同的 h₀ 与 b。

实验结果

研究问题

• RQ1 在各向同性与均匀湍流中，为何纵向与横向速度结构函数的标度指数存在差异？
• RQ2 涡旋纤维中何种物理机制导致纵向与横向速度增量表现出不同的标度行为？
• RQ3 能否在不使用拟合参数或 ESS 假设的前提下，从理论上解析推导出 ζₙ∥ 与 ζₙ⊥ 的差异？
• RQ4 涡旋纤维的几何特性——特别是轴对称性与强弯曲——如何影响标度指数？
• RQ5 何种奇异性谱 D(h) 的函数形式能够准确再现纵向与横向结构函数的观测标度指数？

主要发现

• 该模型预测 ζₙ∥ = 0.367n − 1.12×10⁻²n²（当 n ≤ 16.3 时），且当 n > 16.3 时 ζₙ∥ = 3，与实验和 DNS 数据高度一致。
• 该模型预测 ζₙ⊥ = 0.391n − 1.91×10⁻²n²（当 n ≤ 10.2 时），且当 n > 10.2 时 ζₙ⊥ = 2，与 Benzi 和 Gotoh 的 DNS 结果完全吻合。
• 标度指数的差异源于不同的纤维形态：轴对称纤维主导横向结构函数，而强弯曲纤维主导纵向结构函数。
• 该模型满足精确条件 ζ₂ = ζ₂⊥（当 n=2 时），尽管存在微小偏差（1.6×10⁻²），但可通过引入更高阶 D(h) 项进行修正，而无需改变视觉拟合效果。
• 横向与纵向情况下推导出的 D(h) 函数在物理上具有本质区别：D⊥(0) = 3 与 D∥(0) = 0，反映了在 h=0 处不同的标度行为。
• 该模型基于涡旋纤维动力学与多重分形标度的物理原理，为所有 n 提供了 ζₙ∥ 与 ζₙ⊥ 的完整解析框架，且无需任何拟合参数。