[论文解读] Variation on the Kolmogorov Forcing: Asymptotic Dissipation Rate Driven by Harmonic Forcing

Bertrand Rollin, Yves Dubief+1 more|arXiv (Cornell University)|Mar 11, 2009

Fluid Dynamics and Turbulent Flows参考文献 22被引用 15

一句话总结

本研究探讨了谐和体力的形状如何影响三维周期湍流中渐近能量耗散率。通过直接数值模拟（DNS）和变分优化，结果表明无量纲耗散因子 β 显著依赖于力的形状——在无限雷诺数下定性符合理论边界，但定量预测仍有限。

ABSTRACT

The relation between the shape of the force driving a turbulent flow and the upper bound on the dimensionless dissipation factor $β$ is presented. We are interested in non-trivial (more than two wave numbers) forcing functions in a three dimensional domain periodic in all directions. A comparative analysis between results given by the optimization problem and the results of Direct Numerical Simulations is performed. We report that the bound on the dissipation factor in the case of infinite Reynolds numbers have the same qualitative behavior as for the dissipation factor at finite Reynolds number. As predicted by the analysis, the dissipation factor depends strongly on the force shape. However, the optimization problem does not predict accurately the quantitative behavior. We complete our study by analyzing the mean flow profile in relation to the Stokes flow profile and the optimal multiplier profile shape for different force-shapes. We observe that in our 3D-periodic domain, the mean velocity profile and the Stokes flow profile reproduce all the characteristic features of the force-shape. The optimal multiplier proves to be linked to the intensity of the wave numbers of the forcing function.

研究动机与目标

• 研究非平凡的、多波数的强迫函数对三维周期湍流中无量纲耗散因子 β 上限的影响。
• 将变分优化导出的解析边界与中等雷诺数下的直接数值模拟（DNS）结果进行比较。
• 研究不同强迫形状下平均速度剖面、斯托克斯流剖面与最优乘子剖面之间的关系。
• 评估 Doering 等人（2003）提出的优化框架在有限雷诺数流动中是否能准确预测定量耗散行为。

提出的方法

• 在三维周期域中建立体力求解的不可压缩纳维-斯托克斯系统，采用形式为 f(y) = sin(ky) + A·sin(3ky) 的谐和强迫函数。
• 应用变分优化方法，基于 Doering 与 Constantin 的方法，推导无量纲耗散因子 β = εl/U³ 的上界。
• 对满足周期边界条件的纳维-斯托克斯方程执行直接数值模拟（DNS），以计算 β 和平均速度剖面。
• 将 DNS 结果与解析边界进行比较，重点关注 β 和速度剖面随强迫幅度及波数构成的变化关系。
• 分析最优乘子剖面 ψₘ，以理解其与强迫函数及平均速度场的关系。
• 利用归一化的力、平均速度和最优乘子剖面，可视化不同强迫幅度和波数下的结构关联。

实验结果

研究问题

• RQ1 多谐波体力的形状如何影响三维周期湍流中无量纲耗散因子 β 的上界？
• RQ2 在有限雷诺数的 DNS 中，基于无限雷诺数推导出的 β 解析边界在定性和定量上分别在多大程度上成立？
• RQ3 湍流中平均速度剖面在三维周期域中如何反映强迫函数的结构特征？
• RQ4 最优乘子剖面 ψₘ 与强迫函数或平均速度剖面之间存在何种关系？
• RQ5 在强迫函数中增加更高波数分量是否会引发可预测的、非线性的 β 增加？

主要发现

• 无量纲耗散因子 β 对强迫函数的形状高度敏感，除非振幅足够大，否则高波数分量对耗散的贡献较小。
• 通过变分优化导出的 β 解析上界在无限雷诺数下能定性捕捉强迫形状的依赖关系，且该行为在中等雷诺数的 DNS 中仍保持定性一致。
• 湍流中平均速度剖面再现了强迫函数的所有特征，如斜率变化和局部极值，尤其当强迫包含多个波数时更为明显。
• 最优乘子剖面 ψₘ 主要对强迫函数中的主导波数敏感，当次要分量振幅较小时，对它们不敏感。
• 随着次要波数分量振幅的增加，ψₘ 逐渐反映组合波数结构，其振幅缓慢减小，并在次要项占主导时趋于稳定。
• ψₘ 与平均速度剖面之间的关系不明确——即使平均速度剖面中存在次要特征，ψₘ 也无法反映这些特征，表明乘子与流动响应之间存在解耦。