Modeling of nonlinear ultrasound propagation in tissue from array transducers

A computationally efficient model capable of simulating finite-amplitude ultrasound beam propagation in water and in tissue from phased linear arrays and other transducers of arbitrary quasiplanar geometry is described. It is based on a second-order operator splitting approach [Tavakkoli et al., J. Acoust. Soc. Am. 104, 2061-2072 (1998)], with a fractional step-marching scheme, whereby the effects of diffraction, attenuation, and nonlinearity can be computed independently over incremental steps. This approach is an extension to that of Christopher and Parker [J. Acoust. Soc. Am. 90, 507-521; 90, 488-499 (1991)], wherein linear and nonlinear effects are propagated separately over incremental steps, and the computation of the diffractive substeps are based on an angular spectrum technique with a modified sampling scheme for accurate and efficient implementation of diffractive propagation from nonradially symmetric sources. Results of the model are compared with published data. Predicted field profiles for nonlinear propagation in tissue from realistic array transducers using the pulse inversion method are presented.