[論文レビュー] Galerkin Method for the numerical solution of the RLW equation by using exponential B-splines

Melis Zorşahin Görgülü, İdris Dağ+1 more|arXiv (Cornell University)|Apr 21, 2015

Fractional Differential Equations Solutions参考文献 33被引用数 4

ひとこと要約

本稿では、浅水域波およびプラズマ波をモデル化する非線形偏微分方程式（PDE）としての正則化された長い波（RLW）方程式を解くために、指数型Bスプラインを用いたガレルキン有限要素法を提案する。この手法は高い精度と不変量の保存性を達成し、コロケーション法を上回り、既存のガレルキン手法と同等またはそれ以上の性能を示す。ベンチマークテストでは、ソリトン波の伝播、衝突、境界からの励起による波生成のシナリオで検証されている。

ABSTRACT

In this paper, the exponential B-spline functions are used for the numerical solution of the RLW equation. Three numerical examples related to propagation of single solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.

研究の動機と目的

• 非線形分散波をモデル化するRLW方程式を解くための、指数型Bスプラインに基づく堅牢な数値法の開発。
• 多項式スプラインに比べて急勾配や振動的挙動をより効果的に扱える指数型Bスプラインの特徴を活用し、既存のスプラインベース手法の改善。
• 単一ソリトン波の伝播、二重ソリトン波の衝突、境界からの駆動による波生成という3つのベンチマーク問題を通じて、手法の精度と保存性の妥当性を検証。
• 既存の数値手法と比較し、誤差ノルムおよび不変量の保存性の観点から、優れたまたは同等の性能を示すことを示す。

提案手法

• 指数型Bスプライン関数をガレルキン有限要素法における試験関数として用い、RLW方程式の解を近似する。
• RLW方程式の弱形式を採用し、反復的に解く非線形代数方程式系に帰着させる。
• 形状パrameter p を用いた双曲線関数を用いて指数型Bスプラインを定義し、急激な変化や減衰を示す解の近似を改善する。
• 時間積分にはクランク・ニコルソンスキームを採用し、時間方向に2次精度を達成するとともに、無条件安定性を確保する。
• 領域外の拡張された Knot を用いて境界条件を処理し、境界で同次Neumann条件を適用する。
• シミュレーション全体を通して不変量 $ C_1 $, $ C_2 $, および $ C_3 $ をモニタリングし、長期的な安定性と精度を評価する。

実験結果

リサーチクエスチョン

• RQ1 従来の多項式Bスプラインやコロケーション法と比較して、指数型BスプラインはRLW方程式に対してより高精度かつ安定した数値解を提供できるか？
• RQ2 指数型Bスプラインガレルキン法は、長時間にわたるシミュレーションにおいて不変量 $ C_1 $, $ C_2 $, および $ C_3 $ をどれほど良好に保存するか？
• RQ3 本手法は、ソリトン波の伝播、二重ソリトン波の衝突、境界からの駆動による波生成のシミュレーションにおいて、どの程度の精度を示すか？
• RQ4 特に波形や振幅の捉え方において、粗い空間的・時間的メッシュでも高い精度を維持できるか？
• RQ5 RLW方程式に対する以前のガレルキン法およびコロケーション法と比較して、本手法の性能と誤差ノルムはどのように評価できるか？

主な発見

• 単一ソリトン波テストにおいて、指数型Bスプラインガレルキン法（EBSGM）は $ L_inity $ 誤差ノルムが $ 1.90 imes 10^{-3} $ 未満に抑えられ、高い精度を示した。
• 二重ソリトン波衝突テストでは、時間経過に伴い不変量 $ C_1 $, $ C_2 $, および $ C_3 $ の相対誤差が $ 0.1\text{--}0.3\text{ per cent} $ の範囲に収まり、優れた保存性を示した。
• 波生成テストでは、振幅 $ U_0 = 2 $ のソリトン波の形成が成功裏にシミュレートされ、$ t = 100 $ 時点で振幅が約 3.76 に安定した。これは先行研究と整合的である。
• 計算された $ t = 100 $ 時点の波の振幅（3.76）は、[10, 25, 34] で得られた結果とよく一致しており、本手法の信頼性を裏付けた。
• 特に波伝播中の振幅誤差を最小化する点で、コロケーションに基づく指数型Bスプライン法よりも優れた精度を示した。
• 長時間シミュレーションにおいて $ C_1 $, $ C_2 $, および $ C_3 $ の保存性が安定しており、$ C_1 $ が 8.94 から 78.26 にわずかに増加するにとどまり、堅牢性が確認された。