Modeling of bioheat equation for skin and a preliminary study on a noninvasive diagnostic method for skin burn wounds

Heat transfer in a unit three-dimensional skin tissue with an embedded vascular system of actual histology structure is computed in the present work. The tissue temperature and the blood temperatures in artery and vein vessels are solved with a multi-grid system. The mean temperature of the tissue over the cross-section of the unit skin area is evaluated. The resulting one-dimensional function is regarded as the temperature of healthy tissue (or injured skin but the blood perfusion is still normally working) for large area of skin in view of the symmetric and periodic structure of the paired artery–vein vessels in nature. A three-dimensional bioheat equation then is formulated by the superposition of the skin burn wound effect and the healthy skin temperature with and without thermal radiation exposure. When this bioheat equation is employed to simulate ADT process on burn wounds, the decaying factor of the skin surface temperature is found to be a sharply decreasing function of time in the self-cooling stage after a thermal radiation heating. Nevertheless, the boundary of non-healing (needing surgery) and healing regions in a large burn wound can be estimated by tracking the peak of the gradient of decaying factor within 30 s after the thermal radiation is turned off. Experimental studies on the full ADT procedure are needed to justify the assumptions in the present computation.