[論文レビュー] Continuous-stage Runge-Kutta methods based on weighted orthogonal polynomials

Wensheng Tang|arXiv (Cornell University)|May 25, 2018

Numerical methods for differential equations参考文献 43被引用数 9

ひとこと要約

本稿では、重み付き直交多項式（例えばジャコビ多項式）を用いて、複雑な非線形次数条件を解く必要を回避する、連続段ルンゲ・クッタ（csRK）法を構築するための新規フレームワークを提案する。整備された重み付き補間型求積理論を活用することで、シンプレクティック法や対称法といった幾何的積分法を体系的に構築可能となり、csRK法の適用範囲がルジャンドル多項式に限定されなくなる。

ABSTRACT

We develop continuous-stage Runge-Kutta methods based on weighted orthogonal polynomials in this paper. There are two main highlighted merits for developing such methods: Firstly, we do not need to study the tedious solution of multi-variable nonlinear algebraic equations associated with order conditions; Secondly, the well-known weighted interpolatory quadrature theory appeared in every numerical analysis textbook can be directly and conveniently used. By introducing weight function, various orthogonal polynomials can be used in the construction of Runge-Kutta-type methods. It turns out that new families of Runge-Kutta-type methods with special properties (e.g., symplectic, symmetric etc.) can be constructed in batches, and hopefully it may produce new applications in numerical ordinary differential equations.

研究の動機と目的

• 従来のcsRK法がルジャンドル多項式に依存していたという制限を克服すること。
• 任意の重み付き直交多項式を用いたcsRK法の構築に向けた一般化されたフレームワークを構築すること。
• 多変数非線形代数系の解法を回避することで、次数条件の導出を簡略化すること。
• 求積に基づく設計により、構造保存型積分法（例：シンプレクティック、対称）の体系的構築を可能にすること。
• 直交多項式技術の適用範囲を、数値ODEにおける連続段法にまで拡張すること。

提案手法

• 重み関数を用いてバッハー係数の直交多項式基底を定義するcsRK定式化を導入する。
• 数値解析で標準的である重み付き補間型求積則を用い、係数 $A_{\tau,\sigma}$, $B_\tau$, および $C_\tau$ を直接構築する。
• バッハー核を直交多項式の展開として表現する。例えば $A_{\tau,\sigma} = \frac{1}{2} + \sum_{i+j \text{ odd}} \alpha_{(i,j)} L_i(\tau)L_j(\sigma)$ であり、実数係数 $\alpha_{(i,j)}$ を持つ。
• 重み付き多項式（例：ジャコビ、ルジャンドル）の直交性および求積性の性質を用いて、安定性と精度を確保する。
• ガウス＝クリスティッフェル求積則の理論を応用し、最適な収束性を示す高次法を導出する。
• 本手法が特定の直交多項式族の選択に依存しないことを示し、広範な一般化が可能であることを明らかにする。

実験結果

リサーチクエスチョン

• RQ1 重み付き直交多項式を用いることで、複雑な非線形次数条件を解かずにcsRK法を構築可能か？
• RQ2 既存の求積理論をどのように活用すれば、所望の幾何的性質を持つcsRK法を設計できるか？
• RQ3 本フレームワークを用いて、ルジャンドル多項式を超えてシンプレクティック法や対称法のcsRK法を体系的に生成可能か？
• RQ4 重み関数がcsRK法のクラスを一般化する役割を果たすか？
• RQ5 多項式に基づくアプローチにより、高次幾何的積分法を構築可能か？

主な発見

• 次数条件に関連する多変数非線形代数系の解法を回避するため、直交多項式展開を用いることで、その必要性が排除される。
• 重み付き補間型求積理論により、csRK係数の直接的かつ体系的な構築が可能となり、設計プロセスが著しく簡素化される。
• ジャコビ多項式や他の重み付き直交多項式を用いて、シンプレクティック法および対称法の新規な族を構築可能である。
• 本フレームワークは一般性に富んでおり、ルジャンドル多項式に限定されないため、さまざまな多項式族への応用が可能である。
• 本アプローチは、将来的にエネルギー保存型や共役シンプレクティック法などの他の幾何的積分法の開発を支援する。
• 本手法は、スペクトル法とcsRK法との間に深い関係を示しており、高次精度および超収束の可能性を示唆している。