[論文レビュー] Some infinite matrix analysis, a Trotter product formula for dissipative operators, and an algorithm for the incompressible Navier-Stokes equation

Joerg Kampen|arXiv (Cornell University)|Dec 11, 2012

Advanced Mathematical Physics Problems参考文献 11被引用数 3

ひとこと要約

本稿では、n次元トーラス上での非圧縮性ナビエ＝ストークス方程式に対して、フーリエ空間における無限次元ODE系を用いた画期的なグローバル数値スキームを提示する。分散作用素に対するトロッター積分公式と無限行列解析を活用し、滑らかで多項式的に減衰する初期データに対して、時間拡張に基づく減衰機構を介して解のグローバル存在と一様有界性を確立する。任意の粘性係数 $ u > 0 $ に対して、強い双対ソボレフノルムでの収束を保証する。本手法は事前推定を提供し、収束保証付きのアルゴリズムフレームワークを構築する。

ABSTRACT

We introduce a global scheme on the n-torus of a controlled incompressible Navier-Stokes equation in terms of a coupled controlled infinite ODE-system of Fourier-modes with smooth data. We construct a scheme of global approximations related to linear partial integrodifferential equations in dual space which are uniformly bounded in dual Sobolev spaces with polynomially decaying modes. The scheme is based on some infinite matrix algebra related to weakly singular integrals, and a Trotter-product formula for dissipative operators which leads to rigorous existence results and a uniform bound for the solutions of the successive approximating linearized equations in dual space which may be otherwise represented as formal solutions in the sense of an iterated Dyson formalism.

研究の動機と目的

• 滑らかな初期データを有するn次元トーラス上での非圧縮性ナビエ＝ストークス方程式に対するグローバルに収束する数値スキームの開発。
• ナビエ＝ストークス系から導かれる繰り返し線形偏積分微分方程式の解に対する厳密な事前境界の確立。
• 強い双対ソボレフノルムにおけるグローバル上界を保つ時間拡張に基づく減衰機構の構築。
• 分散作用素および無限次元系の文脈において、オイラー型トロッター積分スキームの収束証明の構築。
• 任意の正の粘性係数 $ u > 0$ に対して強いノルムでの収束を保証する、実および複素フーリエ基底を用いたアルゴリズムフレームワークの定義。

提案手法

• ナビエ＝ストークス方程式を、多項式的減衰を有する列空間における無限ベクトルとして速度モードを表現する無限次元ODE系に変換する。
• 時間拡張変換を用いた制御的・自己制御的スキームを導入し、分散モードによる減衰を誘発することでグローバル有界性を保証する。
• 弱特異積分を含む無限行列代数に依存し、分散作用素に対するトロッター積分公式を適用して解の時間発展を近似する。
• 実数および複素数のフーリエ基底の両方でスキームを定式化し、実空間上での数値実装と誤差制御を可能にする。
• 繰り返し線形方程式に対する事前推定を導出し、強い双対ソボレフ空間におけるコンパクトネスの議論により各時間ステップで境界が保持されることを示す。
• 付録でオイラー型トロッタースキームの収束を厳密に分析し、パラメータ $ r, u, u r^2 = u_0 $ に対する条件が安定性および有界な成長を保証する。

実験結果

リサーチクエスチョン

• RQ1 滑らかな初期データに対して、すべての時間にわたり解の均一有界性を保証する、トーラス上での非圧縮性ナビエ＝ストークス方程式のグローバル数値スキームを構築可能か？
• RQ2 分散作用素に対するトロッター積分公式を無限次元系に厳密に適用することで、存在性および収束結果を得られるか？
• RQ3 時間拡張に基づく減衰機構は、強い双対ソボレフ空間における解ノルムの上界をどのように維持するか？
• RQ4 本スキームにおける線形化・繰り返し方程式に対して、外部制御関数に依存せずにグローバル存在証明が可能な事前推定を確立できるか？
• RQ5 フーリエモードに基づくアルゴリズムフレームワークは、任意の粘性係数 $ u > 0$ に対して強いノルムでの収束をどのように保証するか？

主な発見

• 滑らかでフーリエモードが多項式的に減衰する初期データに対して、n次元トーラス上での非圧縮性ナビエ＝ストークス方程式の解のグローバル存在が証明された。
• $ p > rac{n}{2}+1 $ に対して $ H^p $ 型双対ノルムにおける一様上界が保証され、減衰機構により各時間ステップで境界が維持される。
• スキームにおける線形化・繰り返し方程式に対して事前推定が確立され、強い双対ソボレフ空間におけるコンパクトネスの議論が可能となった。
• オイラー型トロッター積分スキームは $ u r^2 = u_0$ の条件下で収束し、長さ $ a $ の時間区間における元の速度成分の成長が $ u_0 u^2$ で有界となる。この成長は小さな粘性係数によって相殺可能である。
• アルゴリズムは実および複素フーリエ基底の両方で数値実装可能であり、任意の $ u > 0$ に対して強いノルムでの収束が保証される。外部制御関数は選択的であり、ゼロモードにのみ影響を及える。
• 解析により、乱流は系の動的性質として理解可能であり、文献に長年提唱されてきた仮説と整合的であることが示された。