The aim of this study is to evaluate the effect of pulsatile blood flow in thermally significant blood vessels on the thermal lesion region during thermal therapy of tumor. A sinusoidally pulsatile velocity profile for blood flow was employed to simulate the cyclic effect of the heart beat on the blood flow. The evolution of temperature field was governed by the energy transport equation for blood flow together with Pennes' bioheat equation for perfused tissue encircling the blood vessel. The governing equations were numerically solved by a novel multi-block Chebyshev pseudospectral method and the accumulated thermal dose in tissue was computed. Numerical results show that pulsatile velocity profile, with various combinations of pulsatile amplitude and frequency, has little difference in effect on the thermal lesion region of tissue compared with uniform or parabolic velocity profile. However, some minor differences on the thermal lesion region of blood vessel is observed for middle-sized blood vessel. This consequence suggests that, in this kind of problem, we may as well do the simulation simply by a steady uniform velocity profile for blood flow.
