A meshless method for modelling two-phase flows with phase transition is described. The method is based on consideration of three systems: viscous-vortex blobs, thermal-blobs and droplets; and can be applied for numerical simulation of 2D nonisothermal flows of 'gas-evaporating droplets' in the framework of the one-way coupled two-fluid approach. The carrier phase is viscous incompressible gas. The dispersed phase is presented by a cloud of identical spherical droplets, and, due to evaporation, the radius and mass of droplets are time dependent. The carrier phase parameters are calculated using the viscous-vortex and thermal-blob method; the dispersed phase parameters are calculated using the Lagrangian approach. Two applications have been considered: (i) a standard benchmark -Lamb vortex; (ii) a cold spray injected into a hot quiescent gas. In the latter problem three cases corresponding to three droplet sizes were investigated. The smallest droplets (of the three cases considered) are more readily entrained by the carrier phase and form ring-like structures; the flow shows better mixing. Larger droplets evaporate less intensively. The medium sized droplets collect into two narrow bands stretched along the jet axis. The largest droplets form a two-phase jet, which remains close to the jet axis.