The Eulerian stochastic fields (ESF) method, which is based on the transport equation of the
joint subgrid scalar probability density function, is applied to Large Eddy Simulation of a turbulent
dilute spray flame. The approach is coupled with a tabulated chemistry approach to represent
the subgrid turbulence–chemistry interaction. Following a two-way coupled Eulerian–Lagrangian
procedure, the spray is treated as a multitude of computational parcels described in a Lagrangian
manner, each representing a heap of real spray droplets. The present contribution has two objectives:
First, the predictive capabilities of the modeling framework are evaluated by comparing simulation
results using 8, 16, and 32 stochastic fields with available experimental data. At the same time, the
results are compared to previous studies, where the artificially thickened flame (ATF) model was
applied to the investigated configuration. The results suggest that the ESF method can reproduce
the experimental measurements reasonably well. Comparisons with the ATF approach indicate that
the ESF results better describe the flame entrainment into the cold spray core of the flame. Secondly,
the dynamics of the subgrid scalar contributions are investigated and the reconstructed probability
density distributions are compared to common presumed shapes qualitatively and quantitatively in
the context of spray combustion. It is demonstrated that the ESF method can be a valuable tool to
evaluate approaches relying on a pre-integration of the thermochemical lookup-table
