[論文レビュー] A Numerical Study Of Nonlinear Dispersive Wave Models With Spectravvave

Henrik Kalisch, Daulet Moldabayev+1 more|arXiv (Cornell University)|Jun 5, 2016

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

ひとこと要約

本稿では、擬スペクトル法と数値連続法を組み合わせたPythonベースの数値解法、SpectraVVaveを提示する。このフレームワークにより、非線形分散波方程式の移動波解を計算可能となり、Whitham方程式におけるカスプを持つ波、修正されたBenjamin–Ono方程式における波の消失、Benjamin方程式における分岐曲線の交差といった複雑な波動行動の発見が可能になる。

ABSTRACT

In nonlinear dispersive evolution equations, the competing effects of nonlinearity and dispersion make a number of interesting phenomena possible. In the current work, the focus is on the numerical approximation of traveling-wave solutions of such equations. We describe our efforts to write a dedicated Python code which is able to compute traveling-wave solutions of nonlinear dispersive equations of the general form \begin{equation*} u_t + [f(u)]_{x} + \mathcal{L} u_x = 0, \end{equation*} where $\mathcal{L}$ is a self-adjoint operator, and $f$ is a real-valued function with $f(0) = 0$. The SpectraVVave code uses a continuation method coupled with a spectral projection to compute approximations of steady symmetric solutions of this equation. The code is used in a number of situations to gain an understanding of traveling-wave solutions. The first case is the Whitham equation, where numerical evidence points to the conclusion that the main bifurcation branch features three distinct points of interest, namely a turning point, a point of stability inversion, and a terminal point which corresponds to a cusped wave. The second case is the so-called modified Benjamin-Ono equation where the interaction of two solitary waves is investigated. It is found that is possible for two solitary waves to interact in such a way that the smaller wave is annihilated. The third case concerns the Benjamin equation which features two competing dispersive operators. In this case, it is found that bifurcation curves of periodic traveling-wave solutions may cross and connect high up on the branch in the nonlinear regime.

研究の動機と目的

• 非線形分散方程式の定常的かつ対称的な移動波解を計算するための柔軟で使いやすい数値ツールの開発。
• 非線形性と分散性の競合する方程式における移動波の分岐構造と安定性の解明。
• カスプ波のような特異解の存在と性質、および波の消失のような波動相互作用現象の探求。
• AUTO や Wavetrain のような低レベルの連続法ソフトウェアに対するオープンソースで拡張可能な代替手段の提供。
• IPython ノートブックとスペクトル離散化を用いた、インタラクティブで高精度な波形計算の実現。

提案手法

• ソルバーは、自己随伴作用素 $\mathcal{L}$ とフラックス項 $f(u)$ の離散化にフーリエ変換を用いる擬スペクトル法を採用する。
• 波の速度や振幅といったパrameter空間における解分岐を追跡するために、数値連続法が用いられる。
• コードはモジュラーなクラス構造で構成され、Equation, Discretization, BoundaryCondition, Solver, Navigation の各クラスがワークフロー制御を担う。
• SpectraVVaveパッケージは、対称波形の半分のみを計算し、基本波長を $2\pi$ に設定する。
• 残差系はPDEの残差と境界条件を組み合わせ、ニュートン・ラプソン法によって強制される。
• Equationクラスを継承して新しい方程式を定義し、フーリエ乗算子とフラックス関数を定義することで、さまざまな方程式をサポートする。

実験結果

リサーチクエスチョン

• RQ1 Whitham方程式の主要分岐分岐には、折り返し点、安定性の逆転点、カスプ波による終端が存在するか？
• RQ2 修正されたBenjamin–Ono方程式において、2つのソリトン波が相互作用し、小さい方の波が完全に消失するか？
• RQ3 Benjamin方程式の周期的移動波解は、非線形領域において分岐曲線が交差するか？
• RQ4 異なる分散作用素は、移動波解の存在と正則性にどのように影響を与えるか？
• RQ5 スペクトル法と数値連続法は、特異な波形を正確に捉えるために果たす役割は何か？

主な発見

• Whitham方程式の主要分岐分岐には、折り返し点、安定性の逆転点、カスプ波による終端が存在し、波の特異性に関する数値的証拠を確認した。
• 修正されたBenjamin–Ono方程式において、2つのソリトン波が相互作用し、小さい方の波が完全に消失することが確認され、強い非線形波動相互作用の兆候が示された。
• Benjamin方程式において、周期的移動波解の分岐曲線が高振幅領域で交差し、2つの分散作用素が共存する非線形領域における複雑な挙動が明らかになった。
• SpectraVVaveコードは、スペクトルコロケーション法と擬スペクトル法を用いて、対称的移動波解を高精度に計算できた。
• ソルバーはIPythonノートブックを介して波動解のインタラクティブな探索を可能にし、従来のC/Fortranベースの連続法パッケージに比べて使いやすさが著しく向上した。
• フレームワークは、非線形分散系においてカスプ解が自然に分岐分岐の終端点として出現することを示した。