Numerical analysis of thermal response on a non-linear model of multi-layer skin under heating and cooling processes

Based on the non-linear model of multi-layer skin, the present study is performed under heating and cooling processes with linear and non-linear generalized boundary conditions. The numerical outcome is obtained utilizing Runge Kutta (4,5) along with the finite difference scheme and the accuracy of this scheme is shown graphically by comparing it with an accurate analytical outcome in a special case. When the value of gamma increases, the skin temperature is gradually higher in the case of a linear boundary condition and gradually lower in the case of a non-linear boundary condition with respect to linear boundary condition at epidermis-dermis (ED) interface. During cooling, the heat effect has slightly quick vanished in a non-linear boundary condition than in a linear boundary condition. The effect of a non-linear boundary condition gradually decreases as the value of heat transfer coefficient rises, then it reflecting the nature of the first kind of linear boundary condition. The exponential blood perfusion rate has a slightly quick vanishes of heat effect in comparison to constant and linear blood perfusion rate. To investigate the conduct of the temperature distribution in multi-layer skin, all the impacts are depicted graphically.