[论文解读] The role of microscopic fluctuations in transition prediction

Paolo Luchini|ArXiv.org|Apr 13, 2008

Phase Equilibria and Thermodynamics参考文献 4被引用 7

一句话总结

本文表明，源自涨落-耗散定理的微观热涨落，可作为剪切流中湍流转捩的物理解释充分且一致的初始扰动源，其结果与经验观测的 $e^N$ 转捩准则 $N \approx 9\text{--}10$ neper 完全一致。通过随机线性稳定性理论，研究证明热噪声本身即可驱动边界层不稳定性达到破坏阈值，从而在转捩预测模型中无需引入外部扰动。

ABSTRACT

The commonly accepted description of transition to turbulence in shear flows requires the presence of an external source of disturbances that get amplified by an essentially linear mechanism up to the point where breakdown to turbulence occurs. Microscopic fluctuations are shown here to provide just the right amount of initial disturbances to match the predictions of linear stability theory.

研究动机与目标

• 探究微观热涨落是否可作为剪切流中湍流转捩所需初始扰动的物理解释来源。
• 解决长期存在的转捩预测难题：即假设存在外部扰动（如声学、涡量扰动），但其真实特性无法合理指定。
• 检验热噪声作为放大器输入的物理下限，是否可重现经验观测的 $e^N$ 转捩阈值。
• 利用随机线性系统理论，建立热噪声与剪切流中湍流起始之间的理论与定量联系。

提出的方法

• 将随机驱动线性系统的形式化方法应用于包含代表热涨落的随机应力张量的纳维-斯托克斯方程。
• 采用拉普拉斯变换与本征模分解，表达系统对白噪声输入的响应，重点关注主导不稳定模态。
• 利用涨落-耗散定理推导速度扰动的均方根（rms）振幅，其中噪声功率谱与温度及流体物性相关。
• 采用雷诺数 $Re$、参考长度 $h$ 及新引入的材料长度 $\lambda = k_B T / (\rho \nu^2)$ 进行无量纲化，得到响应的无量纲形式。
• 针对准平行流构建空间受激问题，使结果可与 $e^N$ 方法在距离上的放大特性进行比较。
• 对空间受激问题进行数值积分，结果确认仅由热噪声即可产生 $u_{\text{rms}} \approx 1\%$ 的来流速度，且此时 $N \approx 9\text{--}10$。

实验结果

研究问题

• RQ1 仅靠微观热涨落是否足以解释剪切流中湍流转捩所需的初始扰动？
• RQ2 边界层对热噪声的响应是否重现了经验观测的 $e^N$ 转捩准则 $N \approx 9\text{--}10$ neper？
• RQ3 热噪声强度与边界层中托尔明-施利希特林波放大之间的定量关系为何？
• RQ4 由流体物性导出的无量纲参数 $\lambda/h$ 如何影响转捩阈值？
• RQ5 $e^N$ 方法的经验阈值是否与线性放大器系统中热噪声所设定的理论下限一致？

主要发现

• 由热噪声引起的均方根速度涨落，在 $N \approx 9\text{--}10$ 处达到来流速度的 $1\%$，与 $e^N$ 方法预测的转捩位置完全一致。
• 理论估算 $N = 0.5 \log(h/\lambda)$ 在空气温度 300 K、边界层厚度 $h = 1$ mm 时给出 $N \approx 9.00$，与经验阈值完全吻合。
• 标准空气在 300 K 与大气压下，特征长度 $\lambda = 1.508 \times 10^{-11}$ m，为由 $k_B T / (\rho \nu^2)$ 导出的材料属性。
• 该模型表明，热噪声为任何放大器（包括边界层）的输入噪声设定了有限且不可逾越的下限，因此在转捩预测中无需引入外部扰动。
• 数值结果证实，即使在无大气声或涡量扰动的情况下，仅由热涨落即可导致 $N \approx 9\text{--}10$ 处发生转捩。
• 响应的无量纲形式 $\left<|C^{(1)}|^2\right> \propto \frac{\lambda}{h} Re^{-3} \times \text{空间积分}$ 表明，当 $\lambda/h \sim e^{-2N}$ 时，$e^N$ 阈值自然由热噪声产生。