A fast Green’s function method for one-way sound propagation in the atmosphere

Abstract

A Green's function method is used to derive a fast, general algorithm for one-way wave propagation. The algorithm is applied to outdoor sound propagation. The general method is not limited to atmospheric sound propagation, however, and can be applied to other problems, such as sound propagation in the ocean and electromagnetic wave propagation. The new algorithm, called "GF-PE" (Green's function method for the parabolic equation), reduces to the well-known Fourier split-step algorithm for the parabolic equation (PE) when no boundary conditions are imposed (e.g., at a ground surface). With the GF-PE, range steps many wavelengths long are possible, while with a PE algorithm based on a finite-difference range step, such as the Crank-Nicolson method, the range steps are typically limited to a fraction of a wavelength. Because of its longer range step, the new algorithm is 40-450 times faster than PE algorithms that use the Crank-Nicolson method. For outdoor sound propagation over a locally reacting ground surface, the computed GF-PE field is the sum of three terms: a direct wave, a specularly reflected wave, 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 CrankNicolson method. The GF-PE algorithm is shown to be accurate and approximately two orders of magnitude faster than a PE based on the Crank-Nicolson method. Hence, the new algorithm opens the door to some useful new computational capabilities such as real-time predictions on desktop computers, fast pulse calculations, and practical three-dimensional calculations.

Exploiting Supporting Poles to Increase Road Traffic Noise Shielding of Tree Belts

“…In the other plane (plane number 2, see Fig. 1), the sound-soil interaction is predicted (in absence of trees) using the Green's Function Parabolic Equation (GFPE) method [13] [14], which is computationally much faster, basically since stepping in the propagation direction can be performed at multiples of the wavelength. The latter technique accounts for the ground surface impedance discontinuities (from rigid ground to forest floor, and then from forest floor to grassland) when sound travels from the different road segments towards the receiver.…”

A forward-advancing wave expansion method for numerical solution of large-scale sound propagation problems

“…Gilbert and Di [69] developed a "fast PE m ethod" for computing sound prop agation in the atmosphere, the G reen's function parabolic equation method (G FPE). The G FPE m ethod is more efficient than the Crank-Nicholson PE method.…”

Electromagnetic wave propagation through simulated atmospheric refractivity fields

“…The PE calculations presented in this paper were done with an algorithm originally designed for acoustic propagation and modified for EM propagation [Gilbert and Di, 1993]. The model is called the "Green's function parabolic equation model" (GFPE), and for angle-independent impedance boundary conditions at the Earth's surface it is mathematically equivalent to the PE algorithm of Kuttler and Docker), [1991].…”