The availability of reliable modeling tools and input data required for the prediction of surface removal rate from the lithium fluoride targets irradiated by the intense photon beams is essential for many practical aspects. This study is motivated by the practical implementation of soft X-ray (SXR) or extreme ultraviolet (XUV) lasers for the pulsed ablation and thin film deposition. Specifically, it is focused on quantitative description of XUV laser-induced desorption/ablation from lithium fluoride, which is a reference large band-gap dielectric material with ionic crystalline structure. Computational framework was proposed and employed here for the reconstruction of plume expansion dynamics induced by the irradiation of lithium fluoride targets. The morphology of experimentally observed desorption/ablation craters were reproduced using idealized representation (two-zone approximation) of the laser fluence profile. The calculation of desorption/ablation rate was performed using one-dimensional thermomechanic model (XUV-ABLATOR code) taking into account laser heating and surface evaporation of the lithium fluoride target occurring on a nanosecond timescale. This step was followed by the application of two-dimensional hydrodynamic solver for description of laser-produced plasma plume expansion dynamics. The calculated plume lengths determined by numerical simulations were compared with a simple adiabatic expansion (blast-wave) model.