Simulation of thermal ablation by high-intensity focused ultrasound with temperature-dependent properties

Real-Time Reconstruction of HIFU Focal Temperature Field Based on Deep Learning

Objective and Impact Statement: High-intensity focused ultrasound (HIFU) therapy is a promising noninvasive method that induces coagulative necrosis in diseased tissues through thermal and cavitation effects, while avoiding surrounding damage to surrounding normal tissues. Introduction: Accurate and real-time acquisition of the focal region temperature field during HIFU treatment marked enhances therapeutic efficacy, holding paramount scientific and practical value in clinical cancer therapy. Methods: In this paper, we initially designed and assembled an integrated HIFU system incorporating diagnostic, therapeutic, and temperature measurement functionalities to collect ultrasound echo signals and temperature variations during HIFU therapy. Furthermore, we introduced a novel multimodal teacher-student model approach, which utilizes the shared self-expressive coefficients and the deep canonical correlation analysis layer to aggregate each modality data, then through knowledge distillation strategies, transfers the knowledge from the teacher model to the student model. Results: By investigating the relationship between the phantoms, in vitro, and in vivo ultrasound echo signals and temperatures, we successfully achieved real-time reconstruction of the HIFU focal 2D temperature field region with a maximum temperature error of less than 2.5 °C. Conclusion: Our method effectively monitored the distribution of the HIFU temperature field in real time, providing scientifically precise predictive schemes for HIFU therapy, laying a theoretical foundation for subsequent personalized treatment dose planning, and providing efficient guidance for noninvasive, nonionizing cancer treatment.

Evaluation of the dual-frequency transducer for controlling thermal ablation morphology using frequency shift keying signal

PURPOSE

The catheter-based ultrasound (CBUS) can reach the target tissue directly and achieve rapid treatment. The frequency shift keying (FSK) signal is proposed to regulate and evaluate tumor ablation by a miniaturized dual-frequency transducer.

METHODS

A dual-frequency transducer prototype (3 × 7 × 0.4 mm) was designed and fabricated for the CBUS applicator (OD: 3.8 mm) based on the fundamental frequency of 5.21 MHz and the third harmonic frequency of 16.88 MHz. Then, the acoustic fields and temperature field distributions using the FSK signals (with 0, 25, 50, 75, and 100% third harmonic frequency duty ratios) were simulated by finite element analysis. Finally, tissue ablation and temperature monitoring were performed in phantom and ex vivo tissue, respectively.

RESULTS

At the same input electrical power (20 W), the output acoustic power of the fundamental frequency of the transducer was 10.03 W (electroacoustic efficiencies: 50.1%), and that of the third harmonic frequency was 6.19 W (30.6%). As the third harmonic frequency duty ratios increased, the shape of thermal lesions varied from strip to droplet in simulated and phantom experimental results. The same trend was observed in ex vivo tests.

CONCLUSION

Dual-frequency transducers excited by the FSK signal can control the morphology of lesions.

SIGNIFICANCE

The acoustic power deposition of CBUS was optimized to achieve precise ablation.

Thermometry using entropy imaging of ultrasound radio frequency signal time series

Low intensity focused ultrasound (LIFU) is a novel approach that could activate drug release and considerably improve the delivery of anticancer drug. LIFU treatment has some features like is able to penetrate deep into the tissue and being non-invasive, as a consequence LIFU displays great capability for controlling the drug release and improving the chemotherapy treatment efficiency. The goal of this study is to research the feasibility of the entropy parameter of RF time series of ultrasound backscattered signals for measuring the changes in temperature induced by a LIFU device. Entropy Imaging is a technique for reconstructing ultrasound images based on the average uncertainty of time-series in a signal. Furthermore, the Shannon Entropy can quantify the uncertainty of a random process and is usually used as a measure for the information content of probability distributions. In this study, we use the Entropy Imaging method for measuring the LIFU-induced temperature changes in the deep region of ex vivo porcine tissue samples. The results obtained show that the changes of entropy parameter of RF time series signal are proportional to temperature changes recorded by a calibrated thermocouple in the temperature range of 37–47°C. In conclusion, in this study we show that Shannon entropy of RF time series signal possesses promising features like succinctly capturing the available information in a system by considering the uncertainty in a given data that can be used, as a new method, to measure temperature changes non-invasively and quantitatively in the deep region of tissue.

