[論文レビュー] Implementation for blow up of tornado-type solutions for complex version of 3D Navier-Stokes system

Максим Дмитриевич Арнольд, А. В. Хохлов|ArXiv.org|Jun 26, 2008

Navier-Stokes equation solutions参考文献 3被引用数 3

ひとこと要約

本稿では、複素3次元ナビエ＝ストークス系における tornadotype 解の自己相似的爆発シナリオを、純虚数の初期データを用いたフーリエ変換方程式を数値的に実装している。初期条件を (0,0,R) を中心とするスケーリングされたエルミート関数のコンpactサポートを持つ線形結合として構築することで、臨界時刻 $T_{\text{cr}}$ の近くで多項式的爆発挙動を示すエネルギー成長をシミュレートした。観察された成長率は15–20のオーダーであり、SinaとLiが予測する理論的漸近的挙動 $1/(T_{\text{cr}}-t)^5$ に近づいている。

ABSTRACT

We consider Cauchy problem for Fourier transformation of 3-dimensional Navier-Stokes system with zero external force. Using initial data purposed by Dong Li and Ya.G.Sinai we implement self-similar regime producing fast growing behavior of the energy of solution while time tends to critical value.

研究の動機と目的

• SinaとLiが提唱した、複素3次元ナビエ＝ストークス系におけるトナドタイプ解の理論的爆発シナリオを、数値的に実現すること。
• スケーリングされたエルミート関数を用いた線形結合として、(0,0,R) を中心とするフーリエ空間におけるコンパクトサポートを持つ初期条件を構築し、自己相似的爆発ダイナミクスの数値シミュレーションを可能にすること。
• 有限の臨界時刻 $T_{\text{cr}}$ に近づくに従い、解のエネルギー成長率を調査し、理論的 $1/(T_{\text{cr}}-t)^5$ 渐近的挙動を検証すること。
• 離散化パラメータおよび初期条件の形に小さな摂動を加えた場合の爆発挙動のロバストネスを検証すること。

提案手法

• 外部力ゼロ、粘性係数1として、フーリエ変換された3次元ナビエ＝ストークス方程式を用い、振幅パラメータ $A$ におけるべき級数展開により式 (1) の積分方程式を解いた。
• 初期データは $v_A(k,0) = A c_0(k)$ として構築され、$c_0(k)$ は (0,0,R) を中心とする半径 $r$ の球内にコンパクトサポートを持つ関数であり、$D=3$ までのスケーリングされたエルミート関数の線形結合で構成された。
• 3次元構造化グリッド上に、解を数値的に進める。解像度は $nX=nY=160$、$nZ=1440$、$dt=0.001$ であり、3次元畳み込み計算にはIntel MKLを用いた。
• エネルギー $\mathcal{E}(t) = \int |v(k,t)|^2 dk$ は、$L^2$ ノルムの積分の平方根として計算され、連続する時間ステップでモニタリングされた。
• 自己相似的挙動の分析には、$g_p(y,s) = h_p((0,0,R)+\sqrt{R}y,s)$ を用いた正規化が採用され、$p \to \infty$ の極限における不安定モードに注目した。再帰関係のスペクトル解析を基にした。
• 爆発の安定性は、時間ステップ、空間解像度、初期条件の形を最大25%まで摂動させることでテストした。自然なスケーリングに従って解の構造が変化するが、一貫した爆発挙動を示した。

実験結果

リサーチクエスチョン

• RQ1 フーリエ変換方程式とコンパクトサポートを持つ初期データを用いた場合、3次元ナビエ＝ストークス系におけるトナドタイプ解の自己相似的爆発シナリオを数値的に実現できるか？
• RQ2 このようなシミュレーションにおいて、臨界時刻 $T_{\text{cr}}$ に近づくに従い、エネルギー成長率 $\mathcal{E}(t)$ はどの程度の値をとるか？
• RQ3 離散化パラメータおよび初期条件の形に小さな摂動を加えた場合、爆発挙動はどの程度ロバストか？
• RQ4 数値シミュレーションは、SinaとLiが予測する理論的 $1/(T_{\text{cr}}-t)^5$ 渐近的挙動と整合性を示すか？
• RQ5 再帰関係 $h_p(k,s)$ の不安定固有モードは、爆発ダイナミクスの出現にどのような役割を果たすか？

主な発見

• エネルギー $\mathcal{M}(t) = \sqrt{\mathcal{E}(t)}$ は、$t=0.000$ で 1.00000 から $t=0.042$ で 0.3664 に減少したが、その後急激に増加し、$t=0.045$ で 19.972、$t=0.046$ で 101.20、$t=0.047$ で 698.57 に達した。これは $T_{\text{cr}}$ の近くで急速な成長を示している。
• 別の離散化（空間ステップ 0.20328、$dt=0.0008$）でも、$t=0.0424$ で 0.3047 から $t=0.0480$ で 1027.7 にエネルギーが上昇し、パrameterの変更に対しても一貫した爆発挙動を示した。
• 観察されたエネルギー成長率は15から20の間であり、理論的 $1/(T_{\text{cr}}-t)^5$ 率より著しく高いが、計算上の制限を考慮すれば予想される範囲内に収まっている。
• 時間ステップ、空間解像度、初期条件の形に最大25%の摂動を加えても、爆発挙動は安定しており、自然なスケーリングに従って解の構造が変化した。
• 数値解は明確な空間的構造を示し、フーリエ空間に「クラウド」が観察された：最初のクラウドは $z \sim 20$ に、その後のクラウドは $z \sim 40$ と $z \sim 60$ に位置しており、$h_p$ のサポートの自己相似的スケーリング $O(Rp)$ と整合的であった。
• 結果は、有限時間内に爆発するエネルギーを示す10パrameterの初期データ族の存在を支持するが、解像度の制限により正確な $1/(T_{\text{cr}}-t)^5$ 率には達しなかった。