A probability densityjanction (PDF) method is extended to turbulent spray flows. The PDF transport equation of the mixture fraction for turbulent spray flows is deduced, modeled, and solved. The numerical results of the methanol vapor mass fraction for a nonreacting spray, which are obtained using the PDF method, are in good agreement with experimental data and improve the results from the moment closure method Furthermore, the shapes of the probability density function of the mixture fraction at different positions, which are computed by the Monte Carlo method, are presented and analyzed It appears that the spray source changes the value of the mean mixture fraction, but it does not change the shape of its PDF A comparison of the Monte Carlo PDF with the standard beta PDF shows that the beta function fails to describe the shape of the PDF. With the definition of appropriate local maximum and minimum values of the mixture fraction, a modified beta PDF is suitable to reflect the shape of the Monte Carlo PDF very well.