The plane waves method for axisymmetric Helmholtz problems

The plane waves method is employed for the solution of Dirichlet and Neumann boundary value problems for the homogeneous Helmholtz equation in two- and three-dimensional domains possessing radial symmetry. The appropriate selection of collocation points and unitary direction vectors in the method leads to circulant and block circulant coefficient matrices in two and three dimensions, respectively. We propose efficient matrix decomposition algorithms which make use of fast Fourier transforms for the solution of the systems resulting from such a discretization. In conjunction with the method of particular solutions, the method is extended to the solution of inhomogeneous axisymmetric Helmholtz problems. Several numerical examples are presented.