[论文解读] Robust residual-based a posteriori error estimators for mixed finite element methods for fourth order elliptic singularly perturbed problems

Shaohong Du, Runchang Lin+1 more|arXiv (Cornell University)|Sep 15, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 14被引用 6

一句话总结

本论文为小扰动参数ε的四阶奇异摄动椭圆问题的混合有限元方法提出了一种基于残差的鲁棒后验误差估计器。通过提出一种平衡ε相关项与ε无关项的新误差度量，作者推导出在所有ε下均保持一致有效的估计器，数值实例验证了其最优收敛率和边界层的有效分辨率。

ABSTRACT

We consider mixed finite element approximation of a singularly perturbed fourth-order elliptic problem with two different boundary conditions, and present a new measure of the error, whose components are balanced with respect to the perturbation parameter. Robust residual-based a posteriori estimators for the new measure are obtained, which are achieved via a novel analytical technique based on an approximation result. Numerical examples are presented to validate our theory.

研究动机与目标

• 解决在小ε下设计混合有限元方法用于四阶奇异摄动问题的鲁棒后验误差估计器的挑战。
• 提出一种新误差度量，以平衡扰动参数ε与解正则性的影响，避免不平衡范数。
• 开发在ε → 0时仍保持鲁棒（一致有效）的基于残差的估计器，克服标准估计器随ε减小而性能下降的局限性。
• 通过数值方法验证理论估计器，证明其能够捕捉边界层并实现最优收敛率。

提出的方法

• 引入一种新误差度量，结合通量ψ = −Δu的能量范数与u的H¹半范数，具体为(ε²‖∇ψ‖² + ‖ψ‖² + |u|₁²)¹ᐟ²，以平衡ε相关与ε无关分量。
• 基于一种新颖的分析技术（基于逼近结果），推导出针对该新误差度量的基于残差的后验误差估计器。
• 采用由单元级估计量ηψ,T与ηu,T驱动的自适应有限元加密，以解析ψ中的边界层。
• 使用有效性指标验证估计器的可靠性，确保在不同ε值下估计误差均接近真实误差。
• 将估计器应用于两个基准问题，其边界层行为已知或可预期，采用标记参数θ的自适应网格加密。
• 通过证明收敛率与ε无关，特别是能量范数中达到O(h¹ᐟ²)的收敛率，验证估计器的鲁棒性。

实验结果

研究问题

• RQ1 能否构造一种新误差度量，以平衡四阶奇异摄动问题中ε相关与ε无关分量的贡献？
• RQ2 基于该新度量推导的基于残差的后验估计器是否对扰动参数ε具有鲁棒性？
• RQ3 所提出的估计器能否有效捕捉并解析ψ = −Δu中的陡峭边界层，尤其在ε极小时？
• RQ4 估计器在不同ε值（包括ε ≈ 10⁻⁷）下是否保持最优收敛率？
• RQ5 估计器的有效性指标是否在所有ε下一致有界，表明误差估计的可靠性？

主要发现

• 所提出的误差度量(ε²‖∇ψ‖² + ‖ψ‖² + |u|₁²)¹ᐟ²能有效平衡ε与解正则性的影响，避免因不平衡范数而无法捕捉层行为的问题。
• 基于残差的后验估计器对ε具有鲁棒性，在所有测试ε值下均保持一致可靠性，包括ε = 10⁻⁷。
• 数值实验表明，估计器在能量范数中实现约O(h¹ᐟ²)的最优收敛率，且与ε无关。
• 由估计器驱动的自适应网格加密成功解析了ψ中的四个陡峭边界层，经可视化与网格加密模式验证。
• ‖ψ − ψh‖ℰ的估计器有效性指标在所有ε值下均接近1（例如介于0.99至1.13之间），证实了估计器的可靠性。
• 当ε = 10⁻⁵时，经过22次加密迭代后，最小网格尺寸h_min(ε)达到6.10×10⁻⁵，表明在ε尺度上成功解析了边界层。