[论文解读] Resonances and instabilities in symmetric multistep methods

Gerald D. Quinlan|arXiv (Cornell University)|Jan 12, 1999

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

一句话总结

本文識別出在非線性行星軌道積分中應用對稱多步法（特別是高階方法）時，先前被忽略的共振與不穩定現象。雖然這些方法在長期能量與位置誤差表現上優於 Störmer 方法，但在特定步長（例如每軌道 60、90、120 步）下會出現指數型誤差增長，即使處於線性週期性區間內，也可能嚴重影響準確性。

ABSTRACT

The symmetric multistep methods developed by Quinlan and Tremaine (1990) are shown to suffer from resonances and instabilities at special stepsizes when used to integrate nonlinear equations. This property of symmetric multistep methods was missed in previous studies that considered only the linear stability of the methods. The resonances and instabilities are worse for high-order methods than for low-order methods, and the number of bad stepsizes increases with the number frequencies present in the solution. Symmetric methods are still recommended for some problems, including long-term integrations of planetary orbits, but the high-order methods must be used with caution.

研究动机与目标

• 探討在非線性軌道動力學中使用對稱多步法時，意外出現能量誤差尖峰的原因。
• 確定為何先前僅著眼於線性穩定性的研究未能察覺這些不穩定性與共振現象。
• 評估高階對稱多步法在存在此類問題的情況下，是否仍適用於長期行星積分。
• 評估行星之間引力相互作用對共振與不穩定現象發生頻率與嚴重程度的影響。
• 比較對稱方法與 Störmer 方法在最大誤差與平均誤差方面的表現，特別是在真實的多行星系統中。

提出的方法

• 使用 8 階對稱方法 SY8 和 8 階 Störmer 方法對開普勒系統與行星軌道系統進行數值積分。
• 系統性地改變每軌道的步長（例如 54–81 步/軌道），以識別共振與不穩定的步長。
• 分析能量與經度誤差隨時間的變化，以區分指數型增長（不穩定性）與線性增長（共振）。
• 使用多步方程的求和形式計算速度，以最小化速度誤差的人為效應，同時保留不穩定特性。
• 進行擾動分析，以理解非線性系統中不穩定性的來源，並根據線性穩定性理論提供指導。
• 比較單行星系統（僅木星）與雙行星系統（木星–土星）以評估引力耦合對共振行為的影響。

实验结果

研究问题

• RQ1 為何對稱多步法在線性穩定的情況下，於特定步長仍會出現突然且巨大的能量誤差？
• RQ2 當對稱多步法應用於木星軌道等非線性系統時，共振與不穩定性的來源為何？
• RQ3 行星之間的引力相互作用如何影響共振與不穩定現象的發生頻率與嚴重程度？
• RQ4 最大誤差指標在多大程度上高估了對稱多步法在長期行星積分中的實際表現？
• RQ5 儘管存在峰值誤差尖峰，對稱方法的平均誤差表現是否仍使其優於 Störmer 方法？

主要发现

• 對稱多步法在特定步長（如每軌道 60、90、120、150 步）下會出現指數型誤差增長（不穩定性），即使處於線性週期性區間內。
• 在每軌道 60 步時，能量誤差會指數型增長至約 0.25 的峰值，隨後衰減並以約 550 個軌道週期待續振盪。
• 共振（特徵為線性誤差增長）出現在 5、6 或 9 的倍數步長處（例如 54、55、60、63、66、70、72、81），在木星–土星系統中，60、63、72 和 81 步/軌道處出現明顯誤差尖峰。
• 在 h ≈ 80（即 54 步/軌道）附近的不穩定性在雙行星系統中依然存在，且未因引力耦合而消除，顯示其對步長具有本質上的敏感性。
• 儘管在共振點處最大誤差極大，但對稱方法的平均能量誤差遠小於最大誤差，特別是在非共振步長下，可小數個數量級。
• 在木星–土星系統中，對稱方法的平均誤差遠小於其最大誤差，而 Störmer 方法的平均與最大誤差則相近，顯示即使考慮峰值誤差，對稱方法在時間平均下仍表現更優。