A three-level agglomeration model coupled with crystal growth is developed. It accounts for Brownian, laminar, and turbulent agglomeration. The desupersaturation profiles, the particle size distributions, the average sizes, and variances (or standard deviations), as well as the instantaneous agglomeration degrees for each mechanism, can be calculated as functions of time. The model is applied to the crystallization of an amorphous solid into a crystalline polymorph in a batch crystallizer. A runaway phenomenon is detected for agglomeration when crystals are switching over from the Brownian regime to the laminar one: this switchover significantly affects the desupersaturation curve and the crystal shapes.