On the use of higher-order azimuthal schemes in 3-D PE modeling

In this paper, the issue of using higher-order finite difference schemes to handle the azimuthal derivative term in a three-dimensional parabolic equation based model is addressed. The three-dimensional penetrable wedge benchmark problem is chosen to illustrate the accuracy and efficiency of the proposed schemes. Both point source and modal initializations of the pressure field are considered. For each higher-order finite difference scheme used in azimuth, the convergence of the numerical solution with respect to the azimuth is investigated and the CPU times are given. Some comparisons with solutions obtained from another 3-D model [J. A. Fawcett, J. Acoust. Soc. Am. 93, 2627-2632 (1993)] are presented. The numerical simulations show that the use of a higher-order scheme in azimuth allows one to reduce the required number of points in the azimuthal direction while still obtaining accurate solutions. The higher-order schemes have approximately the same efficiency as a FFT-based approach (in fact, may outperform it slightly); however, the finite difference approach has the advantage that it may be more flexible than the FFT approach for various PE approximations.