Large-eddy simulation (LES) is a promising technique for accurate prediction of reacting multiphase flows in practical gas-turbine engines. These combustors involve complex physical phenomena of primary atomization of liquid sheet/jet and secondary breakup, droplet evaporation, turbulent mixing of fuel vapor with oxidizer, and combustion dynamics. This chapter summarizes advances made in modeling spray fields with LES of turbulent reacting flows in realistic combustor configurations. Specifically, details of subgrid models for droplet dynamics including breakup, evaporation, deformation, droplet dispersion, and finite-size droplets are presented in the context of an Eulerian–Lagrangian simulation methodology on unstructured grids. Effectiveness of LES with advanced spray models in predicting spray behavior in a patternation study of realistic Pratt and Whitney injector is described.