We investigate the ultra weak variational formulation (UWVF) of the 2-D Helmholtz equation using a new choice of basis functions. Traditionally the UWVF basis functions are chosen to be plane waves. Here, we instead use first kind Bessel functions. We compare the performance of the two bases. Moreover, we show that it is possible to use coupled plane wave and Bessel bases in the same mesh. As test cases we shall consider propagating plane and evanescent waves in a rectangular domain and a singular 2-D Helmholtz problem in an L-shaped domain.
