Calculation by the Method of Finite Differences of the Temperature Distribution in Layered Tissues

A numerical approach to obtain the temperature distribution in layered tissues with thermal source and cooling inside the tissues is presented in this paper. This approach can be applied to problems with a distributed internal thermal source produced by external radiation such as microwave, ultrasound, and shortwave, or by internal increase of metabolic rate in live subjects. The exact mechanism of cooling is not discussed here, but it is assumed to be due to blood flow in vivo. The calculation of the temperature distribution in layered tissues is based on a linear model of the tissues that consists of several layers of different kinds of tissues separated by parallel boundaries. Some simple mathematical forms of the cooling functions are also assumed in this paper. The results obtained agree well with six available experiments.