Dynamic modeling of photothermal interactions for laser-induced interstitial thermotherapy: parameter sensitivity analysis

A two-dimensional model was developed to model the effects of dynamic changes in the physical properties on tissue temperature and damage to simulate laser-induced interstitial thermotherapy (LITT) treatment procedures with temperature monitoring. A modified Monte Carlo method was used to simulate photon transport in the tissue in the non-uniform optical property field with the finite volume method used to solve the Pennes bioheat equation to calculate the temperature distribution and the Arrhenius equation used to predict the thermal damage extent. The laser light transport and the heat transfer as well as the damage accumulation were calculated iteratively at each time step. The influences of different laser sources, different applicator sizes, and different irradiation modes on the final damage volume were analyzed to optimize the LITT treatment. The numerical results showed that damage volume was the smallest for the 1,064-nm laser, with much larger, similar damage volumes for the 980- and 850-nm lasers at normal blood perfusion rates. The damage volume was the largest for the 1,064-nm laser with significantly smaller, similar damage volumes for the 980- and 850-nm lasers with temporally interrupted blood perfusion. The numerical results also showed that the variations in applicator sizes, laser powers, heating durations and temperature monitoring ranges significantly affected the shapes and sizes of the thermal damage zones. The shapes and sizes of the thermal damage zones can be optimized by selecting different applicator sizes, laser powers, heating duration times, temperature monitoring ranges, etc.