Experimental validation of acoustic and thermal modeling in heterogeneous phantoms using the hybrid angular spectrum method

Evaluation of acoustic-thermal simulations of in vivo magnetic resonance guided focused ultrasound ablative therapy

PURPOSE

To evaluate numerical simulations of focused ultrasound (FUS) with a rabbit model, comparing simulated heating characteristics with magnetic resonance temperature imaging (MRTI) data collected during in vivo treatment.

METHODS

A rabbit model was treated with FUS sonications in the biceps femoris with 3D MRTI collected. Acoustic and thermal properties of the rabbit muscle were determined experimentally. Numerical models of the rabbits were created, and tissue-type-specific properties were assigned. FUS simulations were performed using both the hybrid angular spectrum (HAS) method and k-Wave. Simulated power deposition patterns were converted to temperature maps using a Pennes' bioheat equation-based thermal solver. Agreement of pressure between the simulation techniques and temperature between the simulation and experimental heating was evaluated. Contributions of scattering and absorption attenuation were considered.

RESULTS

Simulated peak pressures derived using the HAS method exceeded the simulated peak pressures from k-Wave by 1.6 ± 2.7%. The location and FWHM of the peak pressure calculated from HAS and k-Wave showed good agreement. When muscle acoustic absorption value in the simulations was adjusted to approximately 54% of the measured attenuation, the average root-mean-squared error between simulated and experimental spatial-average temperature profiles was 0.046 ± 0.019 °C/W. Mean distance between simulated and experimental COTMs was 3.25 ± 1.37 mm. Transverse FWHMs of simulated sonications were smaller than in in vivo sonications. Longitudinal FWHMs were similar.

CONCLUSIONS

Presented results demonstrate agreement between HAS and k-Wave simulations and that FUS simulations can accurately predict focal position and heating for in vivo applications in soft tissue.

The influence of bone model geometries on the determination of skull acoustic properties

In this study, we investigated the impact of various simulated skull bone geometries on the determination of skull speed of sound and acoustic attenuation values via optimization using transmitted pressure amplitudes beyond the bone. Using the hybrid angular spectrum method (HAS), we simulated ultrasound transmission through four model sets of different geometries involving sandwiched layers of diploë and cortical bone in addition to three models generated from CT images of ex-vivo human skull-bones. We characterized cost-function solution spaces for each model and, using optimization, found that when a model possessed appreciable variations in resolvable layer thickness, the predefined attenuation coefficients could be found with low error (RMSE < 0.01 Np/cm). However, we identified a spatial frequency cutoff in the models' geometry beyond which the accuracy of the property determination begins to fail, depending on the frequency of the ultrasound source. There was a large increase in error of the attenuation coefficients determined by the optimization when the variations in layer thickness were above the identified spatial frequency cutoffs, or when the lateral variations across the model were relatively low in amplitude. For our limited sample of three CT-image derived bone models, the attenuation coefficients were determined successfully. The speed of sound values were determined with low error for all models (including the CT-image derived models) that were tested (RMSE < 0.4 m/s). These results illustrate that it is possible to determine the acoustic properties of two-component models when the internal bone structure is taken into account and the structure satisfies the spatial frequency constraints discussed.