Ke. Gilbert et X. Di, A FAST GREEN-FUNCTION METHOD FOR ONE-WAY SOUND-PROPAGATION IN THE ATMOSPHERE, The Journal of the Acoustical Society of America, 94(4), 1993, pp. 2343-2352
A Green's function method is used to derive a fast, general algorithm
for one-way wave propagation. The algorithm is applied to outdoor soun
d propagation. The general method is not limited to atmospheric sound
propagation, however, and can be applied to other problems, such as so
und propagation in the ocean and electromagnetic wave propagation. The
new algorithm, called ''GF-PE'' (Green's function method for the para
bolic equation), reduces to the well-known Fourier split-step algorith
m for the parabolic equation (PE) when no boundary conditions are impo
sed (e.g., at a ground surface). With the GF-PE, range steps many wave
lengths long are possible, while with a PE algorithm based on a finite
-difference range step, such as the Crank-Nicolson method, the range s
teps are typically limited to a fraction of a wavelength. Because of i
ts longer range step, the new algorithm is 40-450 times faster than PE
algorithms that use the Crank-Nicolson method. For outdoor sound prop
agation over a locally reacting ground surface, the computed GF-PE fie
ld is the sum of three terms: a direct wave, a specularly reflected wa
ve, and a surface wave. With the new method, the air-ground impedance
condition is treated exactly and results in an analytic expression for
the surface wave contribution. Numerical results from the GF-PE model
are presented and compared to exact calculations, fast-field program
(FFP) calculations, and PE results computed with the Crank-Nicolson me
thod. The GF-PE algorithm is shown to be accurate and approximately tw
o orders of magnitude faster than a PE based on the Crank-Nicolson met
hod. Hence, the new algorithm opens the door to some useful new comput
ational capabilities such as real-time predictions on desktop computer
s, fast pulse calculations, and practical three-dimensional calculatio
ns.