[論文レビュー] Is Symplectic-Energy-Momentum Integration Well-Posed?

Yosi Shibberu|ArXiv.org|Aug 4, 2006

Numerical methods for differential equations参考文献 19被引用数 3

ひとこと要約

この論文は、拡張された位相空間における離散時間ハミルトニアン（DTH）方程式を分析することで、シンプレクティック・エネルギー・モーメンタム（SEM）積分法の適切な定義（well-posedness）を検討する。非線形ハミルトニアン系、特に非線形振り子を含む系において、特定の位相空間の点で、任意に小さな時間刻みでもDTH方程式の解が存在しないことがあることが示され、SEM積分法がシンプレクティック性およびエネルギー保存性を有するが、常に適切に定義されているとは限らないという仮定に疑問を呈する。

ABSTRACT

We provide new existence and uniqueness results for the discrete-time Hamilton (DTH) equations of a symplectic-energy-momentum (SEM) integrator. In particular, we identify points in extended-phase space where the DTH equations of SEM integration have no solution for arbitrarily small time steps. We use the nonlinear pendulum to illustrate the main ideas.

研究の動機と目的

• シンプレクティック・エネルギー・モーメンタム（SEM）積分法が適切に定義されているかどうかを、離散時間ハミルトニアン（DTH）方程式の可解性を分析することで調査すること。
• DTH方程式が、任意に小さな時間刻みに対しても解を持たない、拡張された位相空間内の特定の点を同定すること。
• 線形系に限らない非線形ハミルトニアン系へと、DTHダイナミクスの存在および一意性に関する結果を拡張すること。
• 変分的定式化が理論的には魅力的であるが、その場合に解が得られない条件を明確にすること。

提案手法

• 時間を一般化座標およびその共共役運動量をもつものとして扱うことで、拡張された位相空間における変分原理を用いて離散時間ハミルトニアン（DTH）ダイナミクスを定式化する。
• 時間ステップの中点におけるハミルトニアンの制約を含む静的作用原理からDTH方程式を導出する。
• 作用とエネルギー制約に関連する関数 $ g(s, z_k) $ の振る舞いを解析することで、得られた非線形代数方程式の可解性を分析する。
• 中間値の定理および単調性の議論を用いて、時間刻みパラメータ $ s_k $ に対する解の存在および一意性を決定する。
• テイラー展開およびリプシッツ型推定から得られる境界を用いて、解が存在しない可能性のある位相空間内の領域を特定する。
• 非線形振り子を具体的な例として用い、特定のエネルギーおよび運動量条件のもとで解の存在が失敗することを示す。

実験結果

リサーチクエスチョン

• RQ1 SEM積分法のDTH方程式が、任意に小さな時間刻みに対しても解を持たない条件は何か？
• RQ2 時間刻みが無限に小さくても、拡張された位相空間内の特定の点でDTH方程式が適切に定義されないことがあるか？
• RQ3 DTH定式化におけるエネルギーおよび運動量制約は、離散方程式の可解性にどのように影響するか？
• RQ4 解の存在が失敗することは、$ \mathcal{H}_k $, $ \psi_k $, および $ \Lambda_k $ などの位相空間変数の観点から定量的に記述できるか？
• RQ5 解の存在が失敗することは、Geの定理を念頭に置いた場合、SEM積分法の理論的裏付けと矛盾するのか？

主な発見

• 特定の拡張された位相空間の点において、非線形ハミルトニアン系を含む系でも、任意に小さな時間刻みに対し、SEM積分法のDTH方程式の解が存在しないことがある。
• 非線形振り子では、$ \mathcal{H}_k / \psi_k > \frac{2}{3}(\psi_k / \psi_k')^2 $ である場合に解が存在しなくなり、$ (0, S_k) $ 内で解が存在しない閾値を示す。
• $ S_k \geq 6 $ かつ $ \mathcal{H}_k / \psi_k > \frac{2}{3}(\psi_k / \psi_k')^2 $ である場合、$ (0, S_k) $ 内に解が存在せず、特定の条件下では $ (0, 6/5) $ および $ (6/5, S_k) $ の両領域に解が共存する可能性がある。
• $ \mathcal{H}_k / \psi_k > \frac{1}{16}\Lambda_k^2(2 - S_k) $ かつ $ S_k < 6/5 $ のとき、$ (-S_k, 0) $ で解が存在しなくなり、非可解領域が生じる。
• 関数 $ g(s, z_k) $ は $ s $ に関して厳密に単調増加であり、解が存在する場合には一意性が保証される（Lemma 10 (ii) および (iii) による）。
• $ \mathcal{H}_k / \psi_k' < -\frac{1}{16}\Lambda_k^3 $ である場合、$ [-\Lambda_k, \Lambda_k] $ 内に解が存在せず、時間刻みパラメータにおける不適切な定義の領域が明確に示された。