[論文レビュー] Solution of the Cauchy problem for the Navier - Stokes and Euler equations

A. Tsionskiy, Mikhail Tsionskiy|arXiv (Cornell University)|Sep 11, 2010

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

ひとこと要約

本稿では、2次元および3次元のナビエ＝ストークス方程式およびオイラー方程式のコーシー問題に対して、フーリエ変換およびラプラス変換を用いた、新しい反復的解析的解法を提示する。非線形項を反復的スキームで再定式化し、周波数領域および時間領域における線形方程式を解くことで、収束する級数解が得られ、物理的に整合性のある初期条件および外力条件下で数値的に検証され、滑らかで減衰する解が得られた。

ABSTRACT

Some known results regarding the Euler and Navier-Stokes equations were obtained by different authors. Existence and smoothness of the Navier-Stokes solutions in two dimensions have been known for a long time. Leray $\cite{jL34}$ showed that the Navier-Stokes equations in three space dimensions have a weak solution. Scheffer and Shnirelman obtained weak solution of the Euler equations with compact support in spacetime. Caffarelli-Kohn-Nirenberg improved Scheffer's results, and F.-H. Lin simplified the proof of the results of J. Leray. Many problems and conjectures about the behavior of solutions of the Euler and Navier-Stokes equations are described in the book of Bertozzi and Majda or Constantin. Solutions of the Navier-Stokes and Euler equations with initial conditions (Cauchy problem) for two and three dimensions are obtained in the convergence series form by the iterative method using the Fourier and Laplace transforms in this paper. For several combinations of problem parameters numerical results were obtained and presented as graphs.

研究の動機と目的

• R^N (N=2,3) におけるナビエ＝ストークス方程式およびオイラー方程式のコーシー問題を解くための体系的で解析的なアプローチの開発。
• フーリエ変換およびラプラス変換に基づく反復的手法を用いて、解の存在および滑らかさを確立すること。
• 無限遠における物理的に現実的な減衰条件の下で、速度場および圧力場の収束する級数解を提供すること。
• さまざまな粘性係数および外力パラメータに対して、解析的解を数値的に検証し、安定性および収束性を示すこと。

提案手法

• 非線形対流項を加算および減算することで、各反復ステップでその項を外力項に分離する形に方程式を再定式化する。
• 反復スキームにより、各ステップ j における線形偏微分方程式の系列が定義され、外力項 f_j は直前のステップの速度場 u_{j-1} に依存する： f_j = f - (u_{j-1}·∇)u_{j-1}。
• 空間変数にフーリエ変換を適用し、各フーリエモードごとに時間に関する常微分方程式（ODE）に変換することで、解法を簡略化する。
• その後、時間変数にラプラス変換を適用し、ODE をラプラス領域における代数方程式に変換する。
• ラプラス＝フーリエ領域における解を解析的および数値的に逆変換し、時間領域における速度場および圧力場を再構築する。
• 初期データおよび外力項の減衰条件（式 (5) および (6) で指定）により、級数解の収束が保証される。

実験結果

リサーチクエスチョン

• RQ1 R^2 および R^3 におけるナビエ＝ストークス方程式のコーシー問題に対して、変換手法を用いて収束する解析的解を構成可能か？
• RQ2 変換された線形方程式を用いた反復的スキームは、ナビエ＝ストークス方程式の非線形性をどのように扱うか？
• RQ3 粘性係数および外部外力パラメータを変化させた場合の解の挙動はいかにか？滑らかで減衰する性質を保つのか？
• RQ4 この手法は、ν=0 のオイラー方程式へと拡張可能か？物理的に意味のある解が得られるか？
• RQ5 さまざまなパラメータの組み合わせに対して、級数解の数値的精度および収束速度はどの程度か？

主な発見

• 反復的手法により、2次元および3次元空間におけるナビエ＝ストークス方程式およびオイラー方程式の収束する級数解が得られた。
• N=2 の場合、2次元フーリエ変換およびラプラス変換を用い、変換領域における速度および圧力の明示的表現が得られた。
• 粘性係数 ν = 0.75, 1, 1.5 および外力振幅 F_k = 0.2 から 1 のさまざまな条件に対して、数値結果が提示され、時間経過に伴い安定で減衰する解が得られた。
• 図 5.1.1 から 5.1.15 までのグラフは、異なる N 値およびパラメータ設定における速度成分の時間発展を示し、滑らかで有界な挙動を確認した。
• 各ステップで非線形項を反復的に再定式化することで、外力項に再定義された線形偏微分方程式が、初期条件が既知の形で解かれるようになり、非線形性が適切に扱われた。
• 解は非圧縮性条件（div u = 0）を満たし、空間無限遠でゼロに減衰するため、物理的期待と整合的である。