Premium
A Petrov–Galerkin formulation for the alpha interpolation of FEM and FDM stencils: Applications to the Helmholtz equation
Author(s) -
Nadukandi Prashanth,
Oñate Eugenio,
Garcia Julio
Publication year - 2011
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.3291
Subject(s) - stencil , finite element method , helmholtz equation , mathematics , interpolation (computer graphics) , petrov–galerkin method , galerkin method , polygon mesh , mathematical analysis , finite difference , boundary value problem , geometry , physics , classical mechanics , motion (physics) , computational science , thermodynamics
SUMMARY A new Petrov–Galerkin (PG) method involving two parameters, namely α 1 and α 2 , is presented, which yields the following schemes on rectangular meshes: (i) a compact stencil obtained by the linear interpolation of the Galerkin FEM and the classical central finite difference method (FDM), should the parameters be equal, that is, α 1 = α 2 = α ; and (ii) the nonstandard compact stencil presented in ( Int. J. Numer. Meth. Engng 2011; 86:18–46) for the Helmholtz equation if the parameters are distinct, that is, α 1 ≠ α 2 . The nonstandard compact stencil is obtained by taking the linear interpolation of the diffusive terms (specified by α 1 ) and the mass terms (specified by α 2 ) that appear in the stencils obtained by the standard Galerkin FEM and the classical central FDM, respectively. On square meshes, these two schemes were shown to provide solutions to the Helmholtz equation that have a dispersion accuracy of fourth and sixth order, respectively ( Int. J. Numer. Meth. Engng 2011; 86:18–46). The objective of this paper is to study the performance of this PG method for the Helmholtz equation using nonuniform meshes and the treatment of natural boundary conditions. Copyright © 2011 John Wiley & Sons, Ltd.