[論文レビュー] 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

ひとこと要約

本稿は、2次元非圧縮性エーラー流体のポisson括弧とハミルトニアンにのみ依存して導かれる、定数なしの標準的分解を導入する。この方法により、自然に大スケールの凝集体（$\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.

研究の動機と目的

• 2次元非圧縮性エーラー流体における自発的スケール分離を説明する長年の課題、特にコherentな渦構造と折れ線型エネルギースペクトルの出現を解明すること。
• ハミルトニアンとポアソン括弧構造にのみ依存する、パrameterフリーで数学的に厳密な渦度の大型成分と小規模成分への分解法を構築すること。
• エーラー方程式の内在的幾何学に基づいて、クライチナンの前向きのエノストロフィー・キャスケードと後ろ向きのエネルギー・キャスケードの動的メカニズムを提供すること。
• 小スケールの力学を乗法的ノイズとして扱うことで、確率的モデル還元の基盤を確立すること。
• シミュレーションや地球物理学的流れで観測される $k^{-3}$ から $k^{-1}$ へのエネルギースペクトルの遷移の起源を、'折れ線'の挙動を捉える標準的分解によって説明すること。

提案手法

• 2次元エーラー方程式のLie–Poisson構造を用いて、球面上でポアソン括弧とハミルトニアンから標準的渦度分解 $\omega = \omega_s + \omega_r$ を導出する。
• ツァイトリンのエーラー方程式の截断モデルにヘルミート行列のスペクトル分解を適用し、渦度場の標準的分解を可能にする。
• $\omega_s$ をストリーム関数 $\psi$ の安定化部分に沿った $\omega$ の直交射影と定義し、$\omega_r$ を補空間と定義することで、エネルギーとエノストロフィーの保存を保証する。
• 射影作用素のLie微分と交換子を用いて $\omega_s$ と $\omega_r$ の結合した時間発展方程式を導出し、$\omega_s$ は $\dot{\omega}_s = [\Pi_\psi, \mathcal{L}_X]\omega$ に従って進化する。
• スペクトル法を用いて球面上で分割を数値的に実装し、時間分割スキームで方程式を解き、エネルギーおよびエノストロフィーのスペクトルを追跡する。
• ベクトル場 $b$ を定義するための暗黙的関係 $\{b, \psi\} = \Pi_\psi^\perp \Delta^{-1} \{\psi, \omega_r\}$ を用い、ポアソン構造と整合性を持つように動的生成を保証する。

実験結果

リサーチクエスチョン

• RQ1 2次元エーラー方程式のハミルトニアンとポアソン括弧から直接、定数なしの標準的渦度分解を導出できるか？
• RQ2 この内在的分割は、数値的および実験的観測と整合する大スケールの渦凝集体と小スケールの揺らぎを自然に生じさせるか？
• RQ3 標準的分解により、クライチナンが予測した前向きのエノストロフィー・キャスケードと後ろ向きのエネルギー・キャスケードの二重キャスケードを再現できるか？
• RQ4 得られた $\omega_s$ と $\omega_r$ のエネルギースペクトルが、実験的に観測された '折れ線' 挙動（$k^{-3}$ と $k^{-5/3}$ から $k^{-1}$ の勾配）を示すか？
• RQ5 この分解が、小スケールの力学を乗法的ノイズとしてモデル化する確率的モデル還元の基盤として機能できるか？

主な発見

• 標準的渦度分解 $\omega = \omega_s + \omega_r$ は、ポアソン括弧とハミルトニアンによって一意に定まり、外部パrameter やフィルタリングを一切必要としない。
• 分解は動的整合性を持つ：$\omega_s$ は大スケールの渦凝集体へ進化し、$\omega_r$ は小スケールの揺らぎを捉えるが、スケール分離を強制的に導入していない。
• 数値的シミュレーションにより、$\omega_s$ は理論的予測に一致する $k^{-3}$ のエネルギースペクトル勾配を示す。
• $\omega_r$ は $k^{-5/3}$ から $k^{-1}$ の間のエネルギースペクトル勾配を示し、2次元乱流における前向きのエノストロフィー・キャスケードと整合する。
• システムはエノストロフィーの前向きキャスケード（$\omega_r$ の平均エノストロフィーの増加）とエネルギーの後ろ向きキャスケード（$\omega_r$ の平均エネルギーの減少）を示し、クライチナンの物理的現象論を裏付けた。
• エネルギーとエノストロフィーは一貫して $H(\omega) = H(\omega_s) - H(\omega_r)$ および $E(\omega) = E(\omega_s) + E(\omega_r)$ に分割され、物理的保存則が保たれている。