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

• Published in: Mathematical biosciences Vol. 281; pp. 82 - 91
• Main Authors: Kumar, P., Kumar, Dinesh, Rai, K.N.
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 01.11.2016; Elsevier Science Ltd
• Subjects: Blood; Boundary conditions; Cartesian coordinates; Computer simulation; Dimensionless analysis; DPL bioheat transfer model; Energy Transfer; Finite difference method; Finite element wavelet Galerkin method; Generalized coordinate system; Heat conductivity; Heat generation; Heat transfer; Humans; Hyperthermia; Hyperthermia, Induced; Mathematical models; Models, Biological; Non-Fourier boundary conditions; Perfusion; Phase lag; Simulation; Spherical coordinates; Temperature dependence; Temperature effects; Therapy; Thermal conductivity; Thermal damage; Thermal therapy; Thermodynamic properties; Wavelet analysis; Non-Fourier boundary conditions; Generalized coordinate system; Finite element wavelet Galerkin method; DPL bioheat transfer model; Thermal therapy; Thermal damage
• ISSN: 0025-5564; 1879-3134
• DOI: 10.1016/j.mbs.2016.08.013

•Surrounding tissues are less affected in 2nd and 3rd kind of B.C. when used with spherical coordinate.•FEWGM localizes small scale variation of solution and fast switching of functional basis.•Wave or diffusion-like behavior in tissues dominates depends whether τq exceeds or precedes τT.•Effect of blood perfusion rate, thermal conductivity and metabolic heat generation rate on tissue temperature has been studied. 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.