Simulation of hemispherical transducers for transcranial HIFU treatments using the hybrid angular spectrum approach
© Almquist et al; licensee BioMed Central Ltd. 2015
Published: 30 June 2015
Transcranial focused ultrasound is emerging as a promising treatment modality for a variety of disorders including essential tremor and glioblastoma. However, traditional full-wave methods of simulating ultrasound pressure fields for these treatments, such as FDTD, are computationally intensive. The situation is exacerbated by the large area insonified by the transcranial transducers, which are typically hemispherical to spread the intensity over a large area of the skull and prevent burns. The Hybrid Angular Spectrum (HAS) approach  has been used to rapidly model ultrasonic beam propagation.
Here we present a method for adapting the HAS technique to hemispherical transducers.
The HAS approach assumes that the pressure pattern incident on the front face of the 3D rectangular model is defined on a plane nominally perpendicular to the propagation axis of the transducer. This is problematic for a hemispherical transducer that wraps around the model. To circumvent this limitation, we divide the transducer into seven sections of spatially clustered elements. For each section, the elements, along with the acoustic model, are rotated into a coordinate system with the direction of propagation away from the center of the section. The HAS method is applied using the specific elements of the section, then the resulting pressure pattern is rotated and interpolated back to the original coordinate frame. The fully insonified field is the superposition of the resulting pressure patterns from each section. The resulting field can be used in temperature simulations or to verify the usefulness of phase corrections. Simulations of a 1024-element 650-kHz InSightec ExAblate transducer were carried out using this method on a selected set of retrospective clinical data. The acoustic speed of sound, density and attenuation values for the skull were derived from a CT scan with 0.43 x 0.43 x 1.0-mm resolution using a previously published method of conversion. The simulated skull/brain model contained 512 x 512 x 195 voxels.
Results and conclusions
There is an increase of approximately 25% in peak intensity when phase correction is employed.
The authors thank Eyal Zadicario and Gilat Schiff of InSightec for providing data used in these simulations. This work was funded by NIH grant 5R01EB013433-03 and by the Focused Ultrasound Foundation.
- Vyas U, Christensen D: Ultrasound beam simulations in inhomogeneous tissue geometries using the hybrid angular spectrum method. Trans. UFFC. 2012, 59 (6): 1093-1100.View ArticleGoogle Scholar
- Pichardo S, Sin V, Hynynen K: Multi-frequency characterization of the speed of sound and attenuation coefficient for longitudinal transmission of freshly excised human skulls. Phys. Med. Biol. 2011, 56 (1): 219-10.1088/0031-9155/56/1/014.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.