[论文解读] A Posteriori Error Analysis of hp-FEM for singularly perturbed problems

Jens Markus Melenk, Thomas P. Wihler|arXiv (Cornell University)|Aug 27, 2016

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

一句话总结

本文针对奇异摄动边值问题提出了一个鲁棒的后验误差估计，适用于hp-有限元方法，即使微分算子非强制时亦成立。该误差界在网格尺寸h和多项式阶p方面完全显式，且对摄动参数保持鲁棒，确保在具有挑战性的奇异摄动区域中实现可靠的自适应方法。

ABSTRACT

We consider the approximation of singularly perturbed linear second-order boundary value problems by hp-finite element methods. In particular, we include the case where the associated differential operator may not be coercive. Within this setting we derive an a posteriori error estimate for a natural residual norm. The error bound is robust with respect to the perturbation parameter and fully explicit with respect to both the local mesh size h and the polynomial degree p.

研究动机与目标

• 为奇异摄动问题中的hp-有限元方法开发一个可靠的后验误差估计器。
• 解决底层微分算子不具强制性的情况，此类情况会增加标准误差分析的复杂性。
• 确保误差界对摄动参数保持鲁棒，这是奇异摄动问题中的关键要求。
• 为自适应加密提供在局部网格尺寸h和多项式阶p方面完全显式的误差估计。

提出的方法

• 该方法基于针对奇异摄动问题结构量身定制的自然残差范数，推导出后验误差估计。
• 采用对偶论证方法，在双线性形式缺乏强制性时，仍能控制能量范数中的误差。
• 该分析在统一框架下同时考虑h-加密（网格尺寸）和p-富集（多项式阶数）。
• 所得误差估计器完全显式，仅依赖于局部单元残差和解数据。
• 通过保持与摄动参数无关的统一有界性，确保方法的鲁棒性。
• 适用于一般线性二阶边值问题，包括边界层现象。

实验结果

研究问题

• RQ1 能否为具有非强制算子的奇异摄动问题中的hp-FEM构造一个可靠的后验误差估计器？
• RQ2 在摄动参数趋近于零的极限下，误差界是否对摄动参数保持鲁棒？
• RQ3 如何使误差估计器在h和p方面完全显式，以支持自适应网格加密？
• RQ4 残差范数在准确捕捉边界层效应方面起到什么作用？
• RQ5 能否在不假设微分算子强制性的前提下推导出该估计器？

主要发现

• 后验误差估计器对摄动参数保持鲁棒，可在所有参数范围内维持统一有界性。
• 误差界在局部网格尺寸h和多项式阶p方面完全显式，可在每个单元上直接计算。
• 该方法适用于微分算子不具强制性的问题，将hp-FEM的应用范围扩展到更一般的奇异摄动问题。
• 基于残差的误差估计器能准确捕捉边界层现象，即使在存在剧烈梯度时亦然。
• 理论框架支持自适应hp-加密策略，其收敛性不依赖于摄动参数。