Initial Value Techniques for the Helmholtz and Maxwell Equations

Abstract

We study the initial value problem of the Helmholtz equation with spatially variable wave number. We show that it can be stabilized by suppressing the evanescent waves. The stabilized Helmholtz equation can be solved numerically by a marching scheme combined with FFT. The resulting algorithm has complexity $n^2 \log n$ on a $n \times n$ grid. We demonstrate the efficacy of the method by numerical examples with caustics. For the Maxwell equation the same treatment is possible after reducing it to a second order system. We show how the method can be used for inverse problems arising in acoustic tomography and microwave imaging.