Analytical solution of the parabolic and hyperbolic heat transfer equations with constant and transient heat flux conditions on skin tissue

In this article, the parabolic (Pennes bioheat equation) and hyperbolic (thermal wave) bioheat transfer models for constant, periodic and pulse train heat flux boundary conditions are solved analytically by applying the Laplace transform method for skin as a semi-infinite and finite domain. The bioheat transfer analysis with transient heat flux on skin tissue has only been studied by Pennes equation for a semi-infinite domain. For modeling heat transfer in short duration of an initial transient, or when the propagation speed of the thermal wave is finite, there are major differences between the results of parabolic and hyperbolic heat transfer equations. The non-Fourier bioheat transfer equation describes the thermal behavior in the biological tissues better than Fourier equation. The outcome of transient heat flux condition shows that by penetrating into the depths beneath the skin subjected to heat, the amplitude of temperature response decreases significantly. The blood perfusion rate can be predicted using the phase shift between the surface temperature and transient surface heat flux. The thermal damage of the skin is studied by applying both the parabolic and hyperbolic bioheat transfer equations.