Equation-based interpolation and incremental unknowns for solving the Helmholtz equation

Pascal Poullet, Amir Boag

Research output: Contribution to journalArticlepeer-review


A multilevel preconditioner for the Helmholtz equation using two types of incremental unknowns (IU) is developed. The transition between the two occurs when the mesh size reaches a predetermined fraction of the wavelength - roughly one quarter wavelength. Conventional IUs based on bilinear interpolation are employed for fine meshes like in an earlier paper by the authors. For coarse mesh sizes, novel IUs are defined using a Helmholtz/wave equation-based interpolation. The interpolation coefficients for the coarse meshes with dimension higher than one are derived numerically for stencils resembling integral representations for interior points. In two dimensions, the IUs are located on crosses surrounded by square contours. Numerical experiments on 1D and 2D examples have demonstrated the efficacy of the proposed approach in reducing the condition numbers and accelerating convergence for coarse grids with mesh sizes exceeding the wavelength.

Original languageEnglish
Pages (from-to)1148-1156
Number of pages9
JournalApplied Numerical Mathematics
Issue number11
StatePublished - Nov 2010


  • Helmholtz equation
  • Incremental unknowns
  • Iterative methods
  • Multilevel methods
  • Preconditioning


Dive into the research topics of 'Equation-based interpolation and incremental unknowns for solving the Helmholtz equation'. Together they form a unique fingerprint.

Cite this