Magnetic fluid hyperthermia simulations in evaluation of SAR calculation methods

• Published in: Physica medica 2020-03, Vol.71, p.39-52
• Main Authors: Papadopoulos, Costas, Efthimiadou, Eleni K., Pissas, Michael, Fuentes, David, Boukos, Nikolaos, Psycharis, Vassilis, Kordas, George, Loukopoulos, Vassilios C., Kagadis, George C.
• Format: Article
• Language: English
• Subjects: Algorithms; Box-Lucas; Calorimetry; Computer Simulation; COMSOL; Corrected slope method; Hot Temperature; Humans; Hyperthermia, Induced - methods; ILP; Initial slope method; Linear Models; Magnetic fluid hyperthermia; Magnetic nanoparticles; Magnetics; Magnetite Nanoparticles; Normal Distribution; Reproducibility of Results; SAR; Simulations; Specific absorption rate
• ISSN: 1120-1797; 1724-191X
• DOI: 10.1016/j.ejmp.2020.02.011

•Different calculation methods were used to estimate SAR values of magnetite MNPs.•A detailed MFH numerical model was developed in COMSOL Multiphysics.•Reference SAR values were generated in silico to evaluate calculation methods.•A novel method was presented to accurately compute the temperature-dependent SAR. The purpose of this study is to employ magnetic fluid hyperthermia simulations in the precise computation of Specific Absorption Rate functions -SAR(T)-, and in the evaluation of the predictive capacity of different SAR calculation methods. Magnetic fluid hyperthermia experiments were carried out using magnetite-based nanofluids. The respective SAR values were estimated through four different calculation methods including the initial slope method, the Box-Lucas method, the corrected slope method and the incremental analysis method (INCAM). A novel numerical model combining the heat transfer equations and the Navier-Stokes equations was developed to reproduce the experimental heating process. To address variations in heating efficiency with temperature, the expression of the power dissipation as a Gaussian function of temperature was introduced and the Levenberg-Marquardt optimization algorithm was employed to compute the function parameters and determine the function’s effective branch within each measurement’s temperature range. The power dissipation function was then reduced to the respective SAR function. The INCAM exhibited the lowest relative errors ranging between 0.62 and 15.03% with respect to the simulations. SAR(T) functions exhibited significant variations, up to 45%, within the MFH-relevant temperature range. The examined calculation methods are not suitable to accurately quantify the heating efficiency of a magnetic fluid. Numerical models can be exploited to effectively compute SAR(T) and contribute to the development of robust hyperthermia treatment planning applications.