[论文解读] Rigorous numerics of finite-time singularities in dynamical systems - methodology and applications

Kaname Matsue|arXiv (Cornell University)|Nov 6, 2017

Nonlinear Dynamics and Pattern Formation参考文献 21被引用 4

一句话总结

本文提出了一种严格的数值方法，用于验证动力系统中的有限时间奇异性，包括有限时间灭绝、紧支集解（compactons）和奇异拐点解（singular canards），通过时间尺度去奇化和稳定流形上的李雅普诺夫函数验证实现。该方法利用区间算法和验证积分，实现了对灭绝时间、支集大小和穿越时间等关键量的带保证误差界的验证计算。

ABSTRACT

This paper aims at providing rigorous numerical computation procedure for finite-time singularities in dynamical systems. Combination of time-scale desingularization as well as Lyapunov functions validation on stable manifolds of invariant sets for desingularized vector fields with standard integration procedure for ordinary differential equations give us validated trajectories of dynamical systems involving finite-time singularities. Our focus includes finite-time extinction, traveling wave solutions with half-line or compact support, and singular canards in fast-slow systems, including rigorous validations of enclosures of extinction, finite-passage times or size of supports for compactons. Such validated solutions lead to a plenty of composite wave solutions for degenerate parabolic equations, for example, with concrete information of profiles and evolutions. The present procedure also provides a universal aspect of finite-time singularities with rigorous numerics, combining with rigorous numerics of blow-up solutions in preceding works.

研究动机与目标

• 开发一种通用的、严格的数值框架，用于验证动力系统中的有限时间奇异性，包括灭绝、紧支集解和拐点解。
• 解决验证此类奇异性是否发生及其发生时间的根本性挑战，并提供保证的误差界。
• 通过结合时间尺度去奇化与验证积分和李雅普诺夫函数分析，将严格数值方法扩展至奇异常微分方程。
• 为灭绝时间、支集大小和行波解中的有限穿越时间等关键量提供具体且经过验证的包围区间。
• 在统一的计算机制下统一处理多种有限时间奇异性，与此前关于爆破解的研究相联系。

提出的方法

• 通过 $ \frac{dt}{d\tau} = T(x(t)) $ 实现时间尺度去奇化，将奇异常微分方程转化为适合标准数值积分的正则系统。
• 应用区间算法和验证积分，在相空间中计算轨迹的严格包围区间，保证数学上的误差界。
• 在去奇化向量场的稳定流形上使用李雅普诺夫函数验证，严格确认收敛至退化状态。
• 通过以参数化平衡点为中心的仿射变换构造隔离块，并利用预测-校正优化减少块的大小。
• 利用锥条件和覆盖关系对连接轨道进行拓扑验证，例如在去奇化系统中存在同宿或异宿轨道。
• 结合基于隐函数定理的参数延续与块构造方法，确保在参数依赖系统中具有鲁棒性和高精度。

实验结果

研究问题

• RQ1 能否在动力系统中严格验证有限时间灭绝，包括精确的灭绝时间？
• RQ2 如何通过数值计算以保证的界验证紧支集解的支集大小？
• RQ3 在快速-慢系统中，使用严格数值方法检测和验证奇异拐点解的机制是什么？
• RQ4 能否建立一个统一的计算框架，用于处理包括爆破、淬灭和灭绝在内的多种有限时间奇异性？
• RQ5 如何将时间尺度去奇化与验证积分和李雅普诺夫函数结合，以确保轨迹包围区间的正确性？

主要发现

• 该方法通过计算灭绝时间的验证包围区间并保证误差界，成功验证了有限时间灭绝。
• 通过在去奇化系统中进行验证轨迹积分，严格计算并包围了紧支集解的支集大小。
• 利用所提出的框架，对靠近退化状态的轨迹的有限穿越时间实现了数学上的确定性计算。
• 预测-校正方法相比标准中心化方法显著减小了隔离块的大小，提升了精度和效率。
• 该框架在单一严格的数值机制下统一处理了有限时间奇异性，包括紧支集解、拐点解和灭绝现象。
• 该方法通过引入去奇化和稳定流形分析，将先前针对爆破解的严格数值方法扩展至非爆破奇异性。