Babich’s Expansion and High-Order Eulerian Asymptotics for Point-Source Helmholtz Equations

The usual geometrical-optics expansion of the solution for the Helmholtz equation of a point source in an inhomogeneous medium yields two equations: an eikonal equation for the traveltime function, and a transport equation for the amplitude function. However, two difficulties arise immediately: one is how to initialize the amplitude at the point source as the wavefield is singular there; the other is that in even-dimension spaces the usual geometrical-optics expansion does not yield a uniform asymptotic approximation close to the source. Babich (USSR Comput Math Math Phys 5(5):247–251, 1965) developed a Hankel-based asymptotic expansion which can overcome these two difficulties with ease. Starting from Babich’s expansion, we develop high-order Eulerian asymptotics for Helmholtz equations in inhomogeneous media. Both the eikonal and transport equations are solved by high-order Lax–Friedrichs weighted non-oscillatory (WENO) schemes. We also prove that fifth-order Lax–Friedrichs WENO schemes for eikonal equations are convergent when the eikonal is smooth. Numerical examples demonstrate that new Eulerian high-order asymptotic methods are uniformly accurate in the neighborhood of the source and away from it.

[1]  Jean-David Benamou,et al.  An Introduction to Eulerian Geometrical Optics (1992–2002) , 2003, J. Sci. Comput..

[2]  Robert J. Plemmons,et al.  Nonnegative Matrices in the Mathematical Sciences , 1979, Classics in Applied Mathematics.

[3]  R. LeVeque.  Finite Difference Methods for Ordinary and Partial Differential Equations: Steady-State and Time-Dependent Problems (Classics in Applied Mathematics Classics in Applied Mathemat) , 2007 .

[4]  Chi-Wang Shu,et al.  Fast Sweeping Fifth Order WENO Scheme for Static Hamilton-Jacobi Equations with Accurate Boundary Treatment , 2010, J. Sci. Comput..

[5]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[6]  Jianliang Qian,et al.  Fast Huygens sweeping methods for Helmholtz equations in inhomogeneous media in the high frequency regime , 2014, J. Comput. Phys..

[7]  Jianliang Qian,et al.  A stopping criterion for higher-order sweeping schemes for static Hamilton-Jacobi equations , 2010 .

[8]  Jianliang Qian,et al.  A Slowness Matching Eulerian Method for Multivalued Solutions of Eikonal Equations , 2003, J. Sci. Comput..

[9]  Gilbert Strang,et al.  Accurate partial difference methods , 1964 .

[10]  Jesse Freeman,et al.  in Morse theory, , 1999 .

[11]  Linbin Zhang,et al.  Eikonal Solver In the Celerity Domain , 2002 .

[12]  Joseph B. Keller,et al.  The high‐frequency asymptotic field of a point source in an inhomogeneous medium , 1963 .

[13]  Hongkai Zhao,et al.  High Order Fast Sweeping Methods for Static Hamilton–Jacobi Equations , 2006, J. Sci. Comput..

[14]  Hongkai Zhao,et al.  Higher-order schemes for 3D first-arrival traveltimes and amplitudes , 2012 .

[15]  J. Ortega.  Matrix Theory: A Second Course , 1987 .

[16]  Randall J. LeVeque,et al.  Finite difference methods for ordinary and partial differential equations - steady-state and time-dependent problems , 2007 .

[17]  V. M. Babich.  The short wave asymptotic form of the solution for the problem of a point source in an inhomogeneous medium , 1965 .

[18]  S. Osher,et al.  Lax-Friedrichs sweeping scheme for static Hamilton-Jacobi equations , 2004 .

[19]  Hongkai Zhao,et al.  Fast sweeping method for the factored eikonal equation , 2009, J. Comput. Phys..

[20]  Jianliang Qian,et al.  Fast Sweeping Methods for Factored Anisotropic Eikonal Equations: Multiplicative and Additive Factors , 2012, J. Sci. Comput..

[21]  A. Pica.  Fast And Accurate Finite-difference Solutions of the 3D Eikonal Equation Parametrized In Celerity , 1997 .

[22]  Jianliang Qian,et al.  Factored singularities and high-order Lax-Friedrichs sweeping schemes for point-source traveltimes and amplitudes , 2011, J. Comput. Phys..

[23]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[24]  Jianliang Qian,et al.  An adaptive finite-difference method for traveltimes and amplitudes , 2002 .

[25]  John E. Vidale,et al.  Rapid calculation of seismic amplitudes , 1990 .

[26]  Jianliang Qian,et al.  High-Order Factorization Based High-Order Hybrid Fast Sweeping Methods for Point-Source Eikonal Equations , 2014, SIAM J. Numer. Anal..

[27]  E. Coddington,et al.  Theory of Ordinary Differential Equations , 1955 .

[28]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..