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

ひとこと要約

本稿は、トーラス上での2次元ナビエ＝ストークス方程式におけるバー状態の力学的不安定性について、動的システムの観点から解説する。バー状態の周囲における線形化作用素を近似し、$\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.

研究の動機と目的

• 2次元非圧縮性ナビエ＝ストークス流れにおける高レイノルズ数下での、長期間にわたる準定常的挙動を示すバー状態のメカニズムを解明すること。
• バー状態の周囲における線形化ダイナミクスを分析し、これらの状態への迅速な収束を引き起こすメカニズムを同定すること。
• 線形化作用素の近似が、粘性時間スケールよりも著しく速い減衰率をもたらすことを示すこと。
• 理論的減衰率が全線形作用素に対して最適であることを数値的証拠で示すこと。
• 非線形ナビエ＝ストークス方程式への線形解析の拡張に向けた基盤を構築すること。

提案手法

• バー状態の周囲における時間に依存する線形化作用素を、高次の非局所項を無視することで近似し、スペクトル特性が解析可能な簡略化された作用素を得る。
• スペクトル解析を用いて、簡略化された作用素が $\mathcal{O}(e^{-\sqrt{\nu}t})$ の割合で減衰することを示し、これは粘性減衰率 $\mathcal{O}(e^{-\nu t})$ よりも著しく速い。
• 変数変換 $w = \sqrt{1 + \Delta^{-1}} u$ を適用して、線形作用素を非対称成分を持つ形に変換し、スペクトル解析を容易にする。
• 全線形化作用素の固有値を数値計算し、近似から得られる理論的減衰率と比較する。
• フーリエ解析とバイオ・サバール則を用いて渦度と速度の関係を定め、バー状態解をコサインモードの形で明示的に定義する。
• バー状態がエーラー方程式の最大エントロピー解であるという事実を活用し、2次元乱流における物理的に関連のある長期間にわたる状態と結びつける。

実験結果

リサーチクエスチョン

• RQ1 2次元ナビエ＝ストークス流れにおけるバー状態は、粘性時間スケールでは不安定であるにもかかわらず、なぜ長期間にわたって持続するのか？
• RQ2 動的システムの観点から、数値シミュレーションで観測されるバー状態への迅速な収束を説明できるか？
• RQ3 線形化近似における減衰率 $\mathcal{O}(e^{-\sqrt{\nu}t})$ が、全線形化作用素に対しても最適であるか？
• RQ4 非自己随伴性は、線形化ダイナミクスにおける高速減衰モードを可能にする役割を果たすか？
• RQ5 バー状態の解析フレームワークを非線形ナビエ＝ストークス方程式へ拡張できるか？

主な発見

• バー状態の周囲における簡略化された線形化作用素は、$\mathcal{O}(e^{-\sqrt{\nu}t})$ の割合で減衰を示し、これは粘性減衰率 $\mathcal{O}(e^{-\nu t})$ よりも著しく速い。
• 全線形化作用素の固有値に対する数値計算により、減衰率 $\mathcal{O}(e^{-\sqrt{\nu}t})$ が、近似だけでなく実際のシステムに対しても最適であることが確認された。
• 近似線形作用素の不変部分空間により、解がバー状態へ迅速に収束することが保証され、準安定性のメカニズムが明らかになった。
• この高速減衰のメカニズムは、バーガース方程式とは異なり、$\sqrt{\nu}$ に比例する小さな固有値に起因するものであり、固有関数展開における大きな係数とは異なる。
• 線形化作用素は非自己随伴的であり、これが高速減衰モードの存在と時間スケールの分離を可能にする。
• 結果から、非線形ナビエ＝ストークス方程式への解析の拡張への道筋が示唆されるが、線形化作用素が時間に依存するため、標準的なスペクトル射影法の適用が複雑になる。