In recent times, the LES-PBE-PDF framework has been developed to couple large eddy simulation (LES) and population balance models (PBE) for the description of soot formation in turbulent hydrocarbon flames. This approach is based on a modeled evolution equation for the LES-filtered probability density function (pdf) associated with the instantaneous gas composition and soot particle size distribution. Here, the interaction of turbulence with chemical reactions and soot formation can be represented without approximations on part of the chemical and soot formation kinetics, while effects due to turbulent transport and molecular diffusion require closure. In view of an efficient numerical solution scheme, we previously proposed to combine a statistically equivalent reformulation of the joint scalar-number density pdf based on Eulerian stochastic fields with a time-explicit adaptive grid discretization in particle size space and a fractional time stepping scheme. In this article, we present algorithmic aspects and relay implementational details for a consistent semi-discrete formulation of the PBE fractional step as well as an effective dynamic load balancing scheme for both the chemical reaction and PBE fractional steps. Considering soot formation in the Delft III turbulent diffusion flame as a test case, we show that the persisting load imbalance is almost negligible on average and give evidence of linear strong scaling on a modern high performance computer for moderate numbers of compute nodes.