Numerical simulation of burn injuries with temperature-dependent thermal conductivity and metabolism under different surface heat sources.

In the present paper, the phenomena of heat transport inside human forearm tissue are studied through a one-dimensional nonlinear bioheat transfer model under the influence of various boundary and interface conditions. In this study, we considered temperature-dependent thermal conductivity and metabolic heat to predict temperature distribution inside the forearm tissue. We have studied the temperature distribution inside inner tissue and bone because it has been found that burn injuries are mostly affected by layer thickness. The temperature distribution inside human forearm tissue is analyzed using the finite difference and bvp4c numerical techniques. To examine the accuracy of present numerical code, we compare the obtained numerical result with the exact analytical result in a specific case and find an excellent agreement with the exact results. We also validated our present numerical code with a hybrid scheme based on Runge-Kutta (4,5) and finite difference technique and found it in good compliance. From the obtained results, we observed that the homogeneous heat flux has a greater impact on the temperature at the outer surface of the skin, but the sinusoidal heat flux has a greater impact on the temperature of the subcutaneous layer and inner tissue. It is found that there is no burn injury in the first type of heat source (T=44°C), but it may occur in the second and third types of heat sources. It has been observed that by raising the blood perfusion rate and reducing the values of reference metabolic heat, coefficient of thermal conductivity, and heat fluxes, we can manage and reduce burn injuries and achieve hyperthermia temperature.