Compact optimal quadratic spline collocation methods for the Helmholtz equation
Date
2011Source
Journal of Computational PhysicsVolume
230Issue
8Pages
2880-2895Google Scholar check
Keyword(s):
Metadata
Show full item recordAbstract
Quadratic spline collocation methods are formulated for the numerical solution of the Helmholtz equation in the unit square subject to non-homogeneous Dirichlet, Neumann and mixed boundary conditions, and also periodic boundary conditions. The methods are constructed so that they are: (a) of optimal accuracy, and (b) compact that is, the collocation equations can be solved using a matrix decomposition algorithm involving only tridiagonal linear systems. Using fast Fourier transforms, the computational cost of such an algorithm is O(N2logN) on an N×N uniform partition of the unit square. The results of numerical experiments demonstrate the optimal global accuracy of the methods as well as superconvergence phenomena. In particular, it is shown that the methods are fourth-order accurate at the nodes of the partition. © 2011 Elsevier Inc.