A velocity-scalar joint filtered density function (VSJFDF) transport equation was derived for the fluid flow field experienced by the particles in the evaporating gas-droplet two-phase flow. In this model, the large eddy simulation (LES) is used to simulate the fluid flow and the Lagrangian method is used to simulate the motion of the particle. In order to consider the sub-grid scale (SGS) effects of the gas velocity and scalar fluctuation transport and evaporation, we utilized the transport equation of the velocity-scalar joint filtered density function (VSJFDF) of the gas flow experienced by the droplets. With the developed LES/VSJFDF model, numerical simulations of an evaporating two-phase spray flow were carried out. The numerical results with the VSJFDF model, the VFDF (velocity Filtered density function) model and non-FDF (Filtered density function) model were compared with the experimental data, simultaneously. The comparison proved that the LES/VSJFDF model can definitely improve the accuracy of the prediction on the droplets mean diameter. (C) 2014 Elsevier Ltd. All rights reserved.