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

In this paper, which is part II in a series of two, the pre-asymptotic error analysis of the continuous interior penalty finite element method (CIP-FEM) and the FEM for the Helmholtz equation in two and three dimensions is continued. While part I contained results on the linear CIP-FEM and FEM, the present part deals with approximation spaces of order $p\ge 1$. By using a modified duality argument, pre-asymptotic error estimates are derived for both methods under the condition of $\frac{kh}{p}\le C_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 pollution errors of both methods in $H^1$-norm are $O(k^{2p+1}{h}^{2p})$ if $p=O(1)$ and are $O\left(\frac{k}{p^2}\right(\frac{kh}{\sigma p}{)}^{2p})$ if the exact solution $u\in H^2(\text{Om})$ which coincide with existent dispersion analyses for the FEM on Cartesian grids. Here $\text{si}$ is a constant independent of $k,h,p$, and the penalty parameters. Moreover, it is proved that the CIP-FEM is stable for any $k,h,p>0$ and penalty parameters with positive imaginary parts. Besides the advantage of the absolute stability of the CIP-FEM compared to the FEM, the penalty parameters may be tuned to reduce the pollution effects.