Numerical approximations of a population balance model for coupled batch preferential crystallizers