[論文レビュー] Resonances and instabilities in symmetric multistep methods

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

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

ひとこと要約

本稿は、非線形惑星軌道積分に適用された対称多段法（特に高次のもの）において、これまでに無視されてきた共鳴および不安定性を特定している。これらの方法はストルマー法と比較して長期的なエネルギー誤差および位置誤差の挙動が優れているが、特定の刻み幅（例：1軌道あたり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.

研究の動機と目的

• 非線形軌道力学に用いられる対称多段法において、予期しないエネルギー誤差の急増の原因を調査すること。
• 線形安定性にのみ注目した先行研究が、なぜこれらの不安定性および共鳴を逃したのかを特定すること。
• これらの問題が存在するにもかかわらず、高次の対称多段法が長期的惑星積分に有効であるかどうかを評価すること。
• 重力的相互作用が共鳴および不安定性の発生頻度および深刻度に与える影響を評価すること。
• 現実的な多惑星系において、対称法とストルマー法の最大誤差対平均誤差の観点からの性能を比較すること。

提案手法

• 8次対称法SY8および8次ストルマー法を用いて、ケプラー的および惑星的軌道系の数値積分を実施。
• 1軌道あたりの刻み幅（例：54～81ステップ/軌道）を体系的に変化させ、共鳴および不安定な刻み幅を同定。
• 時間経過に伴うエネルギー誤差および経度誤差の分析により、指数的増大（不安定性）と線形増大（共鳴）を区別。
• 速度誤差のアーティファクトを最小限に抑えるために、多段方程式の和形式を用い、同時に不安定性の特徴を保持。
• 線形安定性理論を基盤として、非線形系における不安定性の起源を理解するための摂動解析を実施。
• 単一惑星（木星のみ）および二惑星（木星–土星）系の比較により、重力結合が共鳴行動に与える影響を評価。

実験結果

リサーチクエスチョン

• RQ1 線形安定性に属するにもかかわらず、特定の刻み幅で対称多段法が急激に大きなエネルギー誤差を示すのはなぜか？
• RQ2 惑星軌道のような非線形系に適用した際、対称多段法に共鳴および不安定性が生じる原因は何か？
• RQ3 惑星間の重力的相互作用が、共鳴および不安定性の発生頻度および深刻度に与える影響は何か？
• RQ4 最大誤差指標が、長期的惑星積分における対称多段法の実際の性能をどの程度誇張しているのか？
• RQ5 ピーク誤差のスパイクが存在するにもかかわらず、平均誤差の挙動が依然としてストルマー法よりも優れていると言えるか？

主な発見

• 対称多段法は、1軌道あたり60、90、120、150ステップなどの特定の刻み幅で指数的誤差増大（不安定性）を示し、周期性の線形区間内にあっても同様である。
• 1軌道あたり60ステップでは、エネルギー誤差が指数的に増大し、ピークで約0.25に達した後、約550軌道周期の周期で減衰し振動を示す。
• 共鳴（線形誤差増大を特徴とする）は、5、6、9の倍数である刻み幅（例：54、55、60、63、66、70、72、81）で発生し、木星–土星系では1軌道あたり60、63、72、81ステップで誤差スパイクが観測される。
• h ≈ 80（1軌道あたり54ステップ）付近の不安定性は二惑星系でも持続しており、重力結合によって解消されないため、刻み幅に対する本質的感度が示唆される。
• 共鳴付近では最大誤差が著しく大きいが、平均エネルギー誤差は最大誤差に対して数個のオーダー小さく、特に共鳴刻み幅から離れた領域では顕著である。
• 木星–土星系では、対称法の平均誤差は最大誤差に比べて著しく小さく、一方ストルマー法では平均誤差と最大誤差がほぼ同等であるため、時間平均では対称法が依然として優位であることが示された。