Analysis and Simulation of Dynamic Heat Transfer and Thermal Distribution in Burns with Multilayer Models Using Finite Volumes

Burns represent a significant medical challenge, and the development of theoretical models has the potential to contribute to the advancement of new diagnostic tools. This study aimed to perform numerical simulations of the Pennes bioheat transfer equation, incorporating heat generation terms due to the body’s immunological response to thermal injury, as well as changes in skin thermal parameters and blood perfusion for each burn type. We propose the incorporation of specific parameters and boundary conditions related to multilayer perfusion into the Pennes bioheat model. Using the proposed layered skin model, we evaluate temperature differences to establish correlations for determining burn depth. In this investigation, 1D and 3D algorithms based on the finite volume method were applied to capture transient and spatial thermal variations, with the resulting temperature distributions demonstrating the ability of the proposed models to describe the expected thermal variations in healthy and burned tissue. This work demonstrates the potential of the finite volume method to approximate the solution of the Pennes biothermal equation. Overall, this study provides a computational framework for analyzing heat transfer in burn injuries and highlights the relevance of mathematical simulations as a tool for future research on infrared thermography in medicine.