[论文解读] Nodal auxiliary a posteriori error estimates

Yuwen Li, Ludmil Zikatanov|arXiv (Cornell University)|Oct 14, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 68被引用 4

一句话总结

本文提出了一种新型框架，通过将误差控制与有限元外微分计算（FEEC）和节点辅助空间设置中的子空间校正及预条件化技术相联系，实现了有限元方法中后验误差估计的创新。该框架利用连续的正则分解和虚构空间引理，为 H(curl)、H(div)、Hodge 拉普拉斯及线弹性问题构建了鲁棒且与参数无关的误差估计器，从而在具有奇异扰动范数的各类偏微分方程中实现了可靠且高效的误差界。

ABSTRACT

We introduce and explain key relations between a posteriori error estimates and subspace correction methods viewed as preconditioners for problems in infinite dimensional Hilbert spaces. We set the stage using the Finite Element Exterior Calculus and Nodal Auxiliary Space Preconditioning. This framework provides a systematic way to derive explicit residual estimators and estimators based on local problems which are upper and lower bounds of the true error. We show the applications to discretizations of $δd$, curl-curl, grad-div, Hodge Laplacian problems, and linear elasticity with weak symmetry. We also provide a new regular decomposition for singularly perturbed H(d) norms and parameter-independent error estimators. The only ingredients needed are: well-posedness of the problem and the existence of regular decomposition on continuous level.

研究动机与目标

• 开发一个统一框架，将后验误差估计与无限维希尔伯特空间中的子空间校正和预条件化方法相联系。
• 推导出显式、可计算的误差估计器，确保对 de Rham 复形上广泛类别的偏微分方程具有可靠性和高效性。
• 为奇异扰动的 H(d) 范数构造新的正则分解，实现与材料参数无关的误差界。
• 将鲁棒误差估计的适用性扩展至不定和混合变分问题，包括具有弱对称性的线弹性问题。
• 证明对偶空间 V′ 的 Riesz 表示可通过 H¹ 正则分解进行预条件化，从而将误差估计简化为已建立的 H⁻¹ 残差估计。

提出的方法

• 利用虚构空间引理和节点辅助空间预条件化，将误差估计与希尔伯特空间中的算子预条件化相联系。
• 采用向量场和微分形式的连续 H¹ 正则分解，以控制对偶范数中的误差。
• 将对偶空间 V′ 中残差估计的问题，通过 Riesz 表示的谱等价性，转化为 H⁻¹ 范数的估计。
• 应用可交换插值和与外微分运算可交换的正则化插值，以在不同尺度间保持结构。
• 构造一种专为奇异扰动 H(d) 范数设计的新正则分解，确保对材料参数的鲁棒性。
• 利用 Hodge 分解和延拓定理处理混合边界条件和非平凡调和形式。

实验结果

研究问题

• RQ1 能否系统地利用预条件化理论和子空间校正方法推导出后验误差估计器？
• RQ2 如何为 H(curl) 和 H(div) 问题，尤其是在奇异扰动下，构建鲁棒且与参数无关的误差估计器？
• RQ3 连续 H¹ 正则分解在确保误差估计器的可靠性和高效性方面发挥何种作用？
• RQ4 该框架能否扩展至不定系统，如 Hodge 拉普拉斯和具有弱对称性的线弹性问题？
• RQ5 在混合边界条件下，调和形式的结构如何影响误差估计过程？

主要发现

• 所提出的框架为 de Rham 复形上的所有适定问题生成了具有双侧性的后验误差估计器，兼具可靠性和高效性。
• 该方法在 Lipschitz 域上实现了 H(curl) 和 H(div) 问题的参数无关可靠性，即使在奇异扰动范数存在的情况下亦成立。
• 为奇异扰动的 H(d) 范数构造了一种新正则分解，实现了与材料参数无关的鲁棒误差估计。
• 证明了对偶空间 V′ 的 Riesz 表示与 H⁻¹ 范数之和谱等价，从而将误差估计简化为标准残差界。
• 通过将对偶空间分解为理解充分的子空间的笛卡尔积，该方法推广至不定系统，同时保持了鲁棒性。
• 该框架在多个偏微分方程上得到验证，包括 curl-curl、grad-div、Hodge 拉普拉斯及具有弱对称性的线弹性问题，展示了其广泛适用性。