Parameter estimation in high-intensity focused ultrasound therapy

Hyperthermia using High-Intensity Focused Ultrasound (HIFU) is an acoustic therapy for cancer treatment. This technique consists of an increase in the temperature field of the tumor to achieve coagulative necrosis and immediate cell death. Therefore, for having a successful treatment, the physical problem requires to know several properties due to the high variability from individual to individual, or even for the same individual under different physiological conditions. This article presents a numerical simulation of hyperthermia therapy for cancer treatment using HIFU, as well as the estimation of parameters that influence the physical problem. Two mathematical models were considered to solve the forward problem. The acoustic model based on acoustic pressure performs a frequency-domain study, and the bioheat transfer model a time-dependent study. These models were solved using Comsol Multiphysics (R) software in a 2D-axisymmetric rectangular domain to determine the temperature field. Parameter estimation was coded in Matlab Mathworks (R) environment using a Bayesian approach. The Markov Chain Monte Carlo method by the Metropolis-Hastings algorithm was implemented, and the simulated temperature measurements were considered. Results suggest that specific HIFU therapy can be performed for each patient by estimating appropriate parameters for cancer treatment and provides the possibility to define procedures before and during the treatment.