Simulation of three-dimensional nonlinear fields of ultrasound therapeutic arrays

A novel numerical model was developed to simulate three-dimensional nonlinear fields generated by high intensity focused ultrasound (HIFU) arrays. The model is based on the solution to the Westervelt equation; the developed algorithm makes it possible to model nonlinear pressure fields of periodic waves in the presence of shock fronts localized near the focus. The role of nonlinear effects in a focused beam of a two-dimensional array was investigated in a numerical experiment in water. The array consisting of 256 elements and intensity range on the array elements of up to 10 W/cm2 was considered. The results of simulations have shown that for characteristic intensity outputs of modern HIFU arrays, nonlinear effects play an important role and shock fronts develop in the pressure waveforms at the focus.

[1]  L. Gavrilov,et al.  A theoretical assessment of the relative performance of spherical phased arrays for ultrasound surgery , 2000, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[2]  E. Tadmor,et al.  New High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection—Diffusion Equations , 2000 .

[3]  Kullervo Hynynen,et al.  Pre-clinical testing of a phased array ultrasound system for MRI-guided noninvasive surgery of the brain--a primate study. , 2006, European journal of radiology.

[4]  Robin O Cleveland,et al.  Modeling the propagation of nonlinear three-dimensional acoustic beams in inhomogeneous media. , 2007, The Journal of the Acoustical Society of America.

[5]  Lawrence A Crum,et al.  Shock-induced heating and millisecond boiling in gels and tissue due to high intensity focused ultrasound. , 2010, Ultrasound in medicine & biology.

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

[7]  Oleg V. Rudenko,et al.  Theoretical Foundations of Nonlinear Acoustics , 1977 .

[8]  Mathieu Pernot,et al.  3-D real-time motion correction in high-intensity focused ultrasound therapy. , 2004, Ultrasound in medicine & biology.

[9]  Vera A. Khokhlova,et al.  Numerical modeling of finite-amplitude sound beams: Shock formation in the near field of a cw plane piston source , 2001 .

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

[11]  J. Tavakkoli,et al.  Modeling of nonlinear ultrasound propagation in tissue from array transducers. , 2003, The Journal of the Acoustical Society of America.

[12]  M Pernot,et al.  High power transcranial beam steering for ultrasonic brain therapy. , 2003, Physics in medicine and biology.

[13]  Adam Shaw,et al.  Focusing of high-intensity ultrasound through the rib cage using a therapeutic random phased array. , 2010, Ultrasound in medicine & biology.

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

[15]  Algorithm for the nonlinear propagation of acoustic beams from phased arrays and nonplanar sources , 2006, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[16]  C. Moonen,et al.  A method for MRI guidance of intercostal high intensity focused ultrasound ablation in the liver. , 2010, Medical physics.

[17]  J W Hand,et al.  A random phased array device for delivery of high intensity focused ultrasound , 2009, Physics in medicine and biology.

[18]  K. Parker,et al.  New approaches to the linear propagation of acoustic fields. , 1991, The Journal of the Acoustical Society of America.

[19]  Lawrence A Crum,et al.  Acoustic characterization of high intensity focused ultrasound fields: a combined measurement and modeling approach. , 2008, The Journal of the Acoustical Society of America.