Pre-asymptotic Error Analysis of CIP-FEM and FEM for Helmholtz Equation with High Wave Number. Part II: $hp$ Version

Lingxue Zhu,Haijun Wu

DOI: https://doi.org/10.1137/120874643

IF: 3.039

2012-01-01

SIAM Journal on Numerical Analysis

Abstract:In this paper, which is the second in a series of two, thepreasymptotic error analysis of the continuous interior penaltyfinite element method (CIP-FEM) and the FEM for the Helmholtzequation in two and three dimensions is continued. While Part Icontained results on the linear CIP-FEM and FEM, the present partdeals with approximation spaces of order $p \ge 1$. By using amodified duality argument, preasymptotic error estimates arederived for both methods under the condition of $\frac{kh}{p}\leC_0(\frac{p}{k})^{\frac{1}{p+1}}$, where $k$ is the wave number,$h$ is the mesh size, and $C_0$ is a constant independent of $k,h, p$, and the penalty parameters. It is shown that the pollutionerrors of both methods in $H^1$-norm are $O(k^{2p+1}h^{2p})$ if$p=O(1)$ and are $O(\frac{k}{p^2}(\frac{kh}{\sigmap})^{2p})$ if the exact solution $u\in H^2(\Omega)$ whichcoincide with existent dispersion analyses for the FEM onCartesian grids. Here $\sigma$ is a constant independent of $k, h,p$ and the penalty parameters. Moreover, it is proved that theCIP-FEM is stable for any $k, h, p>0$ and penalty parameters withpositive imaginary parts. Besides the advantage of the absolutestability of the CIP-FEM compared to the FEM, the penaltyparameters may be tuned to reduce the pollution effects.