Linear continuous interior penalty finite element method for Helmholtz equation With High Wave Number: One-Dimensional Analysis
Linear continuous interior penalty finite element method for Helmholtz equation With High Wave Number: One-Dimensional Analysis
复制标题
高波数亥姆霍兹方程的线性连续内罚有限元法:一维分析
DOI:
10.1002/num.22054
复制
发表时间:
2016
影响因子:
3.9
通讯作者:
Zhu Lingxue
中科院分区:
文献类型:
--
作者:
Burman Erik;Wu Haijun;Zhu Lingxue
This article addresses the properties of continuous interior penalty (CIP) finite element solutions for the Helmholtz equation. The ‐version of the CIP finite element method with piecewise linear approximation is applied to a one‐dimensional (1D) 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 ‐norm, as the sum of best approximation error plus a pollution term that is the order of the phase difference. It is proved that the pollution effect 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. In particular, we give numerical evidence that the optimal penalty parameter obtained in the 1D case also works very well for the CIP‐FEM on two‐dimensional Cartesian grids.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1378–1410, 2016