[论文解读] Convergence of adaptive FEM for some elliptic obstacle problem with inhomogeneous Dirichlet data

Michael Feischl, Marcus Page+1 more|arXiv (Cornell University)|Jul 13, 2012

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

一句话总结

该论文证明了对于具有非齐次 Dirichlet 边界条件和一般 $ H^2 $ 障碍的椭圆障碍问题，自适应最低阶有限元方法（AFEM）的收敛性。通过将基于残差的误差估计器扩展以控制 Dirichlet 边界条件的振荡，并引入一种改进的能量-估计器耦合方式，作者建立了组合误差量的收缩性质，从而在不依赖嵌套离散空间的前提下确保了收敛性。

ABSTRACT

In this work, we show the convergence of adaptive lowest-order FEM (AFEM) for an elliptic obstacle problem with non-homogeneous Dirichlet data, where the obstacle χ is restricted only by χ in H^2(Ω). The adaptive loop is steered by some residual based error estimator introduced in Braess, Carstensen & Hoppe (2007) that is extended to control oscillations of the Dirichlet data, as well. In the spirit of Cascon et al. (2008), we show that a weighted sum of energy error, estimator, and Dirichlet oscillations satisfies a contraction property up to certain vanishing energy contributions. This result extends the analysis of Braess, Carstensen & Hoppe (2007) and Page & Praetorius (2010) to the case of non-homogeneous Dirichlet data as well as certain non-affine obstacles and introduces some energy estimates to overcome the lack of nestedness of the discrete spaces.

研究动机与目标

• 将自适应 FEM 的收敛性分析扩展至具有非齐次 Dirichlet 条件的椭圆障碍问题。
• 解决由于非协调的 Dirichlet 条件逼近而导致的离散有限元空间缺乏嵌套性的问题。
• 为一般 $ H^2 $ 障碍（包括非仿射情形）建立收敛性。
• 将 Dirichlet 条件振荡纳入误差估计器与分析框架中。
• 证明包含能量误差、估计器和数据振荡的组合误差量的收缩性质。

提出的方法

• 使用扩展以控制非齐次 Dirichlet 条件振荡的基于残差的误差估计器，自适应地对网格进行加密。
• 引入一种改进的误差量 $ \Delta_\ell $，其结合了能量误差、估计器和 Dirichlet 振荡项。
• 通过能量误差、估计器和数据振荡的加权和，建立在趋于消失的扰动下成立的收缩性质。
• 应用弱可靠性估计，以证明当 $ \varrho_\ell \to 0 $ 时，有 $ \|u - U_\ell\|_{H^1} \to 0 $。
• 采用估计器减少而非内部节点性质，使结果对各种加密策略具有鲁棒性。
• 引入具有离散 Dirichlet 数据 $ g_\ell $ 的连续解 $ u\_\ell $，以规避离散空间的非嵌套性问题。

实验结果

研究问题

• RQ1 自适应 FEM 是否能在具有非齐次 Dirichlet 条件和一般 $ H^2 $ 障碍的椭圆障碍问题上收敛？
• RQ2 在障碍问题的后验误差估计中，如何控制 Dirichlet 条件的振荡？
• RQ3 当离散空间非嵌套时，组合误差量的收缩性质是否仍然成立？
• RQ4 收敛性证明能否推广至非仿射障碍和非齐次 Dirichlet 条件的情形？
• RQ5 在非仿射和非齐次情况下，自适应 FEM 的收敛率是否最优？

主要发现

• 即使在离散空间非嵌套的情况下，自适应 FEM 对具有非齐次 Dirichlet 条件和 $ H^2 $ 障碍的椭圆障碍问题仍能收敛。
• 组合误差量 $ \Delta_\ell $ 满足收缩性质 $ \Delta_{\ell+1} \leq \kappa \Delta_\ell + \alpha_\ell $，其中 $ 0 < \kappa < 1 $，且 $ \alpha_\ell \to 0 $。
• 数值实验验证了能量误差的最优收敛率 $ \mathcal{O}(N^{-1/2}) $ 和 Dirichlet 振荡的 $ \mathcal{O}(N^{-3/4}) $ 收敛率。
• 估计器 $ \varrho_\ell $ 具有可靠性和效率，即 $ \|u - U_\ell\|_{H^1} \simeq \varrho_\ell $，数值结果中的平行曲线证实了这一点。
• 自适应加密优于均匀加密，尤其在角落奇点处表现更优，无论自适应参数 $ \theta \in \{0.4, 0.6, 0.8\} $ 取值如何，均能实现最优收敛率。
• 网格加密集中在非活动区域，即离散解未接触障碍的区域，符合预期。