Comparison of Five Numerical Simulation Algorithms in Temperature Prediction for Hollow Microspheres in Magnetic Induction Hyperthermia

Abstract

Magnetic induction hyperthermia (MIH) uses high temperatures beyond 42°C to kill cancer cells. The temperature simulation is an essential part of the hyperthermia planning system. Pre-temperature simulations can help doctors decide the amount and the position of the magnetic agents. Most papers focus on only one numerical heat transfer algorithm to simulate temperature changes in MIH. They usually use software like FLUENT to calculate the temperature variation. Nevertheless, in practice, the core calculation code needs to be written by ourselves, combined with tumor model from CT images. We need to decide which heat transfer algorithm has the practical results. In this paper, we focus on the different forms of FDM. We use five heat transfer algorithms including the explicit format of FDM, Peaceman-Rachford ADI, Brian ADI, Douglas ADI and FVM to simulate. The results of simulations show that the approximations in FDM and FVM lead to relatively large offset compared with other simulations. In addition, the difference between the results of simulations and experiments is caused by the heat power used in simulations, the inaccuracy of physical parameters and measurement error.