Gaussian Beam Summation for Diffraction in Inhomogeneous Media Based on the Grid Based Particle Method
Shingyu Leung 1*, Hongkai Zhao 21 Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong.
2 Department of Mathematics, University of California at Irvine, Irvine, CA 92697-3875, USA.
Received 19 August 2009; Accepted (in revised version) 9 February 2010
Available online 17 May 2010
We develop an efficient numerical method to compute single slit or double slit diffraction patterns from high frequency wave in inhomogeneous media. We approximate the high frequency asymptotic solution to the Helmholtz equation using the Eulerian Gaussian beam summation proposed in [20, 21]. The emitted rays from a slit are embedded in the phase space using an open segment. The evolution of this open curve is accurately computed using the recently developed Grid Based Particle Method  which results in a very efficient computational algorithm. Following the grid based particle method we proposed in [23, 24], we represent the open curve or the open surface by meshless Lagrangian particles sampled according to an underlying fixed Eulerian mesh. The end-points of the open curve are tracked explicitly and consistently with interior particles. To construct the overall wavefield, each of these sampling particles also carry necessary quantities that are obtained by solving advection-reaction equations. Numerical experiments show that the resulting method can model diffraction patterns in inhomogeneous media accurately, even in the occurrence of caustics.AMS subject classifications: 34E05, 78A05, 78A45
Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Dynamic interface, Gaussian beams, diffraction, high frequency asymptotic solution.
Email: email@example.com (S. Leung), firstname.lastname@example.org (H. Zhao)