[论文解读] Metastability and rapid convergence to quasi-stationary bar states for the 2D Navier-Stokes Equations

Margaret Beck, C. Eugene Wayne|arXiv (Cornell University)|Aug 17, 2011

Fluid Dynamics and Turbulent Flows参考文献 17被引用 4

一句话总结

本文为环面上二维Navier-Stokes方程中bar态的亚稳态现象提供了动力系统解释。通过近似bar态附近的线性化算子，并证明其以速率$\mathcal{O}(e^{-\sqrt{\nu}t})$快速衰减，解释了为何这些态尽管在黏性时间尺度上衰减缓慢（$\mathcal{O}(e^{-\nu t})$），却能长时间持续存在，数值证据进一步证实该衰减速率对完整线性算子而言是最优的。

ABSTRACT

Quasi-stationary, or metastable, states play an important role in two-dimensional turbulent fluid flows where they often emerge on time-scales much shorter than the viscous time scale, and then dominate the dynamics for very long time intervals. In this paper we propose a dynamical systems explanation of the metastability of an explicit family of solutions, referred to as bar states, of the two-dimensional incompressible Navier-Stokes equation on the torus. These states are physically relevant because they are associated with certain maximum entropy solutions of the Euler equations, and they have been observed as one type of metastable state in numerical studies of two-dimensional turbulence. For small viscosity (high Reynolds number), these states are quasi-stationary in the sense that they decay on the slow, viscous timescale. Linearization about these states leads to a time-dependent operator. We show that if we approximate this operator by dropping a higher-order, non-local term, it produces a decay rate much faster than the viscous decay rate. We also provide numerical evidence that the same result holds for the full linear operator, and that our theoretical results give the optimal decay rate in this setting.

研究动机与目标

• 解释高雷诺数下二维不可压缩Navier-Stokes流中bar态的长寿命准稳态行为。
• 分析bar态附近的线性化动力学，识别快速收敛至这些态的机制。
• 证明对线性化算子的近似可产生远快于黏性时间尺度的衰减速率。
• 提供数值证据表明，理论衰减速率对完整线性化算子而言是最优的。
• 为将线性分析扩展至完整非线性Navier-Stokes方程奠定基础。

提出的方法

• 通过忽略一个高阶非局部项，对bar态附近的时变线性化算子进行近似，得到一个具有可分析谱性质的简化算子。
• 利用谱分析表明，该简化算子以速率$\mathcal{O}(e^{-\sqrt{\nu}t})$表现出衰减，远快于黏性衰减速率$\mathcal{O}(e^{-\nu t})$。
• 通过变量变换$w = \sqrt{1 + \Delta^{-1}} u$将线性算子转化为具有反对称分量的形式，从而便于谱分析。
• 对完整线性化算子的特征值进行数值计算，以与近似所得的理论衰减速率进行比较。
• 利用傅里叶分析和Biot-Savart定律将涡度与速度关联，并以余弦模态显式定义bar态解。
• 利用bar态是Euler方程最大熵解的事实，将其与二维湍流中物理上相关的长寿命态联系起来。

实验结果

研究问题

• RQ1 为何二维Navier-Stokes流中的bar态尽管在黏性时间尺度上不稳定，却能长时间持续存在？
• RQ2 能否通过动力系统方法解释数值模拟中观察到的快速收敛至bar态的现象？
• RQ3 线性化近似中的衰减速率$\mathcal{O}(e^{-\sqrt{\nu}t})$对完整线性化算子而言是否最优？
• RQ4 非自伴性在实现线性化动力学中快速衰减模态方面起什么作用？
• RQ5 能否将bar态的分析框架扩展至非线性Navier-Stokes方程？

主要发现

• 关于bar态的简化线性化算子表现出速率$\mathcal{O}(e^{-\sqrt{\nu}t})$的衰减，远快于黏性衰减速率$\mathcal{O}(e^{-\nu t})$。
• 对完整线性化算子特征值的数值计算证实，$\mathcal{O}(e^{-\sqrt{\nu}t})$衰减速率不仅对近似成立，对实际系统而言也是最优的。
• 近似线性算子的不变子空间确保了解快速收敛至bar态，从而为亚稳态提供了机制。
• 快速衰减机制与Burgers方程中的机制不同，其源于与$\sqrt{\nu}$成比例的小特征值，而非特征函数展开中的大系数。
• 线性化算子是非自伴的，这支撑了快速衰减模态的存在以及时间尺度的分离。
• 结果表明存在一条路径可将分析扩展至完整非线性Navier-Stokes方程，尽管线性化算子的时间依赖性使得标准谱投影技术变得复杂。