Gilbert Kenneth E
National Center for Physical Acoustics, University of Mississippi, University, Mississippi 38677, USA.
J Acoust Soc Am. 2016 Mar;139(3):1071-80. doi: 10.1121/1.4942591.
Understanding the physics and mathematics underlying a computational algorithm such as the Green's function parabolic equation (GFPE) is both useful and worthwhile. To this end, the present article aims to give a more widely accessible derivation of the GFPE algorithm than was given originally by Gilbert and Di [(1993). J. Acoust. Soc. Am. 94, 2343-2352]. The present derivation, which uses mathematics familiar to most engineers and physicists, begins with the separation of variables method, a basic and well-known approach for solving partial differential equations. The method leads naturally to eigenvalue-eigenfunction equations. A step-by-step analysis arrives at relatively simple, analytic expressions for the horizontal and vertical eigenfunctions, which are sinusoids plus a surface wave. The eigenfunctions are superposed in an eigenfunction expansion to yield a one-way propagation solution. The one-way solution is generalized to obtain the GFPE algorithm. In addition, and equally important, the eigenfunctions are used to give concrete meaning to abstract operator solutions for one-way acoustic propagation. By using an eigenfunction expansion of the acoustic field, together with an operator solution, one can obtain the GFPE algorithm very directly and concisely.