The Key Role of Thermal Relaxation Time on the Improved Generalized Bioheat Equation: Analytical Versus Simulated...

The Pennes bioheat equation is the most widely used model for describing heat transfer in living tissue during thermal exposure. It is derived from the classical Fourier law of heat conduction and assumes energy exchange between blood vessels and surrounding tissues. The literature presents various numerical methods for solving the bioheat equation, with exact solutions developed for different boundary conditions and geometries. However, analytical models based on this framework are rarely reported. This study aims to develop an analytical three-dimensional model using MATHEMATICA software, with subsequent mathematical validation performed through COMSOL simulations, to characterize heat transfer in biological tissues induced by laser irradiation under various therapeutic conditions. The objective is to refine the conventional bioheat equation by introducing three key improvements: (a) incorporating a non-Fourier framework for the Pennes equation, thereby accounting for the relaxation time in thermal response; (b) integrating Dirac functions and the telegraph equation into the bioheat model to simulate localized point heating of diseased tissue; and (c) deriving a closed-form analytical solution for the Pennes equation in both its classical (Fourier-based) and improved (non-Fourier-based) formulations. This paper investigates the nuanced relationship between the relaxation time parameter in the telegraph equation and the thermal relaxation time employed in the bioheat transfer equation. Considering all these aspects, the optimal thermal relaxation time determined for these simulations was 1.16 s, while the investigated thermal exposure time ranged from 0.01 s to 120 s. This study introduces a generalized version of the model, providing a more realistic representation of heat exchange between biological tissue and blood flow by accounting for non-uniform temperature distribution. It is important to note that a reasonable agreement was observed between the two modeling approaches: analytical (MATHEMATICA) and numerical (COMSOL) simulations. As a result, this research paves the way for advancements in laser-based medical treatments and thermal therapies, ultimately contributing to more optimized therapeutic outcomes.