A model for simulating crystallization in imperfectly mixed suspensions was developed based on a population balance equation that takes into account the solid suspension density distributions of different sizes of particles in the crystallizer. A simulation method was developed based on mass balance and population balance. With this method, the particle size distributions in the product and in the crystallizer were predicted. The simulated particle size distributions in the product were verified with experimental data.