A three-dimensional model of the multiphase flow based on the Eulerian-Eulerian approach was implemented using the FlowVision CFD package and, on this basis, a numerical algorithm for study of evaporation of liquid fuel was developed. The created high-performance complex for the carrier and dispersed phases interaction simulation was validated against the well-studied experimental problem of the evaporation and mixing of kerosene emerging from a flat pre-filming airblast atomizer for gas turbine combustors. In this work, the carrier phase is supposed to be air and kerosene vapors, and the dispersed phase is selected as liquid kerosene. Based on the calculated kerosene evaporation drops distributions, an important parameter that characterizes the spray fineness, Sauter mean diameter, is determined. Numerically calculated in the developed model the evaporation rate and Sauter mean diameter of fuel droplets agreed well with the experimental data. In famous works, the air temperature and pressure varied during the experiments. At the same time, in comparison with the calculated data, a stronger influence on the kerosene evaporation was obtained by air temperature than pressure. The dependence on pressure can be seen in the case of taking into account the corresponding changes in the liquid fuel properties. It is also noted that the initial fuel temperature is an important parameter for evaporation. This can be seen in the results of the kerosene evaporation numerical simulation carried out in this study.