New approaches to nonlinear diffractive field propagation.

In many domains of acoustic field propagation, such as medical ultrasound imaging, lithotripsy shock treatment, and underwater sonar, a realistic calculation of beam patterns requires treatment of the effects of diffraction from finite sources. Also, the mechanisms of loss and nonlinear effects within the medium are typically nonnegligible. The combination of diffraction, attenuation, and nonlinear effects has been treated by a number of formulations and numerical techniques. A novel model that incrementally propagates the field of baffled planar sources with substeps that account for the physics of diffraction, attenuation, and nonlinearity is presented. The model accounts for the effect of refraction and reflection (but not multiple reflections) in the case of propagation through multiple, parallel layers of fluid medium. An implementation of the model for axis symmetric sources has been developed. In one substep of the implementation, a new discrete Hankel transform is used with spatial transform techniques to propagate the field over a short distance with diffraction and attenuation. In the other substep, the temporal frequency domain solution to Burgers' equation is implemented to account for the nonlinear accretion and depletion of harmonics. This approach yields a computationally efficient procedure for calculating beam patterns from a baffled planar, axially symmetric source under conditions ranging from quasilinear through shock. The model is not restricted by the usual parabolic wave approximation and the field's directionality is explicitly accounted for at each point. Useage of a harmonic-limiting scheme allows the model to propagate some previously intractable high-intensity nonlinear fields. Results of the model are shown to be in excellent agreement with measurements performed on the nonlinear field of an unfocused 2.25-MHz piston source, even in the near field where the established parabolic wave approximation model fails. Next, the model is used to compare the water path and in situ fields of a medical ultrasound device. Finally, the model is used to calculate the spatial heating rate associated with a nonlinear field and to simulate the phenomenon of saturation-induced beam broadening.