[论文解读] Discrete Compactness for p-Version of Tetrahedral Edge Elements

Ralf Hiptmair|ArXiv.org|Jan 7, 2009

Advanced Numerical Methods in Computational Mathematics参考文献 34被引用 4

一句话总结

本文证明了在 $\boldsymbol{H}(\boldsymbol{curl}, \Omega)$ 中，$p$-版本四面体边有限元的离散紧致性性质，确立了麦克斯韦特征值问题伽辽金逼近的谱正确性。通过利用稳定的向量势构造和基于投影的插值方法，证明了在离散无散子空间中有界序列在 $p \to \infty$ 时于 $L^2(\Omega)$ 中收敛，从而确认了特征值与特征函数的渐近收敛性。

ABSTRACT

We consider the first family of $\Hcurl$-conforming Nedéléc finite elements on tetrahedral meshes. Spectral approximation ($p$-version) is achieved by keeping the mesh fixed and raising the polynomial degree $p$ uniformly in all mesh cells. We prove that the associated subspaces of discretely weakly divergence free piecewise polynomial vector fields enjoy a long conjectured discrete compactness property as $p o\infty$. This permits us to conclude asymptotic spectral correctness of spectral Galerkin finite element approximations of Maxwell eigenvalue problems.

研究动机与目标

• 为四面体网格上的 $p$-版本内德雷列克边元建立离散紧致性性质。
• 解决关于 $p$-版本边有限元方法在麦克斯韦特征值问题中谱正确性的长期猜想。
• 通过证明特征值与特征函数的收敛性，为电磁学中的 $p$-版本谱伽辽金方法提供严格的理论基础。
• 将离散紧致性理论扩展至三维 $hp$-有限元，其中先前结果仅限于 $h$-版本或二维情形。

提出的方法

• 采用固定网格、多项式阶数 $p$ 递增的 $p$-版本有限元方法，以实现谱逼近。
• 应用平滑的庞加莱映射构造方法，生成 $H^1$-稳定的向量势，确保离散亥姆霍兹分解的稳定性。
• 采用基于投影的插值算子，控制旋度算子离散核中逼近误差的控制。
• 依赖于离散向量场在 $L^2$-正交投影下分解为无旋与调和分量的亥姆霍兹分解。
• 通过分别分析无旋分量与调和分量的极限，建立离散解在 $L^2(\Omega)$ 中的收敛性。
• 利用斯特朗格引理与上同调空间估计，通过 $H^s$-正则性（$s > 1/2$）控制离散核逼近误差。

实验结果

研究问题

• RQ1 在一般四面体网格上，$p$-版本四面体边元在 $p \to \infty$ 时是否满足离散紧致性性质？
• RQ2 能否利用 $p$-版本边元证明麦克斯韦特征值问题的谱伽辽金方法具有谱正确性？
• RQ3 在 $L^2(\Omega)$-范数下，$p$-版本边元的离散无散子空间是否在 $L^2(\Omega)$ 中是紧致的？
• RQ4 能否利用离散向量场的 $L^2$-正交分解证明特征函数在 $L^2(\Omega)$ 中的收敛性？

主要发现

• 在 $\boldsymbol{H}(\operatorname{curl}, \Omega)$ 中，$p$-版本四面体边元的离散紧致性性质成立，即离散无散子空间中的任意有界序列在 $L^2(\Omega)$ 中存在收敛子列。
• 离散特征函数的极限属于连续无散空间 $\boldsymbol{H}(\operatorname{curl}, \Omega) \cap \boldsymbol{H}_0(\operatorname{div}_{\boldsymbol{\epsilon}}0, \Omega)$，从而确保谱正确性。
• 通过斯特朗格引理与 $L^2$-误差估计，建立了离散解无旋分量向极限的 $L^2$-投影收敛。
• 由于 $H^s$-正则性与投影算子的 $p$-鲁棒逼近，解的调和分量收敛于连续调和场。
• 离散核逼近误差以 $\mathcal{O}(\log^{3/2}p / p^s)$ 的速率衰减（$s > 1/2$），确认了谱收敛性。
• 结果证实了 $p$-版本边有限元方法在三维麦克斯韦特征值问题中的渐近最优性与谱正确性。