Fast simulation of non-linear pulsed ultrasound fields using an angular spectrum approach.

A fast non-linear pulsed ultrasound field simulation is presented. It is implemented based on an angular spectrum approach (ASA), which analytically solves the non-linear wave equation. The ASA solution to the Westervelt equation is derived in detail. The calculation speed is significantly increased compared to a numerical solution using an operator splitting method (OSM). The ASA has been modified and extended to pulsed non-linear ultrasound fields in combination with Field II, where any array transducer with arbitrary geometry, excitation, focusing and apodization can be simulated. The accuracy of the non-linear ASA is compared to the non-linear simulation program - Abersim, which is a numerical solution to the Burgers equation based on the OSM. Simulations are performed for a linear array transducer with 64 active elements, focus at 40mm, and excitation by a 2-cycle sine wave with a center frequency of 5MHz. The speed is increased approximately by a factor of 140 and the calculation time is 12min with a standard PC, when simulating the second harmonic pulse at the focal point. For the second harmonic point spread function the full width error is 1.5% at -6dB and 6.4% at -12dB compared to Abersim.

[1]  C.E. Shannon,et al.  Communication in the Presence of Noise , 1949, Proceedings of the IRE.

[2]  Jørgen Arendt Jensen,et al.  Medical ultrasound imaging. , 2007, Progress in biophysics and molecular biology.

[3]  F Patat,et al.  Time-domain modeling of nonlinear distortion of pulsed finite amplitude sound beams. , 2000, Ultrasonics.

[4]  Xiang Yan,et al.  Angular Spectrum Decomposition Analysis of Second Harmonic Ultrasound Propagation and its Relation to Tissue Harmonic Imaging , 2006 .

[5]  P. Bloomfield,et al.  Wave envelopes method for description of nonlinear acoustic wave propagation. , 2006, Ultrasonics.

[6]  Lawrence E. Kinsler,et al.  Fundamentals of acoustics , 1950 .

[7]  Peter A. Lewin,et al.  Transducer characterization using the angular spectrum method , 1989 .

[8]  K. Boone,et al.  Effect of skin impedance on image quality and variability in electrical impedance tomography: a model study , 1996, Medical and Biological Engineering and Computing.

[9]  M. Hamilton,et al.  Second-harmonic generation in sound beams reflected from, and transmitted through, immersed elastic solids. , 2001, The Journal of the Acoustical Society of America.

[10]  Jørgen Arendt Jensen,et al.  Ultrasound fields in an attenuating medium , 1993 .

[11]  Xiang Yan,et al.  Statistical model of beam distortion by tissue inhomogeneities in tissue harmonic imaging , 2004 .

[12]  M. Hamilton,et al.  Time‐domain modeling of pulsed finite‐amplitude sound beams , 1995 .

[13]  Prediction of ultrasonic field propagation through layered media using the extended angular spectrum method. , 1994, Ultrasound in medicine & biology.

[14]  T. Varslot,et al.  Computer simulation of forward wave propagation in soft tissue , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[15]  R. Bracewell.  The Fourier Transform and Its Applications , 1966 .

[16]  P. Lewin,et al.  Fast prediction of pulsed nonlinear acoustic fields from clinically relevant sources using time-averaged wave envelope approach: comparison of numerical simulations and experimental results. , 2008, Ultrasonics.

[17]  J. Goodman.  Introduction to Fourier optics , 1969 .

[18]  J. Burgers.  A mathematical model illustrating the theory of turbulence , 1948 .

[19]  R Souchon,et al.  Modeling of pulsed finite-amplitude focused sound beams in time domain. , 1998, The Journal of the Acoustical Society of America.

[20]  M. Averkiou,et al.  Tissue harmonic imaging , 2000, 2000 IEEE Ultrasonics Symposium. Proceedings. An International Symposium (Cat. No.00CH37121).

[21]  P. J. Westervelt.  Parametric Acoustic Array , 1963 .

[22]  E. Zabolotskaya,et al.  Quasi-plane waves in the nonlinear acoustics of confined beams , 1969 .

[23]  C. Chin,et al.  Pulse inversion Doppler: a new method for detecting nonlinear echoes from microbubble contrast agents , 1999, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[24]  Louise Kold Taylor,et al.  Modeling of higher harmonics formation in medical ultrasound systems , 2002, SPIE Medical Imaging.

[25]  Yigang Du,et al.  Feasibility of non-linear simulation for Field II using an angular spectrum approach , 2008, 2008 IEEE Ultrasonics Symposium.

[26]  Jens E. Wilhjelm,et al.  Simulation of non-linear ultrasound fields , 2002 .

[27]  M. Frijlink,et al.  Abersim: A simulation program for 3D nonlinear acoustic wave propagation for arbitrary pulses and arbitrary transducer geometries , 2008, 2008 IEEE Ultrasonics Symposium.

[28]  Jørgen Arendt Jensen,et al.  Users' guide for the Field II program , 2011 .

[29]  J. Jensen,et al.  Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers , 1992, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[30]  A. Polyanin.  Handbook of Linear Partial Differential Equations for Engineers and Scientists , 2001 .