[论文解读] L2-Estimates for a High Order Unfitted Finite Element Method for Elliptic Interface Problems

Christoph Lehrenfeld, Arnold Reusken|arXiv (Cornell University)|Apr 15, 2016

Advanced Numerical Methods in Computational Mathematics被引用 6

一句话总结

本文通过推导最优 L² 误差界，将一种高阶非拟合有限元方法扩展至椭圆界面问题，该方法基于参数化映射技术，利用光滑水平集函数提升几何逼近精度。该方法在 L²-范数下实现了最优收敛率，补充了先前对 H¹-范数的误差分析。

ABSTRACT

In the context of unfitted finite element discretizations the realization of high order methods is challenging due to the fact that the geometry approximation has to be sufficiently accurate. Recently a new unfitted finite element method was introduced which achieves a high order approximation of the geometry for domains which are implicitly described by smooth level set functions. This method is based on a parametric mapping which transforms a piecewise planar interface (or surface) reconstruction to a high order approximation. In the paper [C. Lehrenfeld, A. Reusken, Analysis of a High Order Finite Element Method for Elliptic Interface Problems, arXiv 1602.02970, submitted] an a priori error analysis of the method applied to an interface problem is presented. The analysis reveals optimal order discretization error bounds in the $H^1$-norm. In this paper we extend this analysis and derive optimal $L^2$-error bounds.

研究动机与目标

• 将现有高阶非拟合有限元方法的先验误差分析从 H¹-范数扩展至 L²-范数。
• 为通过隐式定义的界面求解椭圆界面问题，建立 L²-范数下的最优收敛率。
• 通过参数化映射改进几何逼近，确保在非拟合离散化下保持高阶精度。
• 提供一个完整的误差分析框架，支持该方法在 H¹ 和 L² 误差估计两方面的分析。

提出的方法

• 该方法采用参数化映射，将分片平面的界面重构转换为对真实界面的高阶逼近。
• 几何逼近基于光滑水平集函数，从而在表示弯曲界面时实现高阶精度。
• 非拟合有限元方法使用与界面无关的背景网格，并采用稳定化技术处理被切割的单元。
• 分析基于高阶参数化有限元公式，保持最优收敛特性。
• 通过对偶论证和 L²-范数中的插值误差界推导误差估计。
• 即使界面切割背景网格，该方法仍保持最优收敛率。

实验结果

研究问题

• RQ1 能否为应用于椭圆界面问题的高阶非拟合有限元方法推导出最优 L² 误差界？
• RQ2 参数化映射技术如何影响非拟合离散化中的 L²-逼近误差？
• RQ3 当界面由光滑水平集函数隐式定义时，该方法是否在 L²-范数下保持最优收敛率？
• RQ4 几何逼近误差与由此产生的有限元离散化误差在 L²-范数下的相互作用机制是什么？
• RQ5 能否在非拟合方法的背景下，利用对偶技术严格建立 L²-误差分析？

主要发现

• 该方法在椭圆界面问题中实现了最优阶的 L²-误差界，与给定有限元空间所能达到的最佳收敛率一致。
• 参数化映射实现了界面的高阶逼近，这对在 L²-范数下保持最优收敛至关重要。
• 分析结果证实，L²-误差以与 H¹-误差相同的最优速率收敛，验证了该方法的鲁棒性。
• 误差界通过对其偶论证推导得出，确保了估计的紧致性和可靠性。
• 由于高阶几何表示，即使界面切割背景网格，该方法仍保持稳定和高精度。
• 本研究通过在 H¹ 和 L² 范数下完成误差分析，拓展了非拟合有限元方法的理论基础。