[论文解读] Numerical validation of blow-up solutions with quasi-homogeneous compactifications

Kaname Matsue, Akitoshi Takayasu|arXiv (Cornell University)|Jul 19, 2017

Mathematical Biology Tumor Growth参考文献 9被引用 7

一句话总结

本文提出了一种针对在无穷远处具有渐近拟齐次性的有限维ODE中爆破解的严格数值验证方法，通过一种新颖的拟抛物型紧化方法，将发散轨迹转化为终止于紧化无穷远处平衡点的全局轨道。通过结合严格数值积分、李雅普诺夫函数验证以及时间变量去奇化技术，该方法计算了爆破时间的验证上下界及解轮廓的包围区域，并实现了保证误差控制的计算。

ABSTRACT

We provide a numerical validation method of blow-up solutions for finite dimensional vector fields admitting asymptotic quasi-homogeneity at infinity. Our methodology is based on quasi-homogeneous compactifications containing a new compactification, which shall be called a quasi-parabolic compactification. Divergent solutions including blow-up solutions then correspond to global trajectories of associated vector fields with appropriate time-variable transformation tending to equilibria on invariant manifolds representing infinity. We combine standard methodology of rigorous numerical integration of differential equations with Lyapunov function validations around equilibria corresponding to divergent directions, which yields rigorous upper and lower bounds of blow-up times as well as rigorous profile enclosures of blow-up solutions.

研究动机与目标

• 开发一种针对在无穷远处具有渐近拟齐次性的ODE系统中爆破解的严格数值方法。
• 将先前基于紧化的方法扩展至更广泛的向量场类别，超越渐近齐次系统。
• 构建一类新型紧化方法——具体为拟抛物型紧化，以保持拟齐次向量场的尺度结构。
• 利用区间算术与李雅普诺夫函数技术，实现爆破时间与解轮廓的验证计算。
• 提供一个适用于Keller-Segel模型等系统的方法框架，其中爆破在物理与数学上均具有重要意义。

提出的方法

• 该方法采用拟抛物型紧化，将相空间嵌入紧流形中，将无穷远点映射为边界流形。
• 通过时间变量变换（去奇化）重参数化系统，使得爆破解对应于收敛至紧化边界上平衡点的全局轨迹。
• 利用区间算术与仿射算术执行严格数值积分，以在有限时间区间内有界解轨迹。
• 在边界平衡点邻域内验证李雅普诺夫函数，以估计收敛速率并有界化到达无穷远的时间。
• 通过积分李雅普诺夫函数衰减率所得的误差估计，计算爆破时间的验证上下界。
• 使用变换与向量场分量的多项式逼近，以有界紧化系统与原始动力学之间的偏差。

实验结果

研究问题

• RQ1 能否使用基于紧化的方法严格验证拟齐次ODE系统中的爆破解？
• RQ2 如何构造合适的紧化方法，以保持拟齐定向量场的渐近尺度结构？
• RQ3 李雅普诺夫函数在具有验证误差界的情况下估计爆破时间中起什么作用？
• RQ4 该方法能否扩展至齐次系统之外的更一般的拟齐次动力学？
• RQ5 如何将验证数值积分与区间算术结合去奇化技术，以有界爆破轮廓与时间？

主要发现

• 本文提出了一类新型紧化方法——拟抛物型紧化，其推广了适用于拟齐次系统的庞加莱型紧化。
• 爆破解被表征为去奇化向量场的全局轨迹，收敛至紧化边界上的平衡点，从而实现严格验证。
• 通过积分李雅普诺夫函数衰减估计，计算了爆破时间的验证上下界，误差界由区间算术导出。
• 通过在边界平衡点邻域内有界解轨迹，实现了爆破解的验证轮廓包围。
• 当 $k=1$ 时，爆破时间上界被显式计算为 $t_{\text{max}} \leq t_N + \frac{1}{c_{\tilde{N}}} \left\{ 2\|v_1\|c_1^{-1/2}L^{1/2} + \sum_{j=2}^{\max\{2\beta_i\}} \frac{2}{j}\|v_j\|_{\infty}c_1^{j/2-1}L^{j/2} \right\}$，其中 $L$ 为李雅普诺夫函数值。
• 该方法适用于具有非齐次尺度的系统，如 $h(u,v) = u^2 - v$ 的系统，这些系统无法通过标准齐次紧化方法处理。