Numerical Study of Bubble Cloud and Thermal Lesion Evolution During Acoustic Droplet Vaporization Enhanced HIFU Treatment

Acoustic droplet vaporization (ADV) has been proven to enhance high intensity focused ultrasound (HIFU) thermal ablation of tumor. It has also been demonstrated that triggering droplets before HIFU exposure could be a potential way to control both the size and the shape of the thermal lesion. In this paper, a numerical model is proposed to predict the thermal lesion created in ADV enhanced HIFU treatment. Bubble oscillation was coupled into a viscoelastic medium in the model to more closely represent real applications in tissues. Several physical processes caused by continuous wave ultrasound and elevated temperature during the HIFU exposure were considered, including rectified diffusion, gas solubility variation with temperature in the medium, and boiling. Four droplet concentrations spanning two orders of magnitude were calculated. The bubble cloud formed from triggering of the droplets by the pulse wave ultrasound, along with the evolution of the shape and location of the bubble cloud and thermal lesion during the following continuous wave exposure was obtained. The increase of bubble void fraction caused by continuous wave exposure was found to be consistent with the experimental observation. With the increase of droplet concentration, the predicted bubble cloud shapes vary from tadpole to triangular and double triangular, while the thermal lesions move toward the transducer. The results show that the assumptions used in this model increased the accuracy of the results. This model may be used for parametrical study of ADV enhanced HIFU treatment and be further used for treatment planning and optimization in the future.

Influence of temperature-dependent acoustic and thermal parameters and nonlinear harmonics on the prediction of thermal lesion under HIFU ablation

According to the traditional method of high intensity focused ultrasound (HIFU) treatment, the acoustic and thermal characteristic parameters of constant temperature (room temperature or body temperature) are used to predict thermal lesion. Based on the nonlinear spherical beam equation (SBE) and Pennes bio-heat transfer equation, and a new acoustic-thermal coupled model is proposed. The constant and temperature-dependent acoustic and thermal characteristic parameters are used to predict thermal lesion, and the predicted lesion area are compared with each other. Moreover, the relationship between harmonic amplitude ratio (P₂/P₁) and thermal lesion is studied. Combined with the known experimental data of acoustic and thermal characteristic parameters of biological tissue and data fitting method, the relationship between acoustic and thermal characteristic parameters and temperature is obtained; and the thermal lesion simulation calculation is carried out by using the acoustic and thermal characteristic parameters under constant temperature and temperature- dependent acoustic and thermal characteristic parameters, respectively. The simulation results show that under the same irradiation condition, the thermal lesion predicted by temperature-dependent acoustic and thermal characteristic parameters is larger than that predicted by traditional method, and the thermal lesion increases with the decrease of harmonic amplitude ratio.

In Vivo Non-Thermal, Selective Cancer Treatment With High-Frequency Medium-Intensity Focused Ultrasound

Focused ultrasound (FUS) has proven its efficacy in non-invasive, radiation-free cancer treatment. However, the commonly used low-frequency high-intensity focused ultrasound (HIFU) destroys both cancerous and healthy tissues non-specifically through extreme heat and inertial cavitation with low spatial resolution. To address this issue, we evaluate the therapeutic effects of pulsed (60 Hz pulse repetition frequency, 1.45 ms pulse width) high-frequency (20.7 MHz) medium-intensity (spatial-peak pulse-average intensity I_SPPA < 279.1 W/cm², spatial-peak temporal-average intensity I_SPTA < 24.3 W/cm²) focused ultrasound (pHFMIFU) for selective cancer treatment without thermal damage and with low risk of inertial cavitation (mechanical index < 0.66), in an in vivo subcutaneous B16F10 melanoma tumor growth model in mice. The pHFMIFU with 104 μm focal diameter is generated by a microfabricated self-focusing acoustic transducer (SFAT) with a Fresnel acoustic lens. A three-axis positioning system has been developed for automatic scanning of the transducer to cover a larger treatment volume, while a water-cooling system is custom-built for dissipating non-acoustic heat from the transducer surface. Initial testing revealed that pHFMIFU treatment can be applied to a living animal while maintaining skin temperature under 35.6 °C without damaging normal skin and tissue. After eleven days of treatment with pHFMIFU, the treated tumors were significantly smaller with large areas of necrosis and apoptosis in the treatment field compared to untreated controls. Potential mechanisms of this selective, non-thermal killing effect, as well as possible causes of and solutions to the variation in treatment results, have been analyzed and proposed. The pHFMIFU could potentially be used as a new therapeutic modality for safer cancer treatment especially in critical body regions, due to its cancer-specific effects and high spatial resolution.

Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation

BACKGROUND AND OBJECTIVES

During thermal heating surgical procedures such as electrosurgery, thermal ablative treatment and hyperthermia, soft tissue deformation due to surgical tool-tissue interaction and patient movement can affect the distribution of thermal energy induced. Soft tissue temperature must be obtained from the deformed tissue for precise delivery of thermal energy. However, the classical Pennes bio-heat transfer model can handle only the static non-moving state of tissue. In addition, in order to enable a surgeon to visualise the simulated results immediately, the solution procedure must be suitable for real-time thermal applications.

METHODS

This paper presents a formulation of bio-heat transfer under the effect of soft tissue deformation for fast or near real-time tissue temperature prediction, based on fast explicit dynamics finite element algorithm (FED-FEM) for transient heat transfer. The proposed thermal analysis under deformation is achieved by transformation of the unknown deformed tissue state to the known initial static state via a mapping function. The appropriateness and effectiveness of the proposed formulation are evaluated on a realistic virtual human liver model with blood vessels to demonstrate a clinically relevant scenario of thermal ablation of hepatic cancer.

RESULTS

For numerical accuracy, the proposed formulation can achieve a typical 10^-3 level of normalised relative error at nodes and between 10^-4 and 10^-5 level of total errors for the simulation, by comparing solutions against the commercial finite element analysis package. For computation time, the proposed formulation under tissue deformation with anisotropic temperature-dependent properties consumes 2.518 × 10^-4 ms for one element thermal loads computation, compared to 2.237 × 10^-4 ms for the formulation without deformation which is 0.89 times of the former. Comparisons with three other formulations for isotropic and temperature-independent properties are also presented.

CONCLUSIONS

Compared to conventional methods focusing on numerical accuracy, convergence and stability, the proposed formulation focuses on computational performance for fast tissue thermal analysis. Compared to the classical Pennes model that handles only the static state of tissue, the proposed formulation can achieve fast thermal analysis on deformed states of tissue and can be applied in addition to tissue deformable models for non-linear heating analysis at even large deformation of soft tissue, leading to great translational potential in dynamic tissue temperature analysis and thermal dosimetry computation for computer-integrated medical education and personalised treatment.

The Influence of Dynamic Tissue Properties on HIFU Hyperthermia: A Numerical Simulation Study

Accurate temperature and thermal dose prediction are crucial to high-intensity focused ultrasound (HIFU) hyperthermia, which has been used successfully for the non-invasive treatment of solid tumors. For the conventional method of prediction, the tissue properties are usually set as constants. However, the temperature rise induced by HIFU irradiation in tissues will cause changes in the tissue properties that in turn affect the acoustic and temperature field. Herein, an acoustic–thermal coupling model is presented to predict the temperature and thermal damage zone in tissue in terms of the Westervelt equation and Pennes bioheat transfer equation, and the individual influence of each dynamic tissue property and the joint effect of all of the dynamic tissue properties are studied. The simulation results show that the dynamic acoustic absorption coefficient has the greatest influence on the temperature and thermal damage zone among all of the individual dynamic tissue properties. In addition, compared with the conventional method, the dynamic acoustic absorption coefficient leads to a higher focal temperature and a larger thermal damage zone; on the contrary, the dynamic blood perfusion leads to a lower focal temperature and a smaller thermal damage zone. Moreover, the conventional method underestimates the focal temperature and the thermal damage zone, compared with the simulation that was performed using all of the dynamic tissue properties. The results of this study will be helpful to guide the doctors to develop more accurate clinical protocols for HIFU treatment planning.

Effects of an implant on temperature distribution in tissue during ultrasound diathermy

The effects of an implant on temperature distribution in a tissue-mimicking hydrogel phantom during the application of therapeutic ultrasound were investigated. In vitro experiments were conducted to compare the influences of plastic and metal implants on ultrasound diathermy and to calibrate parameters in finite element simulation models. The temperature histories and characteristics of the opaque (denatured) areas in the hydrogel phantoms predicted by the numerical simulations show good correlation with those observed in the in vitro experiments. This study provides an insight into the temperature profile in the vicinity of an implant by therapeutic ultrasound heating typically used for physiotherapy. A parametric study was conducted through numerical simulations to investigate the effects of several factors, such as implant material type, ultrasound operation frequency, implant thickness and tissue thickness on the temperature distribution in the hydrogel phantom. The results indicate that the implant material type and implant thickness are the main parameters influencing the temperature distribution. In addition, once the implant material and ultrasound operation frequency are chosen, an optimal implant thickness can be obtained so as to avoid overheating injuries in tissue.

Development and validation of a MRgHIFU non-invasive tissue acoustic property estimation technique

MR-guided high-intensity focussed ultrasound (MRgHIFU) non-invasive ablative surgeries have advanced into clinical trials for treating many pathologies and cancers. A remaining challenge of these surgeries is accurately planning and monitoring tissue heating in the face of patient-specific and dynamic acoustic properties of tissues. Currently, non-invasive measurements of acoustic properties have not been implemented in MRgHIFU treatment planning and monitoring procedures. This methods-driven study presents a technique using MR temperature imaging (MRTI) during low-temperature HIFU sonications to non-invasively estimate sample-specific acoustic absorption and speed of sound values in tissue-mimicking phantoms. Using measured thermal properties, specific absorption rate (SAR) patterns are calculated from the MRTI data and compared to simulated SAR patterns iteratively generated via the Hybrid Angular Spectrum (HAS) method. Once the error between the simulated and measured patterns is minimised, the estimated acoustic property values are compared to the true phantom values obtained via an independent technique. The estimated values are then used to simulate temperature profiles in the phantoms, and compared to experimental temperature profiles. This study demonstrates that trends in acoustic absorption and speed of sound can be non-invasively estimated with average errors of 21% and 1%, respectively. Additionally, temperature predictions using the estimated properties on average match within 1.2 °C of the experimental peak temperature rises in the phantoms. The positive results achieved in tissue-mimicking phantoms presented in this study indicate that this technique may be extended to in vivo applications, improving HIFU sonication temperature rise predictions and treatment assessment.

Numerical and ex vivo studies of a bioprobe developed for laser-induced thermotherapy (LITT) in contact with liver tissue

This work is based on the production of a bioprobe that is compatible with magnetic resonance imaging (MRI) for laser-induced thermotherapy (LITT) in liver cancer laser therapy. This probe is made of an alumina tube (3-mm diameter) in which an optical fibre is centred and fixed. A shooting window (20mm) is created using a mechanical rectifier. The device is then consolidated by the injection of a transparent and heat-resistant resin. Through numerical modelling, the thermal power damping of the laser source is evaluated as well as the propagation of the heat in the ex vivo liver tissue according to different heating scenarios. These analyses allow for an estimation of the irradiated volume. Ex vivo tests were performed on bovine liver to confirm the adequacy of the bioprobe for LITT and of the irradiated volumes predicted by the numerical model. There was a difference of 8% between the simulations and ex vivo experiments. The pulsed mode heating scenario was the most effective under the experimental conditions.