[论文解读] Canonical scale separation in two-dimensional incompressible hydrodynamics

Klas Modin, Milo Viviani|arXiv (Cornell University)|Feb 2, 2021

Fluid Dynamics and Turbulent Flows参考文献 43被引用 19

一句话总结

本文提出了一种二维不可压缩欧拉流体中涡度的规范、无参数分解，完全基于泊松括号和哈密顿量推导得出。该方法从本质上将涡度分解为大尺度凝聚态（$\omega_s$）与小尺度涨落（$\omega_r$），动态再现了克雷奇南的双重级联：$\omega_s$ 展现 $k^{-3}$ 能量谱，$\omega_r$ 展现从 $k^{-5/3}$ 到 $k^{-1}$ 的标度，分别对应正向涡度通量与反向能量级联。该分解自然源自采时林的截断模型及厄米特矩阵的谱分解。

ABSTRACT

A two-dimensional inviscid incompressible fluid is governed by simple rules. Yet, to characterise its long-time behaviour is a knotty problem. The fluid evolves according to Euler's equations: a non-linear Hamiltonian system with infinitely many conservation laws. In both experiments and numerical simulations, coherent vortex structures, or blobs, emerge after an initial stage. These formations dominate the large-scale dynamics, but small scales also persist. Kraichnan describes in his classical work a forward cascade of enstrophy into smaller scales, and a backward cascade of energy into larger scales. Previous attempts to model Kraichnan's double cascade use filtering techniques that enforce separation from the outset. Here we show that Euler's equations posses an intrinsic, canonical splitting of the vorticity function. The splitting is remarkable in four ways: (i) it is defined solely via the Poisson bracket and the Hamiltonian, (ii) it characterises steady flows, (iii) without imposition it yields a separation of scales, enabling the dynamics behind Kraichnan's qualitative description, and (iv) it accounts for the "broken line" in the power law for the energy spectrum, observed in both experiments and numerical simulations. The splitting originates from Zeitlin's truncated model of Euler's equations in combination with a standard quantum-tool: the spectral decomposition of Hermitian matrices. In addition to theoretical insight, the scale separation dynamics could be used for stochastic model reduction, where small scales are modelled by multiplicative noise.

研究动机与目标

• 为解释二维不可压缩欧拉流体中长期存在的自发尺度分离难题，特别是相干涡结构的出现与折线型能量谱的形成。
• 提出一种基于哈密顿量与泊松括号结构的无参数、数学严谨的涡度分解方法，将涡度精确划分为大尺度与小尺度分量。
• 基于欧拉方程的内在几何结构，建立克雷奇南正向涡度级联与反向能量级联的动力学机制。
• 为随机模型降维奠定基础，将小尺度动力学视为乘法噪声处理。
• 通过一种规范分解解释模拟与大气海洋流动中观测到的 $k^{-3}$ 至 $k^{-1}$ 能量谱转折现象，准确捕捉‘折线’行为。

提出的方法

• 从球面上二维欧拉方程的李–泊松结构出发，仅基于泊松括号与哈密顿量推导出规范涡度分解 $\omega = \omega_s + \omega_r$。
• 将厄米特矩阵的谱分解应用于采时林的欧拉方程截断模型，实现涡度场的规范分解。
• 将 $\omega_s$ 定义为涡度 $\omega$ 在流函数 $\psi$ 的稳定子上的正交投影，$\omega_r$ 为余项，确保能量与涡度守恒。
• 利用李导数与投影算子的交换子，推导 $\omega_s$ 与 $\omega_r$ 的耦合演化方程，其中 $\omega_s$ 按 $\dot{\omega}_s = [\Pi_\psi, \mathcal{L}_X]\omega$ 演化。
• 在球面上采用谱方法实现该分解的数值模拟，使用时间分裂格式求解方程，并追踪能量与涡度谱。
• 通过隐式关系 $\{b, \psi\} = \Pi_\psi^\perp \Delta^{-1} \{\psi, \omega_r\}$ 定义生成动力学的向量场 $b$，确保与泊松结构的一致性。

实验结果

研究问题

• RQ1 能否直接从二维欧拉方程的哈密顿量与泊松括号出发，推导出一种规范的、无参数的涡度分解？
• RQ2 这种内在分解是否自然导致大尺度涡旋凝聚态与小尺度涨落，与数值与实验观测一致？
• RQ3 该规范分解能否再现克雷奇南预测的双重级联——即反向能量级联与正向涡度级联？
• RQ4 所得 $\omega_s$ 与 $\omega_r$ 的能量谱是否表现出实验观测到的‘折线’行为，即分别具有 $k^{-3}$ 与 $k^{-5/3}$ 至 $k^{-1}$ 的斜率？
• RQ5 该分解能否作为随机模型降维的基础，使小尺度动力学可建模为乘法噪声？

主要发现

• 规范涡度分解 $\omega = \omega_s + \omega_r$ 由泊松括号与哈密顿量唯一确定，无需外部参数或滤波处理。
• 该分解具有动力学一致性：$\omega_s$ 演化为大尺度涡旋凝聚态，$\omega_r$ 捕获小尺度涨落，尺度分离未被人为强制引入。
• 数值模拟显示，$\omega_s$ 的能量谱斜率约为 $k^{-3}$，与相干结构的理论预测完全一致。
• 小尺度分量 $\omega_r$ 的能量谱斜率介于 $k^{-5/3}$ 与 $k^{-1}$ 之间，与二维湍流中正向涡度级联的特征相符。
• 系统表现出正向涡度级联（$\omega_r$ 的平均涡度增加）与反向能量级联（$\omega_r$ 的平均能量减少），验证了克雷奇南的物理图像。
• 能量与涡度被一致地划分为 $H(\omega) = H(\omega_s) - H(\omega_r)$ 与 $E(\omega) = E(\omega_s) + E(\omega_r)$，物理守恒律得以保持。