Numerical solution of non-linear dual-phase-lag bioheat transfer equation within skin tissues.

This paper deals with numerical modeling and simulation of heat transfer in skin tissues using non-linear dual-phase-lag (DPL) bioheat transfer model under periodic heat flux boundary condition. The blood perfusion is assumed temperature-dependent which results in non-linear DPL bioheat transfer model in order to predict more accurate results. A numerical method of line which is based on finite difference and Runge-Kutta (4,5) schemes, is used to solve the present non-linear problem. Under specific case, the exact solution has been obtained and compared with the present numerical scheme, and we found that those are in good agreement. A comparison based on model selection criterion (AIC) has been made among non-linear DPL models when the variation of blood perfusion rate with temperature is of constant, linear and exponential type with the experimental data and it has been found that non-linear DPL model with exponential variation of blood perfusion rate is closest to the experimental data. In addition, it is found that due to absence of phase-lag phenomena in Pennes bioheat transfer model, it achieves steady state more quickly and always predict higher temperature than thermal and DPL non-linear models. The effect of coefficient of blood perfusion rate, dimensionless heating frequency and Kirchoff number on dimensionless temperature distribution has also been analyzed. The whole analysis is presented in dimensionless form.

[1]  R. Huggins.  Some practical aspects of a conditional likelihood approach to capture experiments , 1991 .

[2]  H G Klinger,et al.  Relations between heat transfer in perfused biological tissue and the local symmetry components of the vascular system , 1983, Journal of mathematical biology.

[3]  Steven L. Jacques,et al.  Nonlinear finite-element analysis of the role of dynamic changes in blood perfusion and optical properties in laser coagulation of tissue , 1996 .

[4]  Q. Qin,et al.  METHOD OF FUNDAMENTAL SOLUTIONS FOR NONLINEAR SKIN BIOHEAT MODEL , 2014 .

[5]  Keith A. Seffen,et al.  Mathematical Modeling of Skin Bioheat Transfer , 2009 .

[6]  David R. Anderson,et al.  Data-Based Selection of an Appropriate Biological Model: The Key to Modern Data Analysis , 1992 .

[7]  P. Kumar,et al.  Numerical simulation of dual-phase-lag bioheat transfer model during thermal therapy. , 2016, Mathematical biosciences.

[8]  K. N. Rai,et al.  Analysis of classical Fourier, SPL and DPL heat transfer model in biological tissues in presence of metabolic and external heat source , 2016 .

[9]  Da Yu Tzou,et al.  Macro- to Microscale Heat Transfer: The Lagging Behavior , 2014 .

[10]  J. A. J. Stolwijk,et al.  Temperature regulation in man — A theoretical study , 1966, Pflüger's Archiv für die gesamte Physiologie des Menschen und der Tiere.

[11]  T. Lu,et al.  Non-Fourier analysis of skin biothermomechanics , 2008 .

[12]  E. Y. K. Ng,et al.  Prediction and parametric analysis of thermal profiles within heated human skin using the boundary element method , 2010, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[13]  R. M. Cotta,et al.  Hybrid integral transforms analysis of the bioheat equation with variable properties , 2010 .

[14]  Subhash C. Mishra,et al.  Conventional and newly developed bioheat transport models in vascularized tissues: A review , 2013 .

[15]  Jitendra Singh,et al.  Solution of the heat transfer problem in tissues during hyperthermia by finite difference-decomposition method , 2013, Appl. Math. Comput..

[16]  Shanshan Yang,et al.  Analysis of thermal conductivity in living biological tissue with vascular network and convection , 2014 .

[17]  R. Vaillancourt,et al.  Temperature Distributions for Regional Hypothermia Based on Nonlinear Bioheat Equation of Pennes Type: Dermis and Subcutaneous Tissues , 2012 .

[18]  Ruy Freitas Reis,et al.  3D numerical simulations on GPUs of hyperthermia with nanoparticles by a nonlinear bioheat model , 2016, J. Comput. Appl. Math..

[19]  P. Kumar,et al.  A numerical study on dual-phase-lag model of bio-heat transfer during hyperthermia treatment. , 2015, Journal of thermal biology.

[20]  J. Ghazanfarian,et al.  Implementation of Smoothed-Particle Hydrodynamics for non-linear Pennes' bioheat transfer equation , 2015, Appl. Math. Comput..

[21]  Xiaohao Wei,et al.  Heat Conduction: Mathematical Models and Analytical Solutions , 2008 .

[22]  R. Royall.  Statistical Evidence: A Likelihood Paradigm , 1997 .

[23]  Lawrence F. Shampine,et al.  An efficient Runge-Kutta (4,5) pair , 1996 .

[24]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations , 1989 .

[25]  J. Liu,et al.  Parametric studies on the phase shift method to measure the blood perfusion of biological bodies. , 2000, Medical engineering & physics.

[26]  R. C. Mittal,et al.  Numerical simulation of two-dimensional sine-Gordon solitons by differential quadrature method , 2012, Comput. Phys. Commun..

[27]  D. Tzou.  A Unified Field Approach for Heat Conduction From Macro- to Micro-Scales , 1995 .

[28]  Clifford M. Hurvich,et al.  Regression and time series model selection in small samples , 1989 .

[29]  David R. Anderson,et al.  Modeling Survival and Testing Biological Hypotheses Using Marked Animals: A Unified Approach with Case Studies , 1992 .