1993 Fiscal Year Annual Research Report
ハミルトン力学系に対する正準変換型の数値積分公式の開発と応用
| Project/Area Number |
05836036
|
|---|
| Research Institution | National Astronomical Observatory of Japan |
|---|
Principal Investigator |
吉田 春夫 国立天文台, 位置天文・天体力学研究系, 助教授 (70220663)
|
|---|
| Keywords | ハミルトン力学系 / 微分方程式の数値解法 |
|---|
| Research Abstract |
N体問題に代表されるHamilton系H=T(p)+V(q)を数値積分する最も原始的、かつ単純なスキームはオイラー法であるが、オイラー法では相空間の面積要素(symplectic2-form)が一般に保存されず、t=0の(q,p)からt=tauの(q′,p′)への変換が正準変換とはならない。そして、保存すべきエネルギー値の誤差がsecularに増大し、長時間の数値計算の結果の信頼性を失わせる。この事情は局所的な誤差を小さくするRunge‐Kutta型の積分法、および従来の多段型の積分法においても同様である。正準変換型の積分公式(symplectic integrator)は(q,p)から(q′,p′)への変換が厳密に正準変換となるように設計された積分法である。正準変換型の積分法に対して今年度我々は次の点を明らかにした。
1:Symplectic Integratorの応用において一つのネックとなっていた可変時間ステップにおいて、Hut et al(1993)によって示唆されてように、時間ステップをスキーム全体の可逆性を保つように決定すればエネルギーの誤差に長年的な誤差が現れないことを確認した。
2:剛体の自由回転を記述するEulerの方程式などのように保存力学系であるにも係わらず通常の意味でのハミルトン系でない力学系が存在する。そのような系に対するSymplectic Integratorの一般化はLie‐Poisson Integratorと呼ばれるが、この積分法に対してやはりエネルギーの誤差が永年的に増大しないことを示した。
|
