New approaches to the linear propagation of acoustic fields

New algorithms are described that provide insight into linear field propagation and offer significant reductions in computational complexity. The developments presented here include the usage of a recently developed discrete Hankel transform to implement two single step, planar propagation algorithms for baffled, radially symmetric, acoustic pressure or velocity fields; an update on the single step approaches that reduce computational complexity through geometrically determined spatial frequency limitations; and algorithms for extending to multistep propagation. Two equivalent means of introducing arbitrary medium attenuation into the above schemes are presented. Finally, a planar boundary crossing algorithm that accounts for refraction and reflection (but not multiple reflections) is added to one of the multistep propagating algorithms. The resulting algorithm is then used to examine the differences between the corresponding fields of a focused piston source operating in water and in a layered fat/liver (biomedical imaging) medium. The results yield computationally efficient algorithms that can be used for linear propagation of focused or unfocused beams in attenuating, multilayer media, and also provide the basis for a novel nonlinear propagation algorithm.