Mesh-independent prediction of skin burns injury

This paper presents a robust finite element model (FEM) with multiple-layers of varying properties for investigation of burn effects on human skin during a burning process resulting from exposure of skin surface to a contact heat source and a hot moving fluid. Henriques' theory of skin burns is used in conjunction with two-dimensional Pennes bioheat transfer equation for determining the spatial and temporal extent of burn injury. The model developed is a two-dimensional axisymmetric model in cylindrical coordinates. The various tissue layers account for changing thermal properties with respect to skin anatomy. A finite element scheme that uses the backward Euler method is used to solve the problem. The injury processes of skin subsequent to the removal of the heat source (post-burn) will also be inspected. The mesh employed in this model consists of a high density of nodes and elements in which a thorough mesh convergence study was done. A comparison of the transient temperature field computed by this model against Diller's results using the FE technique with a comparatively coarse mesh of 125 elements and experimental data by Orgill et al. has been done in the present study. It concluded that improved accurate solutions have been performed using the robust model developed due to the achievement of a mesh-independent solution.