Numerical simulation for heat transfer in tissues during thermal therapy

Abstract

In this paper, a mathematical model describing the process of heat transfer in biological tissues for different coordinate system during thermal therapy by electromagnetic radiation is studied. The boundary value problem governing this process has been solved using Galerkin's method taking B-polynomial as basis function. The system of ordinary differential equation in unknown time variable, thus obtained, is solved by homotopy perturbation method. The effect of thermal conductivity, antenna power constant, surface temperature, and blood perfusion rate on temperature for different coordinates are discussed. It has been observed that the process is faster in spherical symmetric coordinates in comparison to axisymmetric coordinate and faster in axisymmetric in comparison to Cartesian coordinate. © 2010 Elsevier Ltd.