A fourth‐order compact scheme for the Helmholtz equation: Alpha‐interpolation of FEM and FDM stencils

Published in ''Int. Journal for Numerical Methods in Engineering'' Vol. 86 (1), pp. 18-46, 2011<br />
doi: 10.1002/nme.3043

== Abstract ==

We propose a fourth‐order compact scheme on structured meshes for the Helmholtz equation given by <math> R(\phi):=f(x)+\Delta \phi+\xi 2 \phi =0 </math>. The scheme consists of taking the alpha‐interpolation of the Galerkin finite element method and the classical central finite difference method. In 1D, this scheme is identical to the alpha‐interpolation method (J. Comput. Appl. Math. 1982; 8(1):15–19) and in 2D making the choice <math> \alpha=0.5 </math> we recover the generalized fourth‐order compact Padé approximation (J. Comput. Phys. 1995; 119:252–270; Comput. Meth. Appl. Mech. Engrg 1998; 163:343–358) (therein using the parameter <math> \gamma =2 </math>). We follow (SIAM Rev. 2000; 42(3):451–484; Comput. Meth. Appl. Mech. Engrg 1995; 128:325–359) for the analysis of this scheme and its performance on square meshes is compared with that of the quasi‐stabilized FEM (Comput. Meth. Appl. Mech. Engrg 1995; 128:325–359). In particular, we show that the relative phase error of the numerical solution and the local truncation error of this scheme for plane wave solutions diminish at the rate <math>O((\xi l)^4)</math>, where <math>\xi, l </math> represent the wavenumber and the mesh size, respectively. An expression for the parameter α is given that minimizes the maximum relative phase error in a sense that will be explained in Section 4.5. Convergence studies of the error in the <math>L^2</math> norm, the <math>H^1</math> semi‐norm and the <math>l^\infty</math> Euclidean norm are done and the pollution effect is found to be small. Copyright © 2010 John Wiley & Sons, Ltd.


<pdf>Media:Draft_Samper_192090958_6111_Nadukandi_et_al-2011-International_Journal_for_Numerical_Methods_in_Engineering.pdf</pdf>