[論文レビュー] Spectrally accurate fully discrete schemes for some nonlocal and nonlinear integrable PDEs via explicit formulas

Yvonne Alama Bronsard, Xi Chen+1 more|arXiv (Cornell University)|Jan 1, 2024

Differential Equations and Numerical Methods被引用数 1

ひとこと要約

本稿では、トーラス上でのBenjamin–Ono方程式、Calogero–Sutherland DNLS方程式、および三次Szegő方程式に対して、Laxペア構造から導かれた明示的解公式を用いて、スペクトル的に正確で完全離散化された数値スキームを提示する。これらのスキームは、初期データが $H^s(\mathbb{T})$ に属する $s > 1$ の場合に、空間方向でスペクトル収束を達成し、$L^2$ ノルムにおいて誤差が $K^{-s+1}$ のオーダーで減少する。誤差定数は最終時間に線形に増加するため、計算コストを一定に保ったまま長時間シミュレーションが可能である。

ABSTRACT

We construct fully-discrete schemes for the Benjamin-Ono, Calogero-Sutherland DNLS, and cubic Szegő equations on the torus, which are $ extit{exact in time}$ with $ extit{spectral accuracy}$ in space. We prove spectral convergence for the first two equations, of order $K^{-s+1}$ in $L^2$ norm for initial data in $H^s(\mathbb T)$, $s>1$, with an error constant depending $ extit{linearly}$ on the final time instead of exponentially. These schemes are based on $ extit{explicit formulas}$, which have recently emerged in the theory of nonlinear integrable equations. Numerical simulations show the strength of the newly designed methods both at short and long time scales, thanks to the remarkable fact that the computational cost of the method is independent of the final time. These schemes open doors for the understanding of the long-time dynamics of integrable equations.

研究の動機と目的

• 非局所的かつ非線形な可積分PDEの3つの主要な方程式（Benjamin–Ono、Calogero–Sutherland DNLS、三次Szegő方程式）に対する完全離散化数値スキームの開発。
• Laxペア構造から導かれた明示的解公式を用いて、空間方向でスペクトル的正確性、時間方向で正確性を達成すること。
• 誤差が指数関数的に増加するのではなく、最終時間に線形に増加するという、厳密な収束バウンドの確立。これにより、長時間シミュレーションが可能になる。
• 最近の可積分系に対する解析的解公式の進展と、長時間ダイナミクスを扱う計算手法との橋渡し。

提案手法

• トーラス上でのBenjamin–Ono方程式、CS-DNLS方程式、三次Szegő方程式の明示的解公式（Laxペアに基づく）を活用する。
• スペクトル的空間離散化と、明示的公式による時間に正確な積分を組み合わせて完全離散スキームを構築する。
• 切り捨てられたLax作用素 $L_{u_0,K}$ を用いて離散時間発展作用素 $e^{itA_K}$ を定義し、スペクトル的正確性を保証する。
• 完全なと切り捨てられたLax作用素の差に基づく誤差解析を行い、誤差を $K^{-s}$ のオーダーで評価する。
• ユニタリ時間発展と射影作用素を用いてノルムを制御し、$H^s$ 空間における安定性バウンドを導出する。
• 帰納法と $\|u_0\|_{H^s}$ および時間依存定数を含むノルム推定を用いて収束性を確立する。

実験結果

リサーチクエスチョン

• RQ1 非局所的かつ非線形な可積分PDEに対して、明示的解公式を用いてスペクトル的に正確で完全離散化されたスキームを構築できるか？
• RQ2 初期データが $H^s(\mathbb{T})$ に属する $s > 1$ の場合、$L^2$ ノルムにおける収束速度はどの程度で、誤差は最終時間にどのように依存するか？
• RQ3 誤差定数が指数関数的に増加するのではなく、時間に線形に増加するか？これにより、計算コストを一定に保ったまま長時間シミュレーションが可能になるか？
• RQ4 時間に依存しない計算コストであるにもかかわらず、短時間および長時間スケールでの数値的性能はいかがなものか？

主な発見

• 提案されたスキームは、初期データが $H^s(\mathbb{T})$ に属する $s > 1$ の場合に、空間方向でスペクトル収束を達成し、$L^2$ ノルムにおいて誤差が $K^{-s+1}$ のオーダーで減少する。ここで $K$ はスペクトル切断パラメータである。
• 収束バウンドにおける誤差定数は最終時間 $t$ に線形に増加するが、指数関数的増加ではない。これは、標準的手法に比べて顕著な改善である。
• 本スキームの計算コストは最終時間に依存せず、可積分PDEの効率的かつ長時間のシミュレーションが可能である。
• 数値的シミュレーションにより、短時間および長時間スケールにおいてスキームの頑健さと正確性が確認された。
• 初期データが小さい場合、$\|u - u_K\|_{H^r} \leq C_5(1 + tK)^{2K^{2r-2s}}$ の誤差バウンドにおける誤差定数 $C_5$ は $\|u_0\|_{H^s}$ に比例して線形に減少する。
• 時間 $t = O(K^{-1})$ の初期段階では、最適な減少率 $K^{-s+r}$ を達成しており、初期ダイナミクスにおける高い正確性が示されている。