Large Eddy Simulation (LES) is utilized to investigate soot evolution in a series of turbulent nonpremixed bluff body flames featuring different bluff body diameters. The modeling framework relies on recent development in the soot subfilter Probability Density Function (PDF) model that can correctly account for the distribution of soot with respect to mixture fraction, correcting errors in previous soot subfilter PDF models that significantly overpredict soot oxidation. With the previous soot subfilter PDF model, no soot was observed outside of the recirculation zone in past studies on similar bluff body flames. Results obtained with the current LES modeling approach compare favorably with the experimental measurements of the flow field and the soot volume fraction. Notably, the current LES modeling approach correctly predicts large soot volume fractions in the recirculation zone, a decrease in the soot volume fraction through the high-strain neck region, and then an increase again in the downstream jet-like region. Consistent with the experimental measurements, the larger bluff body diameter, with its larger recirculation zone with longer residence times, generates more soot in the recirculation zone and also more soot in the high-strain neck region. Analysis of the soot volume fraction source terms lead to mechanistic understanding of soot evolution in the entirety of the bluff body flames. Most of the soot generated in the recirculation zone is oxidized but some escapes unoxidized and is passively transported through the neck region. Virtually no new soot forms in the downstream jet-like region, and the increase in the soot volume fraction in the jet-like region is due to acetylene-based surface growth of the soot transported through the neck region. This mechanism could not be predicted with the previous soot subfilter PDF model, with the recent soot subfilter PDF model being critical in the understanding of this basic mechanism.