[論文レビュー] Continuous-stage Runge-Kutta-NystrÖm methods

Wensheng Tang|arXiv (Cornell University)|Jul 9, 2018

Numerical methods for differential equations参考文献 40被引用数 4

ひとこと要約

本稿では、重み関数を組み込み、積分区間を任意の有限または無限区間に一般化することで、連続段階ルンゲ・クッタ・ニーストローム（csRKN）法の新しい枠組みを提案する。直交多項式展開と単純化仮定を活用することで、2階微分方程式に対する高次シンプレクティックかつ対称な積分法の体系的構築が可能となり、数値結果によりハミルトニアン系におけるエネルギーおよび不変量の近似保存が確認された。

ABSTRACT

We develop continuous-stage Runge-Kutta-NystrÖm (csRKN) methods in this paper. By leading weight function into the formalism of csRKN methods and modifying the original pattern of continuous-stage methods, we establish a new and larger framework for csRKN methods and it enables us to derive more effective RKN-type methods. Particularly, a variety of classical weighted orthogonal polynomials can be used in the construction of RKN-type methods. As an important application, new families of symmetric and symplectic integrators can be easily acquired in such framework. Numerical experiments have verified the effectiveness of the new integrators presented in this paper.

研究の動機と目的

• 古典的csRKN枠組みを、重み関数の導入と積分区間の任意の有限または無限領域への一般化によって拡張すること。
• RKN型法の構築に、古典的重み付き直交多項式を用いることで、精度と安定性を向上させること。
• 単純化仮定と直交展開を用いて、高次シンプレクティックかつ対称な積分法を体系的に導出するアプローチを開発すること。
• 数値実験を通じて、新しい積分法がハミルトニアン系における幾何的不変量（エネルギー、角運動量、ランゲ＝レンツ＝ポアソンベクトルなど）をいかに保持するかを示すこと。

提案手法

• csRKN形式に重み関数を導入し、離散段階集合を重み関数で定義される測度を持つ連続区間に変換する。
• 一般区間 $ I $ を用いて新しいcsRKN法を定義し、エルミート、ルジャンドル、アーベル多項式などの古典的重み付き直交多項式の使用を可能にする。
• 次数条件に単純化仮定を適用することで、複雑な非線形代数方程式を解かずに高次法の体系的導出を可能にする。
• 直交多項式展開を用いてブッチャー係数を滑らかな関数として表現し、幾何的性質の解析的取り扱いと統合を容易にする。
• 単純化仮定から導かれる代数的条件を満たす係数関数を確保することで、シンプレクティックかつ対称な積分法を構築する。
• コロケーション法およびガレルキン型アプローチを用いて数値実装を行い、最大ノルムおよびユークリッドノルムに基づく誤差解析を実施する。

実験結果

リサーチクエスチョン

• RQ1 古典的csRKN枠組みは、より広いクラスの直交多項式基底を許容するようにどのように一般化可能か？
• RQ2 csRKN定式化に重み関数を組み込むことで、より効果的かつ柔軟なRKN型積分法が得られるか？
• RQ3 直交多項式展開と単純化仮定を用いて、シンプレクティックかつ対称な積分法をどの程度体系的に導出できるか？
• RQ4 新しい積分法は、ハミルトニアン系における幾何的不変量（エネルギー、角運動量、ランゲ＝レンツ＝ポアソンベクトルなど）をどの程度正確に保存するか？
• RQ5 新しい積分法は、カオス的および非カオス的ハミルトニアン問題における長時間シミュレーションにおいて、どの程度の数値的性能を示すか？

主な発見

• 新しいcsRKN枠組みは、エルミート、ルジャンドル、アーベル多項式などの古典的重み付き直交多項式をRKN型法の構築に効果的に統合できた。
• ヘミート-4法はテストされた積分法の中で優れた精度を示したが、ヘミート-3法は次数が低いため性能が低かった。
• すべてのシンプレクティック積分法は、ハミルトニアンおよびランゲ＝レンツ＝ポアソン不変量を近似的に保存しており、二次的角運動量はマシン精度まで正確に保持された。
• ヘノン＝ハイレス問題において、数値的軌道は期待される領域内に有界に保たれ、正三角形領域からの脱出が確認されなかった。
• 位置および運動量変数の解の誤差は、時間の経過に従い線形に増加し、シンプレクティック積分法に期待される挙動と一致した。
• 単純化仮定と直交展開のおかげで、複雑な非線形次数条件を解かずに高次シンプレクティックおよび対称積分法の導出が可能になった。