Influence of mode conversions in the skull on transcranial focused ultrasound and temperature fields utilizing the wave field separation method: A numerical study*

Transcranial focused ultrasound is a booming noninvasive therapy for brain stimuli. The Kelvin–Voigt equations are employed to calculate the sound field created by focusing a 256-element planar phased array through a monkey skull with the time-reversal method. Mode conversions between compressional and shear waves exist in the skull. Therefore, the wave field separation method is introduced to calculate the contributions of the two waves to the acoustic intensity and the heat source, respectively. The Pennes equation is used to depict the temperature field induced by ultrasound. Five computational models with the same incident angle of 0° and different distances from the focus for the skull and three computational models at different incident angles and the same distance from the focus for the skull are studied. Numerical results indicate that for all computational models, the acoustic intensity at the focus with mode conversions is 12.05% less than that without mode conversions on average. For the temperature rise, this percentage is 12.02%. Besides, an underestimation of both the acoustic intensity and the temperature rise in the skull tends to occur if mode conversions are ignored. However, if the incident angle exceeds 30°, the rules of the over- and under-estimation may be reversed. Moreover, shear waves contribute 20.54% of the acoustic intensity and 20.74% of the temperature rise in the skull on average for all computational models. The percentage of the temperature rise in the skull from shear waves declines with the increase of the duration of the ultrasound.