This article outlines a Monte Carlo approach for simulating crystallization processes both as an alternative to solving the population balance equations (PEE) and as a way of conducting multicoordinate simulations that are too complex to be solved as PBEs. First, simple cases with different combinations of nucleation, growth and aggregation are simulated for batch and continuous systems. The comparison of the results of these test cases with analytical and numerical results shows excellent agreement. New results for two cases with simultaneous growth rate dispersion and aggregation are also discussed. Front these results, a new hypothesis is proposed to explain the observed reduction in the degree of aggregation for larger crystals in the continuous crystallization of potassium sulfate. The simulations show that this experimental observation may be caused by small, slow-growing crystals being swallowed by larger, fast-growing ones.