“…This may, in principle, be obtained directly from MRI thermometry data of the treatments, however, for simulation purposes in Comsol Multiphysics, it is easier to model the temperature rise due to a theoretical acoustic beam than work with noisy data. This may be done using the bioheat equation, that is, we can solve (see O’Neill et al 2012, Wissler 1998):
with boundaries in the +R direction fixed at 37 °C, in the ±z direction set to 21 °C, and an insulation condition at the R=0 axis. For simplicity, we modeled the ultrasound as a Gaussian beam (Wu and Du 1990, O’Neill et al 2012) in cylindrical coordinates, that is:
…”