[论文解读] Is Symplectic-Energy-Momentum Integration Well-Posed?

Yosi Shibberu|ArXiv.org|Aug 4, 2006

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

一句话总结

本文通過分析擴展相空間中的離散時間哈密頓（DTH）方程，研究了辛-能量-動量（SEM）積分器的適定性。結果表明，對於非線性哈密頓系統（包括非線性單擺），在相空間的某些點上，即使時間步長極小，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.

研究动机与目标

• 通過分析離散時間哈密頓（DTH）方程的可解性，研究辛-能量-動量（SEM）積分是否適定。
• 識別擴展相空間中DTH方程即使在時間步長極小時也無解的特定點。
• 將DTH動力學的存在性與唯一性結果從線性系統推廣至非線性哈密頓系統。
• 明確說明變分形式的SEM積分在何種條件下無法產生解，儘管其理論吸引力強。

提出的方法

• 在擴展相空間中，利用變分原理構造離散時間哈密頓（DTH）動力學，將時間視為廣義坐標並引入其共軛動量。
• 從帶有時間步長中點哈密頓約束的駐值作用量原理推導DTH方程。
• 通過研究與作用量和能量約束相關的函數 $ g(s, z_k) $ 的行為，分析所得非線性代數方程的可解性。
• 應用中值定理與單調性論證，確定時間步長參數 $ s_k $ 的解的存在性與唯一性。
• 利用泰勒展開與利普希茨型估計推導的界，確立相空間中解可能不存在的區域。
• 以非線性單擺為具體例子，說明在特定能量與動量條件下，解的存在性如何失效。

实验结果

研究问题

• RQ1 在何種條件下，SEM積分的DTH方程在任意小時間步長下無解？
• RQ2 是否存在擴展相空間中的點，使得DTH方程即使在時間步長無限小時也處於病態？
• RQ3 DTH形式中的能量與動量約束如何影響離散方程的可解性？
• RQ4 能否以相空間變量 $ \mathcal{H}_k $、$ \psi_k $ 和 $ \Lambda_k $ 等形式量化解不存在的情況？
• RQ5 解的不存在是否與SEM積分器的理論依據相矛盾，特別是考慮Ge定理的情況？

主要发现

• 在擴展相空間的某些點上，即使時間步長極小，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 $ 上嚴格單調，根據引理10 (ii) 和 (iii) 可確保解存在時的唯一性。
• 當 $ \mathcal{H}_k / \psi_k' < -\frac{1}{16}\Lambda_k^3 $ 時，$ [-\Lambda_k, \Lambda_k] $ 內無解，顯示時間步長參數中存在一個病態區間。