[论文解读] The Adams-Bashforth-Moulton Integration Methods Generalized to an Adaptive Grid

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

Geophysics and Gravity Measurements被引用 3

一句话总结

该论文将亚当斯-巴什福斯-穆尔顿预测-校正方法推广至自适应网格，实现动态步长调整以维持期望的误差容限。该方法仅需每步两次导数计算即可实现高精度，展示了收敛性，并在仅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.

研究动机与目标

• 将经典的亚当斯-巴什福斯-穆尔顿多步法推广至自适应网格，以提升在导数计算成本较高的问题中的效率。
• 通过预测-校正对的误差估计实现动态步长控制，维持每步的目标相对误差容限。
• 通过每步仅需两次导数计算（无论收敛阶数如何）相比龙格-库塔方法降低计算成本。
• 在物理上具有重要意义的问题上验证该方法：通过托尔曼-奥本海默-沃尔科夫方程确定中子星最大质量。
• 证明即使在自适应、非均匀网格下，高阶收敛性和稳定性仍可保持，通过阶数调节避免类似龙格现象的问题。

提出的方法

• 将N阶显式亚当斯-巴什福斯（AB）预测器与N+1阶隐式亚当斯-穆尔顿（AM）校正器组合，形成预测-校正对。
• 在每个积分步骤中，AB预测器利用先前的导数值提供初始估计，随后由AM方法进行校正以实现更高阶精度。
• AB与AM结果之间的差异提供误差估计，用于动态调整步长，以维持期望的相对误差容限。
• 在每个步骤重新计算AB与AM方法的加权因子，以适应可变步长，但计算成本相比导数计算仍较低。
• 通过隐式特性避免满足柯朗-弗里德里希斯-列维（CFL）条件，从而在光滑区域实现更大的稳定步长。
• 为防止龙格现象，根据问题行为（特别是中子星表面等快速变化区域）动态调节收敛阶数。

实验结果

研究问题

• RQ1 亚当斯-巴什福斯-穆尔顿预测-校正方法能否推广至自适应网格，同时保持其收敛阶数和稳定性？
• RQ2 基于误差估计的动态步长调整是否能在导数计算成本较高的问题中维持精度与效率？
• RQ3 与传统固定网格或龙格-库塔方法相比，自适应ABM方法在求解刚性天体物理问题时，其积分步数与导数计算次数有何差异？
• RQ4 求解托尔曼-奥本海默-沃尔科夫方程时，何种收敛阶数与误差容限可实现最优效率？
• RQ5 步长如何自适应调整以捕捉如中子星表面附近压力急剧下降等陡峭梯度？

主要发现

• 自适应ABM方法在固定网格、固定步长及自适应网格配置下均实现了预期的收敛阶数。
• 使用4-5阶方法与10⁻²容限，仅用23步即以1%精度重现了托尔曼-奥本海默-沃尔科夫中子星质量极限（0.71017188倍太阳质量）。
• 在质量与半径均达到1%精度时，需27步，使用4-5阶方法与10⁻²容限。
• 对于更高精度要求（质量精度10⁻⁸，半径精度10⁻⁵），需131步，使用9-10阶方法与10⁻⁵容限。
• 步长初始呈指数增长，随后在核心区域趋于平稳，而在中子星表面附近因压力急剧下降而显著减小，该现象通过压力分布图可视化得到验证。
• AM校正项在整个积分过程中保持稳定，仅在表面附近出现显著增加，表明该区域具有更高的局部误差敏感性。