[论文解读] Gaussian factor in the distribution arising from the nonextensive statistics approach to fully developed turbulence

A. K. Aringazin, M. I. Mazhitov|arXiv (Cornell University)|Jan 6, 2003

Statistical Mechanics and Entropy参考文献 3被引用 9

一句话总结

该论文针对充分发展湍流中流体粒子加速度的非广延统计力学模型提出了高斯截断幂律修正，通过更准确捕捉加速度概率密度函数的重尾特性，显著提升了与实验数据的拟合度。该修正模型基于Fokker-Planck框架中对伽马分布参数的截断推导得出，与高精度拉格朗日测量数据相比表现出更优的一致性，尤其在 |a|/⟨a²⟩¹ᐟ² = 60 以内的极端尾部区域，优于原始卡方模型和最近提出的对数正态模型。

ABSTRACT

We propose a simple modification, the Gaussian truncation, of the probability density function which was obtained by Beck (2001) to fit the experimental distribution of fluid particle acceleration component from fully developed fluid turbulence inspired by the framework of nonextensive statistical mechanics with the underlying gamma distribution of the model parameter. The modified distribution is given a phenomenological treatment and provides a good fit to new experimental results on the acceleration of fluid particle reported by Crawford, Mordant, and Bodenschatz. We compare the results with the new fit which was recently obtained by Beck (2002) with the use of the underlying log-normal distribution of the model parameter.

研究动机与目标

• 解决原始非广延统计力学模型在捕捉充分发展湍流中流体粒子加速度重尾行为方面的不足。
• 改进理论加速度概率密度函数（PDF）与Crawford、Mordant和Bodenschatz实验数据的拟合度。
• 检验是否通过模型参数分布的表观高斯截断可增强模型与观测到的加速度PDF尾部的一致性。
• 对比修正模型与原始卡方模型及近期提出的对数正态模型在拟合PDF和四阶矩贡献 a⁴P(a) 方面的性能。

提出的方法

• 对控制拉格朗日加速度Langevin型方程中漂移与噪声强度比的模型参数β引入高斯截断。
• 通过在截断的伽马分布上对条件高斯PDF P(a|β)进行积分，推导出修正的稳态概率密度函数，引入新参数 a_c 以控制截断尺度。
• 采用现象学方法通过伽马分布的β建模脉动的能量耗散率，并通过截断以更好地匹配实验尾部分布行为。
• 在Kolmogorov时间尺度与积分时间尺度之间存在时间尺度分离的假设下，使用Stratonovich或Ito形式的Fokker-Planck方程推导稳态PDF。
• 将所得PDF及四阶矩贡献 a⁴P(a) 与实验数据及替代模型（卡方与对数正态）在加速度幅值的全范围内进行比较。

实验结果

研究问题

• RQ1 对模型参数β进行高斯截断是否能提升非广延统计力学模型对充分发展湍流中实验加速度PDF的拟合度？
• RQ2 修正模型在捕捉加速度PDF的极端尾部分布方面表现如何，特别是当 |a|/⟨a²⟩¹ᐟ² > 30 时？
• RQ3 与原始卡方模型和对数正态模型相比，截断模型是否能更准确再现峰结构及 kurtosis 贡献 a⁴P(a)？
• RQ4 引入额外参数 a_c 对模型与实验数据及物理约束的一致性有何影响？

主要发现

• 高斯截断幂律模型在 |a|/⟨a²⟩¹ᐟ² = 60 以内对横向加速度分量的实验数据拟合良好，相较于原始卡方模型，显著改善了尾部分布行为。
• 卡方模型在预测 kurtosis 贡献 a⁴P(a) 的尾部时存在定性失败，而高斯截断模型在 |a|/⟨a²⟩¹ᐟ² < 4 和 |a|/⟨a²⟩¹ᐟ² > 30 的加速度幅值范围内均表现出良好拟合。
• 对数正态模型在 s² = 3 时整体拟合实验PDF表现良好，但在中心区域（|a|/⟨a²⟩¹ᐟ² ≤ 4）与数据存在显著偏差，而该区域主导了测量到的分布。
• 无论是对数正态模型还是高斯截断模型，均未能完全再现实验 a⁴P(a) 曲线中峰的高度与位置，表明其在捕捉完整间歇性结构方面仍存在局限。
• 当 a_c 作为自由参数时，模型能更好地与观测到的低概率尾部分布行为保持一致，表明截断是实现物理一致性的必要修正。
• 参数 a_c 与雷诺数 R_λ 的依赖关系与实验数据中观测到的雷诺数依赖性拉伸指数尾部分布一致。