[论文解读] An Analysis of Galerkin Proper Orthogonal Decomposition for Subdiffusion

Bangti Jin, Zhi Zhou|arXiv (Cornell University)|Aug 25, 2015

Fractional Differential Equations Solutions参考文献 24被引用 3

一句话总结

本文提出了一种用于具有Caputo分数阶导数的亚扩散问题的Galerkin-L1-POD格式，结合有限元法、L1时间半离散化以及本征正交分解（POD），以降低计算成本和存储需求。通过一种新颖的能量分析方法，提供了完整的误差分析，证明了在合理正则性假设下的最优收敛速率，并通过数值实验展示了其高效性。

ABSTRACT

In this work, we develop a novel Galerkin-L1-POD scheme for the subdiffusion model with a Caputo fractional derivative of order $α\in (0,1)$ in time, which is often used to describe anomalous diffusion processes in heterogeneous media. The nonlocality of the fractional derivative requires storing all the solutions from time zero. The proposed scheme is based on continuous piecewise linear finite elements, L1 time stepping, and proper orthogonal decomposition (POD). By constructing an effective reduced-order scheme using problem-adapted basis functions, it can significantly reduce the computational complexity and storage requirement. We shall provide a complete error analysis of the scheme under realistic regularity assumptions by means of a novel energy argument. Extensive numerical experiments are presented to verify the convergence analysis and the efficiency of the proposed scheme.

研究动机与目标

• 为解决由于非局部性导致的具有Caputo分数阶导数的亚扩散问题求解中计算成本和存储需求过高的问题。
• 基于亚扩散问题的解流形，开发一种适用于本征正交分解（POD）的降阶模型。
• 在合理的正则性假设下，对所提出的Galerkin-L1-POD格式提供完整的误差分析。
• 通过广泛的数值实验，验证该格式的效率与收敛性。

提出的方法

• 方法在空间上采用连续分段线性有限元，在时间上采用L1时间步长法处理Caputo分数阶导数。
• 通过全离散解和分数阶差分商构造快照，以生成POD基函数。
• 利用前 $ m $ 个POD基函数构建一个降阶的Galerkin框架，形成低维解空间 $ X_h^m $。
• 误差分析采用一种新颖的能量论证方法，推导出在弱正则性假设下的最优收敛速率。
• 在1D、2D和3D区域上通过广泛的数值实验验证了该格式。

实验结果

研究问题

• RQ1 能否为具有Caputo分数阶导数的亚扩散问题构造一种降阶的Galerkin-L1-POD格式，以显著降低计算成本？
• RQ2 在解的正则性假设合理的情况下，Galerkin-L1-POD格式的收敛行为如何？
• RQ3 与标准全阶方法相比，所提出的格式在效率和精度方面表现如何？
• RQ4 能否开发一种新颖的能量论证方法，以严格分析降阶格式中的误差？

主要发现

• 在合理正则性假设下，包括由于分数阶导数导致的有限光滑性，所提出的Galerkin-L1-POD格式在空间和时间上均实现了最优收敛速率。
• 误差分析基于一种新颖的能量论证，能够有效处理亚扩散解的非局部性和弱正则性。
• 数值实验验证了理论收敛速率，并显著降低了计算成本和存储需求。
• 从解快照中提取的POD基函数能够捕捉解流形的主要模态，从而实现有效的模型降阶。
• 即使解在时间导数上表现出奇异性（如 $ t^{eta-1} $ 行为，其中 $ eta = eta $），该格式仍保持稳定性和精度。