[论文解读] An a posteriori analysis of some inconsistent, nonconforming Galerkin methods approximating elliptic problems

Tristan Pryer|arXiv (Cornell University)|May 16, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 19被引用 5

一句话总结

本文为应用于椭圆问题的不一致、非协调Galerkin方法构建了一个后验误差分析框架，利用重构算子推导出可靠且高效的误差估计量。研究证明，诸如 quadrature 近似或稳定化带来的不一致性可被后验控制，从而即使在非变分或非光滑问题中，也能实现鲁棒的自适应有限元方法。

ABSTRACT

In this work we present an a posteriori analysis for classes of inconsistent, nonconforming schemes approximating elliptic problems. We show the estimates coincide with existing ones for interior penalty type discontinuous Galerkin approximations of the Laplacian and give new estimates for inconsistent discontinuous Galerkin approximation schemes of elliptic problems under quadrature approximation. We also examine the effect of inconsistencies on the a posteriori analysis of schemes applied to an unbalanced problem.

研究动机与目标

• 为应用于椭圆问题的不一致、非协调 Galerkin 格式开发一个统一的后验误差估计框架。
• 分析各种不一致性（特别是 quadrature 近似和稳定化）对后验误差控制的影响。
• 将分析扩展至非变分、非平衡的椭圆问题，推导出 H²-鲁棒的后验误差界。
• 在自适应算法中展示估计量的鲁棒性，特别是在早期细化阶段不一致性主导误差时。
• 证明即使在不一致离散化下，估计量仍能保持可靠性和有效性，其误差界包含数据振荡项。

提出的方法

• 采用受标准后验分析启发的重构算子方法，确保稳定性分析所需的光滑性。
• 构建一个类似于 Strang 引理的抽象后验估计，以量化离散格式中不一致性的效应。
• 将误差估计量分解为三部分：内部残差、跃迁项，以及由 quadrature 或稳定化引起的不一致项。
• 利用重构算子的存在性（无需可计算）推导界，从而在无需显式重构的情况下实现分析。
• 将该框架应用于标准的内罚不连续Galerkin方法及 quadrature 不一致格式，证明其收敛率与理论预测一致。
• 采用由后验估计量驱动的自适应算法，根据误差贡献（包括不一致项）对网格进行细化。

实验结果

研究问题

• RQ1 能否为不一致、非协调 Galerkin 方法构造出即使在格式不一致下仍保持可靠与高效的后验误差估计量？
• RQ2 双线性形式中 quadrature 近似对椭圆问题后验误差估计量的可靠性和有效性有何影响？
• RQ3 在自适应细化中，不一致项在多大程度上可被后验控制？何时其会主导误差估计量？
• RQ4 所提出的框架能否处理具有 H² 正则性的非变分、非平衡椭圆问题？其对应的误差界是什么？
• RQ5 当不一致性显著时，由所提估计量驱动的自适应算法相较于均匀细化的性能如何？

主要发现

• 对于不一致的不连续Galerkin格式，后验估计量与标准内罚方法的现有估计一致，验证了框架的有效性。
• 对于 quadrature 不一致格式，估计量可分解为可计算与不可计算部分，其中不一致项在粗网格阶段占主导地位。
• 在自适应算法中，估计量能有效引导网格在不一致区域和解的非光滑区域进行细化，实现 H² 范数下 O(N⁻¹/²) 的收敛率。
• 当精确解光滑但扩散系数强烈振荡时（如测试3），不一致项最初主导估计量，导致在高振荡区域出现密集细化。
• 对于非光滑解（如属于 H² 但不属于 H³ 的情况），估计量能有效捕捉解的不规则性与不一致效应，引导细化聚焦于奇点区域。
• 数值实验表明，估计量的渐近收敛率与理论预测一致，与定理 4.8 和 4.9 中的理论界相符。