We apply Boltzmann equations for modelling the radiation damage in samples irradiated by photons from free electron laser (FEL). We test this method in a study case of a spherically symmetric xenon cluster irradiated with VUV FEL photons. The results obtained demonstrate the potential of the Boltzmann method for describing the complex and non-equilibrium dynamics of samples exposed to FEL radiation.