Recently developed numerical model, based on the dislocation-driven nucleation of gas bubbles, is used to analyse experimental results on deuterium retention in tungsten under ITER relevant plasma exposure conditions. Focus is put on understanding the relation between exposure temperature and flux on primary features of thermal desorption spectra: peak positions and intensities of the desorption flux. The model allows one to relate the peak positions with the size of plasma induced deuterium bubbles and envisage exposure conditions (temperature and flux) for their formation. Based on the performed analysis, dedicated experimental conditions to validate the model are proposed.