[论文解读] Continuous Interior Penalty Finite Element Method for Helmholtz Equation with High Wave Number: One Dimensional Analysis

Lingxue Zhu, Erik Burman+1 more|arXiv (Cornell University)|Nov 7, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 16被引用 7

一句话总结

本文研究了一维亥姆霍兹方程在高波数下连续内部罚（CIP）有限元方法的性能。证明了当适当地选择实值罚参数时，标准FEM在高k下常见的污染误差被消除，从而在H¹范数下实现最优收敛，误差界为$ kh\|f\| $，并通过严格的误差估计和数值实验加以验证。

ABSTRACT

This paper addresses the properties of Continuous Interior Penalty (CIP) finite element solutions for the Helmholtz equation. The $h$-version of the CIP finite element method with piecewise linear approximation is applied to a one-dimensional model problem. We first show discrete well posedness and convergence results, using the imaginary part of the stabilization operator, for the complex Helmholtz equation. Then we consider a method with real valued penalty parameter and prove an error estimate of the discrete solution in the $H^1$-norm, as the sum of best approximation plus a pollution term that is the order of the phase difference. It is proved that the pollution can be eliminated by selecting the penalty parameter appropriately. As a result of this analysis, thorough and rigorous understanding of the error behavior throughout the range of convergence is gained. Numerical results are presented that show sharpness of the error estimates and highlight some phenomena of the discrete solution behavior.

研究动机与目标

• 分析一维亥姆霍兹方程在高波数下CIP有限元方法的稳定性和收敛性。
• 研究稳定化参数在控制振幅和相位误差（尤其是污染效应）中的作用。
• 推导显式包含波数和网格尺寸依赖关系的预渐近误差估计。
• 确定通过最优罚参数选择消除污染误差的条件。
• 通过数值实验验证理论结果，展示精确的误差估计和离散解的行为。

提出的方法

• 对具有复吸收边界条件的一维亥姆霍兹模型问题应用h-型CIP-FEM，采用分段线性有限元。
• 在罚参数中引入虚部以确保离散适定性，并推导出阶为$ kh + \min(1, k^3 h^2) $的误差估计。
• 对于实值罚参数，通过构造离散格林函数将误差分解为最佳逼近项和相位误差项。
• 推导出与$ |k_h^- - k| $成正比的相位误差估计，表明污染项依赖于离散波数与精确波数的偏差。
• 证明当罚参数$ \gamma $在最优值$ \gamma_o $的$ O(h) $区间内选取时，污染误差消失，从而实现最优$ kh\|f\| $收敛。
• 通过基本解系理论显式求解全局离散系统，以支持稳定性与误差分析。

实验结果

研究问题

• RQ1 CIP-FEM在高波数下一维亥姆霍兹方程中的表现如何，特别是在稳定性和收敛性方面？
• RQ2 稳定化参数的实部或虚部对振幅和相位精度有何影响？
• RQ3 是否可以在CIP-FEM框架中消除污染误差？若可，对罚参数有何条件？
• RQ4 离散波数与精确波数之间的关系如何？这对误差有何影响？
• RQ5 收敛所需的临界自由度数量$ N_c $如何随波数$ k $变化？

主要发现

• 当罚参数具有虚部时，H¹半范数下的误差满足$ \|(u - u_h)'’\lesssim kh + \min(1, k^3 h^2)\|f\| $，确保了适定性和鲁棒性。
• 当罚参数为实数时，误差分解为最佳逼近项和与相位差$ |k_h^- - k| $成正比的污染项。
• 当罚参数$ \gamma $在最优值$ \gamma_o $的$ O(h) $区间内选取时，污染误差被消除，从而获得最优误差估计$ \|(u - u_h)'’\lesssim kh\|f\| $。
• 收敛的临界自由度数$ N_c $约为$ N_c \approx k/\pi $，与最佳逼近和CIP-FEM在$ \gamma = \gamma_o $时的理论预测一致。
• 当$ \gamma = -1/12 $时，预测的$ N_c \sim (k^5/720)^{1/4} $，数值结果证实了大$ k $时的高精度。
• 数值结果验证了误差估计的紧致性，表明当$ \gamma $最优时，相对误差比$ e_c/e_{ba} $保持有界，而当$ \gamma $次优时该比值增长，证实了无污染现象的存在。