Numerical simulation of dual-phase-lag bioheat transfer model during thermal therapy

This paper theoretically investigates the thermal behavior in a living biological tissue under various coordinate systems and different non-Fourier boundary conditions with the dual-phase-lag bioheat transfer model during thermal therapy. The properties of Legendre wavelets together with the finite difference scheme are used to find an approximate analytical solution of the present problem. It has been observed that surrounding healthy tissues are less affected in second and third kind of boundary condition when applied along with spherical symmetric coordinate system. Also greater temperature rise and fast achievement of peak hyperthermia temperature is achieved when second and third kind of boundary conditions are used in combination with Cartesian coordinate system. It is observed that due to the presence of blood perfusion and temperature dependent metabolic heat generation term, the dual-phase-lag bioheat transfer model reduces to Pennes bioheat transfer model only when τ_q=τ_T=0s, not for arbitrary τ_q=τ_T. Further, in case of dual-phase-lag bioheat transfer model wave-like or diffusion-like behavior will dominate depends whether the ratio τ_q/τ_T > 1 or τ_q/τ_T < 1. Effect of temperature dependent metabolic heat generation rate, thermal conductivity and blood perfusion rate on dimensionless temperature are discussed in details. The whole analysis is presented in dimensionless form.