Numerical analysis of local non-equilibrium heat transfer in layered spherical tissue during magnetic hyperthermia

Influence of different heat transfer models on therapeutic temperature prediction and heat-induced damage during magnetic hyperthermia

Magnetic hyperthermia regulates the therapeutic temperature within a specific range to damage malignant cells after exposing the magnetic nanoparticles inside tumor tissue to an alternating magnetic field. The therapeutic temperature of living tissues can be generally predicted using Pennes' bio-heat equation after ignoring both the inhomogeneity of biological structure and the microstructural responses. Although various of the bio-heat transfer models proposed in literature fix these shortages, there is still a lack of a comprehensive report on investigating the discrepancy for different models when applied in the magnetic hyperthermia context. This study compares four different bio-heat equations in terms of the therapeutic temperature distribution and the heat-induced damage situation for a proposed geometric model, which is established based on computed tomography images of a tumor bearing mouse. The therapeutic temperature is also used as an index to evaluate the effect of two key relaxation times for the phase lag behavior on bio-heat transfer. Moreover, this work evaluates the effects of two blood perfusion rates on both the treatment temperature and the cumulative equivalent heating minutes at 43 °C. Numerical analysis results reveal that relaxation times for phase-lag behavior as well as the porosity for living tissues directly affect the therapeutic temperature variation and ultimately the thermal damage for the malignant tissue during magnetic hyperthermia. The dual-phase-lag equation can be converted into Pennes' equation and simple-phase-lag equation when relaxation times meet specific conditions during the process of heat transfer. In addition, different blood perfusion rates can result in an amplitude discrepancy for treatment temperature, but this parameter does not change the characteristics of thermal propagation during therapy.

Thermal therapy induced fluid pressure and stress reductions in a solid tumor

This paper presents an investigation on the interstitial fluid pressure and stress reductions in a vascularized solid tumor using a thermal therapy approach. The solid tumor is modeled as a fluid infiltrated poroelastic medium with a pressure source subjected to spatial heating. The distributions of temperature, interstitial fluid pressure, strains and stresses in a spherical tumor are obtained using a thermoporoelasticity theory in which the extracellular solid matrix and the interstitial fluid have different coefficient of thermal expansion (CTE). The numerical results for a solid tumor subjected to uniform spatial heating indicate that the CTE of the solid matrix of the tumor plays a crucial role in the reductions in the fluid pressure and effective stresses caused by the thermal therapy. The pore pressure and effective stresses are reduced when the CTE of the solid matrix is higher than that of the interstitial fluid. The reductions in fluid pressure and stresses may become significant depending on the difference between the CTEs of the solid matrix and interstitial fluid. The reductions reach the maximum at the tumor center and decrease with increasing radial distance from the tumor center. Finally, the thermally induced fluid flow is directed from the surface towards the center thereby potentially improving the microcirculation in the solid tumor.