Numerical simulations and experimental verification of the extent of HIFU-induced tissue necrosis

Local ultrasonic tissue ablation is induced by a rapid (<3s) rise in temperature in a small ellipsoidal volume (about 13mm(3)) inside the tissue to a cytotoxic level when exposed to a high-intensity focused ultrasound (HIFU) beam. The aim of this study was to develop a numerical tool to predict the location and extent of a necrotic lesion formed locally inside the ex vivo tissue as a result of exposure to a single or multiple HIFU beam, ensuring the efficacy and safety of destroying solid tumors. The proposed tool was based on modelling the non-linear propagation of acoustic waves and heat transfer in heterogeneous media using the k-wave toolbox. The wave propagation equations were solved for two-layer (water/tissue) media. The source of the acoustic waves was a spherical bowl-shaped transducer with a resonance frequency of 1.08 MHz. The distribution of heat sources was determined from the calculated acoustic pressure distribution in the HIFU beam. The obtained temperature distributions during heating and cooling allowed calculation of the thermal dose and prediction of the extent of the necrotic lesion. The obtained results of numerical simulations were compared with the experimental data from previous studies. The mean difference between the calculated and measured length or diameter of a single exposure induced necrotic lesion was approximately 1 mm. In the case of a necrotic lesion induced by multiple exposures, the mean difference between the measured and calculated cross-sectional area of the planned necrotic lesion covered with necrosis was approximately 11.2 %.