LAUSR

We present asymptotic methods for solving high frequency Helmholtz equations in anisotropic media. The methods are motivated by Babich’s expansion that uses Hankel functions of the first kind to approximate the solution of high frequency Helmholtz equation in isotropic media. Within Babich’s expansion, we can derive the anisotropic eikonal equation and a recurrent system of transport equations to determine the phase and amplitude terms of the wave, respectively. In order to reconstruct the wave with the phase and amplitude terms for any high frequencies, they must be computed with high-order accuracy, for which a high-order factorization approach based on power series expansions at the primary source is applied first to resolve the source singularities, after that high-order schemes can be implemented efficiently. Rigorous formulations are derived, and numerical examples are presented to demonstrate the methods.