In comparison with low‐order finite element methods (FEMs), the use of oscillatory basis functions has been shown to reduce the computational complexity associated with the numerical approximation of Helmholtz problems at high wave numbers. We compare two different wave element methods for the 2D Helmholtz problems. The methods chosen for this study are the partition of unity FEM (PUFEM) and the ultra‐weak variational formulation (UWVF). In both methods, the local approximation of wave field is computed using a set of plane waves for constructing the basis functions. However, the methods are based on different variational formulations; the PUFEM basis also includes a polynomial component, whereas the UWVF basis consists purely of plane waves. As model problems we investigate propagating and evanescent wave modes in a duct with rigid walls and singular eigenmodes in an L‐shaped domain. Results show a good performance of both methods for the modes in the duct, but only a satisfactory accuracy was obtained in the case of the singular field. On the other hand, both the methods can suffer from the ill‐conditioning of the resulting matrix system. Copyright © 2008 John Wiley & Sons, Ltd.