[论文解读] Computing Coherent Sets using the Fokker-Planck Equation

Andreas Denner, Oliver Junge+1 more|arXiv (Cornell University)|Dec 11, 2015

Statistical Mechanics and Entropy参考文献 18被引用 4

一句话总结

该论文提出了一种新颖的方法，通过求解福克-普朗克方程，直接在随机扰动下近似转移算子，从而计算动力系统中的相干集。该方法采用谱配置法和指数时间差分法，实现了高精度与高效率。该方法避免了轨迹积分，通过基于福克-普朗克方程的转移算子的奇异向量分析，实现了对相干集的鲁棒检测。

ABSTRACT

We perform a numerical approximation of coherent sets in finite-dimensional smooth dynamical systems by computing singular vectors of the transfer operator for a stochastically perturbed flow. This operator is obtained by solution of a discretized Fokker-Planck equation. For numerical implementation, we employ spectral collocation methods and an exponential time differentiation scheme. We experimentally compare our approach to the more classical method by Ulam that is based on integration of the exact transfer operator.

研究动机与目标

• 开发一种用于识别在小随机扰动下鲁棒的有限维动力系统中相干集的数值方法。
• 用基于福克-普朗克方程的PDE方法替代经典的基于轨迹的方法，以建模随机转移算子。
• 通过直接使用谱配置法求解对流-扩散方程，提高数值效率并减少数值扩散。
• 实现在无需拉格朗日轨迹积分或离格点插值的情况下，准确计算相干集。
• 与经典Ulam方法进行对比验证，证明该方法在检测相干结构方面具有更优性能。

提出的方法

• 该方法通过求解福克-普朗克方程获得随机转移算子，以建模在随机扰动流场下概率密度的演化，从而计算相干集。
• 采用谱配置法对福克-普朗克方程进行离散化，实现在固定且粗略的网格上高阶空间精度。
• 采用指数时间差分（etd）格式处理福克-普朗克方程中拉普拉斯项引起的刚性问题。
• 在高对流区域中使用对流项的反对称形式，以减轻混叠误差。
• 通过时间积分福克-普朗克方程的基函数解，构建最终的转移算子，并计算其奇异向量以识别相干集。
• 该方法避免了显式的轨迹积分与插值，仅依赖于在网格点处对向量场的采样。

实验结果

研究问题

• RQ1 福克-普朗克方程能否作为识别动力系统中相干集的直接计算工具？
• RQ2 基于福克-普朗克方程的转移算子在精度和数值稳定性方面与经典Ulam方法相比如何？
• RQ3 谱配置法与指数时间积分能否为相干集检测提供比基于轨迹的方法更高效、更准确的替代方案？
• RQ4 与Galerkin法或有限差分法相比，求解福克-普朗克方程在减少数值扩散方面能有多大程度的改善？
• RQ5 直接求解福克-普朗克方程能否产生在几何上具有意义且对扰动具有鲁棒性的相干集？

主要发现

• 基于福克-普朗克方程的转移算子是紧算子，将 $L^2$ 中的有界集映射到 $H^1$ 中的有界集，确保了奇异向量的适定性与收敛性。
• 该方法在粗网格下仍能实现高精度，且无需拉格朗日轨迹积分或离格点插值。
• 数值实验表明，福克-普朗克方法在高对流区域中比经典Ulam方法更鲁棒地检测出相干结构。
• 与标准有限差分法或Galerkin方法相比，谱配置法与指数时间积分显著减少了数值扩散与混叠误差。
• 该方法成功识别出四涡旋流场中的相干集，奇异向量揭示了与已知输运屏障对齐的空间相干区域。
• 理论分析证实，基于福克-普朗克方程的算子保持了正性与随机性等关键性质，确保了物理解释的一致性。