[論文レビュー] Numerical solution of one-dimensional Sine--Gordon equation using Reproducing Kernel Hilbert Space Method

Ali Akgül, Inc, Mustafa|arXiv (Cornell University)|Apr 2, 2013

Numerical methods in engineering参考文献 28被引用数 3

ひとこと要約

本稿では、初期条件および境界条件を伴う1次元サインゴルドン方程式を解くための再生核ヒルベルト空間法（RKHSM）を提案する。この手法は、正確な解を再生核ヒルベルト空間内での収束級数として表現し、離散化を回避することで、最小限の計算量で高精度な数値解を得られる。数値例により、一様収束性および10^{-6}未満の誤差を確認した。

ABSTRACT

In this paper, we propose a reproducing kernel Hilbert space method (RKHSM) for solving the sine--Gordon (SG) equation with initial and boundary conditions based on the reproducing kernel theory. Its exact solution is represented in the form of series in the reproducing kernel Hilbert space. Some numerical examples have been studied to demonstrate the accuracy of the present method. The results obtained from the method are compared with the exact solutions and the earlier works. Results of numerical examples show that the presented method is simple and effective.

研究の動機と目的

• 初期条件および境界条件を伴う1次元サインゴルドン方程式を解くためのメッシュフリーな数値法を開発すること。
• 再生核ヒルベルト空間理論を用いて、正確な解の収束級数表現を提供すること。
• 数値例を通じて、本手法の高精度および計算効率を示すこと。
• RKHSMの結果を正確な解および既存の数値法（例：径数基底関数を用いたコロケーション法や有限差分法）と比較し、優れたまたは同等の性能を示すこと。

提案手法

• 本手法は再生核ヒルベルト空間理論を用いて、正確な解が無限級数として表現される解空間を構築する。
• 元の偏微分方程式を解ける演算子方程式に変換するため、再生核空間に有界線形作用素を定義する。
• カーネル関数の切り捨てられた級数を繰り返し用いて解を近似し、正確な解への一様収束を保証する。
• 空間的・時間的離散化を回避することで、大規模な線形方程式系の解法やヤコビ行列の更新の必要性を排除する。
• アルゴリズムはカーネル関数から直接近似解を計算し、計算コストを最小限に抑え、数値的不安定性を回避する。
• 収束解析により、近似解が再生核空間内で正確な解へ一様に収束することを証明する。

実験結果

リサーチクエスチョン

• RQ1 再生核ヒルベルト空間法（RKHSM）は、離散化を伴わず、初期条件および境界条件を伴う1次元サインゴルドン方程式を効果的に解くことができるか？
• RQ2 RKHSMの精度は、正確な解および径数基底関数を用いたコロケーション法や有限差分スキームなどの他の数値法と比較してどうか？
• RQ3 RKHSMの近似解の収束特性はいかにか？ また、正確な解へ一様に近づくか？
• RQ4 ニュートン反復に基づく手法（大規模な線形方程式系の解法を要する）と比較して、RKHSMの計算効率はどの程度か？
• RQ5 RKHSMは、線形および非線形の両方のサインゴルドン方程式に対して、高精度に適用可能か？

主な発見

• RKHSMは高精度を達成しており、テスト例では絶対誤差が常に5 × 10^{-7}未満であり、最良の場合には相対誤差が3.6 × 10^{-7}未満であった。
• 例5.1では、最大絶対誤差が1.756 × 10^{-7}、相対誤差が2.987 × 10^{-7}であった。これは高い精度を示している。
• 例5.2では、対称性に起因する正確な対称性と高い数値的安定性のおかげで、x = 0, ±0.4, ±0.8の点で絶対誤差がゼロとなった。
• 1例あたりのCPU時間は0.546〜1.719秒の間であり、高精度であるにもかかわらず計算コストが低いことを示している。
• 近似解が正確な解へ一様に収束するという事実が確認され、理論的収束結果が裏付けられた。
• 文献[37]の結果と比較したところ、RKHSMは全テスト点で絶対誤差および相対誤差が小さく、優れた精度を示した。