[論文レビュー] Error estimates of a Fourier integrator for the cubic Schr\\"odinger equation at low regularity

Alexander Ostermann, Frédéric Rousset+1 more|arXiv (Cornell University)|Feb 18, 2019

Advanced Mathematical Physics Problems参考文献 34被引用数 52

ひとこと要約

論文はキュービック非線形シュレディンガー方程式のフィルタ付き低正則性フーリエ積分器を導入し、次元1–3でH1データからの global L2 誤差評価を提供し、自然な1/2オーダーの障壁を破る。フィルタ付き離散化と離散ストリッチャーツ型推定を組み合わせ、低正則性での収束速度を向上させる。

ABSTRACT

We present a new filtered low-regularity Fourier integrator for the cubic nonlinear Schr\\"odinger equation based on recent time discretization and filtering techniques. For this new scheme, we perform a rigorous error analysis and establish better convergence rates at low regularity than known for classical schemes in the literature so far. In our error estimates, we combine the better local error properties of the new scheme with a stability analysis based on general discrete Strichartz-type estimates. The latter allow us to handle a much rougher class of solutions as the error analysis can be carried out directly at the level of $L^2$ compared to classical results \\black in dimension $d$, \\black which are limited to higher-order (sufficiently smooth) Sobolev spaces $H^s$ with $s>d/2$. In particular, we are able to establish a global error estimate in $L^2$ for $H^1$ solutions which is roughly of order $\ au^{ {1\\over 2} + { 5-d \\over 12} }$ in dimension $d \\leq 3$ ($\ au$ denoting the time discretization parameter). This breaks the "natural order barrier" of $\ au^{1/2}$ for $H^1$ solutions which holds for classical numerical schemes (even in combination with suitable filter functions).

研究の動機と目的

• キュービック NLS の数値近似における正則性の喪失を動機づけて対処する。
• 解の H1 正則性のみを要求するフィルタリングを備えたフーリエ積分器を開発する。
• 離散ストリッチャーツ型推定を用いて低正則性空間での global L2 誤差境界を確立する。
• 低正則性仮定の下で次元 d=1,2,3 にわたる時間刻みの収束速度を最適化する。

提案手法

• Pi_K という周波数射影演算子とフィルタ関数 phi_1 を用いた1ステップのフーリエ積分器を定義し、非線形相互作用を管理する。
• 局所誤差の構造を分析して有利な形 tau^{1+gamma} |∇|^{gamma} u(t) を利用し、古典的なスキームに典型的な全導関数成長を回避する。
• フィルタ済み線形伝播子 S_K(t) = e^{itΔ} Pi_K に対して、導関数の制御された喪失を含む離散ストリッチャーツ型推定を導入する。
• フィルタパラメータ K = tau^{-alpha/2} をバランスさせ、局所誤差の改善と安定性/ストリッチャーツの損失を天秤にかけて全体のglobal誤差を最適化する。
• 次元 d ≤ 3 に対する解を H1 に保つグローバル L2 誤差界を証明する：τ^{5/6}（d=1）、τ^{3/4}（d=2）、τ^{2/3}（d=3）を対数因子まで。
• 周波数を切り取った連続問題 u^K を用い、ストラッチャーツ推定と u^K - u の差の境界を組み合わせてスキームの安定性と収束性を示す。

実験結果

リサーチクエスチョン

• RQ1 キュービック NLS の初期データが H1 にある場合、低正則性フーリエ積分器はグローバルな収束速度の改善を達成できるか。
• RQ2 フィルタリングされ周波数を切り取った進化に対する離散ストリッチャーツ型推定は、低正則性領域での安定性と誤差境界にどのように影響するか。
• RQ3 d ≤ 3 に対して局所誤差と安定性を均衡させるフィルタパラメータ K（τ に対して）の最適な選択は何か。
• RQ4 低次元での H1 解の自然な τ^{1/2} オーダーの壁をどの程度打破できるか。
• RQ5 実解が H1 に留まる限り、定義解がデフocusing であっても結果は拡張されるか（発散が時間区間内に起こらない場合）？

主な発見

• 提案されたフィルタ付きフーリエ積分器は、次元1–3の H1 データに対し τ^{5/6}、τ^{3/4}、τ^{2/3}（対数因子を含む場合がある）といったグローバル L2 誤差境界を達成する。
• 誤差分析は改善された局所誤差構造と離散ストリッチャーツ型推定を組み合わせ、低正則性を直接 L2 で扱いながら高次ソボ方程式の要件を回避する。
• フィルタ済み伝播子 S_K(t) の離散ストリッチャーツ推定は、導関数の喪失を制御可能な範囲で変換し、一様な導関数喪失境界へと変換できる。
• K = τ^{-α/2} を α ≥ 1 と選ぶと、2つの誤差寄与をバランスさせ、H1 データに対して古典的方法を超える半分未満の収束を L2 で達成する。
• 結果はデフocusing なキュービック NLS に適用可能で、厳密解が H1 に留まる限りフォーカシングの場合にも拡張される（時間区間内に全体的な発散が生じない場合）。
• フィルタ済み解 u^K と真の解 u の周辺での安定性議論を通じて蓄積L2誤差のグローバルな境界を確立し、関連ノルムの時刻依存性を一様に制御する。