[论文解读] A comment on the arguments about the reliability and convergence of chaotic simulations

Shijun Liao|arXiv (Cornell University)|Jan 1, 2014

Quantum chaos and dynamical systems参考文献 23被引用 6

一句话总结

本文展示了通过结合高阶泰勒级数积分与超高压缩精度算术的纯净数值模拟（CNS）方法，能够可靠且收敛地模拟洛伦兹方程在长时间区间内的混沌解。通过使用1200个CPU和4180位精度，作者在[0, 10000]区间内实现了收敛且可靠的混沌解，证实了洛伦兹的观点：在任意有限时间范围内，数值近似可以收敛到真实的混沌解。

ABSTRACT

Yao and Hughes commented (Tellus-A, 60: 803 - 805, 2008) that "all chaotic responses are simply numerical noise and have nothing to do with the solutions of differential equations". However, using 1200 CPUs of the National Supercomputer TH-A1 and a parallel integral algorithm of the so-called "Clean Numerical Simulation" (CNS) based on the 3500th-order Taylor expansion and data in 4180-digit multiple precision, one can gain reliable, convergent chaotic solution of Lorenz equation in a rather long interval [0,10000]. This supports Lorenz's optimistic viewpoint (Tellus-A, 60: 806 - 807, 2008): "numerical approximations can converge to a chaotic true solution throughout any finite range of time".

研究动机与目标

• 解决长期以来关于混沌系统数值模拟是否反映真实解还是仅表现为数值噪声的争议。
• 证明利用先进数值方法可在长时间区间内可靠计算洛伦兹方程的混沌解。
• 验证洛伦兹的假设：在任意有限时间范围内，数值近似可收敛到真实的混沌解。
• 通过使用极端精度和高阶积分，在[0, 10000]区间内实现收敛解，为混沌系统提供数值基准。
• 通过严格且可复现的数值模拟，证实洛伦兹奇异吸引子的存在性和可计算性。

提出的方法

• 采用纯净数值模拟（CNS）方法，结合高阶泰勒级数积分（最高至3500阶）与多精度算术（4180位精度）。
• 在国家超算TH-A1上实现并行积分算法，使用1200个CPU以加速计算并管理内存密集型操作。
• 通过提高泰勒级数阶数（M）和数字精度（N），CNS方法同时减少截断误差与舍入误差，确保在M > M*且N > N*时实现收敛，其中M*和N*为给定时间区间的临界值。
• 通过证明在M和N逐步增加时模拟结果保持不可区分，验证了收敛性，表明数值误差可忽略不计。
• 通过与其他方法（如400位精度的200阶有限元法）交叉比对，验证了CNS方法的可靠性。
• CNS方法应用于洛伦兹方程，固定时间步长Δt = 0.01，经验确定长期模拟的收敛阈值为M* ≈ T/3和N* ≈ 2T/5。

实验结果

研究问题

• RQ1 尽管对初值敏感且存在数值误差，混沌系统（如洛伦兹方程）的数值模拟是否能在长时间有限时间区间内收敛到真实解？
• RQ2 通过使用高阶泰勒级数和超高压缩精度算术，能在多大程度上抑制混沌模拟中的截断误差与舍入误差？
• RQ3 纯净数值模拟（CNS）方法在不同数值方法和硬件平台上是否能产生一致且可靠的结果？
• RQ4 洛伦兹奇异吸引子能否在长时间模拟中被准确计算并验证为稳定且收敛的结构？
• RQ5 是否可以通过使用极端计算资源在[0, 10000]区间内实现收敛解，从而为混沌动力学建立数值基准？

主要发现

• 成功使用1200个CPU和4180位精度，在[0, 10000]区间内计算出洛伦兹方程的可靠且收敛的混沌解。
• CNS方法采用3500阶泰勒级数和4180位精度实现收敛，证实了在任意有限时间范围内，数值近似可收敛到真实的混沌解。
• 通过使用3600阶泰勒级数和4515位精度的独立验证，结果一致，进一步确认了结果的可靠性。
• 通过与400位精度的200阶有限元法进行交叉比对，结果在[0, 1000]区间内保持一致，验证了CNS方法的有效性。
• 本研究证实了洛伦兹的乐观观点：即使在极端敏感条件下，数值近似也能在整个有限时间区间内收敛到真实的混沌解。
• 本工作首次报告在如此长的时间区间内实现收敛且可靠的洛伦兹方程混沌解，为混沌动力学模拟树立了新基准。