A GALERKIN LEAST-SQUARES FINITE-ELEMENT METHOD FOR THE 2-DIMENSIONAL HELMHOLTZ-EQUATION

In this paper a Galerkin least-squares (GLS) finite element method, in which residuals in least-squares form are added to the standard Galerkin variational equation, is developed to solve the Helmholtz equation in two dimensions. An important feature of GLS methods is the introduction of a local mesh parameter that may be designed to provide accurate solutions with relatively coarse meshes. Previous work has accomplished this for the one-dimensional Helmholtz equation using dispersion analysis. In this paper, the selection of the GLS mesh parameter for two dimensions is considered, and leads to elements that exhibit improved phase accuracy. For any given direction of wave propagation, an optimal GLS mesh parameter is determined using two-dimensional Fourier analysis. In general problems, the direction of wave propagation will not be known a priori, In this case, an optimal GLS parameter is found which reduces phase error for all possible wave vector orientations over elements. The optimal GLS parameters are derived for both consistent and lumped mass approximations. Several numerical examples are given and the results compared with those obtained from the Galerkin method. The extension of GLS to higher-order quadratic interpolations is also presented.