A numerical integration scheme for special finite elements for the Helmholtz equation

The theory for integrating the element matrices for rectangular, triangular and quadrilateral finite elements for the solution of the Helmholtz equation for very short waves is presented. A numerical integration scheme is developed. Samples of Maple and Fortran code for the evaluation of integration abscissae and weights are made available. The results are compared with those obtained using large numbers of Gauss-Legendre integration points for a range of testing wave problems. The results demonstrate that the method gives correct results, which gives confidence in the procedures, and show that large savings in computation time can be achieved.

[1]  P. Bettess,et al.  Analysis of the conditioning number of the plane wave approximation for the helmholtz equation , 2000 .

[2]  Darren Redfern,et al.  The Maple handbook , 1994 .

[3]  Roland Bulirsch,et al.  Numerical calculation of the sine, cosine and fresnel integrals , 1967 .

[4]  Masayuki Kikuchi,et al.  Inversion of complex body waves , 1982 .

[5]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[6]  C. Farhat,et al.  The Discontinuous Enrichment Method , 2000 .

[7]  Omar Laghrouche,et al.  Short wave modelling using special finite elements , 2000 .

[8]  Jan Mandel,et al.  The Finite Ray Element Method for the Helmholtz Equation of Scattering: First Numerical Experiments , 1997 .

[9]  P. Bettess,et al.  Diffraction of short waves modelled using new mapped wave envelope finite and infinite elements , 1999 .

[10]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[11]  Armel de La Bourdonnaye,et al.  High frequency approximation of integral equations modeling scattering phenomena , 1994 .

[12]  Edmund Chadwick,et al.  Wave envelope examples for progressive waves , 1995 .

[13]  I. Babuska,et al.  Dispersion and pollution of the FEM solution for the Helmholtz equation in one, two and three dimensions , 1999 .

[14]  R. J. Astley,et al.  Modelling of short wave diffraction problems using approximating systems of plane waves , 2002 .

[15]  P. Ortiz,et al.  An improved partition of unity finite element model for diffraction problems , 2001 .

[16]  P. Bettess,et al.  Coupling of mapped wave infinite elements and plane wave basis finite elements for the Helmholtz equation in exterior domains , 2003 .

[17]  Charbel Farhat,et al.  A discontinuous finite element method for the helmholtz equation , 2000 .

[18]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[19]  J. Avilés,et al.  A boundary method for elastic wave diffraction: Application to scattering of SH waves by surface irregularities , 1982 .

[20]  I. Babuska,et al.  Finite Element Solution of the Helmholtz Equation with High Wave Number Part II: The h - p Version of the FEM , 1997 .

[21]  I. Babuska,et al.  Finite element solution of the Helmholtz equation with high wave number Part I: The h-version of the FEM☆ , 1995 .

[22]  J. Z. Zhu,et al.  The finite element method , 1977 .

[23]  Richard W. Hamming,et al.  Numerical Methods for Scientists and Engineers , 1962 .

[24]  Edmund Chadwick,et al.  Modelling of progressive short waves using wave envelopes , 1997 .

[25]  Ismael Herrera,et al.  Connectivity as an alternative to boundary integral equations: Construction of bases. , 1978, Proceedings of the National Academy of Sciences of the United States of America.

[26]  P. Bettess,et al.  A numerical integration scheme for special quadrilateral finite elements for the Helmholtz equation , 2002 .