[論文レビュー] Research of convergence of the iterative method for solution of the Cauchy problem for the Navier - Stokes equations based on estimated formula

A. Tsionskiy, Mikhail Tsionskiy|arXiv (Cornell University)|Jan 10, 2011

Iterative Methods for Nonlinear Equations被引用数 6

ひとこと要約

本稿では、2次元および3次元のナビエ＝ストークス方程式の初期値問題（カウチ問題）を解くための反復法の収束を、フーリエ変換およびラプラス変換を用いて調査する。収束境界の推定式 $ F_n/( u \mu^4) < 1 $ を確立し、複数の円周モードで数値的に検証した。速度振幅の比較により、すべてのテストパラメータで $ |̄u_2^*| \ll |̄u_1| $ が成立し、収束が著しく速いことが確認された。

ABSTRACT

Solution of the Navier-Stokes equations with initial conditions (Cauchy problem) for 2D and 3D cases was obtained in the convergence series form by iterative method using Fourier and Laplace transforms in paper $\cite{TT02}$. For several combinations of problem parameters numerical results were obtained and presented as graphs. Estimated formula for the border of the parameter area of convergence of the iterative method was obtained in paper $\cite{TT03}$. This paper describes numerical proof of convergence of the iterative method for solution of the Cauchy problem for the Navier - Stokes equations. Usage of estimated formula for the border of the parameter area of convergence of the iterative method is shown for wide ranges of the problem's parameters.

研究の動機と目的

• 2次元および3次元のナビエ＝ストークス方程式のカウチ問題を解く反復法の収束挙動を分析すること。
• 先行研究で導出された反復法の収束境界の推定式を検証すること。
• 粘度、力の振幅、および径方向モードパラメータの変動に対して、反復法の有効性と頑健性を数値的に示すこと。
• 1回目と2回目の反復ステップ間の速度振幅比較を通じて、収束の視覚的および定量的証拠を提供すること。
• 固定された $ F_n $ および $ \nu $ に対して、$ \mu $ の値が高くなるほど収束が改善することを確認すること

提案手法

• 反復法を、初期速度がゼロで極座標系における径方向力が作用するナビエ＝ストークス方程式に適用する。
• 2次元および3次元のカウチ問題に対して、フーリエ変換およびラプラス変換を用いて、級数形の解析解を導出する。
• 収束境界の推定式として $ F_n/( \nu \mu^4) < 1 $ を導出する。この式は、すべての円周モード $ n = 1,2,3,\dots $ に適用可能である。
• $ n = 1 $ から $ 5 $ まで、$ F_n $、$ \mu $、$ \nu $ を変化させた数値シミュレーションを実施。$ F_n $ は $ k = 0,1,2,3 $ に対して $ 10^k/n $ に設定した。
• $ \nu = 0.01 $ および他の値について、$ |\vec{u}_1| $（1回目）および $ |\vec{u}_2^*| $（2回目）の速度振幅を、対数スケールでグラフィカルに比較する。
• $ F_n $-$ \mu $、$ F_n $-速度、$ \nu $-速度の平面にグラフをプロットし、収束領域を $ F_n/( \nu \mu^4) = 1 $ の曲線の上側領域として定義する。

実験結果

リサーチクエスチョン

• RQ1 推定収束境界 $ F_n/( \nu \mu^4) < 1 $ は、ナビエ＝ストークス方程式のカウチ問題における反復法の収束領域を正確に予測できるか？
• RQ2 異なる円周モードおよびパラメータ設定において、反復法の1回目と2回目の反復における速度振幅はどのように変化するか？
• RQ3 粘度 $ \nu $、力の振幅 $ F_n $、径方向モードパラメータ $ \mu $ の広い範囲において、反復法の収束は頑健か？
• RQ4 $ |\vec{u}_1| $ と $ |\vec{u}_2^*| $ の間の定量的関係は何か？ すべてのテスト設定で $ |\vec{u}_2^*| \ll |\vec{u}_1| $ が成立するか？
• RQ5 $ \mu $ を増加させると収束速度にどのような影響があるか？ 固定された $ F_n $ および $ \nu $ に対して、$ \mu $ を大きくすると収束が改善するか？

主な発見

• 反復法は、$ n = 1 $ から $ 5 $ までのすべてのテストされた円周モードで強く収束し、すべてのケースで $ |\vec{u}_2^*| \ll |\vec{u}_1| $ が成立した。
• すべてのテストされたパラメータ組み合わせ（$ F_n $、$ \mu $、$ \nu $）で、$ 0.2 \leq F_n \leq 1000 $ および $ 2.0 \leq \mu \leq 18.0 $ の範囲が、$ F_n/( \nu \mu^4) < 1 $ で定義される収束領域内に位置した。
• $ \nu = 0.01 $ の場合、すべての数値結果が曲線 $ F_n/( \nu \mu^4) = 1 $ の上側に黒丸としてプロットされ、推定境界の有効性が確認された。
• 2回目の反復における速度振幅 $ |\vec{u}_2^*| $ は、1回目の反復における $ |\vec{u}_1| $ より著しく小さく、急激な減衰と強い収束を示した。
• $ \mu $ を大きくするほど収束が改善する傾向があり、固定された $ F_n $ および $ \nu $ に対して、$ \mu $ の値が高くなると収束速度が向上することが示唆された。
• $ 0 \leq r \leq 3 $ の $ \varphi $-平面におけるグラフィカル比較により、すべての半径および角度で $ |\vec{u}_2^*| < |\vec{u}_1| $ であることが確認された。ズームインしたプロットにより、2回目のステップでの低振幅挙動が顕著に可視化された。