Fast multiscale Gaussian beam methods for wave equations in bounded convex domains

Motivated by fast multiscale Gaussian wavepacket transforms and multiscale Gaussian beam methods which were originally designed for pure initial-value problems of wave equations, we develop fast multiscale Gaussian beam methods for initial boundary value problems of wave equations in bounded convex domains in the high frequency regime. To compute the wave propagation in bounded convex domains, we have to take into account reflecting multiscale Gaussian beams, which are accomplished by enforcing reflecting boundary conditions during beam propagation and carrying out suitable reflecting beam summation. To propagate multiscale beams efficiently, we prove that the ratio of the squared magnitude of beam amplitude and the beam width is roughly conserved, and accordingly we propose an effective indicator to identify significant beams. We also prove that the resulting multiscale Gaussian beam methods converge asymptotically. Numerical examples demonstrate the accuracy and efficiency of the method.

[1]  Dongming Wei,et al.  Eulerian Gaussian beam method for high frequency wave propagation in heterogeneous media with discontinuities in one direction , 2012 .

[2]  James Ralston,et al.  Mountain Waves and Gaussian Beams , 2007, Multiscale Model. Simul..

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

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

[5]  Shi Jin,et al.  Gaussian beam methods for the Schrödinger equation with discontinuous potentials , 2014, J. Comput. Appl. Math..

[6]  Jianliang Qian,et al.  The backward phase flow and FBI-transform-based Eulerian Gaussian beams for the Schrödinger equation , 2010, J. Comput. Phys..

[7]  Chunxiong Zheng,et al.  Gaussian beam formulations and interface conditions for the one-dimensional linear Schrödinger equation , 2011 .

[8]  Lexing Ying,et al.  Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation , 2010, J. Comput. Phys..

[9]  Gary Cohen.  Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[10]  Victor P. Maslov,et al.  The complex WKB method for nonlinear equations I , 1994 .

[11]  Qing Huo Liu,et al.  Higher-Order Numerical Methods for Transient Wave Equations , 2003 .

[12]  Björn Engquist,et al.  Gaussian beam decomposition of high frequency wave fields using expectation-maximization , 2011, J. Comput. Phys..

[13]  Hailiang Liu,et al.  Error estimates for Gaussian beam superpositions , 2010, Math. Comput..

[14]  Giovanni Russo,et al.  The Gaussian wave packet transform: Efficient computation of the semi-classical limit of the Schrödinger equation. Part 1 - Formulation and the one dimensional case , 2013, J. Comput. Phys..

[15]  Benjamin S. White,et al.  Gaussian wave packets in inhomogeneous media with curved interfaces , 1987, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[16]  Giovanni Russo,et al.  The Gaussian wave packet transform: Efficient computation of the semi-classical limit of the Schrödinger equation. Part 2. Multidimensional case , 2014, J. Comput. Phys..

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

[18]  HEINZ-OTTO KREISS,et al.  A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data , 2005, SIAM J. Sci. Comput..

[19]  Lexing Ying,et al.  Fast Multiscale Gaussian Wavepacket Transforms and Multiscale Gaussian Beams for the Wave Equation , 2010, Multiscale Model. Simul..

[20]  March,et al.  Gaussian beams and the propagation of singularities , 2014 .

[21]  Lexing Ying,et al.  A Convergent Multiscale Gaussian-Beam Parametrix for the Wave Equation , 2013 .

[22]  Evans M. Harrell.  The Complex WKB Method for Nonlinear Equations 1: Linear Theory (Victor P. Maslov) , 1996, SIAM Rev..

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

[24]  Olof Runborg,et al.  Taylor expansion and discretization errors in Gaussian beam superposition , 2009, 0908.3416.

[25]  Alden Waters.  A Parametrix Construction for the Wave Equation with Low Regularity Coefficients Using a Frame of Gaussians , 2010, 1003.2747.

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

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