Tracing Three-Dimensional Acoustic Wavefronts in Inhomogeneous, Moving Media

Abstract

We present a numerical implementation of an alternative formulation of the geometrical, or ray, acoustics, where wavefronts rather than rays are the primary objects. Rays are recovered as a by-product of wavefront tracing. The alternative formulation of the geometrical acoustics is motivated, first, by the observation that wavefronts are often more stable than rays at long-range sound propagation, and, second, by a need for computationally efficient modeling of high-frequency acoustic fields in three-dimensionally inhomogeneous, moving or motionless fluids. Wavefronts are found as a finite-difference solution to a system of partial differential equations, which is equivalent to the eikonal equation and is a direct implementation of the intuitive Huygens' wavefront construction. The finite-difference algorithm is an extension of the approach originally developed in the framework of an open source Madagascar project. Benchmark problems, which admit exact, analytic solutions of the eikonal equation, are formulated and utilized to verify the finite-difference wavefront tracing algorithm. Huygens' wavefront tracing (HWT) is applied to modeling sound propagation in three-dimensionally inhomogeneous ocean and atmosphere.