A difficulty of the standard Galerkin finite element method has been the ability to accurately resolve oscillating wave solutions at higher frequencies. Many alternative methods have been developed including high-order methods, stabilized Galerkin methods, multi-scale variational methods, and other wave-based discretization methods. In this work, consistent residuals, both in the form of least-squares and gradient least-squares are linearly combined and added to the Galerkin variational Helmholtz equation to form a new generalized Galerkin least-squares method (GGLS). By allowing the stabilization parameters to vary spatially within each element, we are able to select optimal parameters which reduce dispersion error for all wave directions from second-order to fourth-order in nondimensional wavenumber; a substantial improvement over standard Galerkin elements. Furthermore, the stabilization parameters are frequency independent, and thus can be used for both time-harmonic solutions to the Helmholtz equation as well as direct time-integration of the wave equation, and eigenfrequency/eigenmodes analysis. Since the variational framework preserves consistency, high-order accuracy is maintained in the presence of source terms. In the case of homogeneous source terms, we show that our consistent variational framework is equivalent to integrating the underlying stiffness and mass matrices with optimally selected numerical quadrature rules. Optimal GGLS stabilization parameters and equivalent quadrature rules are determined for several element types including: bilinear quadrilateral, linear triangle, and linear tetrahedral elements. Numerical examples on unstructured meshes validate the expected high-order accuracy.
Stable high-order quadrature rules with equidistant points