[论文解读] Continuous Interior Penalty Finite Element Methods for the Helmholtz Equation with Large Wave Number

Haijun Wu|arXiv (Cornell University)|Jun 21, 2011

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

一句话总结

本文提出一种新颖的连续内部罚函数有限元法（CIP-FEM），用于大波数下的Helmholtz方程，采用具有正虚部的复罚参数。该方法实现了无条件稳定性与最优预渐近误差估计，与标准FEM相比显著降低了污染误差，误差界为 $ C_1kh + C_2k^3h^2 $（当 $ k^3h^2 \leq C_0 $ 时），否则为 $ C_1kh + \frac{C_2}{\gamma} $。

ABSTRACT

This paper develops and analyzes some continuous interior penalty finite element methods (CIP-FEMs) using piecewise linear polynomials for the Helmholtz equation with the first order absorbing boundary condition in two and three dimensions. The novelty of the proposed methods is to use complex penalty parameters with positive imaginary parts. It is proved that, if the penalty parameter is a pure imaginary number $ı\ga$ with $0 C_0$ and $kh$ is bounded, where $k$ is the wave number, $h$ is the mesh size, and the $C$'s are positive constants independent of $k$, $h$, and $\ga$. Optimal order $L^2$ error estimates are also derived. The analysis is also applied if the penalty parameter is a complex number with positive imaginary part. By taking $\ga o 0+$, the above estimates are extended to the linear finite element method under the condition $k^3h^2\le C_0$. Numerical results are provided to verify the theoretical findings. It is shown that the penalty parameters may be tuned to greatly reduce the pollution errors.

研究动机与目标

• 解决大波数Helmholtz问题中有限元方法长期存在的污染误差挑战。
• 为高波数问题开发一种无需网格约束的稳定且适定的有限元方法。
• 通过引入复罚参数，改进二维和三维线性有限元的预渐近误差估计。
• 证明通过调节罚参数可显著降低污染误差，提升数值精度。
• 通过令罚参数的虚部趋于零，将理论结果推广至标准FEM。

提出的方法

• 针对具有第一阶吸收边界条件的Helmholtz方程，引入使用分片线性单元的连续内部罚函数有限元法（CIP-FEM）。
• 采用形式为 $ \mathbf{i}\gamma $ 的复罚参数，其中 $ \gamma > 0 $，确保无条件稳定性和适定性。
• 推导 $ H^1 $-范数误差估计：当 $ k^3h^2 \leq C_0 $ 时为 $ C_1kh + C_2k^3h^2 $，当 $ k^3h^2 > C_0 $ 时为 $ C_1kh + \frac{C_2}{\gamma} $，其中常数与 $ k $、$ h $ 和 $ \gamma $ 无关。
• 在相同条件下建立最优的 $ L^2 $-范数误差估计。
• 将分析推广至具有正虚部的一般复罚参数，实现通过参数调节降低误差。
• 在 $ \gamma \to 0^+ $ 的极限下恢复标准FEM的误差估计，该结果在条件 $ k^3h^2 \leq C_0 $ 下成立。

实验结果

研究问题

• RQ1 具有复罚参数的连续内部罚函数FEM能否在大波数下实现Helmholtz方程的无条件稳定性？
• RQ2 CIP-FEM在 $ H^1 $-范数下的预渐近误差估计是什么？其与 $ k $、$ h $ 和 $ \gamma $ 的依赖关系如何？
• RQ3 能否通过使用复罚参数显著减少线性FEM在高波数Helmholtz问题中的污染误差？
• RQ4 罚参数的选择在实际中如何影响收敛行为和相对误差？
• RQ5 CIP-FEM在何种程度上能以更少自由度实现与IPDG等更昂贵方法相当的精度？

主要发现

• 采用复罚参数 $ \mathbf{i}\gamma $ 的CIP-FEM对所有 $ k > 0 $ 均实现无条件稳定与适定，且无需网格约束。
• $ H^1 $-范数误差估计为：当 $ k^3h^2 \leq C_0 $ 时为 $ C_1kh + C_2k^3h^2 $，当 $ k^3h^2 > C_0 $ 时为 $ C_1kh + \frac{C_2}{\gamma} $，其中常数与 $ k $、$ h $ 和 $ \gamma $ 无关。
• 推导出最优的 $ L^2 $-范数误差估计，证实了该方法的鲁棒性。
• 数值结果表明，调节罚参数 $ \gamma $ 可极大降低污染误差，CIP-FEM在低 $ kh $ 值下亦优于标准FEM。
• 当 $ k=100 $ 时，CIP-FEM仅需 35,971 个自由度即可在 $ H^1 $-半范数下达到 30% 的相对误差，而标准FEM需 229,357 个自由度。
• CIP-FEM 的最大波数 $ k_{\text{max}} $ 满足 $ h(k,0.5) \geq 0.001 $ 时为 622，超过标准FEM的两倍，显著降低了污染效应。