Taylor expansion and discretization errors in Gaussian beam superposition

The Gaussian beam superposition method is an asymptotic method for computing high frequency wave fields in smoothly varying inhomogeneous media. In this paper we study the accuracy of the Gaussian beam superposition method and derive error estimates related to the discretization of the superposition integral and the Taylor expansion of the phase and amplitude off the center of the beam. We show that in the case of odd order beams, the error is smaller than a simple analysis would indicate because of error cancellation effects between the beams. Since the cancellation happens only when odd order beams are used, there is no remarkable gain in using even order beams. Moreover, applying the error estimate to the problem with constant speed of propagation, we show that in this case the local beam width is not a good indicator of accuracy, and there is no direct relation between the error and the beam width. We present numerical examples to verify the error estimates.

[1]  L. Hörmander,et al.  The Analysis of Linear Partial Differential Operators I: Distribution Theory and Fourier Analysis , 1983 .

[2]  N. Tanushev.  Superpositions and higher order Gaussian beams , 2008 .

[3]  Shingyu Leung,et al.  Eulerian Gaussian Beams for High Frequency Wave Propagation , 2007 .

[4]  Hailiang Liu,et al.  Recovery of High Frequency Wave Fields from Phase Space-Based Measurements , 2009, Multiscale Model. Simul..

[5]  Jianliang Qian,et al.  Eulerian Gaussian beams for Schrödinger equations in the semi-classical regime , 2009, J. Comput. Phys..

[6]  Mohammad Motamed,et al.  Topics in Analysis and Computation of Linear Wave Propagation , 2008 .

[7]  M. M. Popov,et al.  Computation of wave fields in inhomogeneous media — Gaussian beam approach , 1982 .

[8]  Rolando Magnanini,et al.  On Complex-Valued Solutions to a Two-Dimensional Eikonal Equation. II. Existence Theorems , 2003, SIAM J. Math. Anal..

[9]  M. M. Popov,et al.  Application of the method of summation of Gaussian beams for calculation of high-frequency wave fields , 1981 .

[10]  Shi Jin,et al.  Gaussian beam methods for the Schrodinger equation in the semi-classical regime: Lagrangian and Eulerian formulations , 2008 .

[11]  Luděk Klimeš,et al.  Discretization error for the superposition of Gaussian beams , 1986 .

[12]  Xu Yang,et al.  A numerical study of the Gaussian beam methods for one-dimensional Schrodinger-Poisson equations ⁄ , 2009 .

[13]  Walter Littman,et al.  Studies in partial differential equations , 1982 .

[14]  N. R. Hill,et al.  Prestack Gaussian‐beam depth migration , 2001 .

[15]  Xu Yang,et al.  The Gaussian Beam Methods for Schrodinger-Poisson Equations , 2010 .

[16]  M. Popov.  A new method of computation of wave fields using Gaussian beams , 1982 .

[17]  N. R. Hill,et al.  Gaussian beam migration , 1990 .

[18]  G. Talenti.  On complex-valued solutions to a 2D eikonal equation , 1995 .

[19]  Luděk Klimeš,et al.  Expansion of a high-frequency time-harmonic wavefield given on an initial surface into Gaussian beams , 1984 .

[20]  Radjesvarane Alexandre,et al.  Gaussian beams summation for the wave equation in a convex domain , 2009 .