A Fast Multigrid Approach for Solving the Helmholtz Equation with a Point Source

By Eran Treister

University of British Columbia

Published on


Solving the discretized Helmholtz equations with high wave numbers in large dimensions is a challenging task. However, in many scenarios, the solution of these equations is required for a point source. In this case, the problem can be be reformulated and split into two parts: one in a solution of the eikonal equation for the travel time, and the second is a solution of a complex-valued advection-diffusion-reaction (ADR) equation for the amplitude. The eikonal equation for the travel time can be solved efficiently by several existing solvers in the literature. Given the travel time, the ADR equation for the amplitude can be efficiently solved by standard multigrid methods, since the solution for the amplitude is typically smooth in this case. The only difference between our numerical solution and a solution of the original discretized Helmholtz equations is due to discretization and roundoff errors. Numerical results show that this is an efficient solution algorithm for a wide range of frequencies.

Cite this work

Researchers should cite this work as follows:

  • Eran Treister (2016), "A Fast Multigrid Approach for Solving the Helmholtz Equation with a Point Source," http://nanohub.org/resources/23504.

    BibTex | EndNote


NanoBio Node, Aly Taha

University of Illinois at Urbana-Champaign