[論文レビュー] Rigorous numerics of finite-time singularities in dynamical systems - methodology and applications

Kaname Matsue|arXiv (Cornell University)|Nov 6, 2017

Nonlinear Dynamics and Pattern Formation参考文献 21被引用数 4

ひとこと要約

本稿では、時間スケールの特異性除去と安定多様体上のリャプノフ関数による検証を用いて、有限時間における特異解（消滅、コンパクトン、特異的カナード）を検証する厳密な数値的手法を提示する。この手法により、区間演算と検証済み積分を用いて、誤差の保証された境界を伴う消滅時刻、支持長、通過時刻の検証計算が可能となる。

ABSTRACT

This paper aims at providing rigorous numerical computation procedure for finite-time singularities in dynamical systems. Combination of time-scale desingularization as well as Lyapunov functions validation on stable manifolds of invariant sets for desingularized vector fields with standard integration procedure for ordinary differential equations give us validated trajectories of dynamical systems involving finite-time singularities. Our focus includes finite-time extinction, traveling wave solutions with half-line or compact support, and singular canards in fast-slow systems, including rigorous validations of enclosures of extinction, finite-passage times or size of supports for compactons. Such validated solutions lead to a plenty of composite wave solutions for degenerate parabolic equations, for example, with concrete information of profiles and evolutions. The present procedure also provides a universal aspect of finite-time singularities with rigorous numerics, combining with rigorous numerics of blow-up solutions in preceding works.

研究の動機と目的

• 力学系における有限時間特異解（消滅、コンパクトン、カナードを含む）を検証する普遍的で厳密な数値フレームワークの構築を目的とする。
• このような特異解が発生するかどうか、そしていつ発生するかを、保証された誤差境界を伴って検証するという根本的課題に取り組む。
• 時間スケールの特異性除去と検証済み積分、およびリャプノフ関数解析を組み合わせることで、特異な常微分方程式に対する厳密な数値解析を拡張する。
• 消滅時刻、支持長、および移動波解における有限通過時刻といった主要な量の具体的な検証された包含領域を提供する。
• 爆発解に関する先行研究に接続しながら、多様な有限時間特異解を同一の計算的メカニズムで統一的に取り扱う。

提案手法

• 特異な常微分方程式を標準的な数値積分が可能な正則系に変換するため、$ \frac{dt}{d\tau} = T(x(t)) $ を用いた時間スケールの特異性除去を採用する。
• 区間演算と検証済み積分を適用し、位相空間内での軌道の数学的に保証された誤差境界付き包含領域を計算する。
• 特異性除去されたベクトル場の安定多様体上でのリャプノフ関数による検証を用いて、退化状態への収束を厳密に確認する。
• パラメータ化された平衡点を中心としたアフィン変換を用いて隔離ブロックを構築し、予測子・修正子のリファインメントを用いてブロックサイズを低減する。
• コーン条件とカバー関係を用いて、特異性除去された系における接続軌道（ホモクライニックまたはヘテロクライニック軌道など）を位相的に検証する。
• implicit function theoremに基づくパrametric継続法とブロック構築を組み合わせることで、パラメータ依存系におけるロバスト性と精度を確保する。

実験結果

リサーチクエスチョン

• RQ1 力学系における有限時間の消滅が、正確な消滅時刻とともに厳密に検証可能か？
• RQ2 数値計算を用いて、保証された境界を伴ってコンパクトン解の支持長を検証できるか？
• RQ3 高速・遅速系における特異的カナード解を、厳密な数値解析で検出・検証するメカニズムは何か？
• RQ4 爆発、消滅、クエンチングを含む多様な有限時間特異解を統一した計算フレームワークで扱えるか？
• RQ5 時間スケールの特異性除去を検証済み積分とリャプノフ関数とどのように組み合わせ、軌道の包含領域の正しさを保証できるか？

主な発見

• 本手法は、保証された誤差境界を伴う消滅時刻の検証された包含領域を計算することで、有限時間の消滅を成功裏に検証した。
• コンパクトン解の支持長は、特異性除去された系における検証済み軌道積分を用いて厳密に計算され、包含された。
• 退化状態付近の軌道の有限通過時刻は、提案されたフレームワークを用いて数学的に確実に計算された。
• 予測子・修正子アプローチにより、標準的な中心化手法と比較して隔離ブロックのサイズが低減され、精度と効率が向上した。
• 本フレームワークは、コンパクトン、カナード、消滅といった有限時間特異解の取り扱いを、同一の厳密な数値的メカニズムで統一した。
• 本手法は、特異性除去と安定多様体解析を組み合わせることで、爆発解に関する先行研究の厳密な数値解析を、非爆発的特異解へと拡張した。