In the paper the problem concerning the numerical modeling of thermal processes in domain of biological tissue secured by a layer of protective clothing being in the thermal contact with the environment is discussed. The cross-section of the forearm (2D problem) is treated as the non-homogeneous domain in which the sub-domains of skin tissue, fat, muscle, bone and blood vessels are distinguished. Between skin tissue and protective clothing the air gap is taken into account. This sub-domain is treated as a solid body which thermal conductivity is defined in the special way. The process of external heating is determined by the Neumann boundary condition and the sensitivity analysis with respect to the perturbations of the boundary heat flux is discussed. Both the basic boundary-initial problem and the sensitivity one are solved by means of the control volume method using the Voronoi polygons.