[論文レビュー] The Adams-Bashforth-Moulton Integration Methods Generalized to an Adaptive Grid

Allan Hayes|arXiv (Cornell University)|Apr 16, 2011

Geophysics and Gravity Measurements被引用数 3

ひとこと要約

本稿では、アダムス＝バッシュフォースト＝マウルトン予測補正法を適応的グリッドに一般化し、所望の誤差許容誤差を維持するための動的ステップサイズ調整を可能にした。この手法は、1ステップあたり僅か2回の導出評価で高い精度を達成し、収束性を示し、23ステップで1%の精度でトールマン＝オッペンハイマー＝ヴォルコフの中性子星質量限界（0.71太陽質量）を再現した。

ABSTRACT

We present a generalization of the Adams-Bashforth-Moulton predictor-corrector numerical integration methods to an adaptive grid. The step size may be chosen dynamically in order to maintain a desired relative magnitude of error in each step. We demonstrate that the methods remain convergent to the expected degree, and apply various methods to the famous problem of determining the maximum possible mass of a neutron star supported by pure fermionic exclusion pressure. We reproduce the Tolman-Oppenheimer-Volkoff result of 0.71 solar masses using only 23 integration steps, and reproducing both mass and radius within 1% requires 27. We also present various optimizations and features of our implementation.

研究の動機と目的

• 導出評価が高価な問題において、効率を向上させるために、古典的アダムス＝バッシュフォースト＝マウルトン多段法を適応的グリッドに一般化すること。
• 予測補正ペアからの誤差推定を用いて、各ステップあたりの目標相対誤差許容誤差を維持する動的ステップサイズ制御を可能にすること。
• 収束次数にかかわらず、1ステップあたり2回の導出評価で済むことにより、ルンゲ＝クッタ法と比較して計算コストを低減すること。
• トールマン＝オッペンハイマー＝ヴォルコフ方程式を用いた中性子星最大質量の決定という物理的に有意義な問題に、この手法を実証すること。
• 高次収束性と安定性が、非一様な適応的グリッドでも維持可能であることを示し、特に中性子星表面付近の急変化領域において、ルンゲの現象を避けるために次数チューニングが有効であることを示すこと。

提案手法

• N次の陽的アダムス＝バッシュフォースト（AB）予測子と、N+1次の陰的アダムス＝マウルトン（AM）補正子を組み合わせ、予測補正ペアを形成する。
• 各積分ステップにおいて、AB予測子が以前の導出値を用いて初期推定値を提供し、その後AM法により高次精度に補正する。
• ABとAMの結果の差が誤差推定値として用いられ、所望の相対誤差許容誤差を維持するための動的ステップサイズ調整に使用される。
• ABおよびAM法の重み係数は、各ステップで再計算され、変動するステップサイズに対応するが、導出評価に比べて計算コストは低く抑えられる。
• 陰的であるため、コーラント＝フリードリッヒス＝レヴィ（CFL）条件を回避でき、特に滑らかな領域ではより大きな安定ステップが可能になる。
• ルンゲの現象を回避するため、問題の挙動に応じて収束次数をチューニングする。特に、中性子星表面付近のような急変化領域において有効である。

実験結果

リサーチクエスチョン

• RQ1 アダムス＝バッシュフォースト＝マウルトン予測補正法は、収束次数と安定性を保ちつつ、適応的グリッドに一般化可能か？
• RQ2 誤差推定に基づく動的ステップサイズ調整は、導出評価が高価な問題においても精度と効率を維持できるか？
• RQ3 剛性のある天体物理学的問題において、固定グリッドまたはルンゲ＝クッタ法と比較して、適応的ABM法は、積分ステップ数と導出評価回数の点でどのように優れるか？
• RQ4 トールマン＝オッペンハイマー＝ヴォルコフ方程式を解く際、最適な効率を得るために、どの収束次数と誤差許容誤差が最適か？
• RQ5 ステップサイズは、中性子星表面付近の圧力低下のような急勾配を捉えるためにどのように変化するか？

主な発見

• 適応的グリッド設定下でも、固定グリッド・固定ステップ、および適応的グリッドのすべての設定で、期待される収束次数が達成された。
• 4-5次法と10⁻²の許容誤差を用いて、23ステップでトールマン＝オッペンハイマー＝ヴォルコフの中性子星質量限界（0.71017188太陽質量）を1%の精度で再現した。
• 質量と半径の両方を1%の精度で達成するには、4-5次法と10⁻²の許容誤差を用いて27ステップが必要だった。
• より高い精度（質量で10⁻⁸、半径で10⁻⁵）を達成するには、9-10次法と10⁻⁵の許容誤差を用いて131ステップが必要だった。
• ステップサイズは初期に指数関数的に増加し、核部では plateau に達し、中性子星表面付近で急激に減少する傾向を示した。これは圧力プロファイルの可視化によって確認された。
• AM補正項は、積分の大部分で安定しており、表面付近でのみ顕著に増加しており、局所的な誤差感受性が高まっていることを示している。