Premium
Improvements for the ultra weak variational formulation
Author(s) -
Luostari T.,
Huttunen T.,
Monk P.
Publication year - 2013
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.4469
Subject(s) - basis function , basis (linear algebra) , helmholtz equation , discretization , plane wave , bessel function , mathematics , helmholtz free energy , mathematical analysis , finite element method , wave equation , geometry , physics , boundary value problem , optics , quantum mechanics , thermodynamics
SUMMARY In this paper, we investigate strategies to improve the accuracy and efficiency of the ultra weak variational formulation (UWVF) of the Helmholtz equation. The UWVF is a Trefftz type, nonpolynomial method using basis functions derived from solutions of the adjoint Helmholtz equation. We shall consider three choices of basis function: propagating plane waves (original choice), Bessel basis functions, and evanescent wave basis functions. Traditionally, two‐dimensional triangular elements are used to discretize the computational domain. However, the element shapes affect the conditioning of the UWVF. Hence, we investigate the use of different element shapes aiming to lower the condition number and number of degrees of freedom. Our results include the first tests of a plane wave method on meshes of mixed element types. In many modeling problems, evanescent waves occur naturally and are challenging to model. Therefore, we introduce evanescent wave basis functions for the first time in the UWVF to tackle rapidly decaying wave modes. The advantages of an evanescent wave basis are verified by numerical simulations on domains including curved interfaces.Copyright © 2013 John Wiley & Sons, Ltd.