[论文解读] User guide on Hopf bifurcation and time periodic orbits with pde2path

Hannes Uecker|arXiv (Cornell University)|Aug 2, 2019

Nonlinear Dynamics and Pattern Formation参考文献 29被引用 6

一句话总结

本文为 pde2path（一种使用有限元法求解PDE的MATLAB工具箱）提供了关于霍普夫分支分析与周期轨道延续的全面用户指南。它详细介绍了检测和追踪时间周期解的数值方法，包括对称性感知的延续、通过弗洛quet乘子检测分岔，以及针对具有O(2)对称性、质量守恒和相位约束的系统而设计的专用例程，从而实现了PDE中霍普夫分支及次级分岔的稳健计算。

ABSTRACT

We explain the setup for using the pde2path libraries for Hopf bifurcation and continuation of branches of periodic orbits and give implementation details of the associated demo directories. See [Uecker, Comm. in Comp. Phys., 2019] for a description of the basic algorithms and the mathematical background of the examples. Additionally we explain the treatment of Hopf bifurcations in systems with continuous symmetries, including the continuation of traveling waves and rotating waves in O(2) equivariant systems as relative equilibria, the continuation of Hopf bifurcation points via extended systems, and some simple setups for the bifurcation from periodic orbits associated to critical Floquet multipliers going through +-1.

研究动机与目标

• 为使用pde2path检测和延续具有时间周期解的PDE中的霍普夫分岔，提供实用且以实现为重点的指南。
• 扩展pde2path的功能，以处理具有连续对称性的系统，如O(2)等变性、质量守恒和相位不变性，这些系统需要特殊的相位条件。
• 实现霍普夫分岔点的延续及从周期轨道出发的分支切换，包括周期加倍和霍普夫-劈尖分岔。
• 支持非耗散和不适定PDE（例如带有反向扩散项的最优控制问题），在标准初值求解器失效的情况下，延续方法仍具可行性。
• 记录并标准化用于周期轨道计算的高级例程，包括弧长延续、网格自适应以及事后乘子计算。

提出的方法

• 通过有限元空间半离散化将PDE转化为大规模ODE系统，利用pde2path内置求解器实现稳态解和周期轨道的延续。
• 采用方法线法结合隐式时间积分与弧长延续，追踪周期轨道，使用谐波平衡法和TOM（时间延续法）实现高效计算。
• 通过求解包含特征值条件和相位约束的增广系统，构建扩展系统以定位霍普夫分岔点。
• 应用弗洛quet理论计算单变量矩阵及其特征值（弗洛quet乘子），并提供直接计算与事后计算两种实现方式。
• 引入相位条件（如积分型或配点型）以解决对称系统中的相位模糊性，从而实现相对平衡态和周期轨道的稳定延续。
• 采用预测-校正格式，在扩展系统上进行牛顿迭代，包括使用扩展系统的雅可比矩阵以及时间方向的自适应网格加密。

实验结果

研究问题

• RQ1 如何利用有限元法在具有复杂时空动力学的PDE中可靠地检测和延续霍普夫分岔？
• RQ2 在具有连续对称性的系统（如平移或旋转不变性）中，延续周期轨道所需的数值策略是什么？
• RQ3 如何利用弗洛quet乘子分析检测和定位周期轨道上的分岔（如周期加倍和霍普夫-劈尖分岔）？
• RQ4 在处理非耗散或不适定PDE（如带有反向扩散项的最优控制问题）时，延续框架需要进行哪些修改？
• RQ5 如何正确地制定和实现相位条件，以在对称系统中稳定周期轨道的延续？

主要发现

• pde2path框架成功实现了在1D、2D和3D PDE中检测和延续霍普夫分岔，包括复Ginzburg-Landau方程和Brusselator系统。
• 对于具有O(2)对称性的系统，该方法能够区分行波与驻波，并正确计算相对平衡态和相对周期轨道。
• 相位条件的使用使得在质量守恒和空间平移不变系统（如Kuramoto-Sivashinsky方程和FHN型模型）中实现周期轨道的稳定延续成为可能。
• 通过利用规范型系数的预测-校正格式，实现了从霍普夫轨道出发的分支切换，成功检测到周期加倍和霍普夫-劈尖分岔。
• 事后弗洛quet乘子计算（floqap）与网格自适应（hopftref）的实现显著提升了长时间周期轨道延续的精度与鲁棒性。
• 该方法在不适定PDE中同样有效，例如最优污染控制问题，尽管标准时间积分方法失效，但周期轨道的延续仍具可行性。