[论文解读] Newton's method for computing periodic orbits of the planar three-body problem

Ivan Hristov, Radoslava Hristova+4 more|arXiv (Cornell University)|Nov 21, 2021

Spacecraft Dynamics and Control参考文献 12被引用 5

一句话总结

本文详细实现了牛顿法及其连续类比修正法（CANM），用于计算质量相等且角动量为零的平面三体问题中的周期轨道。通过采用高阶多精度泰勒级数积分并放宽返回接近度阈值，该方法识别出33个先前数据库中未包含的新拓扑家族周期轨道，展示了更强的收敛性与解发现能力。

ABSTRACT

In this paper we present in detail Newton's method and its modification, based on the Continuous analog of Newton's method for computing periodic orbits of the planar three-body problem. The linear system at each step of the method is formed by solving a system of ODEs with the multiple precision Taylor series method. We consider zero angular momentum symmetric initial configuration with parallel velocities, bodies with equal masses and relatively short periods. Taking candidates for the correction method with greater return proximity as usual and correcting with the modified Newton's method, allows us to find some new topological families that are not included in the database in [SCIENCE CHINA Physics, Mechanics & Astronomy 60.12 (2017)]

研究动机与目标

• 开发并详细阐述一种用于计算平面三体问题中周期轨道的稳健数值方法。
• 解决标准双精度方法在捕捉对初值敏感的周期解时的局限性。
• 通过放宽返回接近度阈值并采用改进的牛顿法以提升收敛性，从而改善解的发现能力。
• 通过识别超出现有数据库的新拓扑家族周期轨道来验证该方法的有效性。
• 为研究人员提供一个完整、可实现的框架，以便使用高阶多精度泰勒级数积分复现并拓展结果。

提出的方法

• 该方法采用牛顿迭代，通过高阶多精度泰勒级数积分求解变分方程，构建线性系统。
• 在 $v_x, v_y \in [0,1]$ 范围内采用分辨率为1024×1024的网格搜索算法，生成周期轨道的初始近似值。
• 基于放宽的返回接近度准则（$P(T) < 0.7$ 且在网格上为局部最小值）选择候选解，提高捕获近似周期解的可能性。
• 应用牛顿法的连续类比（CANM）作为校正技术，其收敛域大于经典牛顿法。
• 通过高精度积分数值强制实施周期性条件 $u(T) = u(0)$，以确保精度。
• 通过自由群中的拓扑词对轨道进行分类，其中 $L_f$ 表示词长，实现系统的目录化管理。

实验结果

研究问题

• RQ1 采用高精度积分的改进牛顿法是否能发现超越现有数据库的平面三体问题中新的周期轨道？
• RQ2 在存在数值敏感性的情况下，放宽返回接近度阈值对周期解发现率有何影响？
• RQ3 牛顿法的连续类比（CANM）在改善近似周期初值的收敛性方面发挥何种作用？
• RQ4 高阶多精度泰勒级数积分在多大程度上能提升周期轨道计算的可靠性和精度？
• RQ5 所提出的方法能否系统性地识别并分类对称零角动量三体问题中的新拓扑家族周期轨道？

主要发现

• 该方法成功计算出33个新的周期轨道，周期长达 $T \approx 1279$（归一化周期 $T^* \approx 1.28 \times 10^3$），且未出现在Liao1数据库中。
• 采用放宽的返回接近度阈值（$P(T) < 0.7$）显著增加了可用于校正的可行候选解数量，从而发现了此前遗漏的解。
• 基于CANM的改进牛顿法展现出更大的收敛域，即使在初值与周期性仅有中等接近度时也能实现收敛。
• 计算得到的轨道被划分为33个不同的拓扑家族，词长 $L_f$ 范围从8到52，表明其具有多样的动力学结构。
• 高阶多精度泰勒级数积分能够准确求解变分方程，并实现具有可控误差增长的长期积分。
• 结果证实了清洁数值模拟（CNS）原理的有效性，特别是当与CANM结合时，适用于如三体问题等复杂动力系统的研究。