A fundamental solution-based finite element model for analyzing multi-layer skin burn injury

To understand the physiology of tissue burns for successful clinical treatment, it is important to investigate the thermal behavior of human skin tissue subjected to heat injury. In this paper, a fundamental solution-based hybrid finite element formulation is proposed for numerically simulating steady-state temperature distribution inside a multilayer human skin tissue during burning. In the present approach, since only element boundary integrals are involved, the computational dimension is reduced by one as the fundamental solutions used analytically satisfies the bioheat governing equation. Further, in multi-layer skin modeling, the burn is applied via a heating disk at constant temperature on a part of the epidermal surface of the skin tissue. Numerical results from the proposed approach are firstly verified by comparing them with exact solutions of a simple single-layered model or the results from conventional finite element method. Thereafter, a sensitivity analysis is carried out to reveal the effect of biological and environmental parameters on temperature distribution inside the skin tissue subjected to heat injury.