Antisolvent crystallization is widely used in the production of pharmaceuticals. Although it has been observed experimentally that the crystal size distribution is strongly influenced by the imperfect mixing of the antisolvent with the solution, these effects have not been adequately quantified. In this work, a turbulent computational fluid dynamics (CFD) code was coupled with a multienvironment probability density function (PDF) model, which captures the micromixing in the subgrid scale, and the population balance equation, which models the evolution of the crystal size distribution. The population balance equation (PBE) was discretized along the internal coordinate using a high-resolution central scheme. The presence of solids was addressed by treating the suspension as a pseudo-homogeneous phase with a spatial variation in the effective viscosity. This coupled CFD-PDF-PBE algorithm was applied to an antisolvent crystallization process in an agitated semibatch vessel, where the rising liquid level was modeled by a dynamic mesh. The effects of agitation speed, addition mode, and scale-up on the local primary nucleation and size-dependent growth and dissolution rates, as well as the crystal size distribution, were numerically investigated.