In the paper we consider the problem with moving bounds that models the kinetics of hydrogen desorption from hydrides of metals. Phase change, size reduction, desorption processes on the surface, and heat absorption on the phase bound are taken into consideration. The differences scheme for the corresponding boundary-value problem with moving bounds is presented for the experimental method of thermal desorption spectrometry (TDS). The results of numerical modeling are presented.