[论文解读] Space-time least-squares finite elements for parabolic equations

Thomas Führer, Michael Karkulik|arXiv (Cornell University)|Nov 5, 2019

Advanced Numerical Methods in Computational Mathematics被引用 2

一句话总结

本文提出了一种基于一阶系统弱形式下最小化L²残差泛函的时空最小二乘有限元方法，用于热方程。该方法确保了均匀稳定性、对称正定的线性系统，并支持内置的后验误差估计器的全时空自适应，数值实验在单纯形网格上实现了最优收敛速率。

ABSTRACT

We present a space-time least squares finite element method for the heat equation. It is based on residual minimization in L2 norms in space-time of an equivalent first order system. This implies that (i) the resulting bilinear form is symmetric and coercive and hence any conforming discretization is uniformly stable, (ii) stiffness matrices are symmetric, positive definite, and sparse, (iii) we have a local a-posteriori error estimator for free. In particular, our approach features full space-time adaptivity. We also present a-priori error analysis on simplicial space-time meshes which are highly structured. Numerical results conclude this work.

研究动机与目标

• 开发一种稳定且符合要求的时空有限元方法，用于抛物型PDE，避免时间推进格式的局限性。
• 通过最小二乘变分格式，确保对任意符合要求的离散空间均具有均匀稳定性。
• 通过最小二乘泛函导出的局部误差指标，实现全时空自适应。
• 在结构化的单纯形时空网格上提供先验误差估计。
• 在1D和2D问题的数值实验中，展示最优收敛速率和鲁棒性能。

提出的方法

• 将热方程表述为一阶系统：∂tu − div σ = f，σ − ∇u = 0，初始条件为 u(0) = u₀。
• 定义最小二乘泛函 j(u, σ) = ∫₀ᵀ ‖∂tu − div σ − f‖²_{L²(Ω)} dt + ∫₀ᵀ ‖σ − ∇u‖²_{L²(Ω)} dt + ‖u(0) − u₀‖²_{L²(Ω)}。
• 在 H¹(0,T;L²(Ω)) × L²(0,T;H¹(Ω)) 的乘积空间上最小化 j(u, σ)，并施加适当的迹和正则性约束。
• 从泛函导出对称且强制的双线性形式，确保离散化设置下稳定且稀疏的刚度矩阵。
• 通过将泛函分解为单元贡献，构造局部后验误差估计器。
• 实现由误差估计器驱动的自适应加密，支持局部时空网格丰富化。

实验结果

研究问题

• RQ1 时空最小二乘有限元方法是否能在抛物问题中对任意符合要求的离散空间实现均匀稳定性？
• RQ2 在时空域上对L²范数的残差进行最小化，是否能产生对称、正定且稀疏的代数系统？
• RQ3 最小二乘泛函是否可用于推导可靠且高效的局部误差估计器以支持自适应加密？
• RQ4 该方法在具有不同正则性的问题上，于单纯形时空网格上的收敛速率如何？
• RQ5 在收敛速率和计算效率方面，自适应加密与均匀加密相比表现如何？

主要发现

• 该方法对任意符合要求的离散空间均实现均匀稳定性，确保了对任意符合离散化方案的鲁棒性。
• 所得刚度矩阵为对称、正定且稀疏，可高效地使用预处理共轭梯度等迭代求解器求解。
• 由最小二乘泛函导出的后验误差估计器天然具有局部性，可有效支持时空自适应。
• 数值实验表明，1D问题中自适应加密的收敛速率约为0.45，优于均匀加密的0.25。
• 在具有奇点的2D问题（如钝角）中，自适应加密将收敛速率从约0.2（均匀）提升至约0.24，证明了其在处理低正则性问题中的有效性。
• 对于光滑解，该方法在L²范数下达到N⁻¹/²阶的收敛速率，与最优h-加密行为一致。