[論文レビュー] Solitary wave simulations of the Boussinesq Systems

Özlem Ersoy Hepson, İdiris Dağ+1 more|arXiv (Cornell University)|May 16, 2016

Nonlinear Waves and Solitons参考文献 19被引用数 10

ひとこと要約

本稿では、ソリタリーモードの伝播をモデル化する正則化および古典的ボウシネスク系の初期境界値問題を数値的に解くために、指数型3次Bスプラインを用いた共役法を提案する。特定のパrameter選択において、先行する拡張型3次Bスプライン手法よりも高い精度を達成し、特にU(x,t)の誤差低減において顕著である。t=5におけるL∞ノルム比較により妥当性が検証された。

ABSTRACT

In the study, the collocation method based on exponential cubic B-spline functions is proposed to solve one dimensional Boussinesq systems numerically. Two initial boundary value problems for Regularized and Classical Boussinesq systems modeling motion of traveling waves are considered. The accuracy of the method is validated by measuring the error between the numerical and analytical solutions. The numerical solutions obtained by various values of free parameter $p$ are compared with some solutions in literature.

研究の動機と目的

• 1次元ボウシネスク系におけるソリタリーモード解の高精度な数値シミュレーション手法の開発を目的とする。
• 流体力学における長波をモデル化する非線形で分散性のある偏微分方程式を解くという課題に取り組む。
• 本稿で提案する指数型3次Bスプラインコロケーション法の性能を、文献に登場する既存の拡張型3次Bスプライン手法と比較すること。
• 離散的最大誤差ノルムを用いた誤差評価と、解析解および先行数値結果との比較を通じて、手法の精度を検証すること。

提案手法

• 空間離散化に指数型3次Bスプライン関数に基づくコロケーション法を採用する。
• 時間方向離散化には、クランク・ニコルソンスキームおよび差分法を用いる。
• 支配方程式を、指数型3次Bスプライン基底関数を用いて選択された節点でコロケートすることで、系を解く。
• 指数型3次Bスプラインにおける自由パrameter ζを調整することで、数値的精度を最適化する。
• 本手法は、古典的ボウシネスク系（CBS）と正則化ボウシネスク系（RBS）の2つの初期境界値問題に適用された。
• 数値解は、L∞誤差ノルムを計算し、解析解およびAksoyの拡張型3次Bスプライン手法の結果と比較することで妥当性を検証した。

実験結果

リサーチクエスチョン

• RQ1 指数型3次Bスプラインコロケーション法は、ボウシネスク系におけるソリタリーモードシミュレーションにおいて、既存の拡張型3次Bスプライン手法よりも高い精度を達成できるか？
• RQ2 自由パrameter ζの選択が、U(x,t)およびV(x,t)成分の数値解の精度にどのように影響するか？
• RQ3 長時間にわたるシミュレーションにおいて、本手法はソリタリーモードの形状と速度をどの程度正確に保持できるか？
• RQ4 異なる時間刻みとパrameter値の下で、本手法とAksoyの拡張型3次Bスプライン手法との間で、U(x,t)およびV(x,t)の数値誤差はどのように比較されるか？
• RQ5 古典的および正則化ボウシネスク系に対して、指数型3次Bスプラインコロケーション法の安定性および収束性はどのように振る舞うか？

主な発見

• t=5におけるU(x,t)のL∞誤差は、本手法がAksoyの拡張型3次Bスプライン手法（λ=0）と比較して著しく低く抑えられ、ζ=0.0000086530のとき0.100474×10⁻³、Aksoy法では0.551502×10⁻³にまで低下した。
• V(x,t)に関しては、両手法とも同程度の精度を示し、L∞誤差は約10⁻¹⁴のオーダーに留まり、V成分における高い安定性と精度が裏付けられた。
• 本手法は、形状の歪みなしに移動波をシミュレートでき、Cₛ=1/3で移動する波のピーク位置がt=15でx=5に一致した。
• t=5における誤差分布は、U(x,t)で最大誤差が約1×10⁻⁴、V(x,t)で約1×10⁻¹⁴であり、高精度が確認された。
• ζ=0.0000086530の選択が最適であり、Aksoy法（λ=0）と比較してU(x,t)誤差を80%以上削減した。
• 本手法は長時間シミュレーションにおいても頑健で信頼性が高く、さまざまな時間刻みにおいて波の形状を保持し、低誤差を維持した。