Modeling of high-intensity focused ultrasound-induced lesions in the presence of cavitation bubbles

The classical "Bio Heat Transfer Equation (BHTE)" model is adapted to take into account the effects of oscillating microbubbles that occur naturally in the tissue during high-intensity focused ultrasound (HIFU) treatment. First, the Gilmore-Akulichev model is used to quantify the acoustic pressure scattered by microbubbles submitted to HIFU. Because this scattered pressure is not monochromatic, the concept of harmonic attenuation is introduced and a global attenuation coefficient is estimated for bubble-filled tissues. The first results show that this global attenuation coefficient varies significantly with respect to several parameters such as the frequency and the density of microbubbles in the medium, but also with respect to the incident acoustic pressure which thus becomes a transcendental function. Under these conditions, a layer-by-layer modeling, in the direction of propagation, is proposed to calculate the ultrasonic beam. Finally, the BHTE is solved and the HIFU-induced lesions are estimated by the calculation of the thermal dose. Using this model, it can be observed first that, when the firing power increases, the lesion develops clearly in the direction of the transducer, with a shape agreeing with in vivo experimentation. Next, it is observed that the lesion can be significantly modified in size and position, if an interface (skin or inner wall) is simulated as a zone with multiple cavitation nuclei. With a firing power increase, it is also shown how a secondary lesion can appear at the interface and how, beyond a certain threshold, this lesion develops at the main lesion expense. Finally, a better in-depth homogeneity of lesions is observed when the acoustic frequency of HIFU is increased.