Difference equation technique for determining the evolution of radially symmetric ultrasonic fields

A simplified difference equation technique has been developed to calculate the three-dimensional evolution of radially symmetric, complex ultrasonic fields. The parabolic approximation to the Helmholtz equation is obtained by assuming that localized changes in the magnitude of the field with respect to distance along the axis of propagation are small compared to those with respect to radial displacement. The Helmholtz equation in cylindrical coordinates may then be approximated, after separating variables, using truncated series expansions. The resulting difference equation determines complex field values at discrete locations in a plane using corresponding adjacent locations in the preceding plane. Stability criteria which determine acceptable values for various parameters are discussed. The technique has been applied to the propagation of ultrasonic fields generated by two-dimensional focusing and non-focusing, uniform and Gaussian profile, piezoelectric transducers in both profiled and non-profiled media.