Numerical modeling of biological tissue heating : admissible thermal dose

The cylindrical domain of skin tissue subjected to an external heat flux is considered (shape of domain is determined by form of function describing Neumann boundary condition on external surface of the system). The first version of numerical simulations concerns the heterogeneous multi-layered skin tissue domain (epidermis, dermis, subcutaneous region). The thermophysical parameters of successive layers are assumed to be different, but constant. The second version of computations concerns the homogeneous domain, but the mean values of the thermophysical parameters are temperature-dependent (non-linear task). Knowledge of the spatial, time-dependent temperature field allows one to determine the so-called thermal dose and also the degree of tissue destruction. The algorithm presented can be useful in medical practice, among others, at the stage of the hyperthermia therapy scheme.