Numerical solutions of population balance models in preferential crystallization