[論文レビュー] Numerical solutions of Korteweg de Vries and Korteweg de Vries-Burger's equations using computer programming

Mehri Sajjadian|arXiv (Cornell University)|Sep 9, 2012

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

ひとこと要約

本稿では、Korteweg-de Vries (KdV) 方程式およびKorteweg-de Vries-Burgers (KdVB) 方程式を数値的に解くために、sincコロケーション法とθ重み付き有限差分スキームを組み合わせた手法を提案する。この手法は、Mathematica を用いて得られる行列系を解くことで、最大絶対誤差が 1.85×10⁻¹² にまで低くなる高い精度を達成し、KdV 方程式の3つの保存則の検証も行っている。

ABSTRACT

In this paper, numerical and solitonic solutions of Korteweg de Vries(KdV) and Korteweg de Vries-Burger's (KdVB) equations with initial and boundary conditions are calculated by sinc-collocation method. The basis of method is sinc functions. First, discretizing time derivative of KdV and KdVB's equations using a classic finite difference formula and space derivatives by θ- weighted scheme between successive two time lev- els is applied, then Sinc functions are used to solve these two equations. Mathematica programming is used to solve matrix representation of these equations. KdV equation describes behavior of traveling waves which is a third order non-linear partial differential equation (PDE). Maximum absolute errors are given in Tables. The figures show approximate solutions of these two equations. Three conservation laws for KdV's equation are obtained.

研究の動機と目的

• 初期および境界条件を伴うKorteweg-de Vries (KdV) 方程式およびKorteweg-de Vries-Burgers (KdVB) 方程式を解くための堅牢な数値スキームの開発。
• これらの非線形偏微分方程式に対して、sincコロケーション法を適用し、高精度な数値解を得ること。
• 提案された数値手法を用いて、KdV 方程式の3つの保存則（質量、運動量、エネルギー）の保存性を検証すること。
• さまざまなパラメータを用いた数値実験を通じて、スキームの効率性と精度を評価すること。
• Mathematica を用いた行列解法と誤差解析を通じて、この手法の計算上の利点を示すこと。

提案手法

• KdVおよびKdVB方程式における時間微分を古典的な有限差分公式で離散化する。
• 時間ステップ間の空間微分にθ重み付きスキームを適用し、安定性と精度を確保する。
• 基底関数としてsinc関数を用い、コロケーション法を実装することで、空間近似におけるスペクトル的精度を達成する。
• 半離散化されたシステムを行列形式に構築し、Mathematica を用いて数値的に解く。
• さまざまなパラメータ設定および領域設定下で、KdVおよびKdVB方程式の両方に対してこの手法を実装する。
• KdV 方程式の保存則の計算と解析にこの手法を活用し、質量、運動量、エネルギーの不変量が保持されているかを検証する。

実験結果

リサーチクエスチョン

• RQ1 sincコロケーション法とθ重み付き時間離散化を組み合わせた手法は、KdVおよびKdVB方程式の解法において高精度を達成できるか？
• RQ2 提案されたスキームにおいて、最大絶対誤差は時間刻みと空間グリッドサイズにどのように依存するか？
• RQ3 KdV 方程式の3つの保存則が、数値スキームによってどの程度正確に保存されているか？
• RQ4 分散（μ）、粘性（ν）、非線形性（ε）の変化に伴うさまざまなパrameter領域において、この手法はどのように性能を発揮するか？
• RQ5 他の数値手法と比較して、sincコロケーション法の計算効率と収束特性はいかがなものか？

主な発見

• KdVB方程式に対して、空間および時間分解能を細かくした場合（n=16, δt=0.00001）に、最大絶対誤差が 1.85×10⁻¹² にまで低下する。
• KdV方程式に対しては、質量、運動量、エネルギーの3つの保存則が保持されており、スキームの物理的整合性が確認された。
• 数値結果から、L∞およびL2誤差が時間経過とともに徐々に増加することが示され、1つのテストケースではT=0.0009でL∞誤差が1.66×10⁻¹¹に達した。
• μ=0.001、ν=0.001、ε=1、θ=0.5 などのさまざまなパラメータセットにおいて、本手法は高い精度と安定性を示した。
• 図は、KdVおよびKdVB方程式の両方に対して滑らかで正確な近似解を示しており、ソリトン様の挙動が保持されていることが確認された。
• sincコロケーション法に起因する行列系をMathematicaで効率的に処理でき、実装の実用的妥当性が裏付けられた。