[論文レビュー] On the Rate of Error Growth in Time for Numerical Solutions of Nonlinear Dispersive Wave Equations

Hendrik Ranocha, Manuel Quezada de Luna+1 more|arXiv (Cornell University)|Feb 15, 2021

Advanced Mathematical Physics Problems参考文献 83被引用数 13

ひとこと要約

本稿では、非線形分散波方程式の数値解における時間依存的誤差増大を調査し、広範な数値実験を通じて、保存則を満たす数値スキーム—特に非線形不変量を保存するスキーム—が線形誤差増大を示すのに対し、非保存スキームは2次誤差増大を示すことを示している。本研究では、非分散的双曲型系や2次元浅水域方程式を含む広範な方程式クラスに対してもこの挙動を拡張し、保存スキームが一貫して非保存スキームを上回る精度と長期的安定性を示している。

ABSTRACT

We study the numerical error in solitary wave solutions of nonlinear dispersive wave equations. A number of existing results for discretizations of solitary wave solutions of particular equations indicate that the error grows quadratically in time for numerical methods that do not conserve energy, but grows only linearly for conservative methods. We provide numerical experiments suggesting that this result extends to a very broad class of equations and numerical methods.

研究の動機と目的

• 保存的数値スキームが孤立波解に対して観察された有益な線形誤差増大が、既知の事例を超えてより広範な非線形分散波方程式のクラスにまで拡張されるかどうかを調査すること。
• 保存不変量が非二次的あるいは非線形な場合でも、保存スキームが線形誤差増大を維持するかどうかを検討すること。
• 古典的孤立波解が移動対称性を持たない非分散的1階双曲型系に対しても、この挙動が継続するかどうかを調査すること。
• 変動床深さを有する浅水域方程式のような2次元系における保存スキームの性能を評価すること。
• 理論的枠組みでエネルギーまたは質量保存が厳密に保証されていない状況においても、保存スキームが漸近的にだけでなく短時間においても顕著に小さいグローバル誤差を生じることを示すこと。

提案手法

• 和分による部分積分（SBP）作用素と同時近似項（SAT）を用いて、証明可能に安定かつ保存的な空間離散化を構築すること。
• 離散不変量を保存することができる、明示的または低コストの陰的時間積分を許容するリラクゼーションルンゲ・クッタ法を用いること。
• 高次精度の空間離散化を実現するため、フーリエスペクトル近似に基づくスペクトルコロケーション法を適用すること。
• Fornberg-Whitham方程式、Camassa-Holm方程式、Degasperis-Procesi方程式、Holm-Hone方程式、BBM-BBM方程式、および浅水域方程式など、さまざまな方程式に対して保存的および非保存的数値スキームの変種を実装すること。
• 長時間シミュレーションを実施し、保存スキームと非保存スキームの間での時間的誤差増大率を比較すること。
• 数値実験を通じて、時間経過に伴う数値解と参照解との差のL2ノルムを計算することで、誤差増大を評価すること。

実験結果

リサーチクエスチョン

• RQ1 既知の分散方程式（例：KdV、NLS）において保存スキームで観察された時間的線形誤差増大が、より広範な非線形分散波方程式のクラスにまで拡張されるか。
• RQ2 線形誤差増大は保存不変量の具体的な形式に依存するのか、それとも非二次的不変量に対しても成立するのか。
• RQ3 非分散的1階双曲型系—古典的孤立波解が存在しない系—に対しても、保存的数値スキームが優れた誤差挙動を達成できるか。
• RQ4 変動床深さを有する浅水域方程式のような2次元系に対しても、保存スキームの優位性が継続するか。
• RQ5 空間分解能や数値誤差が、浅水域方程式のような複雑な現実世界の系における観察された誤差増大率に及ぼす影響は何か。

主な発見

• Fornberg-Whitham方程式、Camassa-Holm方程式、Degasperis-Procesi方程式、Holm-Hone方程式、BBM-BBM方程式に対する保存的数値スキームは時間的に線形誤差増大を示す一方、非保存スキームは2次誤差増大を示す。
• 誤差増大の改善効果は二次不変量に限定されない。本研究では、保存量が非線形であっても線形誤差増大が維持されることを確認した。
• 標準的孤立波対称性を欠く1次元の変数係数p系に対しても、保存スキームは著しく減少した誤差増大を示し、この現象が移動不変解に限定されないことを示している。
• 変動床深さを有する2次元浅水域シミュレーションでは、保存スキームは非保存スキームと同等の計算コストで、精度を1桁以上上回った。これは、正確な孤立波解が存在しない状況でも同様に成り立つ。
• 短時間においても、保存スキームは非保存スキームよりも顕著に小さいグローバル誤差を生じ、時間スケールを問わず一貫した優位性を示している。
• 本研究では、質量保存が観察された線形誤差増大のための必要条件であることを計算的証明した。質量保存が破られたケースでは、誤差増大が2次に上昇した。