Impact of Treatment Trajectory on Temperature Field Uniformity in Biological Tissue Irradiated by Ultrasound Pulses with Shocks

High intensity focused ultrasound (HIFU) treatments typically involve the ablation of tissue volumes comprising multiple focal sites. One aspect of treatment planning involves the definition of a sequence of ultrasound pulses and corresponding focal sites as the sonication trajectory. Here, numerical simulations of the thermal effects of different trajectories are performed for HIFU exposures delivered to an ex vivo bovine liver sample by a clinical array (Sonalleve V2 3.0T system, Profound Medical Corp., Canada). Simulations consider boiling histotripsy regime with millisecond-long pulses that include shocks. Focusing of the ultrasound beam in tissue was modeled by the Westervelt equation, and the temperature field was modeled by the bioheat equation. To explore different treatment strategies, trajectories were considered with discrete foci located along two or four concentric circles with radii from 2 to 8 mm. Two approaches for traversing these focal sites were compared: In the first approach each discrete focus was sonicated by a sequence of 15 pulses before moving to the next site in the trajectory. In the second approach, each focus was sonicated once before moving to the next site, with sonications over the whole trajectory repeated 15 times. The influence of the trajectory’s size and the pulsing strategy on the temperature field was analyzed. It is shown that the structure of the temperature field is more uniform with a longer time interval between repeated irradiation of each focus, and the optimal time interval ranges from three to six pulse repetition periods.