[论文解读] Error bounds and analysis of proper orthogonal decomposition model reduction methods using snapshots from the solution and the time derivatives

Tanya Kostova, Geoffrey Oxberry+2 more|arXiv (Cornell University)|Jan 8, 2015

Model Reduction and Neural Networks参考文献 21被引用 5

一句话总结

本文比较了两种用于大型常微分方程组的本征正交分解（POD）模型降阶方法：一种仅使用解快照，另一种同时使用解快照和时间导数快照。文中推导了新的误差界，表明包含时间导数可降低近似误差，尤其在第一个被忽略的奇异值较小时效果更明显，并通过数值实验验证了基于导数的方法在使用更少快照时能显著降低误差。

ABSTRACT

Proper orthogonal decomposition methods for model reduction utilize information about the solution at certain time and parameter points to generate a reduced space basis. In this paper, we compare two proper orthogonal decomposition methods for reducing large systems of ODEs. The first method is based on collecting snapshots from the solutions only; the second method uses snapshots from both the solutions and their time derivatives. To compare the methods, we derive new bounds for the 2-norm of the approximation error induced by the each of the methods. The bounds are represented as a sum of two terms: the first depends on the size of the first neglected singular value while the second depends only on the spacings between the snapshots. We performed numerical experiments to compare the errors from the two model reduction methods applied to the semidiscretized FitzHugh-Nagumo system and investigated the relation between the behavior of the numerically observed error and the error bounds. We find that the error bounds, though not tight, provide insights and justification for using time derivative snapshots in POD model reduction for dynamical systems.

研究动机与目标

• 分析并比较两种POD模型降阶方法的精度：一种仅使用解快照，另一种同时使用解快照和时间导数快照。
• 为两种方法推导显式依赖于快照间距和第一个被忽略的奇异值大小的先验误差界。
• 研究在动力系统中引入时间导数信息是否能提高降阶模型的近似精度。
• 为使用时间导数快照在POD-ROM中的优势提供理论依据和数值证据。

提出的方法

• 通过积分误差方程并应用Gronwall不等式，在系统满足光滑性和Lipschitz条件的假设下，推导POD-ROM的误差界。
• 使用Lagrange插值和Hermite插值，将误差表示为快照之间时间间隔和近似阶数的函数。
• 将误差界表述为两项之和：一项依赖于第一个被忽略的奇异值，另一项依赖于快照间距Δi。
• 将该误差界应用于比较两种方法：方法1（仅使用解快照），方法2（使用解快照和时间导数快照）。
• 通过半离散化的FitzHugh-Nagumo系统的数值实验，验证理论误差界并比较实际误差。
• 分析不同快照间距和快照数量下的奇异值分布及误差行为。

实验结果

研究问题

• RQ1 在POD快照矩阵中包含时间导数快照是否能降低降阶模型的近似误差？
• RQ2 所推导的误差界如何依赖于快照之间的间距和第一个被忽略的奇异值的大小？
• RQ3 理论误差界能否准确预测POD-ROM中观测到的数值误差行为？
• RQ4 使用导数快照带来的精度提升是否在不同快照配置和系统动力学下保持一致？

主要发现

• 误差界显示，误差表达式中的第二项依赖于快照间距，使用时间导数的方法中该项为O(Δ⁴)，而仅使用解的方法中为O(Δ²)。
• 数值实验表明，方法2（包含时间导数）产生的近似误差显著小于方法1，尤其在快照数量较少时。
• 观测到的误差降低与理论误差界一致，尽管误差界并非紧致，表明时间导数快照通过减小与间距相关的误差项提升了精度。
• 时间导数快照不在解快照的张成空间中，这为通过其丰富降阶基空间提供了合理依据。
• 结果表明，使用导数快照可在减少快照数量的同时实现更低误差，从而降低计算成本。
• 误差界中误差项的阶数（α）在实际中高于理论预测，表明存在更紧致的理论误差界的可能性。