The numerical simulation of turbulent dilute spray jet flows with poly-disperse acetone sprays is presented, where the gas-phase model includes a transported joint probability density function (PDF) of the gas velocity and the mixture fraction. The solution of the joint PDF transport equation is achieved through a hybrid Eulerian/Lagrangian Monte Carlo particle method. The simplified Langevin model and the interaction-by-exchange-with-the-mean (IEM) model are used for the velocity and scalar evolution of discrete gas particles, and additional terms account for interaction of the gas and the spray evaporation. The spray dynamics is modeled using a Lagrangian discrete parcel method for the description of droplet motion, heating, and evaporation. The interphase mass, momentum, and energy transfer are considered using the point source approximation for dilute sprays in the gas phase and through appropriate terms in the liquid phase. A convective droplet evaporation model is used, and the infinite conductivity model with consideration of non-equilibrium effects based on the Langmuir-Knudsen law is applied. Numerical results are compared with the experimental data set B of Gounder, Kourmatzis and Masri, Sydney, Australia, in terms of droplet size, liquid volume flux as well as mean and fluctuating axial droplet velocities, where the three different spray flows SP2, SP6, and SP7 with different inlet droplet loadings and turbulence levels are simulated. Generally, good agreement is observed. Moreover, the local joint PDFs of the gas velocity and the mixture fraction are presented and analyzed.