Kinetic Boltzmann equations are used to model the ionization and expansion dynamics of xenon clusters irradiated with short intense VUV pulses. This unified model includes predominant interactions that contribute to the cluster dynamics induced by this radiation. The dependence of the evolution dynamics on cluster size, Natoms = 20 − 90000, and pulse fluence, F = 0.05 − 1.5 J/cm 2 , corresponding to intensities in the range, 10 12 − 10 14 W/cm 2 and irradiation times, ≤ 50 fs, is investigated. The predictions obtained with our model are found to be in good agreement with the experimental data. We find that during the exposure the cluster forms a shell structure consisting of a positively charged outer shell and a core of net charge equal to zero. The width of these shells depends on the cluster size. The charged outer shell is large within small clusters (Natoms = 20, 70), and its Coulomb explosion drives the expansion of these clusters. Within the large clusters (Natoms = 2500, 90000) the neutral core is large, and after the Coulomb explosion of the outer shell it expands hydrodynamically. Highly charged ions within the core recombine efficiently with electrons. As a result, we observe a large fraction of neutral atoms created within the core, its magnitude depending on the cluster size. Atomic clusters are excellent objects to test the dynamics within samples irradiated with radiation from short wavelength free-electron-lasers (FELs) [1,2,3]. Their physical properties put them on the border between the solid state and the gas phase. Cluster studies are important for planned experiments with FELs in solid state physics, materials science and for studies of the extreme states of matter [4]. Accurate predictions on the ionization, thermalization and expansion timescales within irradiated samples that can be obtained with cluster experiments are also needed for exploring the limits of experiments on single particle diffraction imaging [5,6,7,8,9,10].During the first cluster experiments performed at the free-electron-laser facility FLASH at DESY with VUV photons of energy, E = 12.7 eV, and power densities up to a few 10 13 W/cm 2 [11] highly charged Xe ions (up to +8) of high kinetic energies were detected. This unexpectedly strong energy absorption could not be explained using the standard approaches [12,13,14]. More specifically, the energy absorbed was almost an order of magnitude larger than that predicted with classical absorption models, and the ion charge states created were much higher than those observed during the irradiation of * Corresponding author. E-mail: ziaja@mail.desy.de isolated atoms in similar conditions. This indicated that at such radiation wavelengths some processes specific to many-body systems were responsible for an enhanced energy absorption.The physics underlying the dynamics within the irradiated clusters is complex. Several interesting theoretical models have been proposed in order to describe the evolution of clusters exposed to intense VUV pulses [15,16,17,18,19,20,21]. The most exp